Characterization of permeblity variation with depth in compartmentalized deep aquifers, geothermal and hydrocarbon reservoirs is important for prediction of flow and transport in complex subsurface environments and directly affects the development of natural and energy resources. In deep formations, temperature gradient can be significant and temperature data can provide valuable information about fluid displacement and conductivity in the vertical extent of the formation. This paper examines the importance of temperature data in resolving permeability distribution with depth by integrating flow and temperature data jointly and individually. We show that incorporation of temperature data in model calibration of deep aquifers can increase the resolution of permeability distribution profile with depth. To illustrate the importance of temperature measurements, we adopt a coupled heat and fluid flow model as a forward model to predict the heat and fluid transport in an a deep reservoir and perform a series of numerical experiments for integration of flow data alone, temperature data alone, and flow and temperature data jointly. For model calibration, we use the Maximum A-Posteriori (MAP) estimation approach and for uncertainty quantification we apply the Randomized Maximum Likelihood (RML). We develop an adjoint model for the coupled fluid and heat flow system of equations to compute the required gradients for model calibration. Investigation of the sensitivity of temperature and production data to the distribution of permeability shows that while fluid flow rate data can primarily resolve the distribution of permeability in the lateral extent of the reservoir, the fluid temperature data, even when measured at the surface, show sensitivity to permeability variability with depth, allowing for a higher resolution profiling of the permeability map. The results elucidate the value of temperature data in enhancing the resolution of the estimated aquifer permeability maps with depth.