A hybrid data, modeling and deep learning approach for prediction of the small-scale dynamics of the upper troposphere
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
During the 2022 High Ice Water Content (HIWC) Flight Campaign, ice water content (IWC) measurements were made at various locations on the NASA DC-8 Airborne Science Laboratory using the Isokinetic Probe Version 2 (IKP2), Ice Crystal Detectors (ICD), and Robust Probe. Correlations of a wing mounted ICD to IKP2 were made to determine the collision/retention efficiency of the ICD total water content (TWC) element in glaciated conditions. Similarly, correlations of a nose mounted ICD near the pitot probes were made to the IKP2 and to the wing-mounted ICD to estimate the ice concentration factor near the DC-8 pitot probes. Results from the 2022 flight campaign were compared to results from the HIWC RADAR II flight campaign in 2018. There was an apparent reduction in the ICD TWC collision/retention efficiency during the 2022 flight campaign. Although reasons for this discrepancy were explored, more analyses are needed to fully understand the causes. However, it was concluded that an estimated ice concentration factor at the nose of 2.5 was consistent between the two flight campaigns. Additionally, a correlation model was developed to relate measurements from a Nose ICD to freestream IWC as measured by the IKP2 to support IWC measurements from the 2022 Convective Processes Experiment - Cabo Verde (CPEX-CV) flight campaign.
During the 2022 High Ice Water Content (HIWC) Flight Campaign, ice water content (IWC) measurements were made at various locations on the NASA DC-8 Airborne Science Laboratory using the Isokinetic Probe Version 2 (IKP2), Ice Crystal Detectors (ICD), and Robust Probe. Correlations of a wing mounted ICD to IKP2 were made to determine the collision/retention efficiency of the ICD total water content (TWC) element in glaciated conditions. Similarly, correlations of a nose mounted ICD near the pitot probes were made to the IKP2 and to the wing-mounted ICD to estimate the ice concentration factor near the DC-8 pitot probes. Results from the 2022 flight campaign were compared to results from the HIWC RADAR II flight campaign in 2018. There was an apparent reduction in the ICD TWC collision/retention efficiency during the 2022 flight campaign. Although reasons for this discrepancy were explored, more analyses are needed to fully understand the causes. However, it was concluded that an estimated ice concentration factor at the nose of 2.5 was consistent between the two flight campaigns. Additionally, a correlation model was developed to relate measurements from a Nose ICD to freestream IWC as measured by the IKP2 to support IWC measurements from the 2022 Convective Processes Experiment - Cabo Verde (CPEX-CV) flight campaign.
Since 1983 an international group of institutions has collected and analyzed satellite radiance measurements from up to five geostationary and two polar orbiting satellites to infer the global distribution of cloud properties and their diurnal, seasonal and interannual variations. The primary focus of the first phase of the project (1983-1995) was the elucidation of the role of clouds in the radiation budget (top of the atmosphere and surface). In the second phase of the project (1995 onwards) the analysis also concerns improving understanding of clouds in the global hydrological cycle. [Location=GLOBAL] [Temporal_Coverage: Start_Date=1983-07-01; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=112 Km; Longitude_Resolution=112 Km; Temporal_Resolution=5-day].
The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.
The U.S. Department of Energy Atmospheric System Research (ASR)-supported Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in the East River Watershed (ERW) of the Upper Colorado River Basin in southwestern Colorado ran from fall 2021 to spring 2023. Two monitoring sites were deployed in the ERW as part of SAIL. The two sites were the Aerosol Observation System (AOS) located on Crested Butte Ski Mountain, and the second ARM Mobile Facility (AMF2), located at the Rocky Mountain Biological Laboratory in Gothic, Colorado. To gain a more comprehensive understanding of aerosols in complex, mountainous terrain, Handix Scientific deployed SAIL-Net, a distributed network of six measurement nodes spanning the domain of the SAIL research area from October 2021 to July 2023. Each node measured aerosol particles between 140 nm and 3.4 μm in diameter using a small portable optical particle spectrometer (POPS; Gao et al. 2016), cloud condensation nuclei (CNN) using a miniature CCN counter (CloudPuck), and ice nucleating particles (INP) using the time-resolved aerosol filter sampler (TRAPS; Creamean et al. 2018). Our approach was similar to other studies that aimed to better characterize and understand aerosols and gas-phase pollutants using networks of lower-cost sensors (Caubel et al. 2019, Kelly et al. 2021, Asher et al. 2022). Such studies have identified neighborhood-level variations in pollutant concentrations (Schneider et al. 2017, Popoola et al. 2018, Caubel et al. 2019). Small-scale variations such as this are poorly represented in models and poorly measured by a single monitoring system (Caubel et al. 2019). Previous work has shown the representation error (the ability of measurements to represent a larger area) increases with complex orography, leading to decreases in model accuracy (Schutgens et al. 2017). The overall goal of SAIL-Net was to improve our understanding of the variability of aerosol in the ERW, thus increasing our knowledge of aerosol-cloud interactions in this region and informing the usefulness of distributed networks of measurements for future studies. We met this goal by answering the following science questions: 1. What is the aerosol temporal variability, and how does aerosol inhomogeneity vary seasonally? Is there significant seasonal variability in sources, or are short-term meteorological conditions the most important determining factor in sources for cloud nuclei? 2. What is the aerosol spatial variability? What are the aerosol characteristics at cloud base, presumably the particles most representative of those acting as cloud nuclei? 3. How should measurement networks be designed to capture aerosol-cloud interactions, and what do they need to measure? Can a single measurement site accurately represent aerosol properties in regions of complex terrain? SAIL-Net consisted of six measurement nodes spread across the ERW near Crested Butte, Colorado. The primary objective in site placement was to select locations that captured the vertical variation in aerosol properties while also spanning the domain of the SAIL campaign. The elevation of the sites ranged from roughly 2750 m along the valley floor of the ERW to approximately 3500 m near the top of Crested Butte Mountain, which is one of the taller peaks in the ERW. The farthest distance between sites was 14 km, while the closest two sites were approximately 1 km apart. Two of the sites were collocated with the ARM SAIL sites; our instruments sat on top of one of the trailers at AOS and another one of our sites was located in a meadow just above AMF2.
A novel technique is used to estimate derivatives of ice effective radius with respect to height near convective cloud tops (dr(sub e)/dz) from airborne shortwave reflectance measurements and lidar. Values of dr(sub e)/dz are about −6 micrometer/km for cloud tops below the homogeneous freezing level, increasing to near 0 micrometer/km above the estimated level of neutral buoyancy. Retrieved dr(sub e)/dz compares well with previously documented remote sensing and in situ estimates. Effective radii decrease with increasing cloud top height, while cloud top extinction increases. This is consistent with weaker size sorting in high, dense cloud tops above the level of neutral buoyancy where fewer large particles are present and with stronger size sorting in lower cloud tops that are less dense. The results also confirm that cloud top trends of effective radius can generally be used as surrogates for trends with height within convective cloud tops. These results provide valuable observational targets for model evaluation.
An artificial neural network (ANN) algorithm, employing several Aqua MODIS infrared channels, the retrieved total cloud visible optical depth, and vertical humidity profiles is trained to detect multilayer (ML) ice-over-water cloud systems as identified by matched CloudSat and CALIPSO (CC) data. The multilayer ANN, or MLANN, algorithm is also trained to retrieve the optical depth and the top and base heights of the upper-layer ice clouds in ML systems. The trained MLANN was applied to independent MODIS data resulting in a combined ML and single layer hit rate of 80% (77%) for nonpolar regions during the day (night). The results are more accurate than currently available methods and the previous version of the MLANN. Upper-layer cloud top and base heights are accurate to ±1.2 km and ±1.6 km, respectively, while the uncertainty in optical depth is ±0.457 and ±0.556 during day and night, respectively. Areas of further improvement and development are identified and will be addressed in future versions of the MLANN.
This paper is the second in a series in which kilometer-scale-resolving observations from the Atmospheric Radiation Measurement program and a cloud-resolving model (CRM) are used to evaluate the single-column model (SCM) version of the National Centers for Environmental Prediction Global Forecast System model. Part I demonstrated that kilometer-scale cirrus properties simulated by the SCM significantly differ from the cloud radar observations while the CRM simulation reproduced most of the cirrus properties as revealed by the observations. The present study describes an evaluation, through a comparison with the CRM, of the SCM's representation of detrainment from deep cumulus and ice-phase microphysics in an effort to better understand the findings of Part I. It is found that detrainment occurs too infrequently at a single level at a time in the SCM, although the detrainment rate averaged over the entire simulation period is somewhat comparable to that of the CRM simulation. Relatively too much detrained ice is sublimated when first detrained. Snow falls over too deep of a layer due to the assumption that snow source and sink terms exactly balance within one time step in the SCM. These characteristics in the SCM parameterizations may explain many of the differences in the cirrus properties between the SCM and the observations (or between the SCM and the CRM). A possible improvement for the SCM consists of the inclusion of multiple cumulus cloud types as in the original Arakawa-Schubert scheme, prognostically determining the stratiform cloud fraction and snow mixing ratio. This would allow better representation of the detrainment from deep convection, better coupling of the volume of detrained air with cloud fraction, and better representation of snow field.
Saturn atmosphere thermal IR emission observations have uncovered two phenomena that may be accounted for by large ammonia ice particles: the depression of thermal brightness near the equator, and inconsistent IR and radio occultation results. Ammonia ice particles yield the opacity that can account for the contrast between the equatorial region and a brighter area near 15 deg S, and may also reconcile the 45-micron brightness of these two regions with the mean temperature structure of Voyager 2 radio occultation results. A cloud model whose ice particles are distributed in an equal ratio with gas particles up to the 100 mbar pressure level is found to fit the equatorial data, while a thinner cloud, or one that does not extend higher than the 400 mbar limit of the convective region matches the 15 deg S data.
Spectrophotometric observations from 2 to 4 microns and from 8 to 13 microns of several infrared sources associated with molecular clouds are reported. Narrow absorption features at 3.08 microns, attributed to interstellar ices, appear in all sources with a molecular cloud in the intervening line of sight. All sources showing ice absorptions also show broad absorption features, attributed to cold silicates, from 8 to 13 microns. The observed ice absorption profiles are all quite similar; however, they do not fit in detail Mie theory predictions of extinction for pure H2O or NH3 ices. The ratio of ice-to-silicate optical depths is found to vary, with most sources showing a ratio in the range 0.1-0.4. The ratio of visual extinction to ice absorption is found to increase rapidly from inside to outside the molecular cloud in NGC 2024.
To investigate the influence of sea ice openings like leads on wintertime Arctic clouds, the air mass transport is exploited as a heat and humidity feeding mechanism which can modify Arctic cloud properties. Cloud microphysical properties in the central Arctic are analysed as a function of sea ice conditions during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in 2019–2020. The Cloudnet classification algorithm is used to characterize the clouds based on remote sensing observations and the atmospheric thermodynamic state from the observatory on board the research vessel (RV) Polarstern. To link the sea ice conditions around the observational site with the cloud observations, the water vapour transport (WVT) being conveyed towards RV Polarstern has been utilized as a mechanism to associate upwind sea ice conditions with the measured cloud properties. This novel methodology is used to classify the observed clouds as coupled or decoupled to the WVT based on the location of the maximum vertical gradient of WVT height relative to the cloud-driven mixing layer. Only a conical sub-sector of sea ice concentration (SIC) and the lead fraction (LF) centred on the RV Polarstern location and extending up to 50 km in radius and with an azimuth angle governed by the time-dependent wind direction measured at the maximum WVT is related to the observed clouds. We found significant asymmetries for cases when the clouds are coupled or decoupled to the WVT and selected by LF regimes. Liquid water path of low-level clouds is found to increase as a function of LF, while the ice water path does so only for deep precipitating systems. Clouds coupled to WVT are found to generally have a lower cloud base and larger thickness than decoupled clouds. Thermodynamically, for coupled cases the cloud-top temperature is warmer and accompanied by a temperature inversion at the cloud top, whereas the decoupled cases are found to be closely compliant with the moist adiabatic temperature lapse rate. The ice water fraction within the cloud layer has been found to present a noticeable asymmetry when comparing coupled versus decoupled cases. This novel approach of coupling sea ice to cloud properties via the WVT mechanism unfolds a new tool to study Arctic surface–atmosphere processes. With this formulation, long-term observations can be analysed to enforce the statistical significance of the asymmetries. Furthermore, our results serve as an opportunity to better understand the dynamic linkage between clouds and sea ice and to evaluate its representation in numerical climate models for the Arctic system.
Many factors control the ice supersaturation and microphysical properties in cirrus clouds. We explore the effects of dynamic forcing, ice nucleation mechanisms, and ice crystal growth rate on the evolution and distribution of water vapor and cloud properties in nighttime cirrus clouds using a one-dimensional cloud model with bin microphysics and remote sensing measurements obtained at the Department of Energy's Atmospheric Radiation Measurement (ARM) Climate Research Facility located near Lamont, OK. We forced the model using both large-scale vertical ascent and, for the first time, mean mesoscale velocity derived from radar Doppler velocity measurements. Both heterogeneous and homogeneous nucleation processes are explored, where a classical theory heterogeneous scheme is compared with empirical representations. We evaluated model simulations by examining both bulk cloud properties and distributions of measured radar reflectivity, lidar extinction, and water vapor profiles, as well as retrieved cloud microphysical properties. Our results suggest that mesoscale variability is the primary mechanism needed to reproduce observed quantities. Model sensitivity to the ice growth rate is also investigated. The most realistic simulations as compared with observations are forced using mesoscale waves, include fast ice crystal growth, and initiate ice by either homogeneous or heterogeneous nucleation. Simulated ice crystal number concentrations (tens to hundreds particles per liter) are typically two orders of magnitude smaller than previously published results based on aircraft measurements in cirrus clouds, although higher concentrations are possible in isolated pockets within the nucleation zone.
An optical ice particle counter was proposed as a companion instrument to the GSFC laser nephelometer. By counting ice particles and total cloud particles (both ice and liquid water), these two instruments may be used to study the balance between ice and water in clouds.
One of the proposed concepts for mitigating the warming effect of increasing greenhouse gases is seeding cirrus cloud with ice nuclei (IN) in order to reduce the lifetime and coverage of cold cirrus that have a net warming impact on the earth's surface. Global model simulations of the net impact of changing upper tropospheric IN have given widely disparate results, partly as a result of poor understanding of ice nucleation processes in the current atmosphere, and partly as a result of poor representation of these processes in global models. Here, we present detailed process-model simulations of tropical tropopause layer (TTL) transport and cirrus formation with ice nuclei properties based on recent laboratory nucleation experiments and field measurements of aerosol composition. The model is used to assess the sensitivity of TTL cirrus occurrence frequency and microphysical properties to the abundance and efficacy of ice nuclei. The simulated cloud properties compared with recent high-altitude aircraft measurements of TTL cirrus and ice supersaturation. We find that abundant effective IN (either from glassy organic aerosols or crystalline ammonium sulfate with concentrations greater than about 100/L) prevent the occurrences of large ice concentration and large ice supersaturations, both of which are clearly indicated by the in situ observations. We find that concentrations of effective ice nuclei larger than about 50/L can drive significant changes in cirrus microphysical properties and occurrence frequency. However, the cloud occurrence frequency can either increase or decrease, depending on the efficacy and abundance of IN added to the TTL. We suggest that our lack of information about ice nuclei properties in the current atmosphere, as well as uncertainties in ice nucleation processes and their representations in global models, preclude meaningful estimates of climate impacts associated with addition of ice nuclei in the upper troposphere. We will briefly discuss the key field measurements needed to constrain ice nucleation processes.
A comparative analysis between observational data from McMurdo Station, Antarctica and the Community Atmosphere Model version 6 (CAM6) simulation is performed focusing on cloud characteristics and their thermodynamic conditions. Ka-band Zenith Radar (KAZR) and High Spectral Resolution Lidar (HSRL) retrievals are used as the basis of cloud fraction and cloud phase identifications. Radiosondes released at 12-hour increments provide atmospheric profiles for evaluating the simulated thermodynamic conditions. Our findings show that the CAM6 simulation consistently overestimates (underestimates) cloud fraction above (below) 3 km in four seasons of a year. Normalized by total in-cloud samples, ice and mixed phase occurrence frequencies are underestimated and liquid phase frequency is overestimated by the model at cloud fractions above 0.6, while at cloud fractions below 0.6 ice phase frequency is overestimated and liquid-containing phase frequency is underestimated by the model. The cloud fraction biases are closely associated with concurrent biases in relative humidity (RH), i.e., high (low) RH biases above (below) 2 km. Frequencies of correctly simulating ice and liquid-containing phase increase when the absolute biases of RH decrease. Cloud fraction biases also show a positive correlation with RH biases. Water vapor mixing ratio biases are the primary contributor to RH biases, and hence, likely a key factor controlling the cloud biases. Finally, this diagnosis of the evident shortfalls of representations of cloud characteristics in CAM6 simulation at McMurdo Station brings new insight in improving the governing model physics therein.
Globally, on average, the release rate of aerosols from clouds is estimated to be about 6,000 Tg/yr and an aerosol particle is estimated to experience about three cloud cycles. Cloud processing of aerosol produces aerosols that have dramatically different physicochemical properties compared to their original form. Cloud processing may transform particles to efficient cloud condensation nuclei (CCN) or ice nucleating particles (INP) that can further participate in cloud formation. We investigated cloud processing of aerosol during the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) Surface Atmosphere Integrated Field Laboratory (SAIL) field campaign using multi-modal offline analysis techniques available at DOE’s Environmental Molecular Sciences Laboratory (EMSL). The main scientific aim of this project was to understand the processes that transform aerosols (chemically, and physically), and influence aerosol-cloud interactions. Aerosol particles undergo significant modifications due to warm (cloud droplets) and cold cloud (ice crystals) processing. Cloud processing of aerosol has important effects on the atmospheric aerosol number and mass concentration, size distribution, phase separation, and chemical composition.
Abstract. A new empirical parameterization (EP) for multiple groups of primary biological aerosol particles (PBAPs) is implemented in the aerosol–cloud model (AC) to investigate their roles as ice nucleating particles (INPs). The EP describes the heterogeneous ice nucleation by (1) fungal spores, (2) bacteria, (3) pollen, (4) detritus of plants, animals, and viruses, and (5) algae. Each group includes fragments from the originally emitted particles. A high-resolution simulation of a midlatitude mesoscale squall line by AC is validated against airborne and ground observations. Sensitivity tests are carried out by varying the initial vertical profiles of the loadings of individual PBAP groups. The resulting changes in warm and ice cloud microphysical parameters are investigated. The changes in warm microphysical parameters, including liquid water content and cloud droplet number concentration, are minimal (<10 %). Overall, PBAPs have little effect on the ice number concentration (<6 %) in the convective region. In the stratiform region, increasing the initial PBAP loadings by a factor of 1000 resulted in less than 40 % change in ice number concentrations. The total ice concentration is mostly controlled by various mechanisms of secondary ice production (SIP). However, when SIP is intentionally shut down in sensitivity tests, increasing the PBAP loading by a factor of 100 has an effect of less than 3 % on the ice phase. Further sensitivity tests revealed that PBAPs have little effect on surface precipitation and on the shortwave and longwave flux (<4 %) for a 100-fold perturbation in PBAPs.