Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble data assimilation”

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

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

NASA SPoRT Modeling and Data Assimilation Research and Transition Activities Using WRF, LIS and GSI

weather research and forecasting ===== The NASA Short‐term Prediction Research and Transition (SPoRT) program has numerous modeling and data assimilation (DA) activities in which the WRF model is a key component. SPoRT generates realtime, research satellite products from the MODIS and VIIRS instruments, making the data available to NOAA/NWS partners running the WRF/EMS, including: (1) 2‐km northwestern‐hemispheric SST composite, (2) daily, MODIS green vegetation fraction (GVF) over CONUS, and (3) NASA Land Information System (LIS) runs of the Noah LSM over the southeastern CONUS. Each of these datasets have been utilized by specific SPoRT partners in local EMS model runs, with select offices evaluating the impacts using a set of automated scripts developed by SPoRT that manage data acquisition and run the NCAR Model Evaluation Tools verification package. SPoRT is engaged in DA research with the Gridpoint Statistical Interpolation (GSI) and Ensemble Kalman Filter in LIS for soil moisture DA. Ongoing DA projects using GSI include comparing the impacts of assimilating Atmospheric Infrared Sounder (AIRS) radiances versus retrieved profiles, and an analysis of extra‐tropical cyclones with intense non‐convective winds. As part of its Early Adopter activities for the NASA Soil Moisture Active Passive (SMAP) mission, SPoRT is conducting bias correction and soil moisture DA within LIS to improve simulations using the NASA Unified‐WRF (NU‐WRF) for both the European Space Agency's Soil Moisture Ocean Salinity and upcoming SMAP mission data. SPoRT has also incorporated real‐time global GVF data into LIS and WRF from the VIIRS product being developed by NOAA/NESDIS. This poster will highlight the research and transition activities SPoRT conducts using WRF, NU‐WRF, EMS, LIS, and GSI.

Case, Jonathan L.↗
NASA NTRS2017

Assimilation of MODIS Dark Target and Deep Blue Observations in the Dust Aerosol Component of NMMB-MONARCH version 1.0

A data assimilation capability has been built for the NMMB-MONARCH chemical weather prediction system, with a focus on mineral dust, a prominent type of aerosol. An ensemble-based Kalman filter technique (namely the local ensemble transform Kalman filter - LETKF) has been utilized to optimally combine model background and satellite retrievals. Our implementation of the ensemble is based on known uncertainties in the physical parametrizations of the dust emission scheme. Experiments showed that MODIS AOD retrievals using the Dark Target algorithm can help NMMB-MONARCH to better characterize atmospheric dust. This is particularly true for the analysis of the dust outflow in the Sahel region and over the African Atlantic coast. The assimilation of MODIS AOD retrievals based on the Deep Blue algorithm has a further positive impact in the analysis downwind from the strongest dust sources of the Sahara and in the Arabian Peninsula. An analysis-initialized forecast performs better (lower forecast error and higher correlation with observations) than a standard forecast, with the exception of underestimating dust in the long-range Atlantic transport and degradation of the temporal evolution of dust in some regions after day 1. Particularly relevant is the improved forecast over the Sahara throughout the forecast range thanks to the assimilation of Deep Blue retrievals over areas not easily covered by other observational datasets.The present study on mineral dust is a first step towards data assimilation with a complete aerosol prediction system that includes multiple aerosol species.

atmospheric dust↗
DOE OSTI2021

Global Temperature Responses to Large Tropical Volcanic Eruptions in Paleo Data Assimilation Products and Climate Model Simulations Over the Last Millennium

Abstract Large volcanic eruptions are one of the dominant perturbations to global and regional atmospheric temperatures on timescales of years to decades. Discrepancies remain, however, in the estimated magnitude and persistence of the surface temperature cooling caused by volcanic eruptions, as characterized by paleoclimatic proxies and climate models. We investigate these discrepancies in the context of large tropical eruptions over the Last Millennium using two state‐of‐the‐art data assimilation products, the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR), and simulations from the National Center for Atmospheric Research Community Earth System Model‐Last Millennium Ensemble (NCAR CESM‐LME). We find that PHYDA and LMR estimate mean global and hemispheric cooling that is similar in magnitude and persistence once effects from eruptions occurring in short succession are removed. The estimates also compare well to Northern‐Hemisphere reconstructions based solely or partially on tree‐ring density, which have been proposed as the most accurate proxy estimates of surface cooling due to volcanism. All proxy‐based estimates also agree well with the magnitude of the mean cooling simulated by the CESM‐LME. Differences remain, however, in the spatial patterns of the temperature responses in the PHYDA, LMR, and the CESM‐LME. The duration of cooling anomalies also persists for several years longer in the PHYDA and LMR relative to the CESM‐LME. Our results demonstrate progress in resolving discrepancies between proxy‐ and model‐based estimates of temperature responses to volcanism, but also indicate these estimates must be further reconciled to better characterize the risks of future volcanic eruptions.

Tejedor, E.↗
NASA NTRS2015

RACORO Continental Boundary Layer Cloud Investigations: 1. Case Study Development and Ensemble Large-Scale Forcings

Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60 h case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in situ measurements from the Routine AAF (Atmospheric Radiation Measurement (ARM) Aerial Facility) CLOWD (Clouds with Low Optical Water Depth) Optical Radiative Observations (RACORO) field campaign and remote sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, kappa, are derived from observations to be approximately 0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing data sets are derived from the ARM variational analysis, European Centre for Medium-Range Weather Forecasts, and a multiscale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in "trial" large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.

Boundary layers↗
NASA NTRS2016

Assimilation of SMOS Brightness Temperatures or Soil Moisture Retrievals into a Land Surface Model

Three different data products from the Soil Moisture Ocean Salinity (SMOS) mission are assimilated separately into the Goddard Earth Observing System Model, version 5 (GEOS-5) to improve estimates of surface and root-zone soil moisture. The first product consists of multi-angle, dual-polarization brightness temperature (Tb) observations at the bottom of the atmosphere extracted from Level 1 data. The second product is a derived SMOS Tb product that mimics the data at a 40 degree incidence angle from the Soil Moisture Active Passive (SMAP) mission. The third product is the operational SMOS Level 2 surface soil moisture (SM) retrieval product. The assimilation system uses a spatially distributed ensemble Kalman filter (EnKF) with seasonally varying climatological bias mitigation for Tb assimilation, whereas a time-invariant cumulative density function matching is used for SM retrieval assimilation. All assimilation experiments improve the soil moisture estimates compared to model-only simulations in terms of unbiased root-mean-square differences and anomaly correlations during the period from 1 July 2010 to 1 May 2015 and for 187 sites across the US. Especially in areas where the satellite data are most sensitive to surface soil moisture, large skill improvements (e.g., an increase in the anomaly correlation by 0.1) are found in the surface soil moisture. The domain-average surface and root-zone skill metrics are similar among the various assimilation experiments, but large differences in skill are found locally. The observation-minus-forecast residuals and analysis increments reveal large differences in how the observations add value in the Tb and SM retrieval assimilation systems. The distinct patterns of these diagnostics in the two systems reflect observation and model errors patterns that are not well captured in the assigned EnKF error parameters. Consequently, a localized optimization of the EnKF error parameters is needed to further improve Tb or SM retrieval assimilation.

SMOS↗
NASA NTRS2009

Assessing the Uncertainty of Land Surface Fluxes Derived from Operational Analyses

Operational numerical weather prediction systems assimilate tremendous amounts of observations to produce global analyses of the atmospheric states. Since the systems include significant contributions from the model components, physical fields, such as the turbulent heat fluxes, can also be derived. However, the physical fields are indirectly constrained by the data assimilation, and so subject to uncertainty in the model parameterizations. Recently, 10 international operational analyses have been collected and reprocessed through the Coordinated Enhanced. Observing Period (CEOP). Comparisons against independent observations (global precipitation and outgoing longwave radiation) show that while there is significant range of the individual analyses, the ensemble of analyses provide data closest the observations. Here, we evaluate the collection of land surface turbulent heat fluxes. The data are available for Oct 2002 - Dec 2004, at 6 hourly, daily and monthly intervals. An ensemble mean and its variance have been computed. we will analyze the range of analysis results at daily and seasonal time scales, and where possible compare with existing observations.

Bosilovich, Michael G.↗
DOE OSTI2023

National CO 2 budgets (2015–2020) inferred from atmospheric CO 2 observations in support of the global stocktake

Accurate accounting of emissions and removals of CO 2 is critical for the planning and verification of emission reduction targets in support of the Paris Agreement. Here, we present a pilot dataset of country-specific net carbon exchange (NCE; fossil plus terrestrial ecosystem fluxes) and terrestrial carbon stock changes aimed at informing countries' carbon budgets. These estimates are based on “top-down” NCE outputs from the v10 Orbiting Carbon Observatory (OCO-2) modeling intercomparison project (MIP), wherein an ensemble of inverse modeling groups conducted standardized experiments assimilating OCO-2 column-averaged dry-air mole fraction (X CO 2 ) retrievals (ACOS v10), in situ CO 2 measurements or combinations of these data. The v10 OCO-2 MIP NCE estimates are combined with “bottom-up” estimates of fossil fuel emissions and lateral carbon fluxes to estimate changes in terrestrial carbon stocks, which are impacted by anthropogenic and natural drivers. These flux and stock change estimates are reported annually (2015–2020) as both a global 1° x 1° gridded dataset and a country-level dataset and are available for download from the Committee on Earth Observation Satellites' (CEOS) website. Across the v10 OCO-2 MIP experiments, we obtain increases in the ensemble median terrestrial carbon stocks of 3.29–4.58 Pg CO2 yr -1 (0.90–1.25 Pg C yr -1 ). This is a result of broad increases in terrestrial carbon stocks across the northern extratropics, while the tropics generally have stock losses but with considerable regional variability and differences between v10 OCO-2 MIP experiments. We discuss the state of the science for tracking emissions and removals using top-down methods, including current limitations and future developments towards top-down monitoring and verification systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Remote Sensing Data from CLARET: A Prototype Cart Data Set

A data set containing radiation, meteorological, and cloud sensor observations is documented. It was prepared for use by the Department of Energy's Atmospheric Radiation Measurement (ARM) program and other interested scientists. These data are a precursor of the types of data that ARM Cloud And Radiation Testbed (CART) sites will provide. The data are from the Cloud Lidar And Radar Exploratory Test (CLARET) conducted by the Wave Propagation Laboratory during autumn 1989 in the Denver-Boulder area of Colorado primarily for the purpose of developing new cloud-sensing techniques on cirrus. After becoming aware of this experiment, ARM scientists requested archival of subsets or the data to assist in the developing ARM program. Five CLARET cases were selected: two with cirrus, one with stratus, one with mixed-phase clouds, and one with clear skies. The cases range from 2 to 9.5 h in length. A pyranometer, pyrgeometer, pyrheliometer, and an infrared radiometer constituted the ensemble of instruments that provided surface radiation data. A lidar, radar, and ceilometer observed the cloud geometrical structure, and visual reports and all-sky camera observations were assimilated to provide cloud cover data. Radiosondes, wind profiler, RASS (profiling virtual temperature), microwave radiometers (observing column integrated liquid water and water vapor), and standard surface measurements provided meteorological data. Satellite data from the stratus case and one cirrus case were analyzed for statistics on cloud cover and top height. The main body of the selected data are available on diskette from the Wave Propagation Laboratory or Los Alamos National Laboratory. In addition to documenting the data set, this report describes CLARET and gives a bibliography of publications associated with the project. Some preliminary results of CLARET' research are also summarized. Simultaneous CO 2 lidar and radar backscatter measurements were shown to provide estimates of the effective radius of ice particles. Simultaneous radar and infrared radiometer data appear useful for estimating column-integrated numbers and average sizes of ice cloud particles. Ice water content obtained with this method compared favorably with values from another empirical technique using radar data alone. Depolarization of the CO 2 lidar signal from ice clouds was surprisingly small, suggesting that calculation of backscatter from nonspherical particles for this lidar is a tractable problem. Examples are also cited of CO 2 lidar measurements of the effective radius of water cloud drop size distributions and of inference of the size of pristine ice crystals that assume a particular orientation in the air. These parameters are all important to radiative transfer through clouds.

Clouds (Meteorology)↗
NASA NTRS2012

Impact of Satellite Atmospheric Motion Vectors In the GMAO GEOS-5 Global Data Assimilation System

The WMO and THORPEX co-sponsored fifth Workshop on the Impact of Various Observing Systems on Numerical Weather Prediction will be organized by the Expert Team on the Evolution of the Global Observing System in Sedona, Arizona, United States, from 22 to 25 May 2012. Participants are expected to come from all the major NWP centres which are active in the area of impact studies. The workshop will be conducted in English. As for the first four workshops it is planned to produce a workshop report to be published as a WMO Technical Report that will include the papers submitted by the participants. The previous four workshops in this series took place in Geneva {April 1997), Toulouse (March 2000), Alpbach (March 2004) and Geneva (May 2008). Results from Observing System Experiments (OSEs), both with global and regional aspects were presented and conclusions were drawn concerning the contributions of the various components of the observing system to the large scale forecast skill at short and medium range (Workshop Proceedings were published as WMO World Weather Watch Technical Reports TD No. 868, 1034, 1228 and 1450). Since then, some significant changes and developments have affected the global observing system and more efforts have been devoted to meso-scale observing and assimilation systems. There has also been a trend toward using techniques other than OSEs to document data impact, such as adjoint-based sensitivity to observations or ensemble-based sensitivity. Field experiments have been carried out, in particular through the THORPEX project, and the use of targeted data has been assessed.

Gelaro, Ronald↗
NASA NTRS2022

Assimilation of SMAP Products for Improving Streamflow Simulations over Tropical Climate Region—Is Spatial Information More Important Than Temporal Information?

Streamflow is one of the key variables in the hydrological cycle. Simulation and forecasting of streamflow are challenging tasks for hydrologists, especially in sparsely gauged areas. Coarse spatial resolution remote sensing soil moisture products (equal to or larger than 9 km) are often assimilated into hydrological models to improve streamflow simulation in large catchments. This study uses the Ensemble Kalman Filter (EnKF) technique to assimilate SMAP soil moisture products at the coarse spatial resolution of 9 km (SMAP 9 km), and downscaled SMAP soil moisture product at the higher spatial resolution of 1 km (SMAP 1 km), into the Soil and Water Assessment Tool (SWAT) to investigate the usefulness of different spatial and temporal resolutions of remotely sensed soil moisture products in streamflow simulation and forecasting. The experiment was set up for eight catchments across the tropical climate of Vietnam, with varying catchment areas from 267 to 6430 km^2 during the period 2017–2019. We comprehensively evaluated the EnKF-based SWAT model in simulating streamflow at low, average, and high flow. Our results indicated that high-spatial resolution of downscaled SMAP 1 km is more beneficial in the data assimilation framework in aiding the accuracy of streamflow simulation, as compared to that of SMAP 9 km, especially for the small catchments. Our analysis on the impact of observation resolution also indicates that the improvement in the streamflow simulation with data assimilation is more significant at catchments where downscaled SMAP 1 km has fewer missing observations. This study is helpful for adding more understanding of performances of soil moisture data assimilation based hydrological modelling over the tropical climate region, and exhibits the potential use of remote sensing data assimilation in hydrology.

soil moisture↗
NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗
NASA NTRS2011

Estimating Evapotranspiration with Land Data Assimilation Systems

Advancements in both land surface models (LSM) and land surface data assimilation, especially over the last decade, have substantially advanced the ability of land data assimilation systems (LDAS) to estimate evapotranspiration (ET). This article provides a historical perspective on international LSM intercomparison efforts and the development of LDAS systems, both of which have improved LSM ET skill. In addition, an assessment of ET estimates for current LDAS systems is provided along with current research that demonstrates improvement in LSM ET estimates due to assimilating satellite-based soil moisture products. Using the Ensemble Kalman Filter in the Land Information System, we assimilate both NASA and Land Parameter Retrieval Model (LPRM) soil moisture products into the Noah LSM Version 3.2 with the North American LDAS phase 2 (NLDAS-2) forcing to mimic the NLDAS-2 configuration. Through comparisons with two global reference ET products, one based on interpolated flux tower data and one from a new satellite ET algorithm, over the NLDAS2 domain, we demonstrate improvement in ET estimates only when assimilating the LPRM soil moisture product.

Peters-Lidard, C. D.↗
DOE OSTI2022

Scalable deep learning for watershed model calibration

Watershed models such as the Soil and Water Assessment Tool (SWAT) consist of high-dimensional physical and empirical parameters. These parameters often need to be estimated/calibrated through inverse modeling to produce reliable predictions on hydrological fluxes and states. Existing parameter estimation methods can be time consuming, inefficient, and computationally expensive for high-dimensional problems. In this paper, we present an accurate and robust method to calibrate the SWAT model (i.e., 20 parameters) using scalable deep learning (DL). We developed inverse models based on convolutional neural networks (CNN) to assimilate observed streamflow data and estimate the SWAT model parameters. Scalable hyperparameter tuning is performed using high-performance computing resources to identify the top 50 optimal neural network architectures. We used ensemble SWAT simulations to train, validate, and test the CNN models. We estimated the parameters of the SWAT model using observed streamflow data and assessed the impact of measurement errors on SWAT model calibration. We tested and validated the proposed scalable DL methodology on the American River Watershed, located in the Pacific Northwest-based Yakima River basin. Our results show that the CNN-based calibration is better than two popular parameter estimation methods (i.e., the generalized likelihood uncertainty estimation [GLUE] and the dynamically dimensioned search [DDS], which is a global optimization algorithm). For the set of parameters that are sensitive to the observations, our proposed method yields narrower ranges than the GLUE method but broader ranges than values produced using the DDS method within the sampling range even under high relative observational errors. The SWAT model calibration performance using the CNNs, GLUE, and DDS methods are compared using R 2 and a set of efficiency metrics, including Nash-Sutcliffe, logarithmic Nash-Sutcliffe, Kling-Gupta, modified Kling-Gupta, and non-parametric Kling-Gupta scores, computed on the observed and simulated watershed responses. The best CNN-based calibrated set has scores of 0.71, 0.75, 0.85, 0.85, 0.86, and 0.91. The best DDS-based calibrated set has scores of 0.62, 0.69, 0.8, 0.77, 0.79, and 0.82. The best GLUE-based calibrated set has scores of 0.56, 0.58, 0.71, 0.7, 0.71, and 0.8. The scores above show that the CNN-based calibration leads to more accurate low and high streamflow predictions than the GLUE and DDS sets. Our research demonstrates that the proposed method has high potential to improve our current practice in calibrating large-scale integrated hydrologic models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A High‐Resolution Tropical Mesoscale Convective System Reanalysis (TMeCSR)

Abstract Modern global reanalysis products have greatly accelerated meteorological research in synoptic‐to‐planetary‐scale phenomena. However, their use in studying tropical mesoscale convective systems (MCSs) and their regional‐to‐global impact has mostly been limited to supplying initial and boundary conditions for MCS‐resolving simulations and providing information about the large‐scale environments of MCSs. These limitations are due to difficulties in resolving tropical MCS dynamics in the relatively low‐resolution global models and that tropical MCSs often occur over poorly observed regions. In this work, a Tropical MCS‐resolving Reanalysis product (TMeCSR) was created over a region with frequent tropical MCSs. This region spans the tropical Indian Ocean, tropical continental Asia, Maritime Continent, and Western Pacific. TMeCSR is produced by assimilating all‐sky infrared radiances from geostationary satellites and other conventional observations into an MCS‐resolving regional model using the Ensemble Kalman Filter. The resulting observation‐constrained high‐resolution (9‐km grid spacing) data set is available hourly during the boreal summer (June‐August) of 2017, during which widespread severe flooding occurred. Comparisons of TMeCSR and European Center for Medium Range Weather Forecast Reanalysis version 5 (ERA5) against independent satellite retrievals indicate that TMeCSR's cloud and multiscale rain fields are better than those of ERA5. Furthermore, TMeCSR better captured the diurnal variability of rainfall and the statistical characteristics of MCSs. Forecasts initialized from TMeCSR also have more accurate rain and clouds than those initialized from ERA5. The TMeCSR and ERA5 forecasts have similar performances with respect to sounding and surface observations. These results indicate that TMeCSR is a promising MCS‐resolving data set for tropical MCS studies.

Chan, Man‐Yau↗
NASA NTRS2011

Land Surface Verification Toolkit (LVT) - A Generalized Framework for Land Surface Model Evaluation

Model evaluation and verification are key in improving the usage and applicability of simulation models for real-world applications. In this article, the development and capabilities of a formal system for land surface model evaluation called the Land surface Verification Toolkit (LVT) is described. LVT is designed to provide an integrated environment for systematic land model evaluation and facilitates a range of verification approaches and analysis capabilities. LVT operates across multiple temporal and spatial scales and employs a large suite of in-situ, remotely sensed and other model and reanalysis datasets in their native formats. In addition to the traditional accuracy-based measures, LVT also includes uncertainty and ensemble diagnostics, information theory measures, spatial similarity metrics and scale decomposition techniques that provide novel ways for performing diagnostic model evaluations. Though LVT was originally designed to support the land surface modeling and data assimilation framework known as the Land Information System (LIS), it also supports hydrological data products from other, non-LIS environments. In addition, the analysis of diagnostics from various computational subsystems of LIS including data assimilation, optimization and uncertainty estimation are supported within LVT. Together, LIS and LVT provide a robust end-to-end environment for enabling the concepts of model data fusion for hydrological applications. The evolving capabilities of LVT framework are expected to facilitate rapid model evaluation efforts and aid the definition and refinement of formal evaluation procedures for the land surface modeling community.

Kumar, Sujay V.↗
NASA NTRS2018

Insight and Evidence Motivating the Simplification of Dual-Analysis Hybrid Systems into Single-Analysis Hybrid Systems

Many hybrid data assimilation systems currently used for NWP employ some form of dual-analysis system approach. Typically a hybrid variational analysis is responsible for creating initial conditions for high-resolution forecasts, and an ensemble analysis system is responsible for creating sample perturbations used to form the flow-dependent part of the background error covariance required in the hybrid analysis component. In many of these, the two analysis components employ different methodologies, e.g., variational and ensemble Kalman filter. In such cases, it is not uncommon to have observations treated rather differently between the two analyses components; recentering of the ensemble analysis around the hybrid analysis is used to compensated for such differences. Furthermore, in many cases, the hybrid variational high-resolution system implements some type of four-dimensional approach, whereas the underlying ensemble system relies on a three-dimensional approach, which again introduces discrepancies in the overall system. Connected to these is the expectation that one can reliably estimate observation impact on forecasts issued from hybrid analyses by using an ensemble approach based on the underlying ensemble strategy of dual-analysis systems. Just the realization that the ensemble analysis makes substantially different use of observations as compared to their hybrid counterpart should serve as enough evidence of the implausibility of such expectation. This presentation assembles numerous anecdotal evidence to illustrate the fact that hybrid dual-analysis systems must, at the very minimum, strive for consistent use of the observations in both analysis sub-components. Simpler than that, this work suggests that hybrid systems can reliably be constructed without the need to employ a dual-analysis approach. In practice, the idea of relying on a single analysis system is appealing from a cost-maintenance perspective. More generally, single-analysis systems avoid contradictions such as having to choose one sub-component to generate performance diagnostics to another, possibly not fully consistent, component.

Todling, Ricardo↗
NASA NTRS1999

Improving Assimilated Global Data Sets using TMI Rainfall and Columnar Moisture Observations

A global analysis that optimally combine observations from diverse sources with physical models of atmospheric and land processes can provide a comprehensive description of the climate systems. Currently, such data products contain significant errors in primary hydrological fields such as precipitation and evaporation, especially in the tropics. In this study, we show that assimilating precipitation and total precipitable water (TPW) retrievals derived from the TRMM Microwave Imager (TMI) improves not only the hydrological cycle but also key climate parameters such as clouds, radiation, and the large-scale circulation produced by the Goddard Earth Observing System (GEOS) data assimilation system (DAS). In particular, assimilating TMI rain improves clouds and radiation in areas of active convection, as well as the latent heating distribution and the large-scale motion field in the tropics, while assimilating TMI TPW heating distribution and the large-scale motion field in the tropics, while assimilating TMI TPW retrievals leads to reduced moisture biases and improved radiative fluxes in clear-sky regions. The improved analysis also improves short-range forecasts in the tropics. Ensemble forecasts initialized with the GEOS analysis incorporating TMI rain rates and TPW yield smaller biases in tropical precipitation forecasts beyond 1 day and better 500 hPa geopotential height forecasts up to 5 days. Results of this study demonstrate the potential of using high-quality space-borne rainfall and moisture observations to improve the quality of assimilated global data for climate analysis and weather forecasting applications

Hou, Arthur Y.↗
NASA NTRS2024

Integrating State Data Assimilation and Innovative Model Parameterization Reduces Simulated Carbon Uptake in the Arctic and Boreal Region

Model representation of carbon uptake and storage is essential for accurate projection of the response of the arctic‐boreal zone to a rapidly changing climate. Land model estimates of LAI and aboveground biomass that can have a marked influence on model projections of carbon uptake and storage vary substantially in the arctic and boreal zone, making it challenging to correctly evaluate model estimates of Gross Primary Productivity (GPP). To understand and correct bias of LAI and aboveground biomass in the Community Land Model (CLM), we assimilated the 8‐day Moderate Resolution Imaging Spectroradiometer (MODIS) LAI observation and a machine learning product of annual aboveground biomass into CLM using an Ensemble Adjustment Kalman Filter (EAKF) in an experimental region including Alaska and Western Canada. Assimilating LAI and aboveground biomass reduced these model estimates by 58% and 72%, respectively. The change of aboveground biomass was consistent with independent estimates of canopy top height at both regional and site levels. The International Land Model Benchmarking system assessment showed that data assimilation significantly improved CLM's performance in simulating the carbon and hydrological cycles, as well as in representing the functional relationships between LAI and other variables. To further reduce the remaining bias in GPP after LAI bias correction, we re‐parameterized CLM to account for low temperature suppression of photosynthesis. The LAI bias corrected model that included the new parameterization showed the best agreement with model benchmarks. Combining data assimilation with model parameterization provides a useful framework to assess photosynthetic processes in LSMs.

land data assimilation↗