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 307 records · Page 17

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↗

Evaluation and Uncertainty of Radioxenon Transport with a Mesoscale Model after the February 2013 Underground Test in North Korea

Abstract The transport of radioxenon released from the February 2013 underground nuclear weapons test in North Korea was analyzed at two receptors—one at the Comprehensive Nuclear Test Ban Treaty site Rn58 in Russia (400 km downwind) and a second at Rn38 in Japan (1000 km downwind). Transport was modeled with two ensembles of mesoscale simulations, one generated with varying initial and lateral boundary conditions taken from the Global Forecasting System uncertainty ensemble, and a second created from different parameterizations and surface conditions. The wind variability was similar for the two ensembles and consistent with observations at 925 mb (1 mb = 1 hPa) but not at the surface. Biases in calculated surface winds and the radioxenon concentration in the ensembles were attributed mainly to poor simulation of the sea breeze at both locations and mountain lee affects at Rn38 in Japan. These wind regimes affected the timing of the surface radioxenon plume at Rn58 and its duration at Rn38. Surface wind variability induced by terrain and land–sea contrast (the sea breeze) also had a significant effect on the surface winds and plume dynamics, including blocking of flow approaching elevated terrain near Vladivostok and the west side of Japan. Increased plume uncertainty was seen at night because of surface wind variability. Measured surface chemical variability was larger than found in the first European Tracer Experiment in central Europe. The study found that horizontal model resolution contributes to uncertainty but not as much as vertical resolution, boundary layer parameterizations, and assimilation of surface meteorological data near the receptor.

54 ENVIRONMENTAL SCIENCES↗

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.↗

Time-resolved particle image velocimetry

Time-resolved particle image velocimetry (TR-PIV) has become widespread in fluid dynamics. Essentially a velocity field movie, the dynamic content provides temporal as well as spatial information, in contrast to conventional PIV offering only statistical ensembles of flow quantities. From these time series arise further analyses such as accelerometry, space-time correlations, frequency spectra of turbulence including spatial variability, and derivation of pressure fields and forces. The historical development of TR-PIV is chronicled, culminating in an assessment of the current state of technology in high-repetition-rate lasers and high-speed cameras. Commercialization of pulse-burst lasers has expanded TR-PIV into more flows, including the compressible regime, and has achieved MHz rates. Particle response times and peak locking during image interrogation require attention but generally are not impediments to success. Accuracy considerations are discussed, including the risks of noise and aliasing in spectral content. Oversampled TR-PIV measurements allow use of multi-frame image interrogation methods, which improve the precision of the correlation and raise the velocity dynamic range of PIV. In combination with volumetric methods and data assimilation, a full four-dimensional description of a flow is not only achievable but becoming standardized. Lastly, a survey of exemplary applications is followed by a few predictions concerning the future of TR-PIV.

47 OTHER INSTRUMENTATION↗

Assimilation of citizen science data in snowpack modeling using a new snow data set: Community Snow Observations

A physically based snowpack evolution and redistribution model was used to test the effectiveness of assimilating crowd-sourced snow depth measurements collected by citizen scientists. The Community Snow Observations project gathers, stores, and distributes measurements of snow depth recorded by recreational users and snow professionals in high mountain environments. These citizen science measurements are valuable since they come from terrain that is relatively undersampled and can offer in situ snow information in locations where snow information is sparse or nonexistent. The present study investigates (1) the improvements to model performance when citizen science measurements are assimilated, and (2) the number of measurements necessary to obtain those improvements. Model performance is assessed by comparing time series of observed (snow pillow) and modeled snow water equivalent values, by comparing spatially distributed maps of observed (remotely sensed) and modeled snow depth, and by comparing fieldwork results from within the study area. The results demonstrate that few citizen science measurements are needed to obtain improvements in model performance, and these improvements are found in 62 % to 78 % of the ensemble simulations, depending on the model year. Model estimations of total water volume from a subregion of the study area also demonstrate improvements in accuracy after CSO measurements have been assimilated. These results suggest that even modest measurement efforts by citizen scientists have the potential to improve efforts to model snowpack processes in high mountain environments, with implications for water resource management and process-based snow modeling.

54 ENVIRONMENTAL SCIENCES↗

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↗

Updating and Evaluating Anthropogenic Emissions for NOAA’s Global Ensemble Forecast Systems for Aerosols (GEFS-Aerosols): Application of an SO 2 Bias-Scaling Method

We updated the anthropogenic emissions inventory in NOAA’s operational Global Ensemble Forecast for Aerosols (GEFS-Aerosols) to improve the model’s prediction of aerosol optical depth (AOD). We used a methodology to quickly update the pivotal global anthropogenic sulfur dioxide (SO 2 ) emissions using a speciated AOD bias-scaling method. The AOD bias-scaling method is based on the latest model predictions compared to NASA’s Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA2). The model bias was subsequently applied to the CEDS 2019 SO 2 emissions for adjustment. The monthly mean GEFS-Aerosols AOD predictions were evaluated against a suite of satellite observations (e.g., MISR, VIIRS, and MODIS), ground-based AERONET observations, and the International Cooperative for Aerosol Prediction (ICAP) ensemble results. The results show that transitioning from CEDS 2014 to CEDS 2019 emissions data led to a significant improvement in the operational GEFS-Aerosols model performance, and applying the bias-scaled SO 2 emissions could further improve global AOD distributions. The biases of the simulated AODs against the observed AODs varied with observation type and seasons by a factor of 3~13 and 2~10, respectively. The global AOD distributions showed that the differences in the simulations against ICAP, MISR, VIIRS, and MODIS were the largest in March–May (MAM) and the smallest in December–February (DJF). When evaluating against the ground-truth AERONET data, the bias-scaling methods improved the global seasonal correlation (r), Index of Agreement (IOA), and mean biases, except for the MAM season, when the negative regional biases were exacerbated compared to the positive regional biases. The effect of bias-scaling had the most beneficial impact on model performance in the regions dominated by anthropogenic emissions, such as East Asia. However, it showed less improvement in other areas impacted by the greater relative transport of natural emissions sources, such as India. The accuracies of the reference observation or assimilation data for the adjusted inputs and the model physics for outputs, and the selection of regions with less seasonal emissions of natural aerosols determine the success of the bias-scaling methods. A companion study on emission scaling of anthropogenic absorbing aerosols needs further improved aerosol prediction.

54 ENVIRONMENTAL SCIENCES↗

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↗

The Re-Analysis of Ozone Profile Data from a 41-Year Series of SBUV Instruments

In this study we present the validation of ozone profiles from a number of Solar Back Scattered Ultra Violet (SBUV) and SBUV/2 instruments that were recently reprocessed using an updated (Version 8.6) algorithm. The SBUV dataset provides the longest available record of global ozone profiles, spanning a 41-year period from 1970 to 2011 (except a 5-year gap in the 1970s) and includes ozone profile records obtained from the Nimbus-4 BUV and Nimbus-7 SBUV instruments, and a series of SBUV(/2) instruments launched on NOAA operational satellites (NOAA 09, 11, 14, 16, 17, 18, 19). Although modifications in instrument design were made in the evolution from the BUV instrument to the modern SBUV(/2) model, the basic principles of the measurement technique and retrieval algorithm remain the same. The long term SBUV data record allows us to create a consistent, calibrated dataset of ozone profiles that can be used for climate studies and trend analyses. In particular, we focus on estimating the various sources of error in the SBUV profile ozone retrievals using independent observations and analysis of the algorithm itself. For the first time we include in the metadata a quantitative estimate of the smoothing error, defined as the error due to profile variability that the SBUV observing system cannot inherently measure. The magnitude of the smoothing error varies with altitude, latitude, season and solar zenith angle. Between 10 and 1 hPa the smoothing errors for the SBUV monthly zonal mean retrievals are of the order of 1 %, but start to increase above and below this layer. The largest smoothing errors, as large as 15-20%, were detected in in the troposphere. The SBUV averaging kernels, provided with the ozone profiles in version 8.6, help to eliminate the smoothing effect when comparing the SBUV profiles with high vertical resolution measurements, and make it convenient to use the SBUV ozone profiles for data assimilation and model validation purposes. The smoothing error can also be minimized by combining layers of data, and we will discuss recommendations for this approach as well. The SBUV ozone profiles have been intensively validated against satellite profile measurements obtained from the Microwave Limb Sounders (MLS) (on board the UARS and AURA satellites), Stratospheric Aerosol and Gas Experiment (SAGE) and Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Also, we compare coincident and collocated SBUV ozone retrievals with observations made by ground-based instruments, such as microwave spectrometers, lidars, Umkehr instruments and balloon-borne ozonosondes. Finally, we compare the SBUV ozone profiles with output from the NASA GSFC GEOS-CCM model. In the stratosphere between 25 and 1 hPa the mean biases and standard deviations are within 5% for monthly mean ozone profiles. Above and below this layer the vertical resolution of the SBUV algorithm decreases and the effects of vertical smoothing should be taken into account. Though the SBUV algorithm has a coarser vertical resolution in the lower stratosphere and troposphere, it is capable of precisely estimating the integrated ozone column between the surface and 25 hPa. The time series of the tropospheric - lower stratospheric ozone column derived from SBUV agrees within 5% with the corresponding values observed by an ensemble of ozone sonde stations in North Hemisphere. Drift of the ozone time series obtained from each SBUV(/2) instrument relative to ground based and satellite measurements are evaluated and some features of individual SBUV(l2) instruments are discussed. In addition to evaluating individual instruments against independent observations, we also focus on the instrument to instrument consistency in the series. Overall, Version 8.6 ozone profiles obtained from two different SBUV(l2) instruments compare within a couple of percent during overlap periods and are consistently varying in time, with some exceptions. Some of the noted discrepancies might bssociated with ozone diurnal variations, since the difference in the local time of the observations for a pair of SBUV(l2) instruments could be several hours. Other issues include the potential short-term drift in measurements as the instrument orbit drifts, and measurements are obtained at high solar zenith angles (>85 ). Based on the results of the validation, a consistent, calibrated dataset of SBUV ozone profiles has been created based on internal calibration only.

Kramarova, Natalya↗

Satellite-Assisted Particulate Matter (SAPM) for the Models, In situ, and Remote Sensing of Aerosols (MIRA) Working Group

Models, In situ, and Remote sensing of Aerosols (MIRA) is an international working group that provides a forum for collaborations amongst these three atmospheric aerosol communities. MIRA consists of a collection of interdisciplinary and independently funded Topic/Project groups with clear goals and generally characterized by requests for additional scientific data. The Satellite-Assisted Particulate Matter (SAPM) Topic group focuses on the study of fine particulate matter (PM2.5), as it is a major contributor to air pollution and negatively impacts human health. Our group aims to provide intercomparisons of various methods and techniques for estimating surface PM2.5 assisted by satellite remote sensors (passive and active), global aerosol models, and in situ aerosol measurements. The overall motivation of our work is to enhance PM2.5 coverage beyond in situ ground stations, which can be limited. In this presentation, we provide an overview of various PM2.5 estimation approaches by current SAPM team members, as each approach has its own strengths and limitations. One approach uses near-surface aerosol extinction retrievals from the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and an assumed value of mass extinction efficiency to derive PM2.5 concentrations (Toth et al 2019). A study using this method found a promising agreement (R=0.60; slope=0.89) with U.S. Environmental Protection Agency (EPA) in situ PM2.5 measurements over the contiguous United States (CONUS) for a 12-year (2007-2018) period (Toth et al 2022). In a different approach combining spaceborne lidar and a model, the Cloud Aerosol Transport System (CATS) lidar and Goddard Earth Observing System (GEOS) model are blended in an ensemble variational assimilation scheme to retrieve aerosol extinction and convert speciated mass concentrations to total PM2.5. A recent study (Matus et al in prep.) found agreement within 2 μg/m3 over the CONUS on the 2016 annual mean when compared to EPA ground-based measurements. Another study explored trends in city aerosol optical depth (AOD) from a spaceborne passive remote sensor, Moderate Resolution Imaging Spectroradiometer (MODIS), and found that these trends agree well with trends in city surface PM2.5, which provides confidence in the use of city AOD for city PM2.5 trend studies (Vohra et al 2021). In a follow-on study, steep and significant trends (~2.5 to 7.8% a-1) in AOD are found for several South Asian cities (e.g., Bangalore and Hyderabad), suggesting rapid growth in PM2.5 (Vohra et al 2022). Also in this presentation, we report results on our current collective SAPM study of India, a country that exhibits high levels of PM2.5 pollution. We show the temporal (i.e., seasonal) and spatial variability of PM2.5 over India using the CALIOP and CATS+model approaches, comparisons between the estimated PM2.5 from the retrieval schemes, and validation with in situ PM2.5 ground-based observations from the Central Pollution Control Board in India. This region provides us with an excellent case study to test our PM2.5 retrievals in heavily polluted scenes and over a complex and varying topography. The SAPM Topic group actively seeks international participants/collaborators in our group, including those working with in situ aerosol measurements (e.g., ACTRIS). We are interested in aerosol datasets in order to either improve our PM2.5 estimates from various approaches (e.g., using in situ mass scattering/absorption coefficient and aerosol hygroscopic properties for various aerosol species) and/or validate the PM2.5 estimates (e.g., using in situ ground-based PM2.5 concentrations).

Travis D Toth↗