Search NASA⌕ Search

SEARCH · Search NASA

Results for “Observation nudging”

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 19 records

Comparing the interactions between particulate matter and cloud properties over two populated cities in Texas using WRF-Chem fine-resolution modeling

Accurate modeling of aerosol-cloud interactions is essential for reliable weather and air quality simulations, given their significant impact on precipitation patterns, cloud dynamics, and aerosol distributions. This study employed the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) to examine the impact of enhanced meteorological simulations, achieved through advanced microphysics parameterization supported by data assimilation techniques, on air quality across Texas on August 19 and 20, 2022. We tested four distinct configurations: (1) the Morrison two-moment bulk microphysics scheme, (2) Morrison's with observation nudging, (3) the Spectral Bin Microphysics (SBM), and (4) SBM with observation nudging. While the SBM scheme is known for its detailed representation of aerosol-cloud interactions, our focus was on how improvements in meteorological accuracy translate to more precise air quality simulations. Our findings demonstrated a progressive improvement in simulation accuracy, starting with the Morrison's scheme and further enhanced by adopting the SBM scheme, complemented by incorporating observation nudging. Specifically, the combination of the SBM scheme and the nudging substantially enhanced the model's ability to capture convective precipitation events, as shown by better alignment with NEXRAD radar reflectivity, with R increasing from –0.21 to 0.82, IOA from 0.10 to 0.87, and NMB decreasing from 99% to 34% in Houston. The enhanced meteorology translated into more accurate PM 2.5 concentration simulations, particularly through the more accurate representation of aerosol washout during precipitation events. In Houston, the SBM scheme with nudging improved the model's PM 2.5 simulations significantly, with NMB decreasing from –20% to 5% and IOA improving from 0.43 to 0.61. In San Antonio, improvements were also notable, with NMB improved from –27% to –22%, R increased from 0.48 to 0.82, and IOA increased from 0.66 to 0.86. Furthermore, our results underscore the crucial role of accurate meteorological simulations in refining our understanding of aerosol behaviors in relation to precipitation patterns, directly enhancing the reliability and effectiveness of air quality modeling.

54 ENVIRONMENTAL SCIENCES↗

The MJO-QBO Relationship in a GCM with Stratospheric Nudging

Observational studies show a strong connection between the intraseasonal Madden-Julian oscillation (MJO) and the stratospheric quasi-biennial oscillation (QBO): the boreal winter MJO is stronger, more predictable, and has different teleconnections when the QBO in the lower stratosphere is easterly versus westerly. Despite the strength of the observed connection, global climate models do not produce an MJO-QBO link. Here the authors use a current-generation ocean-atmosphere coupled NASA Goddard Institute for Space Studies global climate model (Model E2.1) to examine the MJO-QBO link. To represent the QBO with minimal bias, the model zonal mean stratospheric zonal and meridional winds are relaxed to reanalysis fields from 1980-2017. The model troposphere, including the MJO, is allowed to freely evolve. The model with stratospheric nudging captures QBO signals well, including QBO temperature anomalies. However, an ensemble of nudged simulations still lacks an MJO-QBO connection.

Madden-Julian oscillation (MJO)↗

Biomass burning aerosol radiative effects in the Southeast Atlantic depend strongly on meteorological forcing method

Biomass burning aerosols (BBAs) from African fires may strongly impact Earth’s radiation budget in the southeast Atlantic (SEA), but the sign and magnitude of the overall radiative effect (RE) remain uncertain. Aerosol–climate models are needed to separately quantify direct, indirect, and semi-direct REs. Here, we evaluate improved simulations with the UK Met Office's Unified Model and explore REs resulting from various methods used to match observed meteorology (nudging or running forecasts reinitialized at different frequencies). REs are calculated as differences in radiative fluxes between simulations with and without smoke emissions and with and without aerosol absorption. All model setups agree on net warming for the SEA dominated by the direct effect. Simulated smoke, clouds, and the direct effect agree better with observations than previous studies using the same model, though biases in aerosol extinction and liquid water path remain. Changes in cloud droplet number concentration due to BBA self-lofting influence how cleanly we can separate cloud effects into semi-direct and indirect effects. Total RE, which remains unaffected, ranges from +3.0 to +7.9 W m −2 . The 4.9 W m −2 spread arises mainly from simulated semi-direct effects. Forecasts three days long or less probably do not allow time for plausible differences in boundary layer properties due to semi-direct effects to accumulate. Free running simulations with and without smoke accumulate differences in meteorology that are likely spurious “butterfly effects”. We recommend future research quantifying BBA REs over weeks to months to use meteorological forcing techniques that allow aerosol absorption to affect the boundary layer.

Giuffrida, Eric [Carnegie Mellon University, Pitts↗

Initialization and Predictability of a Coupled ENSO Forecast Model

The skill of a coupled ocean-atmosphere model in predicting ENSO has recently been improved using a new initialization procedure in which initial conditions are obtained from the coupled model, nudged toward observations of wind stress. The previous procedure involved direct insertion of wind stress observations, ignoring model feedback from ocean to atmosphere. The success of the new scheme is attributed to its explicit consideration of ocean-atmosphere coupling and the associated reduction of "initialization shock" and random noise. The so-called spring predictability barrier is eliminated, suggesting that such a barrier is not intrinsic to the real climate system. Initial attempts to generalize the nudging procedure to include SST were not successful; possible explanations are offered. In all experiments forecast skill is found to be much higher for the 1980s than for the 1970s and 1990s, suggesting decadal variations in predictability.

Chen, Dake↗

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗

Changing effects of external forcing on Atlantic–Pacific interactions

Recent studies have highlighted the increasingly dominant role of external forcing in driving Atlantic and Pacific Ocean variability during the second half of the 20th century. This paper provides insights into the underlying mechanisms driving interactions between modes of variability over the two basins. We define a set of possible drivers of these interactions and apply causal discovery to reanalysis data, two ensembles of pacemaker simulations where sea surface temperatures in either the tropical Pacific or the North Atlantic are nudged to observations, and a pre-industrial control run. We also utilize large-ensemble means of historical simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6) to quantify the effect of external forcing and improve the understanding of its impact. A causal analysis of the historical time series between 1950 and 2014 identifies a regime switch in the interactions between major modes of Atlantic and Pacific climate variability in both reanalysis and pacemaker simulations. A sliding window causal analysis reveals a decaying El Niño–Southern Oscillation (ENSO) effect on the Atlantic as the North Atlantic fluctuates towards an anomalously warm state. The causal networks also demonstrate that external forcing contributed to strengthening the Atlantic's negative-sign effect on ENSO since the mid-1980s, where warming tropical Atlantic sea surface temperatures induce a La Niña-like cooling in the equatorial Pacific during the following season through an intensification of the Pacific Walker circulation. The strengthening of this effect is not detected when the historical external forcing signal is removed in the Pacific pacemaker ensemble. The analysis of the pre-industrial control run supports the notion that the Atlantic and Pacific modes of natural climate variability exert contrasting impacts on each other even in the absence of anthropogenic forcing. The interactions are shown to be modulated by the (multi)decadal states of temperature anomalies of both basins with stronger connections when these states are “out of phase”. We show that causal discovery can detect previously documented connections and provides important potential for a deeper understanding of the mechanisms driving changes in regional and global climate variability.

54 ENVIRONMENTAL SCIENCES↗

Extratropical Large-Scale Traveling Weather Systems in the Southern Hemisphere on Mars

From late-autumn through early-spring, the middle- and high-latitudes of both hemispheres of Mars and its predominantly carbon-dioxide atmosphere support mean equator-to-pole thermal contrasts, and then, support a strong mean westerly polar vortex. Observations from orbiting spacecraft indicate that this intense mean baroclinicity-barotropicity supports large-scale eastward traveling weather systems (i.e., transient, traveling synoptic-period waves, on the order of the Rossby deformation scale). On Earth, extratropical weather disturbances arise from wind-shear instabilities, and these are critical components of the terrestrial global circulation. So it is the case for Mars. Large-scale traveling weather systems on Mars serve as agents in the transport of heat, momentum and scalar and tracer quantities (e.g., atmospheric dust, watervapor, ice clouds, chemical species, etc). Such weather systems interact with other large-scale atmospheric circulation components, namely, quasi-stationary (i.e., forced Rossby) modes; global thermal tidal modes; and then, upon large-/continental- geographical scales, upslope/ down-slope flows amongst high relief, low relief, impact basins, and volcanic rises, and more. The character of Mars' traveling extratropical weather disturbances in its southern hemisphere during late winter through early spring is investigated using a high-resolution Mars global climate model (i.e., Mars GCM), and one from the Agency's Mars Climate Modeling Center (MCMC) based at the NASA Ames Research Center. The climate model includes several complex atmospheric physical packages. With such physics modules, our global climate simulations present comparatively well with observations of the planet's current water cycle (Haberle et al.,2019). The climate model is "forced" with an annual dust cycle (i.e., nudged based on MGS/TES observations). Compared to the northern-hemisphere counterparts, the southern synoptic-period weather disturbances and accompanying frontal waves have smaller meridional and zonal scales, and are less intense. Influences of the zonally asymmetric (i.e., east-west varying) topography on southern large-scale weather are investigated, in addition to large-scale up-slope/down-slope flows and the diurnal cycle. A southern storm zone in late winter and early spring presents in the western hemisphere via orographic influences from the Tharsis highlands, and the Argyre and Hellas impact basins. Geographically localized transient-wave activity diagnostics are constructed that illuminate dynamical differences amongst the simulations and these are presented.

Hollingsworth, Jeffery L.↗

The Impact of Winds on AMOC in a Fully-Coupled Climate Model

Here we investigate the role of the atmospheric circulation in AMOC by comparing a fully-coupled large ensemble, a forced-ocean simulation, and new experiments using a fully-coupled global climate model where winds above the boundary layer are nudged towards reanalysis. When winds are nudged north of 45°N, agreement with RAPID array observations of AMOC at 26.5°N improves across several metrics. The phasing of interannual variability is well-captured due to the response of the local Ekman component in both wind-nudging and forced-ocean simulations, however the variance remains underestimated. The mean AMOC strength is substantially reduced relative to the fully-coupled model large ensemble, which is biased high, due to the impact of winds on surface buoyancy fluxes over the subpolar gyre. Nudging winds towards observations also reduces the 1979-2016 trend in AMOC, suggesting that improvement in the representation of the high-latitude atmosphere is important for projecting long-term AMOC changes.

AMOC↗

Seismic Excitation of the Polar Motion, 1977-1993

The mass redistribution in the earth as a result of an earthquake faulting changes the earth's inertia tensor, and hence its rotation. Using the complete formulae developed by CHAO and GROSS (1987) based on the normal mode theory, we calculated the earthquake-induced polar motion excitation for the largest 11,015 earthquakes that occurred during 1977.0-1993.6. The seismic excitations in this period are found to be two orders of magnitude below the detection threshold even with today's high precision earth rotation measurements. However, it was calculated that an earthquake of only one tenth the size of the great 1960 Chile event, if happened today, could be comfortably detected in polar motion observations. Furthermore, collectively these seismic excitations have a strong statistical tendency to nudge the pole towards approximately 140deg E, away from the actual observed polar drift direction. This non-random behavior, similarly found in other earthquake-induced changes in earth rotation and low-degree gravitational field by CHAO and GROSS (1987), manifests some geodynamic behavior yet to be explored.

Chao, Benjamin Fong↗

The Role of Land‐Atmosphere Feedbacks in Midlatitude Wintertime Surface Temperature Variability

Accurately representing synoptic near-surface temperature variability is crucial to predict weather extremes, yet models remain biased. Existing studies primarily attribute wintertime midlatitude near-surface temperature variability to tropospheric large-scale advection, assuming minimal land influence. However, nudging the model's circulation toward observations yields little improvement in wintertime temperature variance over Northern Hemisphere land, suggesting that land-atmosphere interactions also warrant attention. We introduce a new scaling framework for temperature variance that incorporates local land-atmosphere feedbacks. Comparing our framework to the mixing length approach—which links temperature variance to the meridional temperature gradient and air parcel displacement (mixing length)—shows that land-atmosphere feedbacks are inherently embedded in the mixing length, a connection previously overlooked. Roles of land–atmosphere feedbacks are evaluated via model experiments with perturbed circulation, land, or both. We find that longwave radiative damping dominates temperature variance responses over meridional temperature gradient when both land and circulation are perturbed.

atmosheric science↗

Seismic Excitation of the Polar Motion

The mass redistribution in the earth as a result of an earthquake faulting changes the earth's inertia tensor, and hence its rotation. Using the complete formulae developed by Chao and Gross (1987) based on the normal mode theory, we calculated the earthquake-induced polar motion excitation for the largest 11,015 earthquakes that occurred during 1977.0-1993.6. The seismic excitations in this period are found to be two orders of magnitude below the detection threshold even with today's high precision earth rotation measurements. However, it was calculated that an earthquake of only one tenth the size of the great 1960 Chile event, if happened today, could be comfortably detected in polar motion observations. Furthermore, collectively these seismic excitations have a strong statistical tendency to nudge the pole towards approx. 140 deg E, away from the actually observed polar drift direction. This non-random behavior, similarly found in other earthquake-induced changes in earth rotation and low-degree gravitational field by Chao and Gross (1987), manifests some geodynamic behavior yet to be explored.

Chao, Benjamin Fong↗

Validation of the Archived CERES Surface and Atmosphere Radiation Budget (SARB) at SGP

The CERES Surface and Atmosphere Radiation Budget (SARB) product (Charlock et al, 2002) includes the vertical profile of broadband SW, broadband LW, and 8-12 micron window (WN) fluxes; upwelling and downwelling at TOA, 70 hPa, 200 hPa, 500 hPa, and the surface; and for all-sky and clear-sky conditions. We test the archived CERES TRMM record of SARB for January-August 1998 and focus on discrepancies with ground-based measurements at SGP. The CERES SARB is generated by a highly modified Fu-Liou radiative transfer code (Fu and Liou, 1993). The most critical inputs for this application are cloud optical properties (fractional area, optical depth, particle size and phase, height of top, and estimate of geometrical thickness Minnis et al., 2002) from the narrowband VIRS imager. Numerous VIRS pixels (approx. 2km resolution at nadir) are matched to each of the large (approx. 20km) CERES broadband footprints (Wielicki et al, 1996). Other inputs include temperature and humidity from ECMWF (Rabier et al, 1998) , NCEP ozone profiles from SBUV and TOVS (Yang et al, 2001), aerosol optical thickness (AOT) from the Model for Atmospheric Transport and Chemistry (MATCH) aerosol assimilation (Collins et al., 2001) or alternately from the VIRS imager (Ignatov and Stowe, 2000). VIRS AOT is available for clear and partly cloudy ocean footprints during daylight; and only when viewing geometry renders a contribution from sunglint as unlikely. For other footprints, AOT is taken from MATCH. AOT is apportioned into fractions of dust (Tegan and Lacis, 1996), sea salt, sulfate, dust, soluble organic, insoluble organic, and soot (Hess et al., 1996) using the 6-hourly MATCH output. Tuned fluxes are retrieved by adjusting inputs to nudge computed TOA fluxes toward CERES observations (Rose et al, 1997). In clear conditions, the fields of humidity, surface skin temperature, surface albedo and AOT are adjusted to produce a closer match of computed and observed fluxes at TOA. When CERES footprints have clouds, the cloud optical thickness, fractional area within the footprint, and temperature of cloud top are adjusted by the tuning algorithm. Both tuned and untuned fluxes are archived, as are the respective adjustments to any parameters at the surface or within the atmosphere.

Charlock, Thomas P.↗

Improving thermodynamic nudging in the E3SM Atmosphere Model version 2 (EAMv2): strategy and hindcast skills on weather systems

Nudging techniques are commonly employed to constrain atmospheric simulations toward observed states, facilitating model evaluation and sensitivity studies. However, if applied improperly – particularly to thermodynamic variables such as temperature and humidity – nudging can distort physical processes and introduce spurious biases, undermining the credibility of the simulations. This study presents an improved nudging implementation that applies vertically modulated tendencies to reduce adverse impacts on model physics. The framework is tested in version 2 of the Energy Exascale Earth System Model (EAMv2) using a suite of hindcast simulations nudged toward ERA5 reanalysis. We systematically evaluate the individual and combined effects of nudging wind, temperature, and humidity fields on the model's ability to represent large-scale atmospheric states and high-impact weather systems. Results show that the revised strategy – particularly when nudging temperature and humidity at selected levels – enhances hindcast skill by improving agreement with ERA5 without degrading the hydrological cycle or precipitation processes. Additional improvements in surface temperature, outgoing longwave radiation, and precipitation biases are achieved through targeted nudging of land-surface variables. The proposed approach strengthens the representation of large-scale conditions relevant to tropical cyclones, atmospheric rivers, and extratropical cyclones in the low-resolution EAMv2. These findings demonstrate that carefully designed thermodynamic nudging, especially of temperature and humidity, improves the realism of constrained simulations and broadens the utility of nudged EAMv2 for atmospheric modeling, machine learning, and high-impact weather research.

Atmospheric river↗

Hunga Tonga–Hunga Ha′apai Volcano Impact Model Observation Comparison (HTHH-MOC) project: experiment protocol and model descriptions

The 2022 Hunga volcanic eruption injected a significant amount of water vapor and a moderate amount of sulfur dioxide into the stratosphere, causing observable responses in the climate system. We have developed a model–observation comparison project to investigate the evolution of volcanic water and aerosols and their impacts on atmospheric dynamics, chemistry, and climate, using several state-of-the-art chemistry climate models. The project goals are (1) to evaluate the current chemistry–climate models to quantify their performance in comparison to observations and (2) to understand atmospheric responses in the Earth system after this exceptional event and investigate the potential impacts in the projected future. To achieve these goals, we designed specific experiments for direct comparisons to observations, for example from balloons and the Microwave Limb Sounder satellite instrument. Experiment 1 consists of two sets of free-running ensemble experiments from 2022 to 2031: one with fixed sea-surface temperatures and sea ice and one with coupled ocean. These experiments will help to understand the long-term evolution of water vapor and aerosols; quantify HTHH effects on stratospheric and mesospheric temperatures, dynamics, and transport; understand the impact of dynamic changes on ozone chemistry; quantify the net radiative forcings; and evaluate any surface climate impact. Experiment 2 is a nudged-run experiment from 2022 to 2023 using observed meteorology. To allow participation of more climate models with varying complexities of aerosol simulation, we include two sets of simulations in Experiment 2: Experiment 2a is designed for models with internally generated aerosol, while Experiment 2b is designed for models using prescribed aerosol surface area density. This experiment will help to analyze H 2 O and aerosol evolution, quantify the net radiative forcings, understand the impacts on mid-latitude and polar O 3 chemistry, and allow close comparisons with observations.

Zhu, Yunqian [Univ. of Colorado, Boulder, CO (Unit↗

Present and future aerosol impacts on Arctic climate change in the GISS-E2.1 Earth system model

The Arctic is warming 2 to 3 times faster than the global average, partly due to changes in short-lived climate forcers (SLCFs) including aerosols. In order to study the effects of atmospheric aerosols in this warming, recent past (1990–2014) and future (2015–2050) simulations have been carried out using the GISS-E2.1 Earth system model to study the aerosol burdens and their radiative and climate impacts over the Arctic (>60°N), using anthropogenic emissions from the Eclipse V6b and the Coupled Model Intercomparison Project Phase 6 (CMIP6) databases, while global annual mean greenhouse gas concentrations were prescribed and kept fixed in all simulations. Results showed that the simulations have underestimated observed surface aerosol levels, in particular black carbon (BC) and sulfate (SO2−4), by more than 50 %, with the smallest biases calculated for the atmosphere-only simulations, where winds are nudged to reanalysis data. CMIP6 simulations performed slightly better in reproducing the observed surface aerosol concentrations and climate parameters, compared to the Eclipse simulations. In addition, simulations where atmosphere and ocean are fully coupled had slightly smaller biases in aerosol levels compared to atmosphere-only simulations without nudging. Arctic BC, organic aerosol (OA), and SO2−4 burdens decrease significantly in all simulations by 10 %–60 % following the reductions of 7 %–78 % in emission projections, with the Eclipse ensemble showing larger reductions in Arctic aerosol burdens compared to the CMIP6 ensemble. For the 2030–2050 period, the Eclipse ensemble simulated a radiative forcing due to aerosol–radiation interactions (RFARI) of −0.39 ± 0.01 W/sq. m, which is −0.08 W/sq. m larger than the 1990–2010 mean forcing (−0.32 W/sq. m), of which −0.24 ± 0.01 W/sq. m was attributed to the anthropogenic aerosols. The CMIP6 ensemble simulated a RFARI of −0.35 to −0.40 W/sq. m for the same period, which is −0.01 to −0.06 W/sq. m larger than the 1990–2010 mean forcing of −0.35 W/sq. m. The scenarios with little to no mitigation (worst-case scenarios) led to very small changes in the RFARI, while scenarios with medium to large emission mitigations led to increases in the negative RFARI, mainly due to the decrease in the positive BC forcing and the decrease in the negative SO2−4 forcing. The anthropogenic aerosols accounted for −0.24 to −0.26 W/sq. m of the net RFARI in 2030–2050 period, in Eclipse and CMIP6 ensembles, respectively. Finally, all simulations showed an increase in the Arctic surface air temperatures throughout the simulation period. By 2050, surface air temperatures are projected to increase by 2.4 to 2.6 °C in the Eclipse ensemble and 1.9 to 2.6 °C in the CMIP6 ensemble, compared to the 1990–2010 mean. Overall, results show that even the scenarios with largest emission reductions leads to similar impact on the future Arctic surface air temperatures and sea-ice extent compared to scenarios with smaller emission reductions, implying reductions of greenhouse emissions are still necessary to mitigate climate change.

Aerosols↗

Sensitivity of the MAR Regional Climate Model Snowpack to the Parameterization of the Assimilation of Satellite-Derived Wet-Snow Masks on the Antarctic Peninsula

Both regional climate models (RCMs) and remote sensing (RS) data are essential tools in understanding the response of polar regions to climate change. RCMs can simulate how certain climate variables, such as surface melt, runoff and snowfall, are likely to change in response to different climate scenarios but are subject to biases and errors. RS data can assist in reducing and quantifying model uncertainties by providing indirect observations of the modeled variables on the present climate. In this work, we improve on an existing scheme to assimilate RS wet snow occurrence data with the “Modèle Atmosphérique Régional” (MAR) RCM and investigate the sensitivity of the RCM to the parameters of the scheme. The assimilation is performed by nudging the MAR snowpack temperature to match the presence of liquid water observed by satellites. The sensitivity of the assimilation method is tested by modifying parameters such as the depth to which the MAR snowpack is warmed or cooled, the quantity of water required to qualify a MAR pixel as “wet” (0.1 % or 0.2 % of the snowpack mass being water), and assimilating different RS datasets. Data assimilation is carried out on the Antarctic Peninsula for the 2019–2021 period. The results show an increase in meltwater production (+66.7 % on average, or +95 Gt), along with a small decrease in surface mass balance (SMB) (−4.5 % on average, or −20 Gt) for the 2019–2020 melt season after assimilation. The model is sensitive to the tested parameters, albeit with varying orders of magnitude. The prescribed warming depth has a larger impact on the resulting surface melt production than the liquid water content (LWC) threshold due to strong refreezing occurring within the top layers of the snowpack. The values tested for the LWC threshold are lower than the LWC for typical melt days (approximately 1.2 %) and impact results mainly at the beginning and end of the melting period. The assimilation method will allow for the estimation of uncertainty in MAR meltwater production and will enable the identification of potential issues in modeling near-surface snowpack processes, paving the way for more accurate simulations of snow processes in model projections.

regional climate models↗

Evaluation of the GEOS-Chem UCX Stratosphere in the GEOS Composition Forecast System

The NASA GEOS Composition Forecast (GEOS-CF) system provides 3-dimensional atmospheric composition analyses and forecasts to the public in near-real time at the high spatial resolution of 25 km. While the main focus of this new product is on tropospheric air quality information, the GEOS-Chem chemistry model (v12) used in this system includes the UCX stratospheric chemistry mechanism. Here, we describe the GEOS-CF system and provide comparisons against remote-sensed observations for stratospheric composition, including measurements of HCl, ClO, NO2, and O3. The GEOS-CF nudges the stratospheric ozone towards the GEOS Forward Processing (GEOS FP) assimilated ozone product; as a result the stratospheric ozone analysis in theGEOS-CF agrees well with observations. Additionally, with the inclusion of the GEOS-Chem UCX stratospheric chemistry mechanism in GEOS-CF, 5-day forecasts, especially during the abnormal 2020 NH polar spring, capture the chemical and dynamical changes missed by the GEOS FP system, which tends towards climatology. The GEOS-CF is a new tool for the research community providing near-real time 3-dimensional gridded information on atmospheric composition throughout the troposphere and stratosphere.

Stratosphere↗

NASA GEOS Composition Forecast Modeling System GEOS-CF v1.0: Stratospheric Composition

The NASA Goddard Earth Observing System (GEOS) Composition Forecast (GEOS-29CF) provides recent estimates and five-day forecasts of atmospheric composition to the public in near-real time. To do this, the GEOS Earth system model is coupled with the GEOS-Chem tropospheric-stratospheric unified chemistry extension (UCX) to represent composition from the surface to the top of the GEOS atmosphere (0.01 hPa). The GEOS-CF system is described, including updates made to the GEOS-Chem UCX mechanism within GEOS-CF for improved representation of stratospheric chemistry. Comparisons are made against balloon, lidar and satellite observations for stratospheric composition, including measurements of ozone (O3) and important nitrogen and chlorine species related to stratospheric O3recovery. The GEOS-CF nudges the stratospheric O3towards the GEOS Forward Processing (GEOS FP) assimilated O3product; as a result the stratospheric O3in the GEOS-CF historical estimate agrees well with observations. During abnormal dynamical and chemical environments such as the 2020 polar vortexes, the GEOS-CF O3forecasts are more realistic than GEOS FP O3forecasts because of the inclusion of the complex GEOS-Chem UCX stratospheric chemistry. Overall, the spatial patterns of the GEOS-CF simulated concentrations of stratospheric composition agree well with satellite observations. However, there are notable biases – such as low NOx and HNO3 in the polar regions and generally low HCl throughout the stratosphere – and future improvements to the chemistry mechanism and emissions are discussed. GEOS-CF is a new tool for the research community and instrument teams observing trace gases in the stratosphere and troposphere, providing near-real-time three-dimensional gridded information on atmospheric composition.

GEOS-CF↗