Goddard Earth Observing System Research Articles

Contributions of Initial Conditions and Meteorological Forecast to Subseasonal-to-Seasonal Hydrological Forecast Skill in Western Tropical South America

Abstract Hydrological predictions at subseasonal-to-seasonal (S2S) time scales can support improved decision-making in climate-dependent sectors like agriculture and hydropower. Here, we present an S2S hydrological forecasting system (S2S-HFS) for western tropical South America (WTSA). The system uses the global NASA Goddard Earth Observing System S2S meteorological forecast system (GEOS-S2S) in combination with the generalized analog regression downscaling algorithm and the NASA Land Information System (LIS). In this implementation study, we evaluate system performance for 3-month hydrological forecasts for the austral autumn season (March–May) using ensemble hindcasts for 2002–17. Results indicate that the S2S-HFS generally offers skill in predictions of monthly precipitation up to 1-month lead, evapotranspiration up to 2 months lead, and soil moisture content up to 3 months lead. Ecoregions with better hindcast performance are located either in the coastal lowlands or in the Amazon lowland forest. We perform dedicated analysis to understand how two important teleconnections affecting the region are represented in the S2S-HFS: El Niño–Southern Oscillation (ENSO) and the Antarctic Oscillation (AAO). We find that forecast skill for all variables at 1-month lead is enhanced during the positive phase of ENSO and the negative phase of AAO. Overall, this study indicates that there is meaningful skill in the S2S-HFS for many ecoregions in WTSA, particularly for long memory variables such as soil moisture. The skill of the precipitation forecast, however, decays rapidly after forecast initialization, a phenomenon that is consistent with S2S meteorological forecasts over much of the world.

Improved simulations of biomass burning aerosol optical properties and lifetimes in the NASA GEOS Model during the ORACLES-I campaign

Abstract. In order to improve aerosol representation in the NASA Goddard Earth Observing System (GEOS) model, we evaluated simulations of the transport and properties of aerosols from southern African biomass burning sources that were observed during the first deployment of the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) field campaign in September 2016. An example case study of 24 September was analyzed in detail, during which aircraft-based in situ and remote sensing observations showed the presence of a multi-layered smoke plume structure with significant vertical variation in single scattering albedo (SSA). Our baseline GEOS simulations were not able to represent the observed SSA variation or the observed organic aerosol-to-black-carbon ratio (OA : BC). Analyzing the simulated smoke age suggests that the higher-altitude, less absorbing smoke plume was younger (∼4 d), while the lower-altitude and more absorbing smoke plume was older (∼7 d). We hypothesize a chemical or microphysical loss process exists to explain the change in aerosol absorption as the smoke plume ages, and we apply a simple loss rate to the model hydrophilic biomass burning OA to simulate this process. We also utilized the ORACLES airborne observations to better constrain the simulation of aerosol optical properties, adjusting the assumed particle size, hygroscopic growth, and absorption. Our final GEOS model simulation with additional OA loss and updated optics showed better performance in simulating aerosol optical depth (AOD) and SSA compared to independent ground- and space-based retrievals for the entire month of September 2016, including the Ozone Monitoring Instrument (OMI) Aerosol Index. In terms of radiative implications of our model adjustments, the final GEOS simulation suggested a decreased atmospheric warming of about 10 % (∼2 W m−2) over the southeastern Atlantic region and above the stratocumulus cloud decks compared to the model baseline simulations. These results improve the representation of the smoke age, transport, and optical properties in Earth system models.

Tracking Carbon Dioxide with Lagrangian Transport Simulations: Case Study of Canadian Forest Fires in May 2021

The large amounts of greenhouse gases, such as carbon dioxide, produced by severe forest fires not only seriously affect the ecosystems in the area where the fires occur but also cause a greenhouse effect that has a profound impact on the natural environment in other parts of the world. Numerical simulations of greenhouse gas transport processes are often affected by uncertainties in the location and timing of the emission sources and local meteorological conditions, and it is difficult to obtain accurate and credible predictions by combining remote sensing satellite data with given meteorological forecasts or reanalyses. To study the regional transport processes and impacts of greenhouse gases produced by sudden large-scale forest fires, this study applies the Lagrangian particle dispersion model Massive-Parallel Trajectory Calculations (MPTRAC) to conduct forward simulations of the CO2 transport process of greenhouse gases emitted from forest fires in the central region of Saskatchewan, Canada, during the period of 17 May to 25 May 2021. The simulation results are validated with the Orbiting Carbon Observatory-2 Goddard Earth Observing System (OCO-2 GEOS) Level 3 daily gridded CO2 product over the study area. In order to leverage the high computational costs of the numerical simulations of the model, we implement the forward simulations on the Tianhe-2 supercomputer platform and the JUWELS HPC system, which greatly improves the computational efficiency through parallel computation and makes near-real-time predictions of atmospheric transport processes feasible.

Simulating the Volcanic Sulfate Aerosols From the 1991 Eruption of Cerro Hudson and Their Impact on the 1991 Ozone Hole

AbstractThe Chilean volcano Cerro Hudson erupted between August 8th and 15th, 1991, injecting between 1.7 and 2.9 Tg of SO2 into the upper troposphere and lower stratosphere. We simulate this injection using the Goddard Earth Observing System Earth system model with detailed sulfur chemistry and sectional aerosol microphysics, focusing on the resulting aerosols and their contribution to the 1991 Antarctic Austral Springtime ozone hole. The simulations show a column ozone deficit (12 DU) in the Southern Hemisphere vortex collar region. The majority of this effect is between 10 and 20 km and due to heterogeneous chemistry. The model shows a 26% decrease in ozone from background levels at these altitudes, compared with in‐situ observations of a 50% decrease. Above 20 km, the dynamical response to the eruption also causes lower ozone values, a novel modeling result. This experiment highlights potential interactions between proposed solar radiation management geoengineering aerosols and volcanic eruptions.

Investigation of observed dust trends over the Middle East region in NASA Goddard Earth Observing System (GEOS) model simulations

Abstract. Satellite observations and ground-based measurements have indicated a high variability in the aerosol optical depth (AOD) in the Middle East region in recent decades. In the period that extends from 2003 to 2012, observations show a positive AOD trend of 0.01–0.04 per year or a total increase of 0.1–0.4 per decade. This study aimed to investigate if the observed trend was also captured by the NASA Goddard Earth Observing System (GEOS) model. To this end, we examined changes in the simulated dust emissions and dust AOD during this period. To understand the factors driving the increase in AOD in this region we also examined meteorological and surface parameters important for dust emissions, such as wind fields and soil moisture. Two GEOS model simulations were used in this study: the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) reanalysis (with meteorological and aerosol AOD data assimilated) and MERRA-2 Global Modeling Initiative (GMI) Replay (with meteorology constrained by the MERRA-2 reanalysis but without aerosol assimilation). We did not find notable changes in the modeled 10 m wind speed and soil moisture. However, analysis of Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) data did show an average decrease of 8 % per year in the region encompassing Syria and Iraq, which prompted us to quantify the effects of vegetation on dust emissions and AOD in the Middle East region. This was done by performing a sensitivity experiment in which we enhanced dust emissions in grid cells where the NDVI decreased. The simulation results supported our hypothesis that the loss of vegetation cover and the associated increase in dust emissions over Syria and Iraq can partially explain the increase in AOD downwind. The model simulations indicated dust emissions need to be 10-fold larger in those grid cells in order to reproduce the observed AOD and trend in the model.

Benchmarking GOCART-2G in the Goddard Earth Observing System (GEOS)

Abstract. The Goddard Chemistry Aerosol Radiation and Transport (GOCART) model, which controls the sources, sinks, and chemistry of aerosols within the Goddard Earth Observing System (GEOS), recently underwent a major refactoring and update, including a revision of the emissions datasets and the addition of brown carbon. A 4-year benchmark simulation utilizing the new version of the model code, termed GOCART Second Generation (GOCART-2G) and coupled to the Goddard Earth Observing System (GEOS) model, was evaluated using in situ and spaceborne measurements to develop a baseline and prioritize future development. A comparison of simulated aerosol optical depth between GOCART-2G and MODIS retrievals indicates the model captures the overall spatial pattern and seasonal cycle of aerosol optical depth but overestimates aerosol extinction over dusty regions and underestimates aerosol extinction over Northern Hemisphere boreal forests, requiring further investigation and tuning of emissions. This MODIS-based analysis is corroborated by comparisons to MISR and selected AERONET stations; however, discrepancies between the Aqua and Terra satellites indicate there is a diurnal component to biases in aerosol optical depth over southern Asia and northern Africa. Despite the underestimate of aerosol optical depth in biomass burning regions in GEOS, there is an overestimate in the surface mass of organic carbon in the United States, especially during the summer months. Over Europe, GOCART-2G is unable to match the summertime peak in aerosol optical depth, opposing the observed late fall and early spring peaks in surface mass concentration. A comparison of the vertical profile of attenuated backscatter to observations from CALIPSO indicates the GEOS model is capable of capturing the vertical profile of aerosol; however, the mid-troposphere plumes of dust in the North Atlantic and smoke in the southeastern Atlantic are perhaps too low in altitude. The results presented highlight priorities for future development with GOCART-2G, including improvements for dust, biomass burning aerosols, and anthropogenic aerosols.

Connecting an Offshore Dry Air Stratospheric Intrusion with the Outbreak of Soberanes Fire 2016

Abstract This study investigates the connection between the arrival of dry stratospheric air and the Soberanes Fire (2016). The Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) and Goddard Earth Observing System Forward Processing model (GEOS-FP) are used for back-trajectories and offshore deep stratospheric intrusion in conjunction with the ignition and outbreak of the fire. The back-trajectory analysis indicates that most air reaching the vertical column was critically dry, exhibiting relative humidity values below 10%. As the fire ignited, dry air arrived from due west at heights of 1–3 km about 24 h prior. During the overnight fire growth, dry air arrived from the northwest to north-northwest at heights of 3.5–5.5 km 48–72 h prior. The synoptic and GEOS-FP analyses demonstrate offshore mid-to-low stratospheric intrusion. On 21 July 2016, an enclosed upper-level low approached the California–Oregon border along the northwesterly subtropical jet stream hours before the fire outbreak. The GEOS-FP results of potential vorticity, specific humidity, and ozone along the back-trajectories to the west and northwest of the fire suggest a stratospheric intrusion event into the mid-to-low troposphere at the back-trajectory start points, and vertical velocity indicates sinking motion. The specific humidity analyzed at the arrival time shows the transport of the abnormally dry air to the Soberanes Fire. Results suggest a connection between dry stratospheric air transported to the Soberanes Fire at ignition and overnight accelerated growth, supported by a dark bank in satellite water vapor imagery. The prediction of low-level transport of dry stratospheric air to the coastal communities could help to predict the occurrence of wildfire outbreaks, or periods of accelerated fire growth.

Adapting subseasonal-to-seasonal (S2S) precipitation forecast at watersheds for hydrologic ensemble streamflow forecasting with a machine learning-based post-processing approach

Accurate and reliable precipitation predictions made by dynamical forecast models could provide crucial information for human socioeconomic activities by enabling hydrologic forecasts at the Subseasonal-to-Seasonal (S2S) timescale. To utilize available S2S precipitation predictions for hydrologic forecasts, post-processing techniques have been applied to adapt the raw S2S precipitation to local watersheds. However, conventional statistical-based post-processing techniques are more focused on correcting the forecast bias, but rather limited in improving the predictive skill of available S2S precipitation forecasts. In this study, we combine the Random Forest classifiers (RF) with the Bias Correction and Spatial Disaggregation (BCSD) to adapt the 10-member ensemble precipitation forecast from the NASA Goddard Earth Observing System model version 5 (GEOS5) at 4 watersheds located in the NCEI South climate region of the United States. The adapted S2S precipitation is further applied for streamflow forecast by forcing a classical lumped hydrologic model. The performance of S2S precipitation as well as the corresponding streamflow predictions are benchmarked with the randomly resampled precipitation and the corresponding Ensemble Streamflow Prediction (ESP) framework-generated streamflow predictions. Evaluation statistics of Kling-Gupta Efficiency (KGE), Continuous Ranked Probability Skill Score (CRPSS), Reliability, Resolution, and Sharpness are employed to evaluate the predictive skill of precipitation and streamflow both deterministically and probabilistically. Our results indicate that dynamical S2S precipitation after forecast adaptation leads to consistently higher deterministic skill over ESP at all forecast lead times and across study watersheds. However, at longer forecast lead times beyond 10–15 days, S2S precipitation with a limited ensemble size does not present higher probabilistic skill than ESP. Our results shows that the joint application of RF and BCSD improves the predictive skill of the raw S2S precipitation at study watersheds in contrast to BCSD. Further, the added predictive skill of S2S precipitation brought by RF propagates into streamflow predictions, predominantly at longer forecast lead times exceeding 10 days. Overall, our results highlight the potential success of future work to apply other data-driven approaches to adapt the raw precipitation to local watersheds for more accurate and reliable streamflow forecasts at the S2S timescale.

Analysis of the GEFS-Aerosols annual budget to better understand aerosol predictions simulated in the model

Abstract. In September 2020, a global aerosol forecasting model was implemented as an ensemble member of the National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Prediction (NCEP) Global Ensemble Forecasting System (GEFS) v12.0.1 (hereafter referred to as “GEFS-Aerosols”). In this study, GEFS-Aerosols simulation results from 1 September 2019 to 30 September 2020 were evaluated using an aerosol budget analysis. These results were compared with results from other global models as well as reanalysis data. From this analysis, the global average lifetimes of black carbon (BC), organic carbon (OC), dust, sea salt, and sulfate are 4.06, 4.29, 4.59, 0.34, and 3.3 d, respectively, with the annual average loads of 0.14, 1.29, 4.52, 6.80, and 0.51 Tg. Compared with the National Aeronautics and Space Administration (NASA) Goddard Earth Observing System–Goddard Chemistry Aerosol and Radiation Transport (GEOS4-GOCART) model, the aerosols in GEFS-Aerosols have a relatively short lifetime because of the faster removal processes in GEFS-Aerosols. Meanwhile, in GEFS-Aerosols, aerosol emissions are the determining factor for the mass and composition of aerosols in the atmosphere. The size (bin) distribution of aerosol emissions is as important as its total emissions, especially in simulations of dust and sea salt. Moreover, most importantly, the strong monthly and interannual variations in natural sources of aerosols in GEFS-Aerosols suggest that improving the accuracy of the prognostic concentrations of aerosols is important for applying aerosol feedback to weather and climate predictions.

Paired Satellite and NWP Precipitation for Global Flood Forecasting

Abstract Precipitation data are known to be the key driver of hydrological simulations. Hence, reliable quantitative precipitation estimates and forecasts are vital for accurate hydrological forecasting. Satellite-based precipitation estimates from Integrated Multi-satellitE Retrievals for GPM Early Run (IMERG-E) and forecasted precipitation from NASA’s Goddard Earth Observing System Forward Processing (GEOS-FP) have shown values in global flood nowcasting and forecasting. However, few studies have comprehensively evaluated their hydrological performance let alone explored the potential value of combining them. Therefore, this study undertakes a quasi-global evaluation of their utility in real-time hydrological monitoring and 1–5-day forecasting with the Dominant River Tracing-Routing Integrated with Variable Infiltration Capacity (VIC) Environment (DRIVE) model. The gauge-corrected IMERG Final Run precipitation estimates and corresponding hydrological simulation are used as the references. Results showed that the hit bias is the dominant error source of IMERG-E, while the false precipitation is more noticeable in GEOS-FP. In terms of hydrological performance, the GEOS-FP-driven model (DRIVE-FP) performance is close to the IMERG-E-driven model (DRIVE-E) performance on day 1, indicating that GEOS-FP could nicely fill the gap of nowcasting caused by the IMERG-E time latency. For longer lead-time forecasts, the bias tends to diminish in most regions, likely because the under- or overestimation in IMERG-E is generally offset by the distinct types of misestimation in GEOS-FP. The skillful initial hydrological conditions present outperformed forecasts in most regions, except for tropical areas where the accuracy of GEOS-FP prevails. Overall, this study provides a valuable view of the combined use of IMERG-E and GEOS-FP precipitation in the context of hydrological nowcasts and forecasts.

Public Health Data Applications Using the CDC Tracking Network: Augmenting Environmental Hazard Information With Lower‐Latency NASA Data

AbstractExposure to environmental hazards is an important determinant of health, and the frequency and severity of exposures is expected to be impacted by climate change. Through a partnership with the U.S. National Aeronautics and Space Administration, the U.S. Centers for Disease Control and Prevention's National Environmental Public Health Tracking Network is integrating timely observations and model data of priority environmental hazards into its publicly accessible Data Explorer (https://ephtracking.cdc.gov/DataExplorer/). Newly integrated data sets over the contiguous U.S. (CONUS) include: daily 5‐day forecasts of air quality based on the Goddard Earth Observing System Composition Forecast, daily historical (1980‐present) concentrations of speciated PM2.5 based on the modern era retrospective analysis for research and applications, version 2, and Moderate Resolution Imaging Spectroradiometer (MODIS) daily near real‐time maps of flooding (MCDWD). Data integrated into the CDC Tracking Network are broadly intended to improve community health through action by informing both research and early warning activities, including (a) describing temporal and spatial trends in disease and potential environmental exposures, (b) identifying populations most affected, (c) generating hypotheses about associations between health and environmental exposures, and (d) developing, guiding, and assessing environmental public health policies and interventions aimed at reducing or eliminating health outcomes associated with environmental factors.

The capabilities of the adjoint of GEOS-Chem model to support HEMCO emission inventories and MERRA-2 meteorological data

Abstract. The adjoint of the GEOS-Chem (Goddard Earth Observing System with Chemistry) model has been widely used to constrain the sources of atmospheric compositions. Here, we designed a new framework to facilitate emission inventory updates in the adjoint of the GEOS-Chem model. The major advantage of this new framework is good readability and extensibility, which allows us to support Harmonized Emissions Component (HEMCO) emission inventories conveniently and to easily add more emission inventories following future updates in GEOS-Chem forward simulations. Furthermore, we developed new modules to support MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, version 2) meteorological data, which allows us to perform long-term analyses with consistent meteorological data for the period 1979–present. The performances of the developed capabilities were evaluated with the following steps: (1) diagnostic outputs of carbon monoxide (CO) sources and sinks to ensure the correct reading and use of emission inventories, (2) forward simulations to compare the modeled surface and column CO concentrations among various model versions, (3) backward simulations to compare adjoint gradients of global CO concentrations to CO emissions with finite-difference gradients, and (4) observing system simulation experiments (OSSEs) to evaluate the model performance in 4D variational (4D-Var) assimilations. Finally, an example application of 4D-Var assimilation was presented to constrain anthropogenic CO emissions in 2015 by assimilating Measurement of Pollution in the Troposphere (MOPITT) CO observations. The capabilities developed in this work are important for better applications of the adjoint of the GEOS-Chem model in the future. These capabilities will be submitted to the standard GEOS-Chem adjoint code base for better development of the community of the adjoint of the GEOS-Chem model.

Rapid O3 assimilations – Part 1: Background and local contributions to tropospheric O3 changes in China in 2015–2020

Abstract. A single ozone (O3) tracer mode was developed in this work to build the capability of the Goddard Earth Observing System model with Chemistry (GEOS-Chem) for rapid O3 simulation. The single O3 tracer simulation demonstrates consistency with the GEOS-Chem full chemistry simulation, with dramatic reductions in computational costs of approximately 91 %–94 %. The single O3 tracer simulation was combined with surface and Ozone Monitoring Instrument (OMI) O3 observations to investigate the changes in tropospheric O3 over eastern China in 2015–2020. The assimilated O3 concentrations demonstrate good agreement with O3 observations because surface O3 concentrations are 43.2, 41.8, and 42.1 ppb and tropospheric O3 columns are 37.1, 37.9, and 38.0 DU in the simulations, assimilations, and observations, respectively. The assimilations indicate rapid rises in surface O3 concentrations by 1.60 (spring), 1.16 (summer), 1.47 (autumn), and 0.80 ppb yr−1 (winter) over eastern China in 2015–2020, and the increasing trends are underestimated by the a priori simulations. More attention is suggested to the rapid increases in the O3 pollution in spring and autumn. We find stronger rises in tropospheric O3 columns over highly polluted areas due to larger local contributions, for example, 0.12 DU yr−1 (North China Plain) in contrast to −0.29 (Sichuan Basin) and −0.25 DU yr−1 (southern China). Furthermore, our analysis demonstrated noticeable contributions of the interannual variability in background O3 to the trends in surface O3 (particularly in the summer) and tropospheric O3 columns over eastern China in 2015–2020. This work highlights the importance of rapid simulations and assimilations to extend and interpret atmospheric O3 observations.

Response of the ozone-related health burden in Europe to changes in local anthropogenic emissions of ozone precursors

Exposure to ozone (O3) is associated with many human health problems, resulting in tens of thousands of premature deaths annually in Europe. This study quantifies the impact of changes in anthropogenic emissions of O3 precursors on premature deaths from long-term O3 exposure in Europe and the impact of emissions changes during 2005–2015 using the nested-grid chemical transport model Goddard Earth Observing System (GEOS)-Chem and its adjoint. In 2015, it is estimated that a 20% decrease in total anthropogenic emissions in our modeled European domain could prevent 1576 (467–3252) premature deaths from respiratory disease (30 years of age), 70% of which is owing to the decrease in nitrogen oxides (NO x ) emissions. Underlying this aggregate effect is substantial spatial variation. In most of Europe, O3 formation is NO x -limited so that NO x emission reductions help to decrease premature deaths. Yet where O3 formation is NO x -saturated (as in parts of the United Kingdom, Benelux and Germany) emission reductions cause more premature deaths through increased ozone exposure. Despite the overall decreases in anthropogenic emissions, the marginal benefit, expressed as the avoided premature deaths per 1 kg km−2 yr−1 reduction in NO x emissions, is found to generally increase during 2005–2015, with a mean value more than doubling over Europe. This highlights the general trend that O3 formation becomes less sensitive to volatile organic compound emissions and more limited by NO x emissions. An important policy implication of increasing marginal benefits is that more costly regulations of NO x emissions are economically justified even as total anthropogenic emission are declining. NO x contributions from road transport, industry, energy, and residential sectors are most affected by the change in the O3 production regime. Consequently, European regulations of NO x emissions targeted at those sectors will yield the highest health benefits per unit NO x emission of all sources.

Direct assimilation of AVHRR satellite radiance measurements in a reanalysis system

AbstractThe Advanced Very‐High‐Resolution Radiometer (AVHRR) is a broad‐band, five‐channel scanner sensing in the visible, near‐infrared and thermal infrared portions of the electromagnetic spectrum. AVHRR instruments onboard polar‐orbiting satellites have data records spanning 30 years. The radiances from the infrared channels of AVHRR have been directly assimilated over the ice‐free ocean in the NASA Goddard Earth Observing System (GEOS) since 2017 to constrain skin sea surface temperature (SST). The GEOS system already uses an advanced bulk SST and a skin SST model which incorporates diurnal warming and cool‐skin temperature on the surface of the ocean. The positive contribution of this effort to the Numerical Weather Prediction (NWP) makes it desirable to extend the skin SST data assimilation procedure to reanalysis applications. The AVHRR data from platforms that were previously unexplored for NWP applications are assimilated at the Global Modeling and Assimilation Office (GMAO) at NASA, for an upcoming reanalysis project. This study uses results from reanalysis experiments to assess the impact of assimilating AVHRR radiances over open ocean on the atmospheric state variables, focusing on the analyzed skin temperatures that include the active skin SST model. It is demonstrated that including the skin temperature variations in atmospheric analysis is generally beneficial. The addition of radiance measurements from AVHRR provides further improvements to the overall reanalysis system performance. by helping to correct the diurnal heating and reducing errors in the background departure of hyperspectral radiance observations, which are an essential component of atmospheric reanalyses.

Surface Radiative Forcing as a Climate-Change Indicator in North India due to the Combined Effects of Dust and Biomass Burning

This study estimates the radiative forcing by biomass burning and dust aerosols over the Indian subcontinent, with emphasis on the Indo-Gangetic Plains (IGP) during the period from January 2021 to April 2021, based on multiple satellite and reanalysis datasets. In this respect, we used retrievals from the Moderate Resolution Spectroradiometer (MODIS) and the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) system, as well as reanalysis data from the Goddard Earth Observing System, version 5 (GEOS-5), the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2), the Copernicus Atmosphere Monitoring Service (CAMS), and ERA-Interim. According to the MERRA-2 and the CAMS, the highest black carbon (BC) concentrations in January 2021 were 7–8 µg m−3, which were significantly lower than measurements performed in main cities along the IGP, such as Patiala, Delhi, and Kanpur. The meteorological data analysis accompanied by the CALIPSO lidar measurements showed that the vertical distribution of total attenuated backscatter (TAB) could reach altitudes of up to ~4–5 km and could be transported over the central Himalayan region. The spatial-averaged daily aerosol radiative forcing (ARF) values over the Indian subcontinent from January 2021 to April 2021 were found to range from −51.40 to −6.08 W m−2 (mean of −22.02 ± 9.19 W m−2), while on a monthly basis, the ARF values varied widely, from −146.24 to −1.63 W m−2 (mean of −45.56 ± 22.85 W m−2) over different parts of the study region. Furthermore, the spatial-averaged daily BC radiative forcing ranged from −2.23 to −0.35 (−1.01 ± 0.40 W m−2), while it varied from −15.29 to −0.31 W m−2 (−2.46 ± 2.32 W m−2) over different regions of southern Asia, indicating a rather small contribution to the total aerosol radiative effect and a large presence of highly scattering aerosols. Our findings highlight the importance of growing biomass burning, in light of recent climate change and the rapid decline in air quality over North India and the Indian Ocean.

Tropospheric ozone data assimilation in the NASA GEOS Composition Forecast modeling system (GEOS-CF v2.0) using satellite data for ozone vertical profiles (MLS), total ozone columns (OMI), and thermal infrared radiances (AIRS, IASI)

The NASA Goddard Earth Observing System Composition Forecast system (GEOS-CF) provides global near-real-time analyses and forecasts of atmospheric composition. The current version of GEOS-CF builds on the GEOS general circulation model with Forward Processing assimilation of meteorological data (GEOS-FP) and includes detailed GEOS-Chem tropospheric and stratospheric chemistry. Here we add 3D variational data assimilation in GEOS-CF to assimilate satellite observations of ozone including MLS vertical profiles, OMI total columns, and AIRS and IASI hyperspectral 9.6 μm radiances. We focus our evaluations on the troposphere. We find that the detailed tropospheric chemistry in GEOS-CF significantly improves the simulated background ozone fields relative to previous versions of the GEOS model, allowing for specification of smaller background errors in assimilation and resulting in smaller assimilation increments to correct the simulated ozone. Assimilation increments are largest in the upper troposphere and are consistent between satellite data sets. The OMI and MLS ozone data generally provide more information than the AIRS and IASI radiances except at high latitudes where the radiances provide more information. Comparisons to independent ozonesonde and aircraft (ATom-4) observations for 2018 show significant GEOS-CF improvement from the assimilation, particularly in the extratropical upper troposphere.

Interactive Stratospheric Aerosol Microphysics‐Chemistry Simulations of the 1991 Pinatubo Volcanic Aerosols With Newly Coupled Sectional Aerosol and Stratosphere‐Troposphere Chemistry Modules in the NASA GEOS Chemistry‐Climate Model (CCM)

AbstractWe have coupled the GEOS‐Chem tropospheric‐stratospheric chemistry mechanism and the Community Aerosol and Radiation Model for Atmospheres (CARMA), a sectional aerosol microphysics module, within the NASA Goddard Earth Observing System Chemistry‐Climate Model (GEOS CCM) in order to simulate the interactions between stratospheric chemistry and aerosol microphysics. We use observations of the 1991 Mount Pinatubo volcanic cloud to evaluate this new version of the GEOS CCM. The GEOS‐Chem chemistry module is used to simulate the oxidation of sulfur dioxide (SO2) more realistically than assuming hydroxyl radical (OH) fields are constant, as OH concentrations in the plume decrease dramatically in the weeks following the eruption. CARMA simulates sulfate aerosols with dynamic microphysical and optical properties. The CARMA‐calculated aerosol surface area is coupled to the chemistry module from GEOS‐Chem for the calculation of heterogeneous chemistry. We use a set of observational and theoretical constraints for Pinatubo to evaluate the performance of this new version of the GEOS CCM. These simulations are specifically compared with satellite and in‐situ observations and provide insights into the connections between the gas‐phase chemistry and the aerosol microphysics of the early plume and how they impact the climatic and chemical changes following a large volcanic eruption. A second, smaller eruption is also included in these simulations, the 15 August 1991, eruption of Cerro Hudson in Chile, which we find essential in explaining the aerosol optical depth in the Southern Hemisphere in 1991.

Application of Satellite‐Based Detections of Arctic Bromine Explosion Events Within GEOS‐Chem

AbstractDuring polar spring, periods of elevated tropospheric bromine drive near complete removal of surface ozone. These events impact the tropospheric oxidative capacity and are an area of active research with multiple approaches for representing the underlying processes in global models. We present a method for parameterizing emissions of molecular bromine (Br2) over the Arctic using satellite retrievals of bromine monoxide (BrO) from the Ozone Monitoring Instrument (OMI). OMI retrieves column BrO with daily near global coverage, and we use the GEOS‐Chem chemical mechanism, run online within the Goddard Earth Observing System Earth System Model to identify hotspots of BrO likely associated with polar processes. To account for uncertainties in modeling background BrO, hotspots are only identified where the difference between OMI and modeled columns exceeds a statistical threshold. The resulting hotspot columns are a lower‐limit for the portion of OMI BrO attributable to bromine explosion events. While these hotspots are correlated with BrO measured in the lower troposphere over the Arctic Ocean, a case study of missing detections of near‐surface BrO is identified. Daily flux of Br2 is estimated from hotspot columns of BrO using internal model parameters. When the emissions are applied, BrO hotspots are modeled with a 5% low bias. The sensitivity of the resulting ozone simulations to the treatment of background uncertainties in the BrO column is demonstrated. While periods of isolated, large (>50%) decreases in surface ozone are modeled, this technique does not simulate the low ozone observed at coastal stations and consistently underestimates ozone loss during March.

Impacts of Marine Surface Pressure Observations from a Spaceborne Differential Absorption Radar Investigated with an Observing System Simulation Experiment

Abstract A new instrument has been proposed for measuring surface air pressure over the marine surface with a combined active/passive scanning multichannel differential absorption radar to provide an estimate of the total atmospheric column oxygen content. A demonstrator instrument, the Microwave Barometric Radar and Sounder (MBARS), has been funded by the National Aeronautics and Space Administration for airborne test missions. Here, a proof-of-concept study to evaluate the potential impact of spaceborne surface pressure data on numerical weather prediction is performed using the Goddard Modeling and Assimilation Office global observing system simulation experiment (OSSE) framework. This OSSE framework employs the Goddard Earth Observing System model and the hybrid 4D ensemble variational Gridpoint Statistical Interpolation data assimilation system. Multiple flight and scanning configurations of potential spaceborne orbits are examined. Swath width and observation spacing for the surface pressure data are varied to explore a range of sampling strategies. For wider swaths, the addition of surface pressures reduces the root-mean-square surface pressure analysis error by as much as 20% over some ocean regions. The forecast sensitivity observation impact tool estimates impacts on the Pacific Ocean basin boundary layer 24-h forecast temperatures for spaceborne surface pressures that are on par with rawinsondes and aircraft and estimates greater impacts than the current network of ships and buoys. The largest forecast impacts are found in the Southern Hemisphere extratropics.