Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASSIMILATION”

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 1,477 records · Page 82

Using Satellite-derived Ice Concentration to Represent Antarctic Coastal Polynyas in Ocean Climate Models

The focus of this paper is on the representation of Antarctic coastal polynyas in global ice-ocean general circulation models (OGCMs), in particular their local, regional, and high-frequency behavior. This is verified with the aid of daily ice concentration derived from satellite passive microwave data using the NASATeam 2 (NT2) and the bootstrap (BS) algorithms. Large systematic regional and temporal discrepancies arise, some of which are related to the type of convection parameterization used in the model. An attempt is made to improve the fresh-water flux associated with melting and freezing in Antarctic coastal polynyas by ingesting (assimilating) satellite ice concentration where it comes to determining the thermodynamics of the open-water fraction of a model grid cell. Since the NT2 coastal open-water fraction (polynyas) tends to be less extensive than the simulated one in the decisive season and region, assimilating NT2 coastal ice concentration yields overall reduced net freezing rates, smaller formation rates of Antarctic Bottom Water, and a stronger southward flow of North Atlantic Deep Water across 30 S. Enhanced net freezing rates occur regionally when NT2 coastal ice concentration is assimilated, concomitant with a more realistic ice thickness distribution and accumulation of High-Salinity Shelf Water. Assimilating BS rather than NT2 coastal ice concentration, the differences to the non-assimilated simulation are generally smaller and of opposite sign. This suggests that the model reproduces coastal ice concentration in closer agreement with the BS data than with the NT2 data, while more realistic features emerge when NT2 data are assimilated.

Stoessel, Achim↗

ENVISAT Land Surface Processes

This is a progress report of the 2nd phase of the project ENVISAT- Land Surface Processes, which has a 3-year scope. In this project, preparative research is carried out aiming at the retrieval of land surface characteristics from the ENVISAT sensors MERIS and AATSR, for assimilation into a system for Numerical Weather Prediction (NWP). Where in the 1st phase a number of first shot experiments were carried out (aiming at gaining experience with the retrievals and data assimilation procedures), the current 2nd phase has put more emphasis on the assessment and improvement of the quality of the retrieved products. The forthcoming phase will be devoted mainly to the data assimilation experiments and the assessment of the added value of the future ENVISAT products for NWP forecast skill. Referring to the retrieval of albedo, leaf area index and atmospheric corrections, preliminary radiative transfer calculations have been carried out that should enable the retrieval of these parameters once AATSR and MERIS data become available. However, much of this work is still to be carried out. An essential part of work in this area is the design and implementation of software that enables an efficient use of MODTRAN(sub 4) radiative transfer code, and during the current project phase familiarization with these new components has been achieved. Significant progress has been made with the retrieval of component temperatures from directional ATSR-images, and the calculation of surface turbulent heat fluxes from these data. The impact of vegetation cover on the retrieved component temperatures appears manageable, and preliminary comparison of foliage temperature to air temperatures were encouraging. The calculation of surface fluxes using the SEBI concept,which includes a detailed model of the surface roughness ratio, appeared to give results that were in reasonable agreement with local measurements with scintillometer devices. The specification of the atmospheric boundary conditions appears a crucial component, and the use of first guess estimates from the RACMO models partially explains the success. Earlier data assimilation experiments with directional surface temperatures have been analysed a bit further and were also compared to results obtained from directly modeling the surface roughness ratio. Results between these calculations and the data assimilation results appeared well comparable, but a full test in which the surface roughness model is allowed to play a free role during the data assimilation process has yet to be carried out. A considerable number of tasks that have yet to be carried out during Phase 3 has been formulated.

vandenHurk, B. J. J. M.↗

Effect of Cross-Correlation on Geomagnetic Forecast Accuracies

Surface geomagnetic observation can determine up to degree L = 14 time-varying spherical harmonic coefficients of the poloidal magnetic field. Assimilation of these coefficients to numerical dynamo simulation could help us understand better the dynamical processes in the Earth's outer core, and to provide more accurate forecast of geomagnetic secular variations (SV). In our previous assimilation studies, only the poloidal magnetic field in the core is corrected by the observations in the analysis. Unobservable core state variables (the toroidal magnetic field and the core velocity field) are corrected via the dynamical equations of the geodynamo. Our assimilation experiments show that the assimilated core state converges near the CMB, implying that the dynamo state is strongly constrained by surface geomagnetic observations, and is pulled closer to the truth by the data. We are now carrying out an ensemble of assimilation runs with 1000 years of geomagnetic and archeo/paleo magnetic record. In these runs the cross correlation between the toroidal and the poloidal magnetic fields is incorporated into the analysis. This correlation is derived from the physical boundary conditions of the toroidal field at the core-mantle boundary (CMB). The assimilation results are then compared with those of the ensemble runs without the cross-correlation, aiming at understanding two fundamental issues: the effect of the crosscorrelation on (1) the convergence of the core state, and (2) the SV prediction accuracies. The constrained dynamo solutions will provide valuable insights on interpreting the observed SV, e.g. the near-equator magnetic flux patches, the core-mantle interactions, and possibly other geodynamic observables.

Kuang, Weijia↗

Towards the Handing of Cloud-Affected Infrared Radiances in the GSI

In the gridpoint statistical interpolation (GSI) data assimilation algorithm, only thermal infrared measurements determined to be uncontaminated by clouds are assimilated. Using this approach, typically only 19-29% of footprints are deemed to have no cloud affects through the measured spectra. This study will discuss the efforts underway at the Global Modeling and Assimilation Office (GMAO) at NASA Goddard Space Flight Center, in conjunction with the Joint Center for Satellite Data Assimilation (JCSDA), to actively assimilate these morecomplicated observations by using a graybody assumption. In the GSI, cloud top pressure and effective cloud amount are retrieved concurrently using a minimum residual method. This study will address the limitations and advantages of the technique and the modifications underway to the assimilation system to incorporate those two parameters into the radiative transfer forward operators and TL/AD calculations. Furthermore, it will explain the efforts underway to incorporate these parameters into the control vector so that they can be altered variationally as part of the minimization.

McCarty, William↗

File Specification for GEOS-5 FP (Forward Processing)

The GEOS-5 FP Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The assimilation is performed at a horizontal resolution of 0.3125-degree longitude by 0.25- degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational FP, using the most recent version of the GEOS-5 assimilation system. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary. Documentation about the current access methods for products described in this document can be found on the GMAO products page: http://gmao.gsfc.nasa.gov/products/.

GSI↗

Soil Moisture Active Passive Mission L4_SM Data Product Assessment (Version 2 Validated Release)

During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.

SMAP L4_SM↗

Aerosol EnKF at GMAO

In the GEOS near real-time system, as well as in MERRA-2 which is the latest reanalysis produced at NASAs Global Modeling and Assimilation Office(GMAO), the assimilation of aerosol observations is performed by means of a so-called analysis splitting method. In line with the transition of the GEOS meteorological data assimilation system to a hybrid Ensemble-Variational formulation, we are updating the aerosol component of our assimilation system to an ensemble square root filter(EnSRF; Whitaker and Hamill (2002)) type of scheme.We present a summary of our preliminary results of the assimilation of column integrated aerosol observations (Aerosol Optical Depth; AOD) using an Ensemble Square Root Filters (EnSRF) scheme and the ensemble members produced routinely by the meteorological assimilation.

MERRA-2 EnSRF↗

The Impact of Satellite Sea Surface Salinity for Prediction of the Coupled Indo-Pacific System

We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO 3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10 degrees South latitude-5 degrees North latitude. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO 3.4 region. In addition, salting near the ITCZ (Intertropical Convergence Zone) leads to a deepening of the mixed layer and thermocline near 8 degrees North latitude. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.

El Nino↗

The Impact of Satellite Sea Surface Salinity for Prediction of the Coupled Indo-Pacific System

We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO3.4 region. In addition, salting near the ITCZ leads to a deepening of the mixed layer and thermocline near 8N. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.

SST↗

The MERRA-2 Aerosol Reanalysis, 1980 Onward. Part II: Evaluation and Case Studies

The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) is NASA’s latest reanalysis for the satellite era (1980-present) using the Goddard Earth Observing System version 5 (GEOS-5) Earth system model. MERRA-2 provides several improvements over its predecessor (MERRA), including the inclusion of interactive aerosols for the entire period. In addition to assimilating bias-corrected Aerosol Optical Depth (AOD) from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Terra and Aqua satellites, it also includes the assimilation of bias-corrected AOD from Advanced Very High Resolution Spectroradiometer (AVHRR) instruments, Multi-angle Imaging SpectroRadiometer (MISR) AOD over bright surfaces, and ground-based Aerosol Robotic Network (AERONET) AOD. This paper is the second of a pair that summarizes our efforts to assess the quality of the MERRA-2 aerosol assimilation. In this study, we first follow previous work performed with version 1 of the MERRA Aerosol Reanalysis (MERRAero) using independent observations. The evaluation of MERRA-2 Absorption Aerosol Optical Depth (AAOD) and ultra-violet Aerosol Index (UV-AI) against the Ozone Monitoring Instrument (OMI) observations show good agreement, particularly over dusty regions where our previous efforts improved model aerosol optical properties. Next, we find that aerosol assimilation system improves the aerosol vertical structure when compared to estimates from the same version of the model without AOD assimilation. A similar conclusion is found for MERRA-2 aerosol surface fine particulate matter (PM (sub 2.5)). However, deficiencies in the forward model such as missing emissions noted during the MERRAero study still explain the MERRA-2 PM(sub 2.5) bias relative to observations over the United States. Finally, to illustrate successes and weaknesses of the AOD assimilation, we focus on the performance of the MERRA-2 aerosol system during several major aerosol events: the Mount Pinatubo eruption in 1991, a Saharan dust event and transportation over the Atlantic Ocean during April 2010, the Rim Fire of summer 2013 in California, and an extreme pollution event over China in January 2013. We conclude with a summary that points to best practices for utilizing the MERRA-2 aerosol reanalysis in future studies.

V Buchard↗

Progress Toward Implementation of Coupled Ocean/Atmosphere in Unified Forward Operator (UFO)

The (Sea-Ice Ocean and Coupled Assimilation) SOCA project focuses on meshing atmospheric and marine data assimilation efforts to create innovative data assimilation systems. One of the main goals ofSOCA is to make use of surface-sensitive radiances to constrain sea-ice and upper ocean fields (e.g., salinity, temperature, sea-ice fraction, sea-ice temperature, etc.). the focus of this research is to build first elements toward an ocean/sea-ice/atmosphere coupled data assimilation capability, with a focus on supporting and developing the assimilation of observations sensitive to multiple sub-domains. We Setup the existing Joint Effort for Data assimilation Integration (JEDI) coupled UFO (CRTM related) for surface sensitive microwave radiances (SST and SSS) focusing on Global Precipitation Measurement (GPM) Imager (GMI) and Soil Moisture Active Passive (SMAP).

Ebrahimi, Hamideh↗

Impact of Gauge-Based Precipitation Corrections on the Skill of SMAP Level-4 Soil Moisture Estimates

The NASA Soil Moisture Active Passive (SMAP) mission provides observations of L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations at a resolution of ~40 km globally every 2-3 days. These observations are routinely assimilated into the NASA Catchment land surface model to generate the Level-4 Soil Moisture (L4_SM) product, which provides global estimates of surface and root-zone soil moisture, soil temperature, and surface fluxes (among others) at 9-km, 3-hourly resolution with ~2.5-day latency. The Catchment land surface model in the L4_SM algorithm is driven with 0.25°, hourly surface meteorological forcing data from the NASA Goddard Earth Observing System (GEOS) "forward-processing" product. Outside of Africa and the high latitudes, the GEOS precipitation forcing is corrected using the Climate Prediction Center Unified (CPCU) gauge-based, 0.5°, daily precipitation product.Soil moisture estimates from the L4_SM product were previously shown to improve over land model-only estimates that do not benefit from the assimilation of Tb observations, thereby demonstrating the value of assimilating SMAP observations for soil moisture estimation. In this presentation, we further isolate the contribution of the gauge-based precipitation corrections to the skill of the L4_SM soil moisture estimates. Specifically, we compare the skill of the L4_SM soil moisture to that of separate model-only and assimilation estimates obtained without the benefit of the gauge-based precipitation corrections.Preliminary results suggest that the soil moisture skill added by the CPCU-based precipitation corrections primarily depends on the quality of the CPCU precipitation product and is greatest in regions where the CPCU gauge network is dense and reliable. Conversely, in regions where the CPCU product is known to be of poor quality, for example in central Australia, the assimilation of SMAP Tb observations provides the most benefit. The presentation will provide an in-depth evaluation of the soil moisture skill of the model-only and assimilation estimates vs. independent in situ and satellite measurements.

Reichle, Rolf↗

Evaluation of 18 Satellite- and Model-Based Soil Moisture Products Using in Situ Measurements From 826 Sensors

Information about the spatiotemporal variability of soil moisture is critical for many purposes, including monitoring of hydrologic extremes, irrigation scheduling, and prediction of agricultural yields. We evaluated the temporal dynamics of 18 state-of-the-art (quasi-)global near-surface soil moisture products, including six based on satellite retrievals, six based on models without satellite data assimilation (referred to hereafter as “open-loop” models), and six based on models that assimilate satellite soil moisture or brightness temperature data. Seven of the products are introduced for the first time in this study: one multi-sensor merged satellite product called MeMo (Merged soil Moisture) and six estimates from the HBV (Hydrologiska Byråns Vattenbalansavdelning) model with three precipitation inputs (ERA5, IMERG, and MSWEP) with and without assimilation of SMAPL3E satellite retrievals, respectively. As reference, we used in situ soil moisture measurements between 2015 and 2019 at 5 cm depth from 826 sensors, located primarily in the USA and Europe. The 3-hourly Pearson correlation (R) was chosen as the primary performance metric. We found that application of the Soil Wetness Index (SWI) smoothing filter resulted in improved performance for all satellite products. The best-to-worst performance ranking of the four single-sensor satellite products was SMAPL3E SWI , SMOS SWI , AMSR2 SWI , and ASCAT SWI , with the L-band-based SMAPL3E SWI (median R of 0.72) outperforming the others at 50 % of the sites. Among the two multi-sensor satellite products (MeMo and ESA-CCI SWI ), MeMo performed better on average (median R of 0.72 versus 0.67), probably due to the inclusion of SMAPL3E SWI . The best-to-worst performance ranking of the six open-loop models was HBV-MSWEP, HBV-ERA5, ERA5-Land, HBV-IMERG, VIC-PGF, and GLDAS-Noah. This ranking largely reflects the quality of the precipitation forcing. HBV-MSWEP (median R of 0.78) performed best not just among the open-loop models but among all products. The calibration of HBV improved the median R by +0.12 on average compared to random parameters, highlighting the importance of model calibration. The best-to-worst performance ranking of the six models with satellite data assimilation was HBV-MSWEP+SMAPL3E, HBV-ERA5+SMAPL3E, GLEAM, SMAPL4, HBV-IMERG+SMAPL3E, and ERA5. The assimilation of SMAPL3E retrievals into HBV-IMERG improved the median R by +0.06, suggesting that data assimilation yields significant benefits at the global scale.

Hylke E. Beck↗

Assessment of Retrieved GMI Emissivity Over Land, Snow, and Sea Ice in the GEOS System

Directly assimilating microwave radiances over land, snow and sea ice remains a significant challenge for data assimilation systems. These data assimilation systems are critical to the success of global numerical weather prediction systems including the Global Earth Observing System-Atmospheric Data Assimilation System (GEOS-ADAS). Extending more surface sensitive microwave channels over land, snow and ice could provide a needed source of data for Numerical Weather Prediction particularly in the Planetary Boundary Layer (PBL). Unfortunately, the accuracy of emissivity models currently available within the GEOS-ADAS along with other data assimilation systems are insufficient to simulate and assimilate radiances. Recently, Munchak et al. (2020) published a 5-year climatological database for retrieved microwave emissivity from the GPM Microwave Imager (GMI) aboard the Global Precipitation Mission (GPM). In this work the database is utilized by modifying the GEOS-ADAS to use this emissivity database in place of the default emissivity value available in the Community Radiative Transfer Model (CRTM), which is the fast radiative transfer model used by the GEOS-ADAS. As a first step, the GEOS-ADAS is run in a so-called “stand-alone” mode to simulate radiances from GMI using the default CRTM emissivity, and replacing the default CRTM emissivity models with values from Munchak et al, 2020. The simulated observations using Munchak et al., 2020 agree more closely with observations from GMI. These results are presented along with a discussion of the implication for GMI observations within the GEOS-ADAS.

GPM↗

M2-SCREAM: A Stratospheric Composition Reanalysis of Aura MLS

Stratospheric composition and transport patterns are changing because of increasing greenhouse gas concentrations and declining levels of ozone depleting substances. Superimposed on forced trends are: dynamical variability at a range of time scales and unforeseen perturbations such as those resulting from injections of volcanic material and smoke from PyroCb events. Chemical reanalyses can be powerful tools for studying the composition of the stratosphere and disentangling forced responses from internal variability. This presentation introduces a new chemical reanalysis of stratospheric constituents developed and produced at NASA’s Global Modeling and Assimilation Office (GMAO) using the recently developed GMAO Constituent Data Assimilation System. Called the MERRA-2 Stratospheric Composition Reanalysis with Aura MLS (M2-SCREAM), this reanalysis consists of assimilated global three-dimensional fields of stratospheric ozone, water vapor, hydrogen chloride (HCl), nitric acid (HNO3) and nitrous oxide (N2O) mixing ratios and covers the period since the beginning of MLS observations in September 2004 to the present. M2-SCREAM assimilates version 4.2 MLS profiles of the five constituents alongside total ozone column from the Ozone Monitoring Instrument. The dynamics and tropospheric water vapor are constrained by assimilated meteorological fields from MERRA-2. The assimilated constituent fields and meteorology are provided to users at a 50-km horizontal resolution and a three-hourly frequency. Monthly analysis uncertainties and flags are also provided to the users. M2-SCREAM is an accurate and dynamically consistent high-resolution data record of the five constituents, all of which are of primary importance to stratospheric chemistry and transport studies. We will present a description of M2-SCREAM and selected results of a process-based evaluation of this product using independent data. As an illustration, we will discuss the perturbation to stratospheric composition from the eruption of Hunga Tonga–Hunga Ha'apai in January 2022 as seen in the reanalysis. We will also propose potential scientific applications of M2-SCREAM and outline plans for an upcoming comprehensive composition reanalysis that is being developed at NASA GMAO.

MERRA-2↗

Importance of Radiative Transfer Models in Atmospheric Remote Sensing

Radiative transfer models (RTMs) play a significant role in the development of satellite instruments for remote sensing applications. These models simulate electromagnetic radiation's propagation through the atmosphere, providing valuable insights into atmosphere-radiation interactions. RTMs facilitate the optimization of satellite instrument designs, ensuring their ability to measure targeted atmospheric and surface properties accurately. Moreover, they aid in simulating instrument’s measurements under various atmospheric conditions, enabling calibration and validation processes to enhance data quality and reliability. RTMs are extensively used in the Observing System Simulation Experiments (OSSE), to generate synthetic observations. By incorporating RTMs into OSSE, we can assess the potential impact of future satellite missions, sensor configurations, and data assimilation techniques. This approach allows for the optimization of satellite instruments and constellations and the evaluation of their influence on weather forecasting, climate monitoring, and other Earth science applications. Another crucial application area of RT models is data assimilation, where they play a fundamental role in combining satellite observations with numerical models to improve atmospheric and environmental predictions. RTMs provide the link between observed radiances and atmospheric parameters, enhancing the accuracy of numerical models and generating more reliable forecasts for weather events, air quality assessments, and climate projections. Moreover, adapting RT models to capture the intricate radiation interactions within the Planetary Boundary Layer will significantly contribute to improving weather forecasting and climate change projections. Current community radiative transfer (RT) models are primarily developed and optimized for operational data assimilation of satellite observations. These models excel at assimilating satellite data into numerical weather prediction models to improve forecast accuracy. However, their focus on data assimilation limits their suitability for other important applications, such as satellite instrument development, OSSE, and Planetary Boundary Layer (PBL) studies. Moreover, for PBL studies, RT models need to be adapted to capture the intricate radiation interactions within this crucial atmospheric layer. Developing RT models that can represent the PBL's unique characteristics, such as surface interactions, will contribute significantly to understanding and predicting weather phenomena, air quality, and climate dynamics. This abstract provides a comprehensive overview of the current status of RT models and highlights their limitations concerning satellite instrument development, OSSE, and PBL studies. Addressing these shortcomings requires concerted efforts to enhance RT models' capabilities and expand their applications beyond data assimilation. By investing in research and development to improve these

Isaac Moradi↗

Detection of CH 4 Hotspots from NASA’s GEOS Composition Analysis System

Recent research indicates that 8 to 12% of the global oil and gas production methane emissions could be attributed to ultra-emitters, which result in high concentration ‘hotspots’ near point sources. Identifying these emissions in near real time provides useful information to the policy makers and private industry, who are working to reduce their impact. To meet this need, scientists are increasingly analyzing satellite data from the TROPOspheric Monitoring Instrument (TROPOMI) instrument aboard ESA’s Sentinel 5-Precursor mission. While direct analysis of TROPOMI level 2 swath data has been successful in identifying some large emission events, identifying hotpots is challenging because of the imaging noise due to a variety of artifacts and limits in daily coverage. Here we explore possible methodologies to detect methane hotspots using a new, gap-filled, and temporally continuous methane product from NASA’s Goddard Earth Observing System (GEOS) Constituent Data Assimilation System (CoDAS), which assimilates column averaged methane mole fractions from the TROPOMI with capabilities to assimilate other remote sensing measurements. The CoDAS has been expanded from a heritage of stratospheric composition and carbon dioxide assimilation allowing for the support of regional modeling, validation with non-coincident operations, and merging variety of datasets. The current work mainly explores Observing System Simulation Experiments with methane GEOS simulations without assimilation to prepare the groundwork for further experiments with assimilated TROPOMI. First, known hotspots based on the known inventory are identified to demonstrate the capability of the system to point out emission hotspots. In the next step, a variety of machine learning techniques such as Self-Organizing Maps and Deep Learning are explored to automate detection of the plumes. Finally, a few approaches to quantify emissions from the identified hotspots are presented and are evaluated against the inventory. The effort is directed toward a future evaluation of the CoDAS based methane monitoring system’s ability to successfully detect and quantify hotspots.

Nikolay Balashov↗

The Land Surface in Current and Planned MERRA Reanalysis Products

Current global atmospheric reanalysis products such as the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5), the NASA Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2), and the Japanese Reanalysis for Three Quarters of a Century (JRA-3Q) provide estimates of land surface states and fluxes, including soil moisture, soil temperature, snow mass, latent and sensible heat fluxes, and runoff, that are widely used in research and applications. These land surface estimates are based on land surface process models and, depending on the reanalysis product, on precipitation observations or the assimilation of land surface observations of soil moisture, soil temperature, snow conditions, and screen-level air temperature and humidity from satellite observations and in situ measurements. In this presentation, we review the land surface modeling and data assimilation components of the suite of current and planned MERRA reanalysis products. In addition to MERRA-2, we will discuss the latest NASA reanalysis, MERRA for the 21st century (M21C), which is currently under production, as well as the development and planning of the next version of the MERRA reanalysis, tentatively labeled MERRA-3. In MERRA-2, observations-based precipitation data products are used to correct the precipitation falling on the land surface. Outside of the high-latitudes and Africa, the daily, 0.5-degree, gauge-based Climate Prediction Center (CPC) Unified (CPCU) product is used. In Africa, the pentad, 2.5-degree, satellite- and gauge-based CPC Merged Analysis of Precipitation (CMAP) product is used. Poleward of 62.5 degrees latitude, the land surface sees the precipitation generated by the atmospheric model in the cycling data assimilation system. This configuration provides improved soil moisture estimates compared to those of the original (version 1) MERRA estimates, which did not benefit from the use of precipitation observations. Moreover, the use of precipitation observations facilitates a seamless spin-up of the land surface initial conditions across the MERRA-2 production streams. The use of a gauge-only precipitation product in MERRA-2 across much of the globe, however, adversely impacts the quality of the MERRA-2 land surface estimates in regions with poor gauge coverage, including most of South America and Australia. Therefore, the forthcoming M21C reanalysis uses satellite- and gauge-based precipitation from the Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement Mission (IMERG). This change results in significant improvements in the quality of the M21C soil moisture estimates in the Southern Hemisphere compared to those from MERRA-2. Planning for MERRA-3 focuses on the assimilation of soil moisture observations from the Soil Moisture Active Passive (SMAP) mission and the Advanced Scatterometer (ASCAT), along with snow cover area fraction observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) to further improve the quality of the land surface estimates from the reanalysis. As a first step towards the assimilation of land surface observations in MERRA-3, the offline (land-only) M21C-Land reanalysis is currently under development as a supplemental M21C product that includes the assimilation of SMAP, ASCAT, and MODIS observations. Preliminary results from M21C and M21C-Land will be discussed in the context of MERRA-2 and plans for MERRA-3.

Rolf Reichle↗