Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASSIMILATION”

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

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

At least 1,621 records · Page 90

An Introduction to the NASA GMAO Coupled Atmosphere-Ocean System - GEOS-S2S Version 3

Recently NASA's Global Modeling and Assimilation Office (GMAO) has developed a new Subseasonal to Seasonal Prediction system Version 3 (GEOS-S2S-3). This upgrade replaces the GEOS-S2S-2 which is NASA's current contribution to the North American Multi-Model Experiment seasonal prediction project (Kirtman et al., 2014). The main improvements for our S2S-3 system include 1) a higher resolution MOM5 (Griffies et al., 2005) ocean model (now 0.25o x 0.25o x 50 layers), 2) an improved atmospheric/ocean interface layer (Akella and Suarez, 2018), and 3) assimilation of a long-track satellite salinity into the ocean model (Hackert et al, 2019). Atmospheric forcing is provided by the NASA MERRA-2 reanalysis (Gelaro et al., 2017). Initialization for the ocean relies on the GMAO ocean reanalysis system which assimilates all available in situ temperature and salinity, satellite sea surface salinity, and sea level using the Local Ensemble Transform Kalman Filter (LETKF) implementation of (Penny et al., 2013) on a 5 day assimilation cycle with 20 fixed ensemble members.In this presentation, we will authenticate our new S2S-3 ocean reanalysis using standard GODAE validation metrics. For example, we will compare gridded fields of mean and standard deviation of the ocean reanalysis versus observed fields. We will show correlation/RMS of model versus observations and temperature and salinity mean profiles for the various basins and latitude bands. Basin-scale volume transports, such as the Atlantic Meridional Overturning Circulation and the Indonesian Throughflow will be validated. Equatorial ocean waves will be compared by decomposing sea level into Kelvin and Rossby components. For each of these metrics, we plan to validate the results and then compare our new S2S-3 against the current production version, S2S-2. Finally, we will compare 9-month seasonal forecasts initialized from these two systems for the tropical Pacific NINO3.4 region over the period 1981-present.

Molod, Andrea↗

Watershed Modeling with Remotely Sensed Big Data: MODIS Leaf Area Index Improves Hydrology and Water Quality Predictions

Traditional watershed modeling often overlooks the role of vegetation dynamics. There is also little quantitative evidence to suggest that increased physical realism of vegetation dynamics in process-based models improves hydrology and water quality predictions simultaneously. In this study, we applied a modified Soil and Water Assessment Tool (SWAT) to quantify the extent of improvements that the assimilation of remotely sensed Leaf Area Index (LAI) would convey to streamflow, soil moisture, and nitrate load simulations across a 16,860 km2 agricultural watershedin the midwestern United States. We modified the SWAT source code to automatically override the model’s built-in semiempirical LAI with spatially distributed and temporally continuous estimates from Moderate Resolution Imaging Spectroradiometer (MODIS). Compared to a “basic” traditional model with limited spatial information, our LAI assimilation model (i) significantly improved daily streamflow simulations during medium-to-low flow conditions, (ii) provided realistic spatial distributions of growing season soil moisture, and (iii) substantially reproduced the long-term observed variability of daily nitrate loads. Further analysis revealed that the overestimation or underestimation of LAI imparted a proportional cascading effect on how the model partitions hydrologic fluxes and nutrient pools. As such, assimilation of MODIS LAI data corrected the model’sLAI overestimation tendency, which led to a proportionally increased rootzone soil moisture and decreased plant nitrogen uptake. With these new findings, our study fills the existing knowledge gap regarding vegetation dynamics in watershed modeling and confirms that assimilation of MODIS LAI data in watershed models can effectively improve both hydrology and water quality predictions.

Adnan Rajib↗

Observing System Simulation Experiments Investigating Atmospheric Motion Vectors and Radiances from a Constellation of 4-5 μm Infrared Sounders

A set of Observing System Simulation Experiments (OSSEs) was performed to investigate the utility of a constellation of passive infrared spectrometers, strategically designed with the aim of deriving the three-dimensional retrievals of the horizontal wind via atmospheric motion vectors (AMVs) from instruments with the spectral resolution of an infrared sounder. The instrument and constellation designs were performed in the context of the Midwave Infrared Sounding of Temperature and humidity in a Constellation for Winds, or MISTiC Winds. The Global Modeling and Assimilation Office OSSE system, which includes a full suite of operational meteorological observations, served as the control. To illustrate the potential impact of this observing strategy, two experiments were performed by adding the new simulated observations to the control. First, perfect (error-free) simulated AMVs and radiances were assimilated. Second, the data were made imperfect by adding realistic modeled errors to the AMVs and radiances that were assimilated. The experimentation showed beneficial impacts on both the mass and wind fields, as based on analysis verification, forecast verification, and the assessment of the observations using the Forecast Sensitivity to Observation Impact (FSOI) metric. In all variables and metrics, the impacts of the imperfect observations were smaller than those of the perfect observations, though much of the positive benefit was retained. The FSOI metric illustrated two key points. First, the largest impacts were seen in the middle troposphere AMVs, which is a targeted capability of the constellation strategy. Second, the addition of modeled errors showed that the assimilation system was unable to fully exploit the 4.3 μm carbon dioxide absorption radiances.

Will McCarty↗

The role of analysis error in the convergence of reanalysis production streams in MERRA-2

Due to production time constraints, most reanalyses are produced in multiple parallel streams instead of a single continuous one. These streams cover separate segments of the reanalysis time period with short overlaps to allow reconstruction of the official record. A fundamental assumption justifying this approach is that the streams will be assimilating the same observations during the periods where they overlap, and so will eventually converge to a similar atmospheric state, making discontinuities at stream junctions negligible. This assumption is revisited in this work by examining the impact of analysis error on the differences between MERRA-2 overlapping streams in three historical periods. Comparison results are shown in terms of standard deviations of stream differences as well as the spectral decomposition of the variance of their differences. Residual differences were found at the end of each year of overlap, with larger values observed in the earlier segments of the presatellite era. By drawing parallels with analysis error statistics estimated from the GMAO OSSE system, these differences are shown to reflect the varying constraint of data with the varying observing network, and to further carry the imprint of errors that the data assimilation process is not able to mitigate. As such, they are unlikely to be reduced by longer spinup periods. The ability of data assimilation to ensure continuity in the parallel streams is put into question when the observing system coverage is inadequate or simply when the data assimilation system as a whole is suboptimal.

Amal El Akkraoui↗

Near real-time air quality forecasts using the NASA GEOS model

NASA's Global Modeling and Assimilation Office (GMAO) produces high-resolution global forecasts for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model has been expanded to provide global near-real-time5-day forecasts of atmospheric composition at unprecedented horizontal resolution of 0.25 degrees (~25 km). This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module (version 12) to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). Satellite observations are assimilated into the system for improved representation of weather and smoke. The assimilation system is being expanded to include chemically reactive trace gases. We discuss current capabilities of the GEOS Constituent Data Assimilation System (CoDAS) to improve atmospheric composition modeling and possible future directions, notably incorporating new observations (TROPOMI, geostationary satellites) and machine learning techniques. We show how machine learning techniques can be used to correct for sub-grid-scale variability, which further improves model estimates at a given observation site.

Air Quality Forecast↗

Observing System Simulation Experiments Investigating Atmospheric Motion Vectors and Radiances from a Constellation of 4–5-μm Infrared Sounders

A set of observing system simulation experiments (OSSEs) was performed to investigate the utility of a constellation of passive infrared spectrometers, strategically designed with the aim of deriving the three-dimensional retrievals of the horizontal wind via atmospheric motion vectors (AMVs) from instruments with the spectral resolution of an infrared sounder. The instrument and constellation designs were performed in the context of the Midwave Infrared Sounding of Temperature and humidity in a Constellation for Winds (MISTiC Winds). The Global Modeling and Assimilation Office OSSE system, which includes a full suite of operational meteorological observations, served as the control. To illustrate the potential impact of this observing strategy, two experiments were performed by adding the new simulated observations to the control. First, perfect (error free) simulated AMVs and radiances were assimilated. Second, the data were made imperfect by adding realistic modeled errors to the AMVs and radiances that were assimilated. The experimentation showed beneficial impacts on both the mass and wind fields, as based on analysis verification, forecast verification, and the assessment of the observations using the forecast sensitivity to observation impact (FSOI) metric. In all variables and metrics, the impacts of the imperfect observations were smaller than those of the perfect observations, although much of the positive benefit was retained. The FSOI metric illustrated two key points. First, the largest impacts were seen in the middle troposphere AMVs, which is a targeted capability of the constellation strategy. Second, the addition of modeled errors showed that the assimilation system was unable to fully exploit the 4.3-μm carbon dioxide absorption radiances.

Will McCarty↗

Underlying Fundamentals of Kalman Filtering for River Network Modeling

The grand challenge of producing hydrometeorological estimates every time and everywhere has motivated the fusion of sparse observations with dense numerical models, with a particular interest on discharge in river modeling. Ensemble methods are largely preferred as they enable the estimation of error properties, but at the expense of computational load and generally with underestimations. These imperfect stochastic estimates motivate the use of correction methods, that is, error localization and inflation, although the physical justifications for their optimality are limited. The purpose of this study is to use one of the simplest forms of data assimilation when applied to river modeling and reveal the underlying mechanisms impacting its performance. Our framework based on assimilating daily averaged in situ discharge measurements to correct daily averaged runoff was tested over a 4-yr case study of two rivers in Texas. Results show that under optimal conditions of inflation and localization, discharge simulations are consistently improved such that the mean values of Nash–Sutcliffe efficiency are enhanced from211.32 to 0.55 at observed gauges and from212.24 to21.10 at validation gauges. Yet, parameters controlling the inflation and the localization have a large impact on the performance. Further investigations of these sensitivities showed that optimal inflation occurs when compensating exactly for discrepancies in the magnitude of errors while optimal localization matches the distance traveled during one assimilation window. These results may be applicable to more advanced data assimilation methods as well as for larger applications motivated by upcoming river-observing satellite missions, such as NASA’s Surface Water and Ocean Topography mission.

Streamflow↗

Global analysis and forecasts of carbon monoxide on Mars

Spatial and temporal variations in the Martian carbon monoxide (CO) cycle have been investigated through combining Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) retrievals of carbon monoxide mixing ratio and Mars Climate Sounder (MCS) temperature profiles with a Martian global circulation model (GCM) to produce the first global reanalysis of the Martian CO cycle. The reanalysis reduces the root mean square error between the forecast and CRISM CO retrievals by a factor of 2–4, dependent on the time of year. Forecasts initiated from the reanalysis show an improved match to standalone CRISM CO retrievals from northern winter solstice to northern summer solstice, indicating the benefit of data assimilation in simulating the CO abundance. The northern summer solstice CO minimum between 10°S – 50°S in the CRISM CO retrievals is found to be caused by a suppression of CO-enriched air leaking from the Hellas and Argyre basins, and is also strongly influenced by the amount of carbon dioxide condensing at the time. The GCM is the first reported to simulate the local CO minimum however it is weaker in strength, as CO-enriched air is also released northward in the area of Argyre basin in the GCM because of a shifted boundary between the southern polar vortex and mid-latitudes. The reanalysis, as a result of the assimilation of MCS temperature profiles, indicates that the polar vortex boundary is northward of Argyre basin and hence no northerly transport of CO-enriched air should be present over this location. Differences in local CO abundance between the reanalysis and GCM are also evident in spatio-temporal regions where no nearby CRISM CO retrievals are available as a result of alterations in local circulation and the sublimation/condensation of carbon dioxide through the combined assimilation of CRISM CO retrievals alongside MCS temperature profiles, illustrating how constraints can be imposed indirectly on the CO cycle through the powerful technique of data assimilation.

James A. Holmes↗

Improving the Use of Surface-Sensitive Radiances in the GMAO Hybrid-4DEnVar System

The planetary boundary layer (PBL) was designated as an incubation-class targeted observable (TO) in the 2018 Decadal Survey. As no single instrument source will provide enough information to constrain the global PBL, the assimilation of a wide range of observations in data assimilation (DA) systems will play a critical role, and GMAO has put in efforts to enhance surface-sensitive radiance assimilation. Although a vast number of microwave radiance observations are used in the GEOS global DA system, very few surface-sensitive radiances are currently used over land due to large uncertainty of land surface emissivity in the CRTM as well as cloud detection issue. Dynamically varying emissivity is retrieved from observations of window channels in the GEOS system for non-scattering FOVs and applied to sounding channel assimilation. Moreover, the original radiance bias correction is found not to work well over land, and shows drifting bias when the original emissivity sensitivity bias predictor with a dynamically varying emissivity retrieval is used. Hence, the radiance bias correction has been modified, and the quality control procedure has been adapted accordingly. Cycled experiment results show improvement in the temperature forecast at lower model levels; further refinement of this system continues.

Yanqiu Zhu↗

Quantifying the Observational Requirements of a Space-borne LiDAR Snow Mission

This study quantifies the level of observational accuracy required from a spaceborne light detection and ranging (LiDAR) snow depth retrieval mission for enabling beneficial impacts for snow estimation. The study is conducted over a region in Western Colorado using a suite of observing system simulation experiments(OSSEs).The Joint UK Land Environment Simulator, version 5.0 (JULES v5.0) is employed to simulate a suite of idealized LiDAR observations, considering a range of LiDAR snow depth retrieval errors, different hypothetical sensor swath widths, and the impact of cloud cover on observability. These simulated observations are then assimilated into the Noah land surface modelwith multi-parameterization options, version 3.6 (Noah-MP v3.6) model. This data assimilation setup is used to systematically evaluate the potential utility of LiDAR observations for improving modeled snow water equivalent (SWE)estimates and water budget variables such as runoff. Results from the OSSE runsshow that, in general, assimilation of synthetic LiDAR observations provide beneficial impacts when theLiDAR snow depth retrieval error standard deviation (σerror) is below 60 cm.Based on comparisons between the realistic (i.e., swath-limited and cloud-attenuated) case and the idealized (i.e., infinite swath width in the absence of cloud cover) case,this study concludes that observations with a conservative error standard deviation threshold of 40 cm (i.e., upper limit of the snow depth retrieval error that adds value to the SWE estimates viaassimilation) are needed for improving modeled snow estimates. More than a 33% reduction in SWEroot mean square errors and more than a 15% increase in correlation coefficientsare achieved when σ error ≤ 40 cmusing a 170-km sensor swath width in the presence of cloud attenuation effects. Further, the integrated hydrologic response, as represented by total (surface and subsurface) runoff estimates during the snow ablation season, are also enhanced when assimilating synthetic LiDAR snow depth retrievals with errors below this level.

Yonghwan Kwon↗

An Update on the Hydrological Land Surface Data and Services at NASA GES DISC

"The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) is one of twelve NASA Earth Observing System data centers that document, process, archive and distribute data from Earth science missions and related projects. The GES DISC hosts many hydrological land surface data products including North American Land Data Assimilation System (NLDAS), Global Land Data Assimilation System (GLDAS), the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS), and assimilated-Gravity Recovery and Climate Experiment (GRACE). The Land Surface Model (LSM) products contain model output, including heat fluxes, rain, snow, soil temperature, soil moisture and runoff; and observational forcing data, including surface pressure, temperature, precipitation, downward shortwave and longwave radiation, humidity, and wind. The temporal resolution of the hydrology data at the GES DISC ranges from hourly to monthly and the spatial resolution 0.1 degree to 1.0 degree. The variability of Earth’s water and energy cycles can be characterized with these spatially and temporally consistent, and quality controlled hydrological land surface data. These data could support ongoing modeling activities, and help improve our understanding of land-surface-atmosphere interactions and their impact on climate. Over the years, new data is added, and older data is reprocessed and updated. This presentation provides a summary table of the hydrological data and discusses recent updates to the data and data services at the GES DISC. New data products include (1) GRACE Global Version 3.0 drought indicators and (2) FLDAS monthly global and Central Asia daily data. Reprocessed and updated products include (1) NLDAS Version 2.0 and (2) GRACE United States Version 4.0 drought indicators. Lastly, this presentation will discuss the opportunities for accessing some of our hydrology data from the cloud. As the hydrology land surface data is moved to the cloud, users can expect faster service for longer time series, notably the hydrology data rods time series service."

Ashley Heath↗

The GEOS Neural Network Retrieval (NNR) for Multi-spectral AOD

One of the difficulties in data assimilation is the need for multi-sensor data merging that can account for temporal and spatial biases between satellite sensors. In the Goddard Earth Observing System Model Version 5 (GEOS-5) aerosol data assimilation system, a neural network retrieval (NNR) is used as a mapping between satellite observed top of the atmosphere (TOA) reflectance and AOD, which is the target variable that is assimilated in the model. By training observations of TOA reflectance from multiple sensors to map to a common AOD dataset (in this case AOD observed by the ground based Aerosol Robotic Network, AERONET), we are able to create a global, homogenous, satellite data record of AOD from multiple sensors. In this presentation, I will present recent updates to the GEOS-5 NNR for estimation of spectral AOD from MODIS and VIIRS, and the potential for multi-channel AOD assimilation to provide constraints on aerosol composition.

GEOS↗

The GEOS Neural Network Retrieval (NNR) for Multi-spectral AOD

One of the difficulties in data assimilation is the need for multi-sensor data merging that can account for temporal and spatial biases between satellite sensors. In the Goddard Earth Observing System Model Version 5 (GEOS-5) aerosol data assimilation system, a neural network retrieval (NNR) is used as a mapping between satellite observed top of the atmosphere (TOA) reflectance and AOD, which is the target variable that is assimilated in the model. By training observations of TOA reflectance from multiple sensors to map to a common AOD dataset (in this case AOD observed by the ground based Aerosol Robotic Network, AERONET), we are able to create a global, homogenous, satellite data record of AOD from multiple sensors. In this presentation, I will present recent updates to the GEOS-5 NNR for estimation of spectral AOD from MODIS and VIIRS, and the potential for multi-channel AOD assimilation to provide constraints on aerosol composition.

Patricia Castellanos↗

Evaluation of TAO Observation System on ENSO Predictions from the GMAO S2S Forecast System

Since the early 1990’s, the Tropical Atmosphere Ocean (TAO) mooring array has been instrumental for observing Kelvin and Rossby wave propagation which is key for El Nino/Southern Oscillation (ENSO) prediction. However, recent funding and programmatical issues have highlighted the need to rigorously assess the impact of the TAO observing system on ENSO predictions. Therefore, we evaluate the TAO using data assimilation observation denial experiments (also known as Observing System Evaluation or OSE experiments). This presentation will evaluate the TAO observing impact on both reanalyses and prediction of the big 2015 El Nino. We have completed reanalyses for July 2014-Dec 2015 for both the CONTROL experiment (that assimilates all available data) and the NOTAO experiment (that is identical to the CONTROL but withholds all TAO observations) using the best available NASA GEOS-S2S V3 seasonal prediction system. Validation of these reanalyses shows that TAO assimilation generally improves comparisons of temperature and salinity versus gridded in situ observations. Temperature is universally improved above the thermocline near the equator, but the biggest improvement is found in the eastern Pacific, just below the thermocline, where the variation of the thermocline defines ENSO events. In addition, a surprising result is that even with relatively few observations, salinity is improved throughout the equatorial region except near 120oW near the surface. ENSO forecasts were performed that were initialized from these CONTROL and NOTAO reanalyses. For the 9-month forecasts which were initialized in January, July, and October 2015, the NINO3.4 SST shows that the CONTROL forecasts are warmer and closer to observations than the NOTAO forecasts. We will show that upwelling and shoaling of the mixed layer amplifies the ENSO signal due to TAO assimilation. Prior to (after) April 2015, this upwelling is caused by relatively stronger Rossby (Kelvin) waves in the CONTROL than in the NOTAO experiments.

E. Hackert↗

A Comprehensive Forward Model for Spaceborne Radar Instruments

We present the development and validation of a comprehensive forward model designed to enhance remote sensing capabilities of spaceborne radar instruments. To overcome limitations in existing models, we integrated a Discrete Dipole Approximation (DDA) cloud scattering database into our Radiative Transfer Model (RTM), focusing on microwave frequencies. By simulating the optical properties of non-spherical frozen hydrometeors, the DDA technique effectively reduced discrepancies between simulated and observed values, surpassing traditional Mie tables. The evaluation of DDA lookup tables involved comparisons with a collocated dataset comprising short-term forecasts and satellite microwave data, providing evidence of their superiority. Additionally, we address the challenges of assimilating active radar measurements, which offer vertically resolved insights into clouds and precipitation. We explored the assimilation of spaceborne radar measurements in Numerical Weather Prediction (NWP) models by integrating a forward radar model, along with its adjoint and tangent linear, into the data assimilation system. Evaluation using CloudSat measurements demonstrated promising agreement between simulations and observations, particularly when the input hydrometeor profiles aligned with the measured reflectivity profiles, showcasing the potential of the developed forward radar model. Moreover, we discuss other challenges in radar measurement assimilation within NWP models, including potential observation errors and biases.

Isaac Moradi↗

Early Comparison of Ocean Analyses Using JEDI/SOCA 3 DVar vs. GEOS-S2S- 3 ODAS LETKF

At NASA’s Global Modeling and Assimilation Office (GMAO), data assimilation (DA) for the next-generation Goddard Earth Observing System Subseasonal-to-Seasonal (GEOS-S 2 S) coupled-model forecast system will transition to the Joint Effort for DA Integration (JEDI) system, which includes the marine DA component SOCA (Sea ice, Ocean, and Coupled Assimilation). It is envisioned that incorporating SOCA into GMAO Earth System Modeling will allow a more systematic approach to assimilating new data types (e.g., SWOT KaRIn), increasing resolution (e.g., 1/12°), and facilitating weakly (and eventually strongly) coupled DA (air-sea-ice, etc.). To prepare for this transition, testing is underway to compare SOCA results at high ocean resolution against the current GEOS-S 2 S Version 3 Ocean DA System (ODAS) 1 using similar initial conditions and observations.

David Russell↗

Smoothing-based compressed state Kalman filter for joint state-parameter estimation: Applications in reservoir characterization and CO 2 storage monitoring

The operation of most engineered hydrogeological systems relies on simulating physical processes using numerical models with uncertain parameters and initial conditions. Predictions by such uncertain models can be greatly improved by Kalman-filter techniques that sequentially assimilate monitoring data. Each assimilation constitutes a nonlinear optimization, which is solved by linearizing an objective function about the model prediction and applying a linear correction to this prediction. However, if model parameters and initial conditions are uncertain, the optimization problem becomes strongly nonlinear and a linear correction may yield unphysical results. In this paper, we investigate the utility of one-step ahead smoothing, a variant of the traditional filtering process, to eliminate nonphysical results and reduce estimation artifacts caused by nonlinearities. Here, we present the smoothing-based compressed state Kalman filter (sCSKF), an algorithm that combines one step ahead smoothing, in which current observations are used to correct the state and parameters one step back in time, with a nonensemble covariance compression scheme, that reduces the computational cost by efficiently exploring the high-dimensional state and parameter space. Numerical experiments show that when model parameters are uncertain and the states exhibit hyperbolic behavior with sharp fronts, as in CO 2 storage applications, one-step ahead smoothing reduces overshooting errors and, by design, gives physically consistent state and parameter estimates. We compared sCSKF with commonly used data assimilation methods and showed that for the same computational cost, combining one step ahead smoothing and nonensemble compression is advantageous for real-time characterization and monitoring of large-scale hydrogeological systems with sharp moving fronts.

42 ENGINEERING↗

Forecasting Responses of a Northern Peatland Carbon Cycle to Elevated CO 2 and a Gradient of Experimental Warming

The ability to forecast ecological carbon cycling is imperative to land management in a world where past carbon fluxes are no longer a clear guide in the Anthropocene. However, carbon–flux forecasting has not been practiced routinely like numerical weather prediction. This study explored (1) the relative contributions of model forcing data and parameters to uncertainty in forecasting flux– versus pool–based carbon cycle variables and (2) the time points when temperature and CO 2 treatments may cause statistically detectable differences in those variables. We developed an online forecasting workflow (Ecological Platform for Assimilation of Data (EcoPAD)), which facilitates iterative data–model integration. EcoPAD automates data transfer from sensor networks, data assimilation, and ecological forecasting. We used the Spruce and Peatland Responses Under Changing Experiments data collected from 2011 to 2014 to constrain the parameters in the Terrestrial Ecosystem Model, forecast carbon cycle responses to elevated CO 2 and a gradient of warming from 2015 to 2024, and specify uncertainties in the model output. Our results showed that data assimilation substantially reduces forecasting uncertainties. Interestingly, we found that the stochasticity of future external forcing contributed more to the uncertainty of forecasting future dynamics of C flux–related variables than model parameters. However, the parameter uncertainty primarily contributes to the uncertainty in forecasting C pool–related response variables. Given the uncertainties in forecasting carbon fluxes and pools, our analysis showed that statistically different responses of fast–turnover pools to various CO 2 and warming treatments were observed sooner than slow–turnover pools. In conclusion, our study has identified the sources of uncertainties in model prediction and thus leads to improve ecological carbon cycling forecasts in the future.

54 ENVIRONMENTAL SCIENCES↗