Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shallow cumulus clouds”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Macro-physical Properties of Shallow Cumulus from Integrated ARM Observations (Final Report)

Fair-weather shallow cumuli (ShCu) play an important role in many climate-related processes. Irregular geometry of ShCu and their strong temporal and spatial variability make it challenging to observe ShCu holistically and to represent them correctly in climate models. To improve ShCu parameterizations, information on both vertically and horizontally resolved cloud properties is required. Commonly, the vertically resolved cloud properties are provided by zenith pointing lidar-radar observations with a very narrow field of view (FOV). Thus, these “pencil-beam” properties may not be representative of a larger surrounding area. Limited number of areal-averaged cloud properties, such as fractional sky cover (FSC), are offered typically by wide-FOV observations. The main goal of our project was to integrate advantages of the narrow-FOV (vertical structure of clouds) and wide-FOV (spatial arrangement of clouds) observations for an improved characterization of single-layer ShCu observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site for an 18-yr period (2000-2017). There are four major accomplishments of our project, First, an updated operational cloud classification for days with ShCu has been suggested and evaluated through a detailed comparison with the manually curated records. Our classification extends successfully the latest ARM cloud type Value Added Product (VAP) based on the Active Remote Sensing of Clouds (ARSCL) cloud product by incorporating both cloud fraction (CF) provided by narrow-FOV ceilometer data and FSC from wide-FOV images offered by a Total Sky Imager (TSI). Moreover, our classification allows one to identify impact of instrumentation changes at the SGP site, namely the transition to KAZRARSCL with the updated cloud radar, on the identification of periods with single-layer ShCu. Second, a new approach that resolves cloud area distributions for a given region (up to 4x4 km 2 ) has been suggested and cloud equivalent diameters (CEDs) have been estimated for the first time. These estimations have been performed over a wide range of cloud sizes (about 0.01–3.5 km) with high temporal resolution (30s) using wide-FOV TSI images and cloud base height (CBH) provided by complementary narrow-FOV lidar measurements. Our simple and computationally inexpensive approach offers a previously unavailable dataset for process studies in the convective boundary layer and evaluation of ShCu parameterizations in cloud-resolving models. Third, a long-term integrated record of ShCu macrophysical properties has been developed. The developed record represents the longest available compilation of events with ShCu and includes (i) a novel visualization of the spatial variability in cloud cover both along- and across-wind directions, (ii) updated estimates of narrow-FOV CF and wide-FOV FSC, (iii) updated narrow-FOV CBH, and (iv) complementary data, such as wind speed and direction from the 915-MHz Radar Wind Profiler (RWP) data. The developed record has been used successfully to assess conventional observational estimates of cloud cover and their sensitivity to the following two factors: (i) instrument-dependent cloud detection and data merging criteria and (ii) FOV configuration. Fourth, co-variability of the ShCu macrophysical properties and environmental parameters has been analyzed for a 3-yr period (2016-2018). Our initial analysis includes diurnal changes of FSCs obtained for clouds with small, moderate and large CEDs and several environmental parameters, such as lifted condensation level (LCL) and mixed layer height (zi). Preliminary results of our analysis suggest that the horizontal extent of ShCu is controlled substantially by the sign and magnitude of difference between these two parameters (zi-LCL): the CED tends to grow with increase of this difference (zi exceeds LCL). We have initiated relationships between the ShCu and key atmospheric parameters that control both the development and evolution of ShCu using our new data product, which combines effectively the advantages of narrow-FOV data offered by zenith pointing cloud radars and lidars and wide-FOV TSI images. While the latest instrumentation at the ARM sites may address these challenging relationships in the future, we believe that the historical ARM data at the SGP site has not yet been fully utilized. Overall, our data product can be used by researchers working on a wide range of climate-related projects. These projects may include (i) a comprehensive evaluation of outputs from the Large-Eddy Simulation (LES) and single-column models for their future improvement, (ii) the representativeness of “short-period” results obtained from the previous model and observational studies and (iii) the planning of future field campaigns with focus on improved understanding of the diurnal cycle of cumulus convection.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Automated Tracking of Shallow Maritime Clouds on Geostationary Imagery to Extract Lifecycle Characteristics

Shallow moist convection is ubiquitous throughout the tropics and represents a key player in boundary layer processes. Satellites have provided valuable insight on shallow clouds, such as size, structure, and geographical coverage, from static views of recurring cloud fields. But determining why certain cloud features appear and persist for different periods requires a time-evolving view of their behaviors. Geostationary satellites provide a unique opportunity to follow the time evolution of individual convective features, given their enhanced spatial and temporal sampling. A cloud-tracking tool was developed to identify properties of cloud lifecycle from the NASA Cloud, Aerosol, and Monsoon Processes Philippines Experiment (CAMP2EX) field campaign of 2019. The mission conducted intensive sampling of shallow cumulus in the West Pacific Ocean, in tandem with Rapid Scan imagery from the Advanced Himawari Imager (AHI) on the Japan Meteorological Agency’s (JMA) Himawari-8 satellite. Shallow cumulus were segmented according to thresholds in 0.5-km visible reflectance and with blurring techniques. Despite being limited to daytime hours, the segmentations yielded the best resolution available for capturing cloud initiation, growth, and decay. The tracking procedure is based on a computer vision package that includes Kalman filters for motion prediction, object overlap search, and the Hungarian (or Kuhn-Munkres) matching algorithm for track designation. AHI radiances available within the tracked cloud boundaries are assembled to form individual spectral histories. The resulting catalog provides thousands of cloud histories for domains measuring only a few degrees in latitude and longitude. We present an overview of the cloud-tracking tool and its results for cloud fields sampled throughout CAMP2EX by the airborne P-3. Cloud tracks were selected from about 10 flights to form ensembles, specific groups of tracks occurring in a region with airborne sampling. Cloud lifecycle properties, including duration and maximum area, are calculated for all ensemble members, and analyzed for cloud behavior and P-3 coincidences. By following this strategy, we quantitatively assess the degree of airborne sampling for specific cloud classes defined by the lifecycle calculations. We can summarize which cloud classes had more sampling, the stage of development during sampling, and general differences in character (e.g., isolated congestus vs. cold-pool producer).

Cloud Tracking↗
DOE OSTI2021

Utilizing a Storm-Generating Hotspot to Study Convective Cloud Transitions: The CACTI Experiment

The Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign was designed to improve understanding of orographic cloud life cycles in relation to surrounding atmospheric thermodynamic, flow, and aerosol conditions. The deployment to the Sierras de Córdoba range in north-central Argentina was chosen because of very frequent cumulus congestus, deep convection initiation, and mesoscale convective organization uniquely observable from a fixed site. The C-band Scanning Atmospheric Radiation Measurement (ARM) Precipitation Radar was deployed for the first time with over 50 ARM Mobile Facility atmospheric state, surface, aerosol, radiation, cloud, and precipitation instruments between October 2018 and April 2019. An intensive observing period (IOP) coincident with the RELAMPAGO field campaign was held between 1 November and 15 December during which 22 flights were performed by the ARM Gulfstream-1 aircraft. A multitude of atmospheric processes and cloud conditions were observed over the 7-month campaign, including numerous orographic cumulus and stratocumulus events; new particle formation and growth producing high aerosol concentrations; drizzle formation in fog and shallow liquid clouds; very low aerosol conditions following wet deposition in heavy rainfall; initiation of ice in congestus clouds across a range of temperatures; extreme deep convection reaching 21-km altitudes; and organization of intense, hail-containing supercells and mesoscale convective systems. These comprehensive datasets include many of the first ever collected in this region and provide new opportunities to study orographic cloud evolution and interactions with meteorological conditions, aerosols, surface conditions, and radiation in mountainous terrain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Observations Pertaining to Precipitation within the Northeast Pacific Stratocumulus-to-Cumulus Transition

Three genuine stratocumulus-to-cumulus transitions sampled during the Cloud System Evolution over the Trades (CSET) campaign are documented. The focus is on Lagrangian evolution of in situ precipitation, thought to exceed radar/lidar retrieved values because of Mie scattering. Two of the three initial stratocumulus cases are pristine [cloud droplet number concentrations ($N_d$) of ~22 cm –3 ] but occupied boundary layers of different depths, while the third is polluted ($N_d$ ~ 225 cm –3 ). Hourly satellite-derived cloud fraction along Lagrangian trajectories indicate that more quickly deepening boundary layers tend to transition faster, into more intense but more occasional precipitation. These transitions begin either in the morning or late afternoon, suggesting that preceding night processes can precondition or delay the inevitable transition. The precipitation shifts toward larger drop sizes throughout the transition as the boundary layers deepen, with aerosol concentrations only diminishing in two of the three cases. Ultraclean ($N_d$ < 1 cm –3 ) cumulus clouds evolved from pristine stratocumulus cloud with unusually high precipitation rates occupying a shallow, well-mixed boundary layer. Results from a simple one-dimensional evaporation model and from radar/lidar retrievals suggest subcloud evaporation likely increases throughout the transition. This, coupled with larger drop sizes capable of lowering the latent cooling profile, facilitates the transition to more surface-driven convection. The coassociation between boundary layer depth and precipitation does not provide definitive conclusions on the isolated effect of precipitation on the pace of the transition. Differences between the initial conditions of the three examples provide opportunities for further modeling studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Marine Shallow Cloud Adjustments to the Presence of Shortwave-Absorbing Aerosols: Advancing Understanding Through a Combined Analysis of Lasic Datasets and Process Modeling (Final Report)

This is a final report for this award, whose purpose was to support the analysis of data emanating from the LASIC (Layered Atlantic Smoke Interaction with Clouds) campaign. (Layered Atlantic Smoke Interactions with Clouds) is a strategy to improve our understanding of aged carbonaceous aerosol, its seasonal evolution, and the mechanisms by which clouds adjust to the presence of the aerosol. The observational strategy centers on deploying the AMF1 cloud, aerosol, and atmospheric profiling instrumentation to Ascension Island, located within the trade-wind shallow cumulus regime (14.50W, 80S) 3000 km offshore of continental Africa. The location is within the latitude zone of the maximum outflow of aerosol, with the deepening boundary layer known to entrain free-tropospheric smoke. The primary activities for LASIC are: 1) to improve current knowledge on aged biomass burning aerosol and its radiative properties as a function of the seasonal cycle; 2) to use surface-based remote sensing to sensitively interrogate the atmosphere for the relative vertical location of aerosol and clouds; 3) to improve our understanding of the cloud adjustments to the presence of shortwave-absorbing aerosol within the vertical column, both through aerosol-radiation and through aerosol-cloud interactions; 4) to aid low cloud parameterization efforts for climate models. This award serves to further these activities and includes a sub-award to Dr. Pablo Saide at UCLA for activities 3) and 4).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗
NASA NTRS2023

Dominant Cloud Controlling Factors for Low-Level Cloud Fraction: Subtropical Versus Extratropical Oceans

To improve cloud feedback understanding and simulation, observations have been used to quantify the rate of change of cloud radiative properties as a function of specific environmental metrics (or cloud controlling factors; CCFs). The study focuses on low-level cloud dominated regions during 2006–2010. For each ocean gridpoint, Spearman's rank correlation coefficients of daily mean observed cloud fraction versus (a) 10-m wind, (b) sensible heat flux (SHF), (c) sea surface temperature, (d) estimated inversion strength (EIS), (e) 850 hPa vertical velocity, and (f) the M parameter (Δ𝜃 800hPa skin ) are sorted to identify the dominant CCFs in both extratropics (30°–60°N/S) and subtropics (30°S–30°N). A novel map for visualizing dominant CCFs for low-level cloud fraction reveals that: Δ𝜃 800hPa skin dominates in the subtropical stratocumulus regions while 10-m winds dominate in shallow cumulus regions but in the extratropics, a different inversion structure diagnostic (EIS) dominates, while SHF dominates in western boundary current areas.

cloud controlling factor↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

First Assessment of Cloud‐Land Coupling in LASSO Large‐Eddy Simulations

Abstract To enhance our understanding of cloud simulations over land, this study provides the first assessment of coupling between cloud and land surface in the Large‐Eddy Simulation (LES) Atmospheric Radiation Measurement Symbiotic Simulation and Observation (LASSO) activity for the shallow convection scenario. The analysis of observation data reveals a diurnal cycle of cloud‐land coupling, which co‐varies with surface fluxes. However, coupled (or decoupled) cumulus clouds are inadequately simulated, manifesting as a too‐high (or low) occurrence frequency during the afternoon. This discrepancy is mirrored by the overestimated cloud liquid water path and cloud‐top height. These overestimations are linked to the overpredicted boundary‐layer development and the easier trigger of shallow convection misrepresented in LES runs. Our study underscores the need to improve the representations of boundary‐layer processes and cloud‐land interactions within LES to better simulate shallow clouds in the future.

58 GEOSCIENCES↗
DOE OSTI2021

Clarifying remotely-retrieved precipitation of shallow marine clouds from the NSF/NCAR Gulfstream V

Precipitation is a key process within the shallow cloud lifecycle. The Cloud System Evolution in the Trades (CSET) campaign included the first deployment of a 94 GHz Doppler radar and 532 nm lidar. Despite a larger sampling volume, initial mean radar/lidar retrieved rain rates (Schwartz et al. 2019) based on the upward-pointing remote sensor datasets are systematically less than those measured by in-situ precipitation probes in the cumulus regime. Subsequent retrieval improvements produce rainrates that compare better to in-situ values, but still underestimate. Retrieved shallow cumulus drop sizes can remain too small and too few, with an overestimated shape parameter narrowing the raindrop size distribution too much. Three potential causes for the discrepancy are explored: the gamma functional fit to the dropsize distribution, attenuation by rain and cloud water, and an underaccounting of Mie dampening of the reflectivity. A truncated exponential fit may represent the dropsizes below a showering cumulus cloud more realistically, although further work would be needed to fully evaluate the impact of a different dropsize representation upon the retrieval. The rain attenuation is within the measurement uncertainty of the radar. Mie dampening of the reflectivity is shown to be significant, in contrast to previous stratocumulus campaigns with lighter rain rates, and may be difficult to constrain well with the remote measurements. An alternative approach combines an a priori determination of the dropsize distribution width based on the in-situ data with the mean radar Doppler velocity and reflectivity. This can produce realistic retrievals, although a more comprehensive assessment is needed to better characterize the retrieval errors.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Exploring the Continuum of Stratiform and Cumulus Cloudiness: Observational Insights Related to Marine Boundary Layer Clouds

Low-level marine clouds take on a wide range of morphological states owing to external forcing and internal feedbacks. Rarely are shallow marine clouds adequately described by a single simplified framework but instead exhibit features that partially conform to several states that could be considered “limiting”. Perhaps the most apparent dichotomy is whether to consider shallow clouds as existing within a surface-driven mixed layer or as convective “plume” elements positioned on top. Do the rules that govern the limiting cases transition smoothly through the hybrid scenarios found in nature, or are there particular regions of attraction or repulsion?We have identified 3 ”axes” along which to investigate how the energetics, thermodynamics and cloud macro structure relate to the intermediate spaces between idealizations.

Ewan Crosbie↗
DOE OSTI2025

Using Machine Learning to Generate a GISS ModelE Calibrated Physics Ensemble (CPE)

A neural network (NN) surrogate of the NASA GISS ModelE atmosphere (version E3) is trained on a perturbed parameter ensemble (PPE) spanning 45 physics parameters and 36 outputs. The NN is leveraged in a Markov Chain Monte Carlo (MCMC) Bayesian parameter inference framework to generate a second posterior constrained ensemble coined a “calibrated physics ensemble,” or CPE. The CPE members are characterized by diverse parameter combinations and are, by definition, close to top-of-atmosphere radiative balance, and must broadly agree with numerous hydrologic, energy cycle and radiative forcing metrics simultaneously. Global observations of numerous cloud, environment, and radiation properties (provided by global satellite products) are crucial for CPE generation. The inference framework explicitly accounts for discrepancies (or biases) in satellite products during CPE generation. We demonstrate that product discrepancies strongly impact calibration of important model parameter settings (e.g., convective plume entrainment rates; fall speed for cloud ice). Structural improvements new to E3 are retained across CPE members (e.g., stratocumulus simulation). Notably, the framework improved the simulation of shallow cumulus and Amazon rainfall while not degrading radiation fields, an upgrade that neither default parameters nor Latin Hypercube parameter searching achieved. Analyses of the initial PPE suggested several parameters were unimportant for output variation. However, many “unimportant” parameters were needed for CPE generation, a result that brings to the forefront how parameter importance should be determined in PPEs. From the CPE, two diverse 45-dimensional parameter configurations are retained to generate radiatively-balanced, auto-tuned atmospheres that were used in two E3 submissions to CMIP6.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Adaptive Grid Redistribution for a 1D Model of Turbulence and Clouds

In global atmospheric models, resolving stratocumulus (Sc) in the vertical is computationally expensive. However, Sc appear only under special meteorological conditions. Therefore, there is motivation to refine the vertical grid levels adaptively. In order to facilitate the possibility of parallelization on graphical processing units, our grid adaptation method prescribes the number of vertical levels a priori. Then grid levels are relocated toward altitude ranges in need of refinement. Because the method relocates existing grid levels, rather than adding extra levels, there is a risk of creating regions with overly coarse grid spacing, that is, voids in the grid mesh. To prevent such voids from forming, a simple method is developed to impose a maximum grid spacing. To decide where to place enhanced resolution, the authors develop an empirical mesh refinement criterion. It refines grid spacing near the ground, near strong temperature gradients, and within clouds. Our grid adaptation method is implemented in a single-column model and evaluated on four test cases: decaying stratocumulus, developing shallow cumulus, a quasi-stationary stratocumulus deck, and the diurnal cycle of a dry boundary layer. In the stratocumulus cases, mesh refinement leads to improvements in both the time evolution of fields and their time averages. The other two cases show smaller differences.

Carstensen, Steffen [Univ. of Wisconsin, Milwaukee↗
DOE OSTI2023

Collaborative Proposal: Improving Understanding of the Internal Structure and Dynamics of Deep Convection Using ARM Observations and Large Eddy Simulations (Final progress report)

This is the final technical report for the DOE Atmospheric System Research funded project entitled "Collaborative Proposal: Improving Understanding of the Internal Structure and Dynamics of Deep Convection Using ARM Observations and Large Eddy Simulations". This project addressed several science questions related to the structure and behavior of deep convective clouds in the atmosphere. It applied knowledge gained from this work to improve the deep cumulus convection parameterization in the Community Earth System Model (CESM). Specific outcomes related to the project include an improved understanding of: (1) how convective cloud environments regulate the size of thermals reaching the upper troposphere, and hence the shallow-to-deep convective transition; (2) the role of dynamic pressure forcing in thermal and convective cloud evolution; and (3) the role of vertical wind shear on the behavior of moist thermals and the transition from shallow to deep convection. Although the primary motivation of this project was to improve basic understanding of thermal dynamics in atmospheric convection, an important component of this work was using this improved understanding to modify the representation of entrainment and vertical velocity profiles in the Zhang-McFarlane cumulus parameterization used in CESM. This work also had a strong observational component, with analyses of Atmospheric Radiation Measurement (ARM) observations to characterize thermal characteristics and behavior. Observations collected during the ARM-supported CACTI field project were used with high-resolution simulations to explore environmental controls on deep convection initiation. Data from the ARM-supported MC3E field project were also used to evaluate the representation of updraft vertical velocities in the Zhang-McFarlane cumulus parameterization. This project directly supported the training of two post-doctoral researchers and led to 16 peer-reviewed published papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Influences of Cloud Microphysics on the Components of Solar Irradiance in the WRF-Solar Model

An accurate forecast of Global Horizontal solar Irradiance (GHI) and Direct Normal Irradiance (DNI) in cloudy conditions remains a major challenge in the solar energy industry. This study focuses on the impact of cloud microphysics on GHI and its partition into DNI and Diffuse Horizontal Irradiance (DHI) using the Weather Research and Forecasting model specifically designed for solar radiation applications (WRF-Solar) and seven microphysical schemes. Three stratocumulus (Sc) and five shallow cumulus (Cu) cases are simulated and evaluated against measurements at the US Department of Energy’s Atmospheric Radiation Measurement (ARM) user facility, Southern Great Plains (SGP) site. Results show that different microphysical schemes lead to spreads in simulated solar irradiance components up to 75% and 350% from their ensemble means in the Cu and Sc cases, respectively. The Cu cases have smaller microphysical sensitivity due to a limited cloud fraction and smaller domain-averaged cloud water mixing ratio compared to Sc cases. Cloud properties also influence the partition of GHI into DNI and DHI, and the model simulates better GHI than DNI and DHI due to a non-physical error compensation between DNI and DHI. The microphysical schemes that produce more accurate liquid water paths and effective radii of cloud droplets have a better overall performance.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

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↗
NASA NTRS2021

Identifying Meteorological Influences on Marine Low Cloud Mesoscale Morphology Using Satellite Classifications

Marine low cloud mesoscale morphology in the southeastern Pacific Ocean is analyzed using a large dataset of machine-learning generated classifications spanning three years. Meteorological variables and cloud properties are composited 10by mesoscale cloud type, showing distinct meteorological regimes of marine low cloud organization from the tropics to the midlatitudes. The presentation of mesoscale cellular convection, with respect to geographic distribution, boundary layer structure, and large-scale environmental conditions, agrees with prior knowledge. Two tropical and subtropical cumuliform boundary layer regimes, suppressed cumulus and clustered cumulus, are studied in detail. The patterns in precipitation, circulation, column water vapor, and cloudiness are consistent with the representation of marine shallow mesoscale convective 15 self-aggregation by large eddy simulations of the boundary layer. Although they occur under similar large-scale conditions, the suppressed and clustered low cloud types are found to be well-separated by variables associated with low-level mesoscale circulation, with surface wind divergence being the clearest discriminator between them, whether reanalysis or satellite observations are used. Clustered regimes are associated with surface convergence and suppressed regimes are associated with surface divergence.

Johannes Mohrmann↗
DOE OSTI2025

Shallow cumulus entrainment: observational retrieval, physical interpretation, and climate impacts (Final grant report)

This document constitutes a final technical report for DE-SC0020083, a research grant awarded by the US Department of Energy Office of Science, to Prof. Daniel Kirshbaum at McGill University. The original research proposal, entitled "Shallow cumulus entrainment: observational retrieval, physical interpretation, and climate impacts", defined four scientific objectives that synthesized ARM observations with turbulence-permitting large-eddy simulations (LES) to quantify and understand the process of cumulus entrainment, whereby a cumulus cloud ingests surrounding air through its lateral edges. These objectives included developing observational climatologies of cumulus bulk entrainment rate at two ARM observatories, physical interpretation of these findings using LES, quantification of the impacts of subcloud circulations on entrainment, and the climate impacts of flow-dependent entrainment variations. In the six years that have elapsed since the grant was awarded (three funded years and three no-cost annual extensions), all four of the objectives were successfully addressed. Four personnel were trained (one B.Sc., one M.Sc. to completion, five years of Ph.D. training, and 2.5 years of postdoctoral training), eight peer-reviewed articles were published in high-impact journals, and 20 presentations were delivered (three invited).

54 ENVIRONMENTAL SCIENCES↗