Trends in vertical wind velocity variability reveal cloud microphysical feedback
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Barahona, Donifan.
Explore the source record for details and available documents.
Abstract Human activities affect the Earth’s climate through modifying the composition of the atmosphere, which then creates radiative forcing that drives climate change. The warming effect of anthropogenic greenhouse gases has been partially balanced by the cooling effect of anthropogenic aerosols. In 2020, fuel regulations abruptly reduced the emission of sulfur dioxide from international shipping by about 80% and created an inadvertent geoengineering termination shock with global impact. Here we estimate the regulation leads to a radiative forcing of $$+0.2\pm 0.11$$ + 0.2 ± 0.11 Wm −2 averaged over the global ocean. The amount of radiative forcing could lead to a doubling (or more) of the warming rate in the 2020 s compared with the rate since 1980 with strong spatiotemporal heterogeneity. The warming effect is consistent with the recent observed strong warming in 2023 and expected to make the 2020 s anomalously warm. The forcing is equivalent in magnitude to 80% of the measured increase in planetary heat uptake since 2020. The radiative forcing also has strong hemispheric contrast, which has important implications for precipitation pattern changes. Our result suggests marine cloud brightening may be a viable geoengineering method in temporarily cooling the climate that has its unique challenges due to inherent spatiotemporal heterogeneity.
Atmospheric models with typical resolution in the tenths of kilometers cannot resolve the dynamics of air parcel ascent, which varies on scales ranging from tens to hundreds of meters. Small-scale wind fluctuations are thus characterized by a subgrid distribution of vertical wind velocity W with standard deviation σ W . The parameterization of σ W is fundamental to the representation of aerosol–cloud interactions, yet it is poorly constrained. Using a novel deep learning technique, this work develops a new parameterization for σ W merging data from global storm-resolving model simulations, high-frequency retrievals of W , and climate reanalysis products. The parameterization reproduces the observed statistics of σ W and leverages learned physical relations from the model simulations to guide extrapolation beyond the observed domain. Incorporating observational data during the training phase was found to be critical for its performance. The parameterization can be applied online within large-scale atmospheric models, or offline using output from weather forecasting and reanalysis products.
Recently NASA's Global Modeling and Assimilation Office (GMAO) has developed a new Subseasonal to Seasonal Prediction system Version 3 (GEOS-S2S-3). This upgrade replaces the GEOS-S2S-2 which is NASA's current contribution to the North American Multi-Model Experiment seasonal prediction project (Kirtman et al., 2014). The main improvements for our S2S-3 system include 1) a higher resolution MOM5 (Griffies et al., 2005) ocean model (now 0.25o x 0.25o x 50 layers), 2) an improved atmospheric/ocean interface layer (Akella and Suarez, 2018), and 3) assimilation of a long-track satellite salinity into the ocean model (Hackert et al, 2019). Atmospheric forcing is provided by the NASA MERRA-2 reanalysis (Gelaro et al., 2017). Initialization for the ocean relies on the GMAO ocean reanalysis system which assimilates all available in situ temperature and salinity, satellite sea surface salinity, and sea level using the Local Ensemble Transform Kalman Filter (LETKF) implementation of (Penny et al., 2013) on a 5 day assimilation cycle with 20 fixed ensemble members.In this presentation, we will authenticate our new S2S-3 ocean reanalysis using standard GODAE validation metrics. For example, we will compare gridded fields of mean and standard deviation of the ocean reanalysis versus observed fields. We will show correlation/RMS of model versus observations and temperature and salinity mean profiles for the various basins and latitude bands. Basin-scale volume transports, such as the Atlantic Meridional Overturning Circulation and the Indonesian Throughflow will be validated. Equatorial ocean waves will be compared by decomposing sea level into Kelvin and Rossby components. For each of these metrics, we plan to validate the results and then compare our new S2S-3 against the current production version, S2S-2. Finally, we will compare 9-month seasonal forecasts initialized from these two systems for the tropical Pacific NINO3.4 region over the period 1981-present.
The NASA/Goddard Global Modeling and Assimilation Office (GMAO) released Version 2 of the Subseasonal to Seasonal (GEOS-S2S) forecast system in the fall of 2017, and it has been producing near-real time subseasonal to seasonal forecasts and a weakly coupled atmosphere-ocean data assimilation record since then. A new version of the coupled modeling and analysis system (Version 3) was released by the GMAO at the end of 2019. The new version runs at higher oceanic resolution than the previous (approximately 1/2 degree for the atmosphere, 1/4 degree for the ocean), and includes interactive earth system model components not typically present in seasonal prediction systems (two moment cloud microphysics for aerosol indirect effect and an interactive aerosol model). The weakly coupled atmosphere-ocean data assimilation system now includes assimilation of sea surface salinity, that has been shown to result in improved ocean mixed layer simulation and ENSO prediction skill.
Ice nucleation active site densities and nucleation rates are commonly determined by measuring the fraction of droplets that freeze at different temperatures. In analysing these experiments it is commonly assumed that variability in the immersed particle surface area and composition within a sample can be neglected, allowing a simple parameterization of ice nucleation efficiencies. Here it is shown that these approximations may result in large error, and translate into a bias in estimated active site densities of about six orders of magnitude. When applied in cloud models it may represent up to 5 K error in cloud glaciation temperature. A solution to this enduring problem is proposed for the first time. It is shown that by performing droplet freezing experiments varying mean particle surface area along with the temperature such a bias can be quantitatively removed. This amounts to differentiating the "freezing rate'' from the ice nucleation coefficient, or, the apparent and the actual active site density of a material. The proposed method will contribute to reduce existing uncertainties in the relative strength of the ice nucleation ability of different aerosol species, and the modeling of mixed-phase clouds.
Observational paucity is a reality in the Arctic Ocean. This is especially true for near-surface variables such as temperature, moisture, heat fluxes and BL clouds. As a result, modeling has become one of the major avenues for understanding current and future Arctic trends. Reanalyses are frequently used to force global ocean circulation and sea-ice models. But in northern high latitudes, model integrations and reanalyses are known to have large uncertainties in temperature and humidity profiles, and in boundary layer cloudiness. These are common sources of error in the surface radiative budget terms.An important way to diagnose these biases spatially and temporally is by using satellite remote sensing data. However, remotely-sensed observations also have large uncertainties, especially in near-surface temperature and relative humidity profiles.In situ observational studies are important in bridging our knowledge gap in regions such as the Arctic Ocean. Here, we utilize airborne and ship observations from the ARISE, ACME-IV, and ASCOS campaigns to construct critical relative humidity (RH) profiles over the Beaufort Sea. Such profiles are used as parameterization inputs in the NASA GOES global model to derive the total water condensate in a model grid-box, which determines the cloud fraction. Currently, the critical RH profiles are derived by global AIRS data, relaying mostly on mid-latitude regions, which are not necessarily relevant to the Arctic.We derive campaign-wide mean, standard deviation, and critical RH values, for grid size of 50x50 km and altitude bins between 50 to 400 m, covering both open-ocean and sea-ice covered regions. We compare profiles over open ocean and sea-ice, and look at correlations between the observed critical RH values and water condensate (by cloud number concentration) from observations versus the modeled ones. We then input our calculated values of minimal critical RH values into a set of GEOS single column model (SCM) simulations over the ARISE and ASCOS regions and compare the differences between the predicted values of cloud liquid water path (LWP), ice water path (IWP) and surface fluxes with the observed ones under the range of input parameterizations. Finally, we discuss the implications on surface radiative budget predictions in this region.
Aerosol emissions modify the properties of clouds hence impacting climate. The aerosol indirect effect may have offset part of the global warming caused by anthropogenic greenhouse gas emissions during the industrial era. It however remains unclear whether the same effect is significant over time scales relevant for seasonal and weather climate prediction. Answering such a question has been difficult since most weather prediction systems lack a proper representation of the aerosol evolution and transport and their interaction with clouds. Even in advanced systems it is not clear to what extent cloud microphysical properties are predictable over subseasonal to seasonal time scales. Such an issue is addressed in this study. We use a set of 30 year, four ensemble member, 9 month lead hindcast simulations of the NASA GEOS seasonal prediction system (GEOS-S2S) to study the predictability of cloud droplet number concentration in warm stratocumulus clouds. The latest version GEOS-S2S system implements interactive aerosol as well as a two moment cloud microphysics scheme therefore it is suitable for studying the aerosol indirect effect on climate. Long term retrievals from the MODIS (Moderate Resolution Imaging Spectroradiometer) are used to validate the model predictions and assess its skill in predicting cloud droplet number concentration.
A new version of the coupled modeling and analysis system used to produce near real time subseasonal to seasonal forecasts was recently released by the NASA/Goddard Global Modeling and Assimilation Office. The new version runs at higher atmospheric resolution than the previous, (approximately 1/2 degree globally), contains a substantially improved model description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter, and now includes the assimilation of along-track sea surface height. Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will also present results from a series of retrospective seasonal forecasts. Results show significant improvements in surface temperatures over much of the northern hemisphere and a much improved prediction of sea ice extent in both hemispheres. Analysis of the ensemble spread shows improvements relative to the previous system, including generally better reliability. The precipitation forecast skill is comparable to previous S2S systems, and the only tradeoff is an increased "double ITCZ", which is expected as we go to higher atmospheric resolution.
A comprehensive ice nucleation parameterization has been implemented in the global chemistry-climate model EMAC to realistically represent ice crystal number concentrations. The parameterization of Barahona and Nenes (2009, hereafter BN09) allows the treatment of ice nucleation, taking into account the competition for water vapour between homogeneous and heterogeneous nucleation and pre-existing ice crystals in cold clouds. Furthermore, the influence of chemically-heterogeneous, polydisperse aerosols is considered via multiple ice nucleating particle spectra, which are included in the parameterization to compute the heterogeneously formed ice crystals. BN09 has been implemented to operate both in the cirrus and in the mixed-phase regimes. Compared to the standard EMAC results, BN09 produces fewer ice crystals in the upper troposphere but higher ice crystal number concentrations in the middle troposphere, especially in the Northern Hemisphere where ice nucleating mineral dust particles are relatively abundant. The comparison with a climatological data set of aircraft measurements shows that BN09 used in the cirrus regime improves the model results and, therefore, is recommended for future EMAC simulations.
In this presentation we present an overview of the GMAO Sub-Seasonal and Seasonal Prediction System, current users and products, and methods for validation and evaluation of the system. Methods for evaluation include baseline evaluations metrics, the ability to simulate key modes of variability, and evaluation of new development areas.
In this presentation we present an overview of the GMAO Sub-Seasonal and Seasonal Prediction System with a focus on the computing time and resources and actual time it takes to complete a full set of hindcasts. The goal is to come up with some solutions to allow us to run more ensemble members for the next version of the system which will be higher resolution and take many more resources.
It is known that ice nucleating particles (INP) immersed within supercooled droplets promote the formation of ice. Common theoretical models used to represent this process assume that the immersed particle lowers the work of ice nucleation without significantly affecting the dynamics of water in the vicinity of the particle. This is contrary to evidence showing that immersed surfaces significantly affect the viscosity and diffusivity of vicinal water. To study how this may affect ice formation this work introduces a model linking the ice nucleation rate to the modification of the dynamics and thermodynamics of vicinal water by immersed particles. It is shown that INP that significantly reduce the work of ice nucleation also pose strong limitations to the growth of the nascent ice germs. This leads to the onset of a new ice nucleation regime, called spinodal ice nucleation, where the dynamics of ice germ growth instead of the ice germ size determines the nucleation rate. Nucleation in this regime is characterized by an enhanced sensitivity to particle area and cooling rate. Comparison of the predicted ice nucleation rate against experimental measurements for a diverse set of species relevant to cloud formation suggests that spinodal ice nucleation may be common in nature.
Weather and climate models have long struggled to realistically simulate the Madden-Julian Oscillation (MJO). Here we present a significant improvement in MJO simulation in NASA's GEOS atmospheric model with the implementation of 2-moment microphysics and the UW shallow cumulus parameterization. Comparing ten-year runs (2007-2016) with the old (1mom) and updated (2mom+shlw) model physics, the updated model has increased intra-seasonal variance with increased coherence. Surface fluxes and OLR are found to vary more realistically with precipitation, and a moisture budget suggests that changes in rain reevaporation and the cloud longwave feedback help support heavy precipitation. Preliminary results also show improved MJO hindcast skill.
Immersion freezing is likely involved in the initiation of precipitation and determines to large extent the phase partitioning in convective clouds. Theoretical models commonly used to describe immersion freezing in atmospheric models are based on the classical nucleation theory which however neglects important interactions near the immersed particle that may affect nucleation rates. This work introduces a new theory of immersion freezing based on two premises. First, immersion ice nucleation is mediated by the modification of the properties of water near the particle-liquid interface, rather than by the geometry of the ice germ. Second, the same mechanism that leads to the decrease in the work of germ formation also decreases the mobility of water molecules near the immersed particle. These two premises allow establishing general thermodynamic constraints to the ice nucleation rate. Analysis of the new theory shows that active sites likely trigger ice nucleation, but they do not control the overall nucleation rate nor the probability of freezing. It also suggests that materials with different ice nucleation efficiency may exhibit similar freezing temperatures under similar conditions but differ in their sensitivity to particle surface area and cooling rate. Predicted nucleation rates show good agreement with observations for a diverse set of materials including dust, black carbon and bacterial ice nucleating particles. The application of the new theory within the NASA Global Earth System Model (GEOS-5) is also discussed.
A new version of the modeling and analysis system used to produce sub-seasonal to seasonal forecasts has just been released by the NASA Goddard Global Modeling and Assimilation Office. The new version runs at higher atmospheric resolution (approximately 12 degree globally), contains a substantially improved model description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter. Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will also present results from a free-running coupled simulation with the new system and results from a series of retrospective seasonal forecasts. Results from retrospective forecasts show significant improvements in surface temperatures over much of the northern hemisphere and a much improved prediction of sea ice extent in both hemispheres. The precipitation forecast skill is comparable to previous S2S systems, and the only trade off is an increased double ITCZ, which is expected as we go to higher atmospheric resolution.
Along with minimizing parameter uncertainty, understanding the cause of temporal and spatial variability of the nucleated ice crystal number, Ni, is key to improving the representation of cirrus clouds in climate models. To this end, sensitivities of Ni to input variables like aerosol number and diameter provide valuable information about nucleation regime and efficiency for a given model formulation. Here we use the adjoint model of the adjoint of a cirrus formation parameterization (Barahona and Nenes, 2009b) to understand Ni variability for various ice-nucleating particle (INP) spectra. Inputs are generated with the Community Atmosphere Model version 5, and simulations are done with a theoretically derived spectrum, an empirical lab-based spectrum and two field-based empirical spectra that differ in the nucleation threshold for black carbon particles and in the active site density for dust. The magnitude and sign of Ni sensitivity to insoluble aerosol number can be directly linked to nucleation regime and efficiency of various INP. The lab-based spectrum calculates much higher INP efficiencies than field-based ones, which reveals a disparity in aerosol surface properties. Ni sensitivity to temperature tends to be low, due to the compensating effects of temperature on INP spectrum parameters; this low temperature sensitivity regime has been experimentally reported before but never deconstructed as done here.
One of the overall scientific goals of the NASA Hurricane and Severe Storm Sentinel (HS-3) field campaign is to better understand the role of the Saharan Air Layer (SAL) in tropical storm development. During the 2012 HS-3 deployment, the Cloud Physics Lidar (CPL) observed dust within SAL air in close proximity to a developing Nadine (September 11, 2012). Throughout the mission, the NASA GEOS-5 modeling system supported HS-3 by providing 0.25 degrees resolution 5-day global forecasts of aerosols, which were used to support mission planning. The aerosol module was radiatively interactive within the GEOS-5 model, but aerosols were not directly coupled to cloud and precipitation processes. In this study we revisit the aerosol forecasts with an updated version of the GEOS-5 model. For the duration of Hurricane Nadine, we run multiday climate simulations leading up to each respective Global Hawk flight with and without aerosol direct interaction. For each set of simulations, we compare simulated dust mass fluxes to identify differences in SAL entrainment related to the interaction between dust aerosols and the atmosphere. We find that the direct effects of dust induce a low level anticyclonic circulation that temporarily shields Nadine from the intrusion of dry air, leading to a more intense storm.