Search NASA⌕ Search

SEARCH · Search NASA

Results for “Convective mixing layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

A Method for Obtaining High Time and Spatial Resolution Convective Cloud Top Data for the TTL

A method for obtaining high time and spatial resolution convective cloud top data for the TTL Leonhard Pfister, Eric Jensen, Rei Ueyama, Eliot Atlas, and Maria Navarro Convective systems in the tropics have a maximum in the cloud top altitude distribution of about 13.5 km. However, there is a significant tail to this distribution -- a few percent reach the cold point tropopause (CPT) at 16.5 km, and there has been clear evidence of convective mass deposited as high as 19 km in the tropics. The region between 13.5 km and the cold point tropopause is transitional, between the free tropical troposphere where convective mixing dominates, and the stratosphere where slow upward ascent dominates. In this region (the Tropical Tropopause Layer), convective injection, slow ascent, and mixing from midlatitudes all have similar time scales. So, even though only a few percent of convective systems reach the CPT, convection is important. Space Based Lidar and cloud radar measurements have yielded information about long term average statistical distributions of cloud altitude as a function of location. However, we also need time-dependent cloud top altitude and cloud top potential temperature information, primarily to understand the water vapor and TTL cloud distributions. This is because the effect of convection depends on the local temperature, and on the subsequent temperature history. Time dependent cloud top information is also needed to understand short-lived tracers because cross-isentropic flow is time and space dependent. This paper presents a method of obtaining time and space dependent convective cloud top theta (and altitude) information using 3-hourly geostationary brightness temperature data, coupled with global 3 -hourly rainfall estimates and temperature analyses. We explore different mixing algorithms to obtain the most reasonable agreement with near-simultaneous observations by cloudsat and calipso. Observations of short-lived tracers from ATTREX, coupled with short-term trajectories are used to test the method's accuracy. An important caveat is the ambiguity of evaluating convective cloud top altitudes under from combined cloudsat and calipso measurements.

cold point tropopause (CPT)↗

An Automated Detection Methodology for Dry Well-Mixed Layers

The intense surface heating over arid land surfaces produces dry well-mixed layers (WML) via dry convection. These layers are characterized by nearly constant potential temperature and low, nearly constant water vapor mixing ratio. To further the study of dry WMLs, we created a detection methodology and supporting software to automate the identification and characterization of dry WMLs from multiple data sources including rawinsondes, remote sensing platforms, and model products. The software is a modular code written in Python, an open source language. Radiosondes from a network of synoptic stations in North Africa were used to develop and test the WML detection process. The detection involves an iterative decision tree that ingests a vertical profile from an input data file, performs a quality check for sufficient data density, and then searches upward through the column for successive points where the simultaneous changes in water vapor mixing ratio and potential temperature are less than the specified maxima. If points in the vertical profile meet the dry WML identification criteria, statistics are generated detailing the characteristics of each layer in the profile. At the end of the vertical profile analysis, there is an option to plot analyzed profiles in a variety of file formats. Initial results show that the detection methodology can be successfully applied across a wide variety of input data and North African environments and for all seasons. It is sensitive enough to identify dry WMLs from other types of isentropic phenomena such as subsidence layers and distinguish the current day’s dry WML from previous days.

Stephen D. Nicholls↗

An Investigation of High Frequency Motions in the Tropical Tropopause Layer near Convection

Indirect evidence indicates a role for vertical mixing in the Tropical Tropopause Layer (TTL). In the past 20 years, high altitude NASA aircraft such as the ER-2, WB-57, and GLobal Hawk have been making 20hz measurements of vertical velocity and other meteorological parameters in the Upper Tropospere-Lower Stratosphere region, many in the tropics, most recently in connection with the Airborne Tropical TRopopause EXperiment (ATTREX). In the stable environment of the UTLS, high frequency activity occurs in bursts, presumably in connection with nearby convection or strong vertical shear associated with larger scale gravity waves. This paper examines tropical high frequency aircraft data to obtain some basic information about the distribution and character of high frequency activity in vertical velocity in the TTL. In particular, we focus on relating the high frequency activity to nearby tropical convection.

tropopause↗

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗

Aerosol backscatter measurements at 10.6 microns with airborne and ground-based CO2 Doppler lidars over the Colorado High Plains. II - Backscatter structure

Measurements of tropospheric aerosol volume backscatter coefficients at 10.6-microns wavelength were obtained with airborne continuous wave and ground-based pulsed CO2 Doppler lidars over the Colorado High Plains during a 20-day period in summer 1982. A persistent 'background' layer was found between 6- and 10-km altitude, with a generally uniform backscatter mixing ratio of about 10 to the -10th sq m/kg per sr. The upper boundary of this background layer varied with the tropopause height; the lower boundary varied with the strength and diurnal cycle of convective mixing in the planetary boundary layer (PBL). For quiescent meteorological conditions, the transition from the PBL to the background layer was usually very sharp, with backscatter decreases sometimes as large as 3 decades in about 70 m. Sharp gradients were also found at the boundaries of shallow (tens of meters) subvisible cirrus clouds. For less stable conditions, associated with vertical aerosol transport by deep cumuliform clouds, backscatter tended to decrease exponentially with altitude.

Bowdle, David A.↗

Small Scale Motions Observed by Aircraft in the Tropical Tropopause Layer - Convective and Non-Convective Environments

Indirect evidence indicates a role for vertical mixing in the Tropical Tropopause Layer (TTL). In particular, detailed model studies suggest that such vertical mixing may be required to explain the value of the water vapor minimum in the TTL. There have been previous observations during the STEP Tropical aircraft campaign (1987) of bursts of high frequency activity associated with convectively generated gravity waves in the tropical western Pacific. Higher frequency, higher quality measurements from NASA high altitude aircraft (ER-2, WB-57, and Global Hawk) have been made available in the last 20 years. These include measurements of vertical velocity and other meteorological parameters. Most recently, during the ATTREX Global Hawk aircraft mission (Airborne Tropical TRopopause EXperiment), there have been extensive measurements at all altitudes of the TTL in both convective (winter western Pacific) and less convective (winter eastern Pacific) regions. This presentation represents an initial analysis of high frequency small scale (a few km max) meteorological measurements from the ATTREX dataset. We obtain some basic information about the distribution and character of high frequency activity in vertical velocity in the TTL. In particular, we focus on relating the high frequency activity to nearby tropical convection and to vertical shears associated with gravity and inertia-gravity waves.

convection↗

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

Integration of satellite and surface data using a radiative-convective oceanic boundary-layer model

A mixing-line boundary-layer model is used to retrieve cloud-top height from satellite-derived cloud-top temperatures, using 700-hPa National Meteorological Center (NMC) analyses and the Comprehensive Ocean and Atmosphere Data Set (COADS) surface data as supporting datasets. Results are compared with the fixed-lapse-rate method of retrieving boundary-layer depth from sea surface temperatures (SST) and cloud-top temperatures. A radiative-convective equilibrium boundary-layer model is used to retrieve boundary-layer structure given SST and surface wind, satellite cloud-top temperatures and cloud fraction, and the 700-hPa NMC thermodynamic analyses. Good agreement is found between the COADS data and the model solutions for low-level temperature and moisture. This suggests that equilibrium boundary-layer models may be of use over remote oceans in the retrieval of boundary-layer structure.

Betts, Alan K.↗

Simulation of TRMM Microwave Imager Brightness Temperature using Precipitation Radar Reflectivity for Convective and Stratiform Rain Areas over Land

Rain is highly variable in space and time. In order to measure rainfall over global land with satellites, we need observations with very high spatial resolution and frequency in time. On board the Tropical Rainfall Measuring Mission (TRMM) satellite, the Precipitation Radar (PR) and Microwave Imager (TMI) are flown together for the purpose of estimating rain rate. The basic method to estimate rain from PR has been developed over the past several decades. On the other hand, the TMI method of rain estimation is still in the state development, particularly over land. The objective of this technical memorandum is to develop a theoretical framework that helps relate the observations made by these two instruments. The principle result of this study is that in order to match the PR observations with the TMI observations in convective rain areas, a mixed layer of graupel and supercooled water drops above the freezing level is needed. On the other hand, to match these observations in the stratiform region, a layer of snowflakes with appropriate densities above the freezing level, and a melting layer below the freezing level, are needed. This understanding can lead to a robust rainfall estimation technique from the microwave radiometer observations.

Prabhakara, C.↗

LES of Temporally Evolving Mixing Layers by an Eighth-Order Filter Scheme

An eighth-order filter method for a wide range of compressible flow speeds (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) are employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) and Reynolds numbers. The high order filter method is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The value of Mc considered is for the TML range from the quasi-incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self similarity property, compressibility effects and the presence of large-scale structure with shocklets for high Mc) are considered for the LES study. The LES results using the same scheme parameters for all studied cases agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002).

Hadjadj, A↗

Comparative Study of Three High Order Schemes for LES of Temporally Evolving Mixing Layers

Three high order shock-capturing schemes are compared for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) ranging from the quasi-incompressible regime to highly compressible supersonic regime. The considered high order schemes are fifth-order WENO (WENO5), seventh-order WENO (WENO7) and the associated eighth-order central spatial base scheme with the dissipative portion of WENO7 as a nonlinear post-processing filter step (WENO7fi). This high order nonlinear filter method (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The LES results by WENO7fi using the same scheme parameter agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002), whereas results by WENO5 and WENO7 compare poorly with experimental data and DNS computations.

Yee, Helen M. C.↗

The Role of Oceanic Processes in the Initiation of Boreal Winter Intraseasonal Oscillations over the Indian Ocean

Observational analyses and a hierarchy of ocean general circulation model (OGCM) experiments were performed to understand the influence of oceanic processes on the warm sea surface temperature anomalies (SSTAs) prior to the convection initiation of boreal winter intraseasonal oscillations (ISOs), including the Madden-Julian Oscillation (MJO), in the equatorial Indian Ocean. We found 39 strong ISOs that passed over the Indian Ocean Warm Pool region during the November-April season of the 2001-2012 period. 17/39 ISO events initiated in the Seychelles-Chagos Thermocline Ridge (SCTR) before propagating eastward; the remaining events initiated in the southern Arabian Sea (6) or Warm Pool (16) regions. The SCTR event set was notable in that it contained more global-scale MJOs (71-76%), as defined by the RMM and OMI indices, than the WP events (25-44%). Additionally, ~24% (44%) of the SCTR (Warm Pool) events were preceded by strong oceanic process-induced SSTAs of similar magnitude to those of shortwave radiative and turbulent heat fluxes. The Arabian Sea events, however, were not associated with statistically significant SSTA signals prior to convection. Based on a mixed layer heat budget analysis, entrainment and upwelling reduction were the dominant oceanic processes contributing to the warming, in contrast with boreal summer, when horizontal advection dominated. We examined several case studies, including primary MJO events, where oceanic Rossby waves were associated with the entrainment and upwelling reduction. Two simple atmospheric boundary layer convergence models revealed that the SSTAs contributed at least half of the total convergence and suggested that the ocean dynamical effect was responsible for the majority of SSTA-forced convergence for those case studies. These results underscore the need for climate prediction models to accurately represent the ocean structure and processes to include the effects of oceanic predictors.

West, B. Jason↗

Turbulence modeling for high speed compressible flows

The following grant objectives were delineated in the proposal to NASA: to offer course work in computational fluid dynamics (CFD) and related areas to enable mechanical engineering students at North Carolina A&T State University (N.C. A&TSU) to pursue M.S. studies in CFD, and to enable students and faculty to engage in research in high speed compressible flows. Since no CFD-related activity existed at N.C. A&TSU before the start of the NASA grant period, training of students in the CFD area and initiation of research in high speed compressible flows were proposed as the key aspects of the project. To that end, graduate level courses in CFD, boundary layer theory, and fluid dynamics were offered. This effort included initiating a CFD course for graduate students. Also, research work was performed on studying compressibility effects in high speed flows. Specifically, a modified compressible dissipation model, which included a fourth order turbulent Mach number term, was incorporated into the SPARK code and verified for the air-air mixing layer case. The results obtained for this case were compared with a wide variety of experimental data to discern the trends in the mixing layer growth rates with varying convective Mach numbers. Comparison of the predictions of the study with the results of several analytical models was also carried out. The details of the research study are described in the publication entitled 'Compressibility Effects in Modeling Turbulent High Speed Mixing Layers,' which is attached to this report.

Chandra, Suresh↗

Lidar observations of mixed layer dynamics - Tests of parameterized entrainment models of mixed layer growth rate

Ground based lidar measurements of the atmospheric mixed layer depth, the entrainment zone depth and the wind speed and wind direction were used to test various parameterized entrainment models of mixed layer growth rate. Six case studies under clear air convective conditions over flat terrain in central Illinois are presented. It is shown that surface heating alone accounts for a major portion of the rise of the mixed layer on all days. A new set of entrainment model constants was determined which optimized height predictions for the dataset. Under convective conditions, the shape of the mixed layer height prediction curves closely resembled the observed shapes. Under conditions when significant wind shear was present, the shape of the height prediction curve departed from the data suggesting deficiencies in the parameterization of shear production. Development of small cumulus clouds on top of the layer is shown to affect mixed layer depths in the afternoon growth phase.

Boers, R.↗

Vertical Transport Processes for Inert and Scavenged Species: TRACE-A Measurements

The TRACE-A mission of the NASA DC-8 aircraft made a large-scale survey of the tropical and subtropical atmosphere in September and October of 1992. Both In-situ measurements of CO (G. Sachsen NASA Langley) and aerosol size (J. Browell group, NASA Langley) provide excellent data sets with which to constrain vertical transport by planetary boundary layer mixing and deep-cloud cumulus convection. Lidar profiles of aerosol-induced scattering and ozone (also by Bremen) are somewhat require more subtle interpretation as tracers, but the vertical information on layering largely compensates for these complexities. The reason this DC-8 dataset is so useful is that very large areas of biomass burning over Africa and South America provide surface sources of appropriate sizes with which to characterize vertical and horizontal motions; the major limitation of our source description is that biomass burning patterns move considerably every few days, and daily burning inventories are a matter of concurrent, intensive research. We use the Penn State / NCAR MM5 model in an assimilation mode on the synoptic and intercontinental scale, and assess the success it shows in vertical transport descriptions. We find that the general level of emissions suggested by the climatological approach (Will. Has, U. of Montana) appears to be approximately correct, possibly a bit low, for this October, 1992, time period. Vertical transport in planetary boundary layer mixing to 5.5 kin was observed and reproduced in our simulations. Furthermore we find evidence that Blackader "transilient" or matrix-transport scheme is needed, but may require some adaptation in our tracer model: CO seems to exhibit very high values at the top of the planetary boundary layer, a process that stretches the eddy-diffusion parameterization. We will report on progress in improving the deep convective transport of carbon monoxide: the Grail scheme as we used it at 100 kin resolution did not transport enough material to the upper troposphere. We expect to be able to attribute this to either parameterization reasons (inadequacy of this parameterization at the large 100km scale) or other reasons. Nevertheless, the qualitative nature of deep transport by clouds shows up well in the simulations. As for scavengable species, the simulations predict tens of micrograms per standard cubic meter of smoke aerosol in the boundary layer. In a straightforward illustration of our simple bulk-mass scavenging parameterization, to one or two micrograms per standard cubic meter of smoke aerosol in the free troposphere just above the source regions: very high concentrations for the free troposphere. We expect to report on comparisons of these predictions to a variety of observations.

Chatfield, Robert B.↗

Stack plume characterization and model assessment with lidar data

A mobile lidar system was used to obtain plume rise and dispersion measurements for validating plume models. Predictions of rise and spread from the two-thirds law were found to agree well with observations. These predictions were used in combination with Lamb's results for neutrally buoyant particle diffusion in the mixed layer to estimate dispersion when convective turbulence became important.

Weil, J. C.↗

Interface shapes during vertical Bridgman growth of (Pb, Sn)Te crystals

Melt-solid interfaces obtained during vertical Bridgman growth of (Pb, Sn)Te crystals were investigated with a quenching technique. The shapes of these interfaces, revealed by etching longitudinally cut sections, were correlated with the composition variations determined by EMPA. These experiments demonstrated that the interface shape can be changed from concave to convex by moving its location from the edge of the cold zone into the hot zone. The metallography and microsegregation near the melt-solid interface were analyzed in detail. A sharp change in composition above the interface indicated the existence of a diffusion boundary layer 40-90 microns thick. This small diffusion boundary layer is consistent with strong convective mixing in the (Pb, Sn)Te melt.

Huang, YU↗

Distribution of tropical tropospheric water vapor

Utilizing a conceptual model for tropical convection and observational data for water vapor, the maintenance of the vertical distribution of the tropical tropospheric water vapor is discussed. While deep convection induces large-scale subsidence that constrains the turbulent downgradient mixing to within the convective boundary layer and effectively dries the troposphere through downward advection, it also pumps hydrometeors into the upper troposphere, whose subsequent evaporation appears to be the major source of moisture for the large-scale subsiding motion. The development of upper-level clouds and precipitation from these clouds may also act to dry the outflow, thus explaining the low relative humidity near the tropopause. A one-dimensional model is developed to simulate the mean vertical structure of water vapor in the tropical troposphere. It is also shown that the horizontal variation of water vapor in the tropical troposphere above the trade-wind boundary layer can be explained by the variation of a moisture source that is proportional to the amount of upper-level clouds. Implications for the nature of water vapor feedback in global warming are discussed.

Sun, De-Zheng↗