Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Assimilation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Implementation of Multidomain Unified Forward Operators (UFO) Within the Joint Effort for Data Assimilation Integration (JEDI): Ocean Applications

The Joint Effort for Data assimilation Integration (JEDI) is a collaborative development led by the Joint Center for Satellite Data Assimilation (JCSDA) in conjunction with NASA, NOAA and the Department of Defense (NAVY and Air Force). The (Sea-Ice Ocean and Coupled Assimilation) SOCA as one of the JCSDA projects, focuses on the application of JEDI to marine data assimilation. One of the goals of SOCA is to make use of surface-sensitive radiances to constrain sea-ice and upper ocean fields (e.g., salinity, temperature, sea-ice fraction, sea-ice temperature, etc.). The first elements toward an ocean/atmosphere coupled data assimilation capability within JEDI, with a focus on supporting and developing the assimilation of radiance observations sensitive to the ocean and atmosphere has been implemented. The direct radiance assimilation of surface sensitive microwave radiances focusing on Global Precipitation Measurement (GPM) Imager (GMI) for the SST Constraint and Soil Moisture Active Passive (SMAP) for the Sea Surface Salinity (SSS) has been the main focus. Also, in UFO the capability to calculate the cool skin layer depth and skin temperature has been implemented similar to the GEOS-5. It has been tested with GMI sea surface temperature retrievals. This is important because Satellite and in-situ observations of the Sea-Surface Temperature (SST) show high variability, including a diurnal cycle and very thin, cool skin layer in contact with the atmosphere, and Incorporating a realistic skin SST is essential for atmosphere-ocean coupled data assimilation.

Unified Forward Operators (UFO)↗
NASA NTRS2020

Assimilation of Vegetation Optical Depth Retrievals from Passive Microwave Radiometry

Vegetation optical depth (VOD) retrievals from passive microwave sensors provide analog estimates of above-ground canopy biomass. This study presents the development and analysis of assimilating VOD retrievals from X-, C-, and L-band passive microwave instruments within the Noah-MP land surface model over the Continental U.S. The results from this study demonstrate that the assimilation of VOD retrievals have a significant beneficial impact on the simulation of evapotranspiration and GPP, particularly over the agricultural areas of the U.S. The improvements in the water and carbon fluxes from the assimilation of VOD from X- and C-band sensors are found to be comparable to those obtained from the assimilation of vegetation indices from optical sensors.The study also quantifies the relative and joint impact of assimilating surface soil moisture and VOD from the Soil Moisture Active Passive (SMAP) mission. The utility of soil moisture assimilation for improving ET is more significant over water-limited regions, whereas VOD DA is more impactful over areas where soil moisture is not the primary controlling factor on ET. The results also indicate that the information on moisture and vegetation states from SMAP can be simultaneously exploited through the joint assimilation of surface soil moisture and VOD. Since passive microwave-based VOD retrievals are available in nearly all weather conditions, their use within data assimilation systems offers the ability to extend and improve the utility obtained from the use of optical/infrared based vegetation retrievals

Sujay V Kumar↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

The Transition of Satellite Observations Assimilated in GEOS to JEDI

The Goddard Earth Observing System (GEOS) is used by NASA’s Global Modeling and Assimilation Office (GMAO) to produce weather, climate, and air quality forecasts and reanalysis datasets. In order to incorporate the Joint Effort for Data assimilation Integration (JEDI) system into GEOS it is necessary to provide a validated set of input observations in JEDI. The GMAO, in collaboration with the Joint Center for Satellite Data Assimilation (JCSDA), is developing the Unified Forward Operator (UFO) in JEDI and adding all the necessary features to replicate the existing capability of the Gridpoint Statistical Interpolation (GSI)–based GEOS atmospheric data assimilations system. GMAO has been adding, validating, and updating procedures to assimilate the vast array of satellite and conventional observations assimilated in GEOS, including the GEOS all-sky microwave radiance assimilation framework to assimilate those observations in the UFO. Robust tests have been conducted to ensure correct configurations of observational data bias correction, quality control, and observation error in UFO and good agreements between UFO and GSI results. Our work on satellite observations will be reported in this presentation.

Jianjun Jin↗
NASA NTRS2013

Assessing the Impact of Pre-gpm Microwave Precipitation Observations in the Goddard WRF Ensemble Data Assimilation System

The forthcoming Global Precipitation Measurement (GPM) Mission will provide next generation precipitation observations from a constellation of satellites. Since precipitation by nature has large variability and low predictability at cloud-resolving scales, the impact of precipitation data on the skills of mesoscale numerical weather prediction (NWP) is largely affected by the characterization of background and observation errors and the representation of nonlinear cloud/precipitation physics in an NWP data assimilation system. We present a data impact study on the assimilation of precipitation-affected microwave (MW) radiances from a pre-GPM satellite constellation using the Goddard WRF Ensemble Data Assimilation System (Goddard WRF-EDAS). A series of assimilation experiments are carried out in a Weather Research Forecast (WRF) model domain of 9 km resolution in western Europe. Sensitivities to observation error specifications, background error covariance estimated from ensemble forecasts with different ensemble sizes, and MW channel selections are examined through single-observation assimilation experiments. An empirical bias correction for precipitation-affected MW radiances is developed based on the statistics of radiance innovations in rainy areas. The data impact is assessed by full data assimilation cycling experiments for a storm event that occurred in France in September 2010. Results show that the assimilation of MW precipitation observations from a satellite constellation mimicking GPM has a positive impact on the accumulated rain forecasts verified with surface radar rain estimates. The case-study on a convective storm also reveals that the accuracy of ensemble-based background error covariance is limited by sampling errors and model errors such as precipitation displacement and unresolved convective scale instability.

microwave observations↗
NASA NTRS2017

Role of Forcing Uncertainty and Background Model Error Characterization in Snow Data Assimilation

Accurate specification of the model error covariances in data assimilation systems is a challenging issue. Ensemble land data assimilation methods rely on stochastic perturbations of input forcing and model prognostic fields for developing representations of input model error covariances. This article examines the limitations of using a single forcing dataset for specifying forcing uncertainty inputs for assimilating snow depth retrievals. Using an idealized data assimilation experiment, the article demonstrates that the use of hybrid forcing input strategies (either through the use of an ensemble of forcing products or through the added use of the forcing climatology) provide a better characterization of the background model error, which leads to improved data assimilation results, especially during the snow accumulation and melt-time periods. The use of hybrid forcing ensembles is then employed for assimilating snow depth retrievals from the AMSR2 (Advanced Microwave Scanning Radiometer 2) instrument over two domains in the continental USA with different snow evolution characteristics. Over a region near the Great Lakes, where the snow evolution tends to be ephemeral, the use of hybrid forcing ensembles provides significant improvements relative to the use of a single forcing dataset. Over the Colorado headwaters characterized by large snow accumulation, the impact of using the forcing ensemble is less prominent and is largely limited to the snow transition time periods. The results of the article demonstrate that improving the background model error through the use of a forcing ensemble enables the assimilation system to better incorporate the observational information.

assimilation↗
NASA NTRS2022

The Joint Assimilation of Remotely Sensed Leaf Area Index and Surface Soil Moisture into a Land Surface Model

This work tests the hypothesis that jointly assimilating satellite observations of leaf area index and surface soil moisture into a land surface model improves the estimation of land vegetation and water variables. An Ensemble Kalman Filter is used to test such hypothesis across the Contiguous United States during April 2015 – December 2018. The performance of the proposed methodology is assessed for several modeled vegetation and water variables (evapotranspiration, net ecosystem exchange, and soil moisture) in terms of random errors and anomaly correlation coefficients against a set of independent validation datasets (i.e., Global Land Evaporation Amsterdam Model, FLUXCOM, and International Soil Moisture Network). Results show that the assimilation of leaf area index mostly improves the estimation of evapotranspiration and net ecosystem exchange, whereas the assimilation of surface soil moisture alone improves surface soil moisture content, especially in the western US, in terms of both root mean squared error and anomaly correlation coefficient. The joint assimilation of vegetation and soil moisture information combines the results of individual vegetation and soil moisture assimilations and reduces errors (and increases correlations with the reference datasets) in evapotranspiration, net ecosystem exchange, and surface soil moisture simulated by the land surface model. However, because soil moisture satellite observations only provide information on the water content in the top 5 cm of the soil column, the impact of the proposed data assimilation technique on root zone soil moisture is limited. This work moves one step forward in the direction of improving our estimation and understanding of land surface interactions using a multi-variate data assimilation approach, which can be particularly useful in regions of the world where ground observations are sparse or missing altogether.

Data assimilation↗
NASA NTRS2024

Influence of SMAP Soil Moisture Retrieval Assimilation on Runoff Estimation Across South Asia

This study was designed to characterize and quantify the influence of surface soil moisture assimilation on estimated runoff (surface flow and baseflow) and hydraulically-routed streamflow across three large river basins in South Asia that are at risk of impending water stress. Soil Moisture Active Passive (SMAP) surface soil moisture retrievals were assimilated into the Noah-MP land surface model. The gridded runoff was hydraulically routed to obtain volumetric streamflow values for a river network using a runoff routing module (Hydrological Modeling and Analysis Platform, HyMAP). The open loop (OL, model-only) and data assimilated (DA, includes soil moisture retrieval assimilation) Noah-MP runoff estimates highlighted the improvements in estimated total runoff across irrigated areas. Soil moisture assimilation impacted baseflow more relative to surface runoff. The HyMAP-based OL and DA streamflow generally underestimated the streamflow at upstream stations and overestimated the streamflow at downstream stations within the Indus basin due to missing physics related to reservoir operations. For stations located in the Ganges–Brahmaputra basins, the OL and DA estimation performance varied. The OL and DA results showed that the assimilation of soil moisture retrievals improves gridded runoff and volumetric streamflow across irrigated areas. Considerable relative change in streamflow (>70% increase in magnitude relative to OL) is noted after assimilation across the highly irrigated lower Indus basin and high precipitation regions in Bangladesh. Improving the modeling system’s representativeness of ground conditions via inclusion of water management information could improve large-scale streamflow modeling across South Asia.

soil moisture↗
NASA NTRS1985

Variational mesoscale satellite data assimilation and initialization

The problems of mesoscale satellite data assimilation were examined. Assimilation of satellite data to improve the forecasts made by mesoscale forecast models was undertaken. Assimilation of high resolution satellite derived temperature data into a mesoscale model with horizontal resolution of 50 to 60 km is reported. Unlike global assimilation, in which a small portion of the forecast model domain is subject to data insertion at virtually every time step, the mesoscale assimilation virtually all of the forecast model domain is subject to data insertion at one time step. The mesoscale problem lends itself naturally to intermittent data assimilation and the forecast model is reinitialized whenever a new satellite pass covers its domain with data. The satellite data assimilation as an initialization problem are discussed.

Sasaki, Y. K.↗
NASA NTRS1998

Stratospheric Assimilation of Chemical Tracer Observations Using a Kalman Filter: Chi-Square Validated Results and Analysis of Variance and Correlation Dynamics - Pt. 2

A Kalman filter system designed for the assimilation of limb-sounding observations of stratospheric chemical tracers, which has four tunable covariance parameters, was developed in Part I (Menard et al. 1998) The assimilation results of CH4 observations from the Cryogenic Limb Array Etalon Sounder instrument (CLAES) and the Halogen Observation Experiment instrument (HALOE) on board of the Upper Atmosphere Research Satellite are described in this paper. A robust (chi)(sup 2) criterion, which provides a statistical validation of the forecast and observational error covariances, was used to estimate the tunable variance parameters of the system. In particular, an estimate of the model error variance was obtained. The effect of model error on the forecast error variance became critical after only three days of assimilation of CLAES observations, although it took 14 days of forecast to double the initial error variance. We further found that the model error due to numerical discretization as arising in the standard Kalman filter algorithm, is comparable in size to the physical model error due to wind and transport modeling errors together. Separate assimilations of CLAES and HALOE observations were compared to validate the state estimate away from the observed locations. A wave-breaking event that took place several thousands of kilometers away from the HALOE observation locations was well captured by the Kalman filter due to highly anisotropic forecast error correlations. The forecast error correlation in the assimilation of the CLAES observations was found to have a structure similar to that in pure forecast mode except for smaller length scales. Finally, we have conducted an analysis of the variance and correlation dynamics to determine their relative importance in chemical tracer assimilation problems. Results show that the optimality of a tracer assimilation system depends, for the most part, on having flow-dependent error correlation rather than on evolving the error variance.

Menard, Richard↗
NASA NTRS1998

Comparison of Surface Ground Temperature from Satellite Observations and the Off-Line Land Surface GEOS Assimilation System

The surface ground temperature (Tg) is an important meteorological variable, because it represents an integrated thermal state of the land surface determined by a complex surface energy budget. Furthermore, Tg affects both the surface sensible and latent heat fluxes. Through these fluxes. the surface budget is coupled with the atmosphere above. Accurate Tg data are useful for estimating the surface radiation budget and fluxes, as well as soil moisture. Tg is not included in conventional synoptical weather station reports. Currently, satellites provide Tg estimates globally. It is necessary to carefully consider appropriate methods of using these satellite data in a data assimilation system. Recently, an Off-line Land surface GEOS Assimilation (OLGA) system was implemented at the Data Assimilation Office at NASA-GSFC. One of the goals of OLGA is to assimilate satellite-derived Tg data. Prior to the Tg assimilation, a thorough investigation of satellite- and model-derived Tg, including error estimates, is required. In this study we examine the Tg from the n Project (ISCCP DI) data and the OLGA simulations. The ISCCP data used here are 3-hourly DI data (2.5x2.5 degree resolution) for 1992 summer months (June, July, and August) and winter months (January and February). The model Tg for the same periods were generated by OLGA. The forcing data for this OLGA 1992 simulation were generated from the GEOS-1 Data Assimilation System (DAS) at Data Assimilation Office NASA-GSFC. We examine the discrepancies between ISCCP and OLGA Tg with a focus on its spatial and temporal characteristics, particularly on the diurnal cycle. The error statistics in both data sets, including bias, will be estimated. The impact of surface properties, including vegetation cover and type, topography, etc, on the discrepancies will be addressed.

Yang, R.↗
NASA NTRS1999

1D Variational Assimilation of TOVS/ATOVS Level 1b Data

Recently, it has been shown that the use of observations from satellite-borne microwave and infrared radiometers in data assimilation systems consistently increases forecast skill. Considerable effort has been expended over the past two decades, particularly with the TIROS Operational Vertical Sounder (TOVS), to achieve this result. The positive impact on forecast skill has resulted from improvements in quality control algorithms, systematic error correction schemes, and data assimilation systems. Despite these recent advances, there are still many issues regarding the use of satellite data in data assimilation systems that remain unresolved. For example, at several centers, much of the TOVS data are not assimilated over land as a result of the difficulty in specifying the surface emissivity. In this study, we add the microwave surface emissivity to the state vector so that it is retrieved rather than specified by a model. The LDVAR temperature and humidity information is then assimilated interactively into the Goddard Earth Observing System - Data Assimilation System (GEOS-DAS). Results of 1DVAR assimilations with both TOVS and ATOVS will be shown. We will compare the results obtained with HIRS/MSU/SSU and HIRS/AMSU-A2.

Joiner, J.↗
NASA NTRS2001

Assimilation of GOES Land Surface Data into a Mesoscale Models

A technique has been developed for assimilating Geostationary Operational Environmental Satellite (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 assimilation technique has been applied to the Oklahoma-Kansas region during the spring-summer 2000 time period when dynamic changes in vegetation cover occur. In April, central Oklahoma is characterized by large NDVI associated with winter wheat while surrounding areas are primarily rangeland with lower NDVI. In July the vegetation pattern reverses as the central wheat area changes to low NDVI due to harvesting and the surrounding rangeland is greener than it was in April. The goal of this study is to determine if assimilating satellite land surface data can improve simulation of the complex spatial distribution of surface energy and water fluxes across this region. The PSU/NCAR NM5 V3 system is used in this study. The grid configuration consists of a 36-km CONUS domain and a 12-km nest over the area of interest. Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint indicates that assimilation of the satellite data results reduces both the bias and RMSE for both state variables. In addition, comparison of model data with ARM/CART EBBR flux observations reveals that the assimilation technique adjusts the bowen ratio in a realistic fashion.

Lapenta, William M.↗
NASA NTRS2002

Land Surface Data Assimilation and the Gulf Coast Sea-Breeze

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 focus of this paper is to examine how the satellite assimilation technique impacts simulations of near-surface meteorology on the 4- to 48-hour time scale. The technique is implemented within the PSU/NCAR MM5 V34 and applied on a 36-km CONUS domain and a 12-km nest centered over Oklahoma. Two Friday periods during 2000 were selected for study, one during the warm season (May-June) and the other in the cold season (Nov-Dec). Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint show that the assimilation technique can improve numerical simulations on both regional and continental scales. Comparison of the simulated surface energy and water fluxes with observations from Energy Balance Bowen Ratio measurements taken at CART/ARM sites in the Southern Great Plains showed the assimilation technique produces a realistic argument of the Bowen ratio in both seasons. Additional runs with the Oregon State University land surface scheme as implemented within MM5 will be made and comparisons with assimilation results will be presented at the meeting.

Lapenta, William M.↗
NASA NTRS2002

Sensitivity of Lower Stratospheric Assimilated Ozone on Error Covariance Modeling and Data Selection

Assimilated ozone is produced at the NASA/Goddard Data Assimilation Office by blending 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 transport model. The current system tends to overestimate the amount of lower stratospheric ozone. This is a region where ozone 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. Since these are all important applications of assimilated ozone products, effort is being directed towards reducing this bias. The SBUV 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 in the lower stratosphere, and its vertical distribution in particular, 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. Accordingly, the sensitivity of the assimilated lower stratospheric ozone fields to changes in the TOMS error-covariance modeling and the SBUV data selection has been investigated. The use of a spatially correlated TOMS error covariance model led to improvements in the product. However, withholding the SBUV/2 data for the layer between 63 and 126 hPa typically degraded the product, a result which vindicates the use of this layer ozone product, despite its known errors. These efforts to improve the lower stratospheric distribution will be extended to include a more advanced forecast error covariance model, and by assimilating ozone products from new instruments on Envisat and EOS Aura.

Stajner, Ivanka↗
NASA NTRS2003

Comparison of Flow-Dependent and Static Error Correlation Models in the DAO Ozone Data Assimilation System

In a data assimilation system the forecast error covariance matrix governs the way in which the data information is spread throughout the model grid. Implementation of a correct method of assigning covariances is expected to have an impact on the analysis results. The simplest models assume that correlations are constant in time and isotropic or nearly isotropic. In such models the analysis depends on the dynamics only through assumed error standard deviations. In applications to atmospheric tracer data assimilation this may lead to inaccuracies, especially in regions with strong wind shears or high gradient of potential vorticity, as well as in areas where no data are available. In order to overcome this problem we have developed a flow-dependent covariance model that is based on short term evolution of error correlations. The presentation compares performance of a static and a flow-dependent model applied to a global three- dimensional ozone data assimilation system developed at NASA s Data Assimilation Office. We will present some results of validation against WMO balloon-borne sondes and the Polar Ozone and Aerosol Measurement (POAM) III instrument. Experiments show that allowing forecast error correlations to evolve with the flow results in positive impact on assimilated ozone within the regions where data were not assimilated, particularly at high latitudes in both hemispheres and in the troposphere. We will also discuss statistical characteristics of both models; in particular we will argue that including evolution of error correlations leads to stronger internal consistency of a data assimilation ,

Wargan, K.↗
NASA NTRS2003

Impact of Flow-Dependent Error Correlations and Tropospheric Chemistry on Assimilated Ozone

The presentation compares different versions of a global three-dimensional ozone data assimilation system developed at NASA's Data Assimilation Office. The Solar Backscatter Ultraviolet/2 (SBUV/2) total and partial ozone column retrievals are the sole data assimilated in all of the experiments presented. We study the impact of changing the forecast error covariance model from a version assuming static correlations with a one that captures a short-term Lagrangian evolution of those correlations. This is further combined with a study of the impact of neglecting the tropospheric ozone production, loss and dry deposition rates, which are obtained from the Harvard GEOS-CHEM model. We compare statistical characteristics of the assimilated data and the results of validation against independent observations, obtained from WMO balloon-borne sondes and the Polar Ozone and Aerosol Measurement (POAM) III instrument. Experiments show that allowing forecast error correlations to evolve with the flow results in positive impact on assimilated ozone within the regions where data were not assimilated, particularly at high latitudes in both hemispheres. On the other hand, the main sensitivity to tropospheric chemistry is in the Tropics and sub-Tropics. The best agreement between the assimilated ozone and the in-situ sonde data is in the experiment using both flow-dependent error covariances and tropospheric chemistry.

Wargan, K.↗
NASA NTRS2005

An Improved Analysis of Forest Carbon Dynamics using Data Assimilation

There are two broad approaches to quantifying landscape C dynamics - by measuring changes in C stocks over time, or by measuring fluxes of C directly. However, these data may be patchy, and have gaps or biases. An alternative approach to generating C budgets has been to use process-based models, constructed to simulate the key processes involved in C exchange. However, the process of model building is arguably subjective, and parameters may be poorly defined. This paper demonstrates why data assimilation (DA) techniques - which combine stock and flux observations with a dynamic model - improve estimates of, and provide insights into, ecosystem carbon (C) exchanges. We use an ensemble Kalman filter (EnKF) to link a series of measurements with a simple box model of C transformations. Measurements were collected at a young ponderosa pine stand in central Oregon over a 3-year period, and include eddy flux and soil C02 efflux data, litterfall collections, stem surveys, root and soil cores, and leaf area index data. The simple C model is a mass balance model with nine unknown parameters, tracking changes in C storage among five pools; foliar, wood and fine root pools in vegetation, and also fresh litter and soil organic matter (SOM) plus coarse woody debris pools. We nested the EnKF within an optimization routine to generate estimates from the data of the unknown parameters and the five initial conditions for the pools. The efficacy of the DA process can be judged by comparing the probability distributions of estimates produced with the EnKF analysis vs. those produced with reduced data or model alone. Using the model alone, estimated net ecosystem exchange of C (NEE)= -251 f 197g Cm-2 over the 3 years, compared with an estimate of -419 f 29gCm-2 when all observations were assimilated into the model. The uncertainty on daily measurements of NEE via eddy fluxes was estimated at 0.5gCm-2 day-1, but the uncertainty on assimilated estimates averaged 0.47 g Cm-2 day-1, and only exceeded 0.5gC m-2 day-1 on days where neither eddy flux nor soil efflux data were available. In generating C budgets, the assimilation process reduced the uncertainties associated with using data or model alone and the forecasts of NEE were statistically unbiased estimates. The results of the analysis emphasize the importance of time series as constraints. Occasional, rare measurements of stocks have limited use in constraining the estimates of other components of the C cycle. Long time series are particularly crucial for improving the analysis of pools with long time constants, such as SOM, woody biomass, and woody debris. Long-running forest stem surveys, and tree ring data, offer a rich resource that could be assimilated to provide an important constraint on C cycling of slow pools. For extending estimates of NEE across regions, DA can play a further important role, by assimilating remote-sensing data into the analysis of C cycles. We show, via sensitivity analysis, how assimilating an estimate of photosynthesis - which might be provided indirectly by remotely sensed data - improves the analysis of NEE.

Williams, Mathew↗