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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Improving Regional Air Quality Forecasting Through Chemical Data Assimilation and Dynamic Emissions Adjustment

Poor air quality (AQ) is one of the most important human-health and environmental problems facing the United States (US). In addition to the detrimental impacts on human- and environmental-health, poor AQ has an economic cost of ~5% of the US gross domestic product (~$790 billion). AQ managers use AQ analyses and modeling to better understand, anticipate, and avoid poor AQ events. Our research focuses on improving AQ analysis/forecast skill, predictability, and emission estimates through improved and more efficient: (i) modeling and data assimilation strategies; (ii) dynamic emissions adjustment strategies; and (iii) use of satellite remote-sensing Earth observations (e.g., MOPITT, IASI, MODIS, OMI, TROPOMI, TEMPO, etc.). This seminar will review: (i) regional chemical weather forecasting/data assimilation with dynamic emissions adjustment with WRF-Chem/DART; (ii) strategies for efficiently assimilating satellite retrieval profiles with ‘compact phase space retrievals’ (CPSRs); (iii) results from joint assimilation of multiple satellite retrievals at medium (12 km × 12 km) and high (4 km × 4 km) spatial resolutions; and (iv) results from observing system simulation experiments (OSSEs) to investigate whether we can recover COVID-period anthropogenic emissions by assimilating synthetic TEMPO NO2 tropospheric column retrievals with dynamic emissions adjustment. Biographical Sketch: Dr. Mizzi is a Senior Research Fellow working and Dr. Johnson at the NASA Ames Research Center. He holds BA and MS degrees in Environmental Science from the University of Virginia, MS and PhD degrees in Applied Mathematics from the University of Colorado at Boulder (CUB), and a JD degree (with an emphasis in Environmental Law) from the University of Colorado School of Law. He worked at the National Center for Atmospheric Research for nearly 25 years on global atmospheric modeling, dynamic and physical initialization, regional hybrid data assimilation, and most recently on regional chemical data assimilation. He also worked as an environmental attorney and consultant for nearly 15 years. He is an expert in numerical modeling and is recognized internationally as a leading expert in regional, chemical data assimilation with dynamic emissions adjustment. Dr. Mizzi became affiliated with NASA Ames in March 2020 to work on improving AQ analysis/forecast skill, predictability, and ‘top-down’ emissions adjustment though the assimilation of Earth observations. An emphasis of his current work is developing methods for assimilating synthetic TEMPO retrievals to quantify the expected benefits of TEMPO relative to existing AQ observations.

Arthur P. Mizzi↗

NASA NTRS2022

Soil Moisture Estimation in South Asia via Assimilation of SMAP Retrievals

A soil moisture retrieval assimilation framework is implemented across South Asia in an attempt to improve regional soil moisture estimation as well as to provide a consistent regional soil moisture dataset. This study aims to improve the spatiotemporal variability of soil moisture estimates by assimilating Soil Moisture Active Passive (SMAP) near-surface soil moisture retrievals into a land surface model. The Noah-MP (v4.0.1) land surface model is run within the NASA Land Information System software framework to model regional land surface processes. NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA2) and Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals (IMERG) provide the meteorological boundary conditions to the land surface model. Assimilation is carried out using both cumulative distribution function (CDF)-corrected (DA-CDF) and uncorrected SMAP retrievals (DA-NoCDF). CDF matching is applied to correct the statistical moments of the SMAP soil moisture retrieval relative to the land surface model. Comparison of assimilated and model-only soil moisture estimates with publicly available in situ measurements highlights the relative improvement in soil moisture estimates by assimilating SMAP retrievals. Across the Tibetan Plateau, DA-NoCDF reduced the mean bias and RMSE by 8.4 % and 9.4 %, even though assimilation only occurred during less than 10 % of the study period due to frozen (or partially frozen) soil conditions. The best goodness-of-fit statistics were achieved for the IMERG DA-NoCDF soil moisture experiment. The general lack of publicly available in situ measurements across irrigated areas limited a domain-wide direct model validation. However, comparison with regional irrigation patterns suggested correction of biases associated with an unmodeled hydrologic phenomenon (i.e., anthropogenic influence via irrigation) as a result of SMAP soil moisture retrieval assimilation. The greatest sensitivity to assimilation was observed in cropland areas. Improvements in soil moisture potentially translate into improved spatiotemporal patterns of modeled evapotranspiration, although limited influence from soil moisture assimilation was observed on modeled processes within the carbon cycle such as gross primary production. Improvement in fine-scale modeled estimates by assimilating coarse-scale retrievals highlights the potential of this approach for soil moisture estimation over data-scarce regions.

Jawairia Ahmad↗

NASA NTRSDate not supplied

Joint Assimilation of Surface, GEO, and LEO CO, NO2, and O3 In Situ and Partial Column/Profile Retrievals with WRF-Chem/DART

Low earth orbiting (LEO) satellites like OMI or TROPOMI and geostationary satellites like TEMPO (to be launched January 2023) provide ozone (O3) retrieval profiles that extend into the stratosphere. Regional models like WRF-Chem generally do not include the stratosphere and have vertical grids that end near the tropopause. This causes difficulties when assimilating O3 retrievals (both column and profile retrievals) because the respective observation forward operators require model O3 profiles that extend into the stratosphere. In this presentation, we look at results from assimilating O3 observation from surface (in situ measurements), multiple LEO (profile retrievals), and GEO (proxy profile retrievals) platforms in as regional model where we use a global chemical transport model (GCTM) to provide O3 upper boundary conditions (BCs) to enable application of the O3 observation forward operators. We will also present results from assimilating O3 retrievals together with carbon monoxide (CO) observations from surface (in situ measurements) and multiple LEO platforms (column and profile retrievals), and nitrogen dioxide (NO2) observations from surface (in situ measurements), multiple LEO (column retrievals), and GEO (proxy column retrievals) platforms. Our joint assimilation experiments show that sometimes the observations complement each, e.g., when assimilating surface and satellite CO or NO2 observations, and sometimes the observations conflict with each other, e.g., when assimilating observations that are biased with respect to other observations (e.g., CO retrievals), or when the assimilation increments conflict due to their chemical interaction (e.g., surface O3 and NO2 retrievals). Results from assimilation of O3 retrieval profiles in regional model are preliminary and will also be discussed in our presentation. Our expectation is that using GCTM O3 as upper BCs in the regional model for evaluation of the O3 forward operator will improve our O3 forecasts. In summary, we will present results from jointly assimilating surface, LEO, and GEO CO, NO2, and O3 in situ and satellite observations in a regional model with dynamic emissions estimation. This study uses proxy GEO observations as a prelude to actual observations because the GEO platform will not be launched until after our presentation.

WRF-Chem/DART↗

DOE OSTI2025

Comparing field and lab quantitative stable isotope probing for nitrogen assimilation in soil microbes

ABSTRACT Soil microbial communities play crucial roles in nutrient cycling and can help retain nitrogen in agricultural soils. Quantitative stable isotope probing (qSIP) is a useful method for investigating taxon-specific microbial growth and utilization of specific nutrients, such as nitrogen (N). Typically, qSIP is performed in a highly controlled lab setting, so the field relevance of lab qSIP studies remains unknown. We conducted and compared tandem lab and field qSIP to quantify the assimilation of 15 N by maize-associated soil prokaryotic communities at two agricultural sites. Here, we show that field qSIP with 15 N can be used to measure taxon-specific microbial N assimilation. Relative 15 N assimilation rates were generally lower in the field, and the magnitude of this difference varied by site. Rates differed by method (lab vs field) for 19% of the top N assimilating genera. The field and lab measures were more comparable when relative assimilation rates were weighted by relative abundance to estimate the proportion of N assimilated by each genus with only ~10% of taxa differing by method. Of those that differed, the taxa consistently higher in the lab were inclined to have opportunistic lifestyle strategies, whereas those higher in the field had niches reliant on plant roots or in-tact soil structure (biofilms, mycelia). This study demonstrates that 15 N-qSIP can be successfully performed using field-incubated soils to identify microbial allies in N retention and highlights the strengths and limitations of field and lab qSIP approaches. IMPORTANCE Soil microbes are responsible for critical biogeochemical processes in natural and agricultural ecosystems. Despite their importance, the functional traits of most soil organisms remain woefully under-characterized, limiting our ability to understand how microbial populations influence the transformation of elements such as nitrogen (N) in soil. Quantitative stable isotope probing (qSIP) is a powerful tool to measure the traits of individual taxa. This method has rarely been applied in the field or with 15 N to measure nitrogen assimilation. In this study, we measured genus-specific microbial nitrogen assimilation in two agricultural soils and compared field and lab 15 N qSIP methods. Our results identify taxa important for nitrogen assimilation in agricultural soils, shed light on the field relevance of lab qSIP studies, and provide guidance for the future application of qSIP to measure microbial traits in the field.

Reed, Kinsey (ORCID:0000000155178664)↗

NASA NTRS1995

Sea surface topography fields of the tropical Pacific from data assimilation

Time series of maps of monthly tropical Pacific dynamic topography anomalies from 1979 through 1985 were constructed by means of assimilation of the tide gauge and expendable bathythermograph (XBT) data into a linear model driven by observed winds. Estimates of error statistics were calculated and compared to actual differences between hindcasts and observations. Four experiements were performed as follows: one with no assimilation, one with assimiation of sea level anomaly data from eight selected island tide gauge stations, one with assimilation of dynamic height anomalies derived from XBT data, and one with both XBT and tide gauge data assimilated. Data from seven additional tide gauge stations were withheld from the assimilation process and used for verification in all four experiements. Statistical objective maps based on data alone were also constructed for comparison purposes. The dynamic response of the model without assimilation was, in general, weaker than the observed response. Assimilation resulted in enhanced signal amplitude in all three assimilation experiments. With few exceptions the error estimated generated by the Kalman filter appeared quite reasonable. Since the error processes cannot be assumed to be white or stationary, we could find no straightforward way to test the formal statistical hypothesis that the time series of differences between the filter ouput and the actual observations were drawn from a population with statistics given by the Kalman filter estimaes. We found that the best filter we could devise was still short of the goal of producing a white innovation sequence. Extensive changes in the assumed error statistics make only marginal differences. The same is true for long time and space scale behavior of different models with richer physics and finer resolution. Better data assimilation results will probably require relaxation of the assumptions of stationarity and serial independence of the errors. Formation of such detailed noise models will require longer time series, with the attendant problems of matching very different data sets.

Miller, Robert N.↗

NASA NTRS1999

Current and Future Plans of the NASA Data Assimilation Office (DAO)

The mission of the Data Assimilation Office (DAO) is to advance the state of the art of data assimilation and produce research-quality assimilated data sets which make optimal use of space-based observations. Development efforts over the last few years have focused on delivering a production data assimilation system in support of NASA's Terra platform. That system, called the Goddard Earth Observing System - version 2 or GEOS-2, represents a major upgrade to the baseline GEOS-1 system employed in NASA's first reanalysis effort. GEOS-2 includes a physical-space three dimensional variational analysis algorithm (the Physical-space Statistical Analysis System or PSAS) and numerous improvements to the general circulation model. The latter include a Soil-Vegetation-Atmosphere Transfer (SVAT) land surface scheme, a level 2.5 moist turbulence scheme and new Short Wave (SW) and Long Wave (LW) radiation code. The system also includes an off-line ozone assimilation system, and the capability to assimilate scatterometer surface winds, and TIROS Operational Vertical Sounder (TOVS) and Special Sensor Microwave Imager (SSM/I) moisture data. GEOS-2 is currently run at 1 degree horizontal resolution and 48 levels extending to O.Olmb. Experimental versions of GEOS-2 are run with a global stretched grid allowing enhanced (e.g. 1/4 deg) regional resolution. Other capabilities being developed include, the assimilation of Tropical Rainfall Measuring Mission (TRMM) precipitation and Global Positioning System (GPS) data, an off-line land surface assimilation system, and a retrospective analysis scheme. The DAO is also engaged in a number of collaborative efforts to help accelerate the development of the next generation data assimilation system. These include, a joint modeling effort between the DAO and NCAR/CGDD to develop a new Global Circulation Model (GCM), and a Department of Energy Lawrence Livermore National Laboratory (DOE/LLNL) collaboration on model parallelization. Plans for the next reanalysis will be discussed in the context of current and near term system quality and computing capabilities, and the need for multiple reanalysis products.

Atlas, Robert↗

NASA NTRS2000

Operational Assimilation of GOES Data into a Mesoscale Model

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite- observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The technique has been employed on a semi-operational basis at the Global Hydrology and Climate Center (GHCC) within the Penn State/National Center for Atmospheric Research (PSU/NCAR) Mesoscale Model (MM5) since 1 November 1998. We performed the assimilation on a model grid centered over the Southeastern US. In addition, a control run without assimilation was performed to provide insight into the performance of the assimilation technique. Bulk verification statistics (BIAS and RMSE) of surface air temperature and relative humidity of more than 250 case days has been performed to date. Results show that assimilation of the satellite data results reduces both the bias and RMSE for simulations of surface air temperature and relative humidity. We are working with forecasters at the National Weather Service Forecast Office located in Birmingham, AL to evaluate the impact of the assimilation on precipitation forecasts. In addition, work is currently underway to determine the ability of both the assimilation technique and Early Eta to sense the significant changes that occur in both leaf expansion and soil moisture conditions during the Spring-Summer of 1999. This will be accomplished via examination of the simulated surface energy budgets from both modeling systems.

Lapenta, William↗

NASA NTRS2000

The Impact of Withholding Observations from TOMS or SBUV Instruments on the GEOS Ozone Data Assimilation System

In a data assimilation system (DAS), model forecast atmospheric fields, observations and their respective statistics are combined in an attempt to produce the best estimate of these fields. Ozone observations from two instruments are assimilated in the Goddard Earth Observing System (GEOS) ozone DAS: the Total Ozone Mapping Spectrometer (TOMS) and the Solar Backscatter Ultraviolet (SBUV) instrument. The assimilated observations are complementary; TOMS provides a global daily coverage of total column ozone, without profile information, while SBUV measures ozone profiles and total column ozone at nadir only. The purpose of this paper is to examine the performance of the ozone assimilation system in the absence of observations from one of the instruments as it can happen in the event of a failure of an instrument or when there are problems with an instrument for a limited time. Our primary concern is for the performance of the GEOS ozone DAS when it is used in the operational mode to provide near real time analyzed ozone fields in support of instruments on the Terra satellite. In addition, we are planning to produce a longer term ozone record by assimilating historical data. We want to quantify the differences in the assimilated ozone fields that are caused by the changes in the TOMS or SBUV observing network. Our primary interest is in long term and large scale features visible in global statistics of analysis fields, such as differences in the zonal mean of assimilated ozone fields or comparisons with independent observations, While some drifts in assimilated fields occur immediately, after assimilating just one day of different observations, the others develop slowly over several months. Thus, we are also interested in the length of time, which is determined from time series, that is needed for significant changes to take place.

Stajner, Ovanka↗

NASA NTRS2000

Results from NASA's Near Real Time Ozone Data Assimilation Systems

Organizations in Europe, Australia, and the United States have recently broadened constituent assimilation activities beyond water vapor, which has been assimilated for years in the numerical weather prediction applications. Many of these activities have focused on ozone, with some efforts focused on the entire suite of reactive constituents that control the ozone distribution. This talk will draw from results from the near real-time ozone data assimilation system being run by NASA's Data Assimilation Office. This system utilizes ozone observations from both the TOMS and the SBUV instrument to generate global synoptic maps of ozone. The initial application of this product is to provide ozone fields to assist in the atmospheric corrections' that are necessary for the retrieval of information from other NASA instruments. The validation of the ozone assimilation system shows that the assimilated product agrees well with independent HALOE and ozonesonde observations. Aside from providing a global synoptic map, there is verifiable geophysical information at higher vertical resolution than either of the date types input into the system. This talk will establish the validation results and enumerate applications of the ozone data assimilation system. Results from exploratory research will be presented. The applications being considered include estimates of tropospheric ozone, provision of ozone fields for interactive retrievals, use of analysis increments from the assimilation to evaluate model performance, and development of long-term consistent three-dimensional global ozone fields. The results from the exploratory studies are promising, and help demonstrate how assumptions made in the development of the ozone assimilation impact the other applications. For instance, RMS errors in the current product are large near the tropopause, which is sensitive to the specification of vertical correlation functions, which in turns impacts the amount of ozone analyzed to be in the troposphere. How these sensitivities impact the different applications will also be discussed.

Rood, Richard B.↗

NASA NTRS2002

Sensitivity Studies for Assimilated Ozone Profiles

An ozone data assimilation system at the NASA/Goddard Data Assimilation Office (DAO) produces three-dimensional global ozone fields. They are obtained by assimilating ozone retrieved from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument and the Earth Probe Total Ozone Mapping Spectrometer (EP TOMS) measurements into an off-line parameterized chemistry and transport model. In this talk we focus on the quality of lower stratospheric assimilated ozone profiles. Ozone in the lower stratosphere plays a key role in the forcing of climate. A biased ozone field in this region will adversely impact calculations of the stratosphere-troposphere exchange and, when used as a first guess in retrievals, the values determined from satellite observations. The SBUV/2 ozone data have a coarse vertical resolution with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the lower stratosphere are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. We investigate the sensitivity of assimilated ozone profiles to changes in a variety of system inputs: TOMS and SBUV/2 data selection, forecast and observations error covariance models, inclusion or omission of a parameterized chemistry model, and different versions of DAO assimilated wind fields used to drive the transport model. Comparisons of assimilated ozone fields with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes.

Stajner, Ivanka↗

NASA NTRS2004

Low Frequency Oscillations in Assimilated Global Datasets Using TRMM Rainfall Observations

Global datasets for the period May-August 1998 from the Goddard Earth Observing System (GEOS) data assimilation system (DAS) with/without assimilated Tropical Rainfall Measuring Mission (TRMM) precipitation are analyzed against European Center for Medium-Range Weather Forecast (ECMWF) output, NOAA observed outgoing longwave radiation (OLR) data, and TRMM measured rainfall. The purpose of this study is to investigate the representation of the Madden-Julian Oscillation (MJO) in GEOS assimilated global datasets, noting the impact of TRMM observed rainfall on the MJO in GEOS data assimilations. A space-time analysis of the OLR data indicates that the observed OLR exhibits a spectral maximum for eastward-propagating wavenumber 1-3 disturbances with periods of 20-60 days in the 0deg-30degN latitude band. The assimilated OLR has a similar feature but with a smaller magnitude. However, OLR spectra from assimilations including TRMM rainfall data show better agreement with observed OLR spectra than spectra from assimilations without TRMM rainfall. Similar results are found for wavenumber 4-6 disturbances. There is a spectral peak for eastward-propagating wavenumber 4-6 disturbances with periods of 20-40 days near the equator, while for westward-moving disturbances, a spectral peak is noted for periods of 30-50 days near 25degN. To isolate the MJO, a 30-50 day band filter is selected for this study. It was found that the eastward-propagating waves from the band-filtered observed OLR between 10degs- 10degN are located in the eastern hemisphere. Similar patterns are evident in surface rainfall and the 850 hPa wind field. Assimilation of TRMM-observed rainfall reveals more distinct MJO features in the analysis than without rainfall assimilation. Similar analyses are also conducted over the Indian summer monsoon and East Asia summer monsoon regions, where the MJO is strongly related to the summer monsoon active-break patterns.

Tao, Li↗

NASA NTRS2004

Assimilation of SBUV Version 8 Radiances into the GEOS Ozone DAS

In operational weather forecasting, the assimilation of brightness temperatures from satellite sounders, instead of assimilation of 1D-retrievals has become increasingly common practice over the last two decades. Compared to these systems, assimilation of trace gases is still at a relatively early stage of development, and efforts to directly assimilate radiances instead of retrieved products have just begun a few years ago, partially because it requires much more computation power due to the employment of a radiative transport forward model (FM). This paper will focus on a method to assimilate SBUV/2 radiances (albedos) into the Global Earth Observation System Ozone Data Assimilation Scheme (GEOS-03DAS). While SBUV-type instruments cannot compete with newer sensors in terms of spectral and horizontal resolution, they feature a continuous data record back to 1978, which makes them very valuable for trend studies. Assimilation can help spreading their ground coverage over the whole globe, as has been previously demonstrated with the GEOS-03DAS using SBUV Version 6 ozone profiles. Now, the DAS has been updated to use the newly released SBUV Version 8 data. We will compare pre]lmlnarv results of SBUV radiance assimilation with the assimilation of retrieved ozone profiles, discuss methods to deal with the increased computational load, and try to assess the error characteristics and future potential of the new approach.

Mueller, Martin D.↗

NASA NTRS2004

Torque Balances on the Taylor Cylinders in the Geomagnetic Data Assimilation

In this presentation we report on our continuing effort in geomagnetic data assimilation, aiming at understanding and predicting geomagnetic secular variation on decadal time scales. In particular, we focus on the effect of the torque balances on the cylindrical surfaces in the core co-axial with the Earth's rotation axis (the Taylor cylinders) on the time evolution of assimilated solutions. We use our MoSST core dynamics,model and observed geomagnetic field at the Earth's surface derived via Comprehensive Field Model (CFM) for the geomagnetic data assimilation. In our earlier studies, a model solution is selected randomly from our numerical database. It is then assimilated with the observations such that the poloidal field possesses the same field tomography on the core-mantel boundary (CMB) continued downward from surface observations. This tomography change is assumed to be effective through out the outer core. While this approach allows rapid convergence between model solutions and the observations, it also generates sevee numerical instabilities: the delicate balance between weak fluid inertia and the magnetic torques on the Taylor cylinders are completely altered. Consequently, the assimilated solution diverges quickly (in approximately 10% of the magnetic free-decay time in the core). To improve the assimilation, we propose a partial penetration of the assimilation from the CMB: The full-scale modification at the CMB decreases linearly and vanish at an interior radius r(sub a). We shall examine from our assimilation tests possible relationships between the convergence rate of the model solutions to observations and the cut-off radius r(sub a). A better assimilation shall serve our nudging tests in near future.

Kuang, Weijia↗

NASA NTRS2009

Satellite-Scale Snow Water Equivalent Assimilation into a High-Resolution Land Surface Model

An ensemble Kalman filter (EnKF) is used in a suite of synthetic experiments to assimilate coarse-scale (25 km) snow water equivalent (SWE) observations (typical of satellite retrievals) into fine-scale (1 km) model simulations. Coarse-scale observations are assimilated directly using an observation operator for mapping between the coarse and fine scales or, alternatively, after disaggregation (re-gridding) to the fine-scale model resolution prior to data assimilation. In either case observations are assimilated either simultaneously or independently for each location. Results indicate that assimilating disaggregated fine-scale observations independently (method 1D-F1) is less efficient than assimilating a collection of neighboring disaggregated observations (method 3D-Fm). Direct assimilation of coarse-scale observations is superior to a priori disaggregation. Independent assimilation of individual coarse-scale observations (method 3D-C1) can bring the overall mean analyzed field close to the truth, but does not necessarily improve estimates of the fine-scale structure. There is a clear benefit to simultaneously assimilating multiple coarse-scale observations (method 3D-Cm) even as the entire domain is observed, indicating that underlying spatial error correlations can be exploited to improve SWE estimates. Method 3D-Cm avoids artificial transitions at the coarse observation pixel boundaries and can reduce the RMSE by 60% when compared to the open loop in this study.

De Lannoy, Gabrielle J.M.↗

NASA NTRS2012

Does Ocean Color Data Assimilation Improve Estimates of Global Ocean Inorganic Carbon?

Ocean color data assimilation has been shown to dramatically improve chlorophyll abundances and distributions globally and regionally in the oceans. Chlorophyll is a proxy for phytoplankton biomass (which is explicitly defined in a model), and is related to the inorganic carbon cycle through the interactions of the organic carbon (particulate and dissolved) and through primary production where inorganic carbon is directly taken out of the system. Does ocean color data assimilation, whose effects on estimates of chlorophyll are demonstrable, trickle through the simulated ocean carbon system to produce improved estimates of inorganic carbon? Our emphasis here is dissolved inorganic carbon, pC02, and the air-sea flux. We use a sequential data assimilation method that assimilates chlorophyll directly and indirectly changes nutrient concentrations in a multi-variate approach. The results are decidedly mixed. Dissolved organic carbon estimates from the assimilation model are not meaningfully different from free-run, or unassimilated results, and comparisons with in situ data are similar. pC02 estimates are generally worse after data assimilation, with global estimates diverging 6.4% from in situ data, while free-run estimates are only 4.7% higher. Basin correlations are, however, slightly improved: r increase from 0.78 to 0.79, and slope closer to unity at 0.94 compared to 0.86. In contrast, air-sea flux of C02 is noticeably improved after data assimilation. Global differences decline from -0.635 mol/m2/y (stronger model sink from the atmosphere) to -0.202 mol/m2/y. Basin correlations are slightly improved from r=O.77 to r=0.78, with slope closer to unity (from 0.93 to 0.99). The Equatorial Atlantic appears as a slight sink in the free-run, but is correctly represented as a moderate source in the assimilation model. However, the assimilation model shows the Antarctic to be a source, rather than a modest sink and the North Indian basin is represented incorrectly as a sink rather than the source indicated by the free-run model and data estimates.

Gregg, Watson↗

NASA NTRS2012

Quantifying Interannual Variability of the UTLS Ozone Using Assimilation of Satellite Data

An accurate representation of spatial and temporal variability of the Upper Troposphere Lower Stratosphere (UTLS) ozone is essential for understanding both the tropospheric ozone budget and ozone s contribution to radiative forcing. The complex, dynamically driven structure of trace gas fields in the UTLS presents a challenge to data-based and modelling studies. Small features are not fully resolved in data from limb-sounding instruments such as the Microwave Limb Sounder on EOS-Aura (the EOS-MLS), but are captured in assimilation of those data as vertical structure is added from the assimilated meteorology. This will be demonstrated using a multi-year assimilation of EOS-MLS observations in the Goddard Earth Observing System, Version 5 (GEOS-5) data assimilation system. The results demonstrate the realism of the seasonal and year to year variability of laminar structures in the mid-latitudinal ozone field between years 2005-2007, for which independent validation data are available from the HIRDLS instrument. The analysis is done in the context of the underlying large scale dynamics. The lifetimes of most research instruments are too short for them to be used throughout the duration of long-term (at least 3 decades) reanalyses. For example, the EOS-MLS instrument has operated since mid-2004 until present. By contrast, Solar Backscatter Ultra Violet (SBUV) measurements provide continuous data since late 1978, but their vertical resolution is insufficient to represent the profile shape in the UTLS. Assimilation of these SBUV/2 observations in the GEOS-5 data assimilation system has hitherto not captured a realistic ozone structure in the UTLS, even though transport studies using GEOS-5 wind fields do show such structures. We show that careful construction of the background error covariance structure in GEOS-5 can lead to more realistic UTLS ozone fields when assimilating SBUV/2 observations. The reasoning behind this will be discussed, emphasizing the need to retain the sharp gradient of ozone concentrations across the tropopause. We analyze the UTLS ozone distributions in multi-year SBUV/2 assimilation experiments, comparing the results against the independent HIRDLSdataset and, for a longer period, with the MLS assimilation and discuss the consequences for tropospheric ozone and radiative forcing.

Wargan, K.↗

NASA NTRS2016

Assimilation of Gridded Terrestrial Water Storage Observations from GRACE into a Land Surface Model

Observations of terrestrial water storage (TWS) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission have a coarse resolution in time (monthly) and space (roughly 150,000 km(sup 2) at midlatitudes) and vertically integrate all water storage components over land, including soil moisture and groundwater. Data assimilation can be used to horizontally downscale and vertically partition GRACE-TWS observations. This work proposes a variant of existing ensemble-based GRACE-TWS data assimilation schemes. The new algorithm differs in how the analysis increments are computed and applied. Existing schemes correlate the uncertainty in the modeled monthly TWS estimates with errors in the soil moisture profile state variables at a single instant in the month and then apply the increment either at the end of the month or gradually throughout the month. The proposed new scheme first computes increments for each day of the month and then applies the average of those increments at the beginning of the month. The new scheme therefore better reflects submonthly variations in TWS errors. The new and existing schemes are investigated here using gridded GRACE-TWS observations. The assimilation results are validated at the monthly time scale, using in situ measurements of groundwater depth and soil moisture across the U.S. The new assimilation scheme yields improved (although not in a statistically significant sense) skill metrics for groundwater compared to the open-loop (no assimilation) simulations and compared to the existing assimilation schemes. A smaller impact is seen for surface and root-zone soil moisture, which have a shorter memory and receive smaller increments from TWS assimilation than groundwater. These results motivate future efforts to combine GRACE-TWS observations with observations that are more sensitive to surface soil moisture, such as L-band brightness temperature observations from Soil Moisture Ocean Salinity (SMOS) or Soil Moisture Active Passive (SMAP). Finally, we demonstrate that the scaling parameters that are applied to the GRACE observations prior to assimilation should be consistent with the land surface model that is used within the assimilation system.

GRACE↗

NASA NTRS2017

Assimilation of Global Radar Backscatter and Radiometer Brightness Temperature Observations to Improve Soil Moisture and Land Evaporation Estimates

Active radar backscatter (s◦) observations from the Advanced Scatterometer (ASCAT) and passive radiometer brightness temperature (TB) observations from the Soil Moisture Ocean Salinity (SMOS) mission are assimilated either individually or jointly into the Global Land Evaporation Amsterdam Model (GLEAM) to improve its simulations of soil moisture and land evaporation. To enable s◦ and TB assimilation, GLEAM is coupled to the Water Cloud Model and the L-band Microwave Emission from the Biosphere (L-MEB) model. The innovations, i.e. differences between observations and simulations, are mapped onto the model soil moisture states through an Ensemble Kalman Filter. The validation of surface (0-10 cm) soil moisture simulations over the period 2010-2014 against in situ measurements from the International Soil Moisture Network (ISMN) shows that assimilating s◦ or TB alone improves the average correlation of seasonal anomalies (Ran) from 0.514 to 0.547 and 0.548, respectively. The joint assimilation further improves Ran to 0.559. Associated enhancements in daily evaporative flux simulations by GLEAM are validated based on measurements from 22 FLUXNET stations. Again, the singular assimilation improves Ran from 0.502 to 0.536 and 0.533, respectively for s◦ and TB, whereas the best performance is observed for the joint assimilation (Ran = 0.546). These results demonstrate the complementary value of assimilating radar backscatter observations together with brightness temperatures for improving estimates of hydrological variables, as their joint assimilation outperforms the assimilation of each observation type separately.

GLEAM↗