Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble data 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 217 records · Page 12

Models, In situ, and Remote sensing of Aerosols (MIRA)

There is a natural partitioning of scientific interest amongst three specialties of aerosol research: modeling, in situ measurements, and remote sensing. The community sees enhanced measurement capabilities when these groups interact, and this strengthens the overall scientific impact on climate and air quality. The Models, In situ, and Remote sensing of Aerosols (MIRA) working group connects members of the different aerosol communities through collaborative projects. What is MIRA? MIRA is a forum that fosters international collaborations amongst the aerosol specialties. MIRA is also a collection of interdisciplinary projects with clear goals that are pursued by small working groups. Finally, MIRA projects are generally characterized by requests for additional scientific data (both observational and modeled). Why? The purpose of MIRA is to contextualize both observations and model results through the encouragement of holistic projects and collaborations. How does MIRA differ from other projects? MIRA focuses on interdisciplinarity to improve measurements and their utility, so MIRA complements the activities of other groups. For example, ensemble model runs of AeroCom could be used in a MIRA project with greater robustness than a similar effort that uses single-model analyses. We present a description of MIRA and brief descriptions of some of the current MIRA topics, which include Satellite-Assisted Particulate Matter (SAPM), Mapping of Aerosol lidar ratios for CALIPSO (MAC), Tables of Aerosol Optics (TAO), and the Harmonization of aerosol Assimilation Models and Retrievals (HAMR). Finally, we discuss the immediate science goals and organization of MIRA. See https://science.larc.nasa.gov/mira-wg/ for more information.

Travis Toth↗

High resolution assimilation of multiple satellite retrievals with emissions adjustment to improve air quality forecasting with WRF-Chem/DART

We will present results from medium (15km, 6hr cycling) and high (4 km, 6 hr cycling) spatiotemporal resolution applications of the WRF-Chem/DART ensemble, regional, air quality (AQ) forecast/assimilation system.The medium-resolution setup is applied to the Discover AQ/Front Range Air Pollution and Photochemistry Experiment (FRAPPE) domain from July 14 to July 29, 2014. The high-resolution setup is applied to a Colorado domain from July 14 to July 29, 2020. For the FRAPPE application, we assimilate MOPITT CO; IASI CO;MODIS AOD; OMI O3, NO2; and AirNow CO, O3, NO2, SO2, PM10, and PM2.5. For the Colorado application, we assimilate the same MOPITT, MODIS, and AirNow constituents as in the FRAPPE application and TROPOMI CO, O3, NO2; and synthetic TEMPO O3and NO2. WRF-Chem/DART integrates the Weather Research and Forecast (WRF) model with on-line chemistry (WRF-Chem) into the Data Assimilation Research Testbed (DART). It assimilates AirNow CO, O3, NO2, SO2, PM10, and PM2.5 measurements, MOPITT CO; IASI CO, O3; OMI O3, NO2, SO2; TROPOMI CO, O3, NO2, SO2; MODIS AOD; and synthetic TEMPO O3 and NO2 total/partial column and/or profile retrievals.WRF-Chem/DART uses: (i) the state augmentation method for adjusting emissions; (ii) state-space localization; and (iii) a near-real time scripting system. We use the medium-resolution FRAPPE application to demonstrate the incremental benefits from assimilating OMI observations with emissions adjustment and the high-resolution Colorado application to demonstrate the incremental benefits from assimilating syntheticTEMPO observations with emissions adjustment. For both applications, we expect that: (i) assimilating chemical observations will increaseAQ forecast skill; (ii) including emissions adjustment will increase forecast skill/predictability; and (iii) including assimilation of synthetic TEMPO observations will further increase forecast skill/predictability.

High resolution↗

Benchmarking a Soil Moisture Data Assimilation System for Agricultural Drought Monitoring

Despite considerable interest in the application of land surface data assimilation systems (LDAS) for agricultural drought applications, relatively little is known about the large-scale performance of such systems and, thus, the optimal methodological approach for implementing them. To address this need, this paper evaluates an LDAS for agricultural drought monitoring by benchmarking individual components of the system (i.e., a satellite soil moisture retrieval algorithm, a soil water balance model and a sequential data assimilation filter) against a series of linear models which perform the same function (i.e., have the same basic inputoutput structure) as the full system component. Benchmarking is based on the calculation of the lagged rank cross-correlation between the normalized difference vegetation index (NDVI) and soil moisture estimates acquired for various components of the system. Lagged soil moistureNDVI correlations obtained using individual LDAS components versus their linear analogs reveal the degree to which non-linearities andor complexities contained within each component actually contribute to the performance of the LDAS system as a whole. Here, a particular system based on surface soil moisture retrievals from the Land Parameter Retrieval Model (LPRM), a two-layer Palmer soil water balance model and an Ensemble Kalman filter (EnKF) is benchmarked. Results suggest significant room for improvement in each component of the system.

soil moisture↗

How Accurately Do Maize Crop Models Simulate the Interactions of Atmospheric CO2 Concentration Levels With Limited Water Supply on Water Use and Yield?

This study assesses the ability of 21 crop models to capture the impact of elevated CO2 concentration [CO2] on maize yield and water use as measured in a 2-year Free Air Carbon dioxide Enrichment experiment conducted at the Thunen Institute in Braunschweig, Germany (Manderscheid et al. 2014). Data for ambient [CO2] and irrigated treatments were provided to the 21 models for calibrating plant traits, including weather, soil and management data as well as yield, grain number, above ground biomass, leaf area index, nitrogen concentration in biomass and grain, water use and soil water content. Models differed in their representation of carbon assimilation and evapotranspiration processes. The models reproduced the absence of yield response to elevated [CO2] under well-watered conditions, as well as the impact of water deficit at ambient [CO2], with 50 percent of models within a range of plus/minus 1 Mg ha(exp. -1) around the mean. The bias of the median of the 21 models was less than 1 Mg ha(exp. -1). However under water deficit in one of the two years, the models captured only 30 percent of the exceptionally high [CO2] enhancement on yield observed. Furthermore the ensemble of models was unable to simulate the very low soil water content at anthesis and the increase of soil water and grain number brought about by the elevated [CO2] under dry conditions. Overall, we found models with explicit stomatal control on transpiration tended to perform better. Our results highlight the need for model improvement with respect to simulating transpirational water use and its impact on water status during the kernel-set phase.

grain number↗

Quantifying the Aerosol Semi-Direct Effect in the NASA GEOS-5 AGCM

Aerosols such as black carbon, dust, and some organic carbon species both scatter and absorb incoming solar radiation. This direct aerosol radiative forcing (DARF) redistributes solar energy both by cooling the surface and warming the atmosphere. As a result, these aerosols affect atmospheric stability and cloud cover (the semi-direct effect, or SDE). Furthermore, in regions with persistent high loadings of absorbing aerosols (e.g. Asia), regional circulation patterns may be altered, potentially resulting in changes in precipitation patterns. Here we investigate aerosol-climate coupling using the NASA Goddard Earth Observing System model version 5 (GEOS-5) atmospheric general circulation model (AGCM), in which we have implemented an online version of the Goddard Chemistry, Aerosol, Radiation and Transport (GOCART) model. GOCART includes representations of the sources, sinks, and chemical transformation of externally mixed dust, sea salt, sulfate, and carbonaceous aerosols. We examine a series of free-running ensemble climate simulations of the present-day period (2000-2009) forced by observed sea surface temperatures to determine the impact of aerosols on the model climate. The SDE and response of each simulation is determined by differencing with respect to the control simulation (no aerosol forcing). In a free-running model, any estimate of the SDE includes changes in clouds due both to atmospheric heating from aerosols and changes in circulation. To try and quantify the SDE without these circulation changes we then examine the DARF and SDE in GEOS-5 with prescribed meteorological analyses introduced by the MERRA analysis. By doing so, we are able to examine changes in model clouds that occur on shorter scales (six hours). In the GEOS-5 data assimilation system (DAS), the analysis is defined as the best estimate of the atmospheric state at any given time, and it is determined by optimally combining a first-guess short-term GCM forecast with all available observations. The Incremental Analysis Update (IAU) is added to the model forecast tendencies to align them with the analysis every six hours, thus preventing longer timescale feedbacks due to the aerosol forcing. We calculate the SDE by comparing model runs with and without aerosols, and the difference in the IAU between these runs is a useful metric with which to evaluate the impact of the SDE on the model atmosphere and clouds. Decreasing the IAU indicates that the aerosol direct and semi-direct effects act to reduce the bias between the model and observations and vice versa.

Randles, Cynthia A.↗

Preliminary Results Cycling GEOS-JEDI with GSI-based Background Errors

The first phase of transitioning the NASA GMAO GEOS atmospheric data assimilation capabilities to JEDI involves the replacement of the Grid-point Statistical Interpolation (GSI) with a corresponding JEDI analysis. This includes taking JEDI's Unified Observation Operator (UFO), its underlying dependencies, and the JEDI solver that enables a hybrid 4DEnVar strategy similar to what is used in the current GEOS-GSI system. Variational analysis involves at least two main components associated with the observation and background cost function terms. The first is directly related to the UFO, which is being carefully validated in a joint collaboration between GMAO and NCEP to demonstrate consistency with corresponding observations usage in GSI. The second component is the background term, which in a hybrid system involves the ability to set up both a climatologically-based term and an ensemble-based term. JEDI provides the means to implement both terms through its BUMP component. Use of BUMP would require a complete re-tune of both climatological and ensemble, which is a non-trivial exercise we would prefer to avoid. As an alternative, the work here studies the results of interfacing the GSI-background error capability (GSIBEC) into JEDI through SABER. With this, the exact same background error covariance formulation used in GSI can be employed in JEDI without need for re-tuning. This brief summary covers the work done to interface GSIBEC into JEDI and shows preliminary results where the background error covariances of the control (GEOS-GSI) and experiment (GEOS-JEDI) are identical in corresponding cycling experiments. The cycling exercise is obviously preliminary and so much can be expected from GEOS-JEDI when compared to GEOS-GSI. There is still a number of features that need closer attention and although in some cases in principle ready to cycle have been intentionally either turned off or not fully exercised (e.g., VarBC is applied but not cycled). Other features are still pending implementation, one such example is the implementation of the Tangent Linear Normal Mode Constraint. Still, results are quite encouraging as hopefully the discussion here illustrates.

Ricardo Todling↗

Dynamically stratified Monte Carlo forecasting

A new method for performing Monte Carlo forecasts is introduced. The method, called dynamic stratification, selects initial perturbations based on a stratification of the error distribution. A simple implementation is presented in which the error distribution used for the stratification is estimated from a linear model derived from a large ensemble of 12-h forecasts with the full dynamic model. The stratification thus obtained is used to choose a small subsample of initial states with which to perform the dynamical Monte Carlo forecasts. Several test cases are studied using a simple two-level general circulation model with uncertain initial conditions. It is found that the method provides substantial reductions in the sampling error of the forecast mean and variance when compared to the more traditional approach of choosing the initial perturbations at random. The degree of improvement, however, is sensitive to the nature of the initial error distribution and to the base state. In practice the method may be viable only if the computational burden involved in obtaining an adequate estimate of the error distribution is shared with the data-assimilation procedure.

Schubert, Siegfried↗

Models, In situ, and Remote sensing of Aerosols (MIRA): Formation of an International Working Group

There is a natural partitioning of scientific interest amongst three specialties of aerosol research: modeling, in situ measurements, and remote sensing. The community sees enhanced measurement capabilities when these groups interact, and this strengthens the overall scientific impact on climate and air quality. The Models, In situ, and Remote sensing of Aerosols (MIRA) working group connects members of the different aerosol communities through collaborative projects. What is MIRA? MIRA is forum that fosters international collaborations amongst the aerosol specialties. MIRA is also a collection of interdisciplinary projects with clear goals that are pursued by small working groups. Finally, MIRA projects are generally characterized by requests for additional scientific data (both observational and modeled). Why? The purpose of MIRA is to contextualize both observations and model results through the encouragement of holistic projects and collaborations. How does MIRA differ from other projects? MIRA focuses on interdisciplinarity to improve measurements and their utility, so MIRA complements the activities of other groups. For example, ensemble model runs of AeroCom could be used in a MIRA project with greater robustness than a similar effort that uses single-model analyses. We first present a description of MIRA and how your project can become a part of the MIRA community. This is followed by brief descriptions of some of the current MIRA projects, which include the Mapping of Aerosol lidar ratios for CALIPSO (MAC), the Tables of Aerosol Optics (TAO), and the Harmonization of aerosol Assimilation Models and Retrievals (HAMR). Finally, we discuss the immediate science goals and organization of MIRA.

Gregory Schuster↗

An Assessment of Actual and Potential Building Climate Zone Change and Variability From the Last 30 Years Through 2100 Using NASA's MERRA and CMIP5 Simulations

Background: In the US, residential and commercial building infrastructure combined consumes about 40% of total energy usage and emits about 39% of total CO2 emission (DOE/EIA "Annual Energy Outlook 2013"). Building codes, as used by local and state enforcement entities are typically tied to the dominant climate within an enforcement jurisdiction classified according to various climate zones. These climate zones are based upon a 30-year average of local surface observations and are developed by DOE and ASHRAE. Establishing the current variability and potential changes to future building climate zones is very important for increasing the energy efficiency of buildings and reducing energy costs and emissions in the future. Objectives: This paper demonstrates the usefulness of using NASA's Modern Era Retrospective-analysis for Research and Applications (MERRA) atmospheric data assimilation to derive the DOE/ASHRAE building climate zone maps and then using MERRA to define the last 30 years of variability in climate zones for the Continental US. An atmospheric assimilation is a global atmospheric model optimized to satellite, atmospheric and surface in situ measurements. Using MERRA as a baseline, we then evaluate the latest Climate Model Inter-comparison Project (CMIP) climate model Version 5 runs to assess potential variability in future climate zones under various assumptions. Methods: We derive DOE/ASHRAE building climate zones using surface and temperature data products from MERRA. We assess these zones using the uncertainties derived by comparison to surface measurements. Using statistical tests, we evaluate variability of the climate zones in time and assess areas in the continental US for statistically significant trends by region. CMIP 5 produced a data base of over two dozen detailed climate model runs under various greenhouse gas forcing assumptions. We evaluate the variation in building climate zones for 3 different decades using an ensemble and quartile statistics to provide an assessment of potential building climate zone changes relative to the uncertainties demonstrated using MERRA. Findings and Conclusions: These results show that there is a statistically significant increase in the area covered by warmer climate zones and a tendency for a reduction of area in colder climate zones in some limited regions. The CMIP analysis shows that models vary from relatively little building climate zone change for the least sensitive and conservation assumptions to a warming of at most 3 zones for certain areas, particularly the north central US by the end of the 21st century.

Stackhouse, Paul W., Jr.↗

Excessive Subsidence of the Southeast Atlantic Biomass Burning Aerosol Plume in the NASA Goddard Earth Observing System (GEOS) Model

Global models struggle to accurately represent the meteorological environment over the Southeast Atlantic Ocean, as evidenced by the poor representation of persistent marine stratocumulus decks off the west coast of southern Africa. This complex regime also experiences significant seasonal biomass burning aerosol plumes that are transported from southern Africa over the Atlantic Ocean, which gradually descend via large scale subsidence and interact radiatively and microphysically with the cloudy environment of the region. In comparisons to space borne lidar and aircraft observations we find that the NASA Goddard Earth Observing System (GEOS) model’s MERRA-2 aerosol and meteorology reanalysis shows that the aerosol plume descends too rapidly due to exaggerated large-scale downward vertical motion. This complicates the evaluation of aerosol properties utilizing aircraft observations collected throughout the ORACLES campaign as well as the model’s representation of the heating rate profiles due to the aerosol itself and its radiative interaction with the semi-permanent cloud deck typically located lower in altitude than the smoke. Two types of modeling experiments using GEOS will be presented to diagnose deficiencies in the model related to the excessive subsidence. Ensembles of free-running model simulations with and without radiatively active biomass burning aerosol will be used to quantify the impact of smoke on the atmospheric circulation. It will be demonstrated that heating within the aerosol layer results in self lofting through a reduction of the subsidence but did not alter the vertical placement or speed of the African Easterly jet. The amount of aerosol absorption required to eliminate the bias in the height of the aerosol plume through self-lofting will then be quantified. Finally, a second type of modeling experiment will be discussed in which the meteorology is constrained to observations. This will allow for biases in the height of the aerosol plume due to the general circulation model to be distinguished from biases that may be induced through the data assimilation process.

Allison Collow↗

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗

Regional Climate Studies with Variable-Resolution Stretched-Grid GCMs

A variable resolution GCM using a global stretched grid with fine resolution over the area(s) of interest, is a viable new approach to regional and subregional climate studies and applications. It is an alternative to the widely used nested grid approach introduced a decade ago as a pioneering step in regional climate modeling. The first version of the SG-GCM based on the GEOS (Goddard Earth Observing System) GCM using a finite-difference approximation, has been developed and thoroughly tested during the last few years. Successful simulations have been performed with the SG-GCM for the anomalous regional climate events of the U.S. 1988 summer drought and 1993 summer flood. They have shown the practical feasibility of the SG-approach for regional climate modeling. The GEOS SG-DAS (Data Assimilation System) incorporating the SG-GCM has also been developed and tested. The assimilated regional fields and diagnostics are used for validating the SG-GCM regional simulations. Two new SG-GCMs are being developed. The first is the SG-version of the new NASANCAR FV-GCM (with the finite-volume (FV) dynamics), and the second is the SG-version of the new GCM with spectral-element dynamics. Both GCMs use the WAR CCM4 physics. Using these advanced numerics will provide increased computational efficiency for the new the SG-GCMs, and will allow us to employ more flexible stretching strategies beneficial for the efficient regional down-scaling. The major current developments are focused on: simulating the 1997-1999 (and beyond) ENSO cycle and related monsoonal circulations, with enhanced regional resolution; studying intraseasonal and interannual regional climate variability for the extended multiyear (AMIP-type) SG-GCM simulations; and studying the impact of ensemble integrations.

Fox-Rabinovitz, Michael↗

Joint Assimilation of CO, O3, NO2, SO2, PM, AOD, NH4, PAN, and HNO3 in Support of the Tropospheric Regional Atmospheric Composition and Emissions Reanalysis (2005 - 2024) (TRACER-I)

NASA Ames Research Center is collaborating with the NOAA Chemical Systems Laboratory (NOAA/CSL), the NASA Jet Propulsion Laboratory (JPL), the NASA Goddard Modeling and Assimilation Office, and the National Center for Atmospheric Research to prepare a 20-year Tropospheric Regional Atmospheric Composition and Emissions Reanalysis (2005 – 2024) (TRACER-I) for the continental United States during the summer ozone (O3) and wildfire seasons (April to October). TRACER-I will be a regional complement to JPL’s global Tropospheric Chemistry Reanalyses II and III. We will use WRF-Chem/DART with NOAA/CSL’s WRF-Chem setup at 12 km x 12 km horizontal resolution, 51 vertical levels, and 30 ensemble members. We will assimilate: (i) conventional meteorological observations; (ii) EPA’s Air Quality System in situ measurements of carbon monoxide (CO), O3, nitrogen dioxide (NO2), sulfur dioxide (SO2), particulate matter (PM) with diameters less than 10 µm (PM10), and PM with diameters less than 2.5 µm (PM2.5); and (iii) MOPITT, MODIS, OMI, TROPOMI, GOME-2a, MLS, TES, SCIAMACHY, CrIS, and TEMPO total/partial column and/or profile retrievals of CO, O3, NO2, SO2, aerosol optical depth (AOD), formaldehyde (HCHO), ammonia (NH4), peroxyacetyl nitrate (PAN), and/or nitric acid (HNO3) with 3-hr cycling. We will present an overview of this project, its status, and available results (likely the analysis of sensitivity experiment results from the assimilation/emissions estimation system).

data assimilation↗

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale of post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of the impacts of fire on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify the impacts of wildfire on runoff generation. We consider the disturbances in the vegetation and soil as the two main factors contributing to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We consider the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration for two case studies. These case studies demonstrate the general applicability of hydrophobicity formulation to serve as a guideline for exploring the range of hydrologic responses post-fire.

Fire Hydrology↗

Clouds in GEOS-5

The GEOS-5 atmospheric model is being developed as a weather-and-climate capable model. It must perform well in assimilation mode as well as in weather and climate simulations and forecasts and in coupled chemistry-climate simulations. In developing GEOS-5, attention has focused on the representation of moist processes. The moist physics package uses a single phase prognostic condensate and a prognostic cloud fraction. Two separate cloud types are distinguished by their source: "anvil" cloud originates in detraining convection, and large-scale cloud originates in a PDF-based condensation calculation. Ice and liquid phases for each cloud type are considered. Once created, condensate and fraction from the anvil and statistical cloud types experience the same loss processes: evaporation of condensate and fraction, auto-conversion of liquid or mixed phase condensate, sedimentation of frozen condensate, and accretion of condensate by falling precipitation. The convective parameterization scheme is the Relaxed Arakawa-Schubert, or RAS, scheme. Satellite data are used to evaluate the performance of the moist physics packages and help in their tuning. In addition, analysis of and comparisons to cloud-resolving models such as the Goddard Cumulus Ensemble model are used to help improve the PDFs used in the moist physics. The presentation will show some of our evaluations including precipitation diagnostics.

Bacmeister, Julio↗

Status and Progress of All-Sky Hyperspectral IR Radiance Assimilation in GEOS

The majority of hyperspectral infrared radiance observations that peak in the lower and mid-troposphere have been excluded from the assimilation in the operational GEOS by the cloud detection scheme. The significance and challenges of assimilating these cloud-affected observations have attracted the attention of the research community, encouraged by the progress in utilizing all-sky microwave radiances affected by clouds and precipitation in many data assimilation systems, including GEOS. While the Global Modeling and Assimilation Office (GMAO) has made significant progress on the assimilation of cloud-cleared infrared radiances, this parallel study marks the first effort to directly assimilate cloud-affected hyperspectral infrared radiances in the GEOS. This initial study is focused on CrIS-FSR wa:wqter-vapor channels. The capability of the GEOS model and CRTM in simulating infrared radiances affected by clouds has been evaluated, and the all-sky brightness temperature Jacobians with respect to temperature, specific humidity, and hydrometeors are examined closely in varying cloud conditions. With hydrometeor control variables in the GEOS, several important aspects of all-sky IR radiance assimilation are investigated. The symmetric cloud effects are assessed for their robustness when used as the cloud proxy in observation error modeling that incorporates the inter-channel correlations with error situation dependency on the amount of cloud. The effective radius, bias correction and quality control processes will also be adapted to accommodate modifications necessary for all-sky infrared assimilation before comprehensive four-dimensional ensemble-variational(4dEnVar) experiments are conducted to assess their impacts on the analysis and forecast performance of GEOS.

Wei Gu↗

An Overview of the GEOS-5 Aerosol Reanalysis

GEOS-5 is the latest version of the NASA Global Modeling and Assimilation Office (GMAO) earth system model. GEOS-5 contains components for atmospheric circulation and composition (including data assimilation), ocean circulation and biogeochemistry, and land surface processes. In addition to traditional meteorological parameters, GEOS-5 includes modules representing the atmospheric composition, most notably aerosols and tropospheric/stratospheric chemical constituents, taking explicit account of the impact of these constituents on the radiative processes of the atmosphere. MERRA is a NASA meteorological reanalysis for the satellite era (1979-present) using GEOS-5. This project focuses on historical analyses of the hydrological cycle on a broad range of weather and climate time scales. As a first step towards an integrated Earth System Analysis (IESA), the GMAO is extending MERRA with reanalyses for other components of the earth system: land, ocean, bio-geochemistry and atmospheric constituents. In this talk we will present results from the MERRA-driven aerosol reanalysis covering the Aqua period (2003-present). The assimilation of Aerosol Optical Depth (AOD) in GEOS-5 involves very careful cloud screening and homogenization of the observing system by means of a Neural Net scheme that translates MODIS radiances into AERONET calibrated AOD. These measurements are further quality controlled using an adaptive buddy check scheme, and assimilated using the Local Displacement Ensemble (LDE) methodology. For this reanalysis, GEOS-5 runs at a nominal 50km horizontal resolution with 72 vertical layers (top at approx. 8Skm). GEOS-5 is driven by daily biomass burning emissions derived from MODIS fire radiative power retrievals. We will present a summary of our efforts to validate such dataset. The GEOS-5 assimilated aerosol fields are first validated by comparison to independent in-situ measurements (AERONET and PM2.5 surface concentrations). In order to asses aerosol absorption on a global scale, we perform a detailed radiative transfer calculation to simulate the UV aerosol index, comparing our results to OMI measurements. By simulating aerosol attenuated backscatter, we use CALIPSO measurements to evaluate the vertical structure of our aerosol estimates, in particular in regions where we have larger discrepancies with OMI. Finally, the consistency of our AOD estimates with estimates from MISR, MODIS/Deep Blue, OMI and PARASOL will be briefly discussed.

da Silva, Arlindo↗

Evaluation of long-term Northern Hemisphere snow water equivalent products

Nine gridded Northern Hemisphere snow water equivalent (SWE) products were evaluated as part of the European Space Agency (ESA) Satellite Snow Product Intercomparison and Evaluation Exercise (SnowPEx). Three categories of datasets were assessed: (1) those utilizing some form of reanalysis (the NASA Global Land Data Assimilation System version 2 – GLDAS-2; the European Centre for Medium-Range Weather Forecasts (ECMWF) interim land surface reanalysis – ERA-Interim/Land and ERA5; the NASA Modern-Era Retrospective Analysis for Research and Applications version 1 (MERRA) and version 2 (MERRA-2); the Crocus snow model driven by ERA-Interim meteorology – Crocus); (2) passive microwave remote sensing combined with daily surface snow depth observations (ESA GlobSnow v2.0); and (3) stand-alone passive microwave retrievals (NASA AMSR-E SWE versions 1.0 and 2.0) which do not utilize surface snow observations. Evaluation included validation against independent snow course measurements from Russia, Finland, and Canada and product intercomparison through the calculation of spatial and temporal correlations in SWE anomalies. The stand-alone passive microwave SWE products (AMSR-E v1.0 and v2.0 SWE) exhibit low spatial and temporal correlations to other products and RMSE nearly double the best performing product. Constraining passive microwave retrievals with surface observations (GlobSnow) provides performance comparable to the reanalysis-based products; RMSE over Finland and Russia for all but the AMSR-E products is ∼50 mm or less, with the exception of ERA-Interim/Land over Russia. Using a seven-dataset ensemble that excluded the stand-alone passive microwave products reduced the RMSE by 10 mm (20 %) and increased the correlation from 0.67 to 0.78 compared to any individual product. The overall performance of the best multiproduct combinations is still at the margins of acceptable uncertainty for scientific and operational requirements; only through combined and integrated improvements in remote sensing, modeling, and observations will real progress in SWE product development be achieved.

Northern Hemisphere snow water equivalent (SWE) pr↗