SEARCH · Engineering Papers

Results for “forecast error”

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

SourceResearch connectionFrom yearThrough year

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

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

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 NTRS2004

Current results from AlRS/AMSU/HSB

AIRS was launched on EOS Aqua on May 4,2002, together with AMSU A and HSB, to form a next generation polar orbiting infrared and microwave atmospheric sounding system. The primary products of AIRS/AMSU/HSB are twice daily global fields of atmospheric temperature-humidity profiles, ozone profiles, sea/land surface skin temperature, and cloud related parameters including OLR. The sounding goals of AIRS are to produce 1 km tropospheric layer mean temperatures with an rms error of 1K, and layer precipitable water with an rms error of 20%, in cases with up to 80% effective cloud cover. Pre-launch simulation studies indicated that these results should be achievable. Minor modifications have been made to the pre-launch retrieval algorithm as alluded to in this paper. Sample fields of parameters retrieved from AIRS/AMSU/HSB data are presented and temperature profiles are validated as a function of retrieved fractional cloud cover. As in simulation, the degradation of retrieval accuracy with increasing cloud cover is small. Select fields are also compared to those contained in the ECMWF analysis, done without the benefit of AIRS data, to demonstrate information that AIRS can add to that already contained in the ECMWF analysis. Assimilation of AIRS temperature soundings in up to 80% cloud cover for the month of January 2003 into the GSFC FVSSI data assimilation system resulted in improved 5 day forecasts globally, both with regard to anomaly correction coefficients and the prediction of location and intensity of cyclones.

Susskind, Joel↗

NASA NTRS2015

Effects of Resolution and Spectral Nudging in Simulation the Effects of Wintertime Atmospheric River Landfalls in the Western US

Landfalling atmospheric rivers (ARs) play a crucial role in the climate of the US Pacific coast region as they are frequently related with heavy precipitation and flash flooding events. Thus, the capability of climate models to accurately simulate AR landfalls and their key hydrologic effects is an important practical concern for WUS, from flood forecasting to future water resources projections. In order to examine the effects of model configuration, including the resolution and spectral nudging, in simulating the climatology of key weather events in the conterminous US, a NASA team has performed a hindcast experiment using the GEOS5 global and the NU-WRF regional models for Nov 1999 - Oct 2010. This study examines the skill of these hindcasts, with different models and their configurations, in simulating key footprints of landfalling ARs in the WUS region. Using an AR-landfall chronology based on the vertically-integrated water vapor flux calculated from the MERRA2 reanalysis, we have analyzed the observed and simulated precipitation and temperature anomalies associated with wintertime AR landfalls along the US Pacific coast. Model skill is measured using metrics including regional means, a skill score based on correlations and mean-square errors, and Taylor diagrams in four WUS Bukovsky regions. Results show that the AR-related anomalies of precipitation is more reliable than of surface temperatures. Model skill also varies according to regions. The AR temperature anomalies are well simulated in most of the WUS region except PNW. For precipitation, simulations with finer spatial resolution tend to generate larger spatial variability and agree better with the PRISM data in most regions. Such a resolution dependence of spatial variability is not found for temperatures; e.g., the MERRA2 reanalysis often outperforms, with similar spatial variability and higher pattern correlations with the PRISM data, finer-resolution NU-WRF runs in simulating temperature variations within subregions. Results from this study will be summarized to assist future (regional) climate experiments for climate change impact assessments and developing adaptationmitigation strategies, the key elements of the National Climate Assessment.

Simulation↗

NASA NTRS1994

Analyses and forecasts with LAWS winds

Horizontal fluxes of atmospheric water vapor are studied for summer months during 1989 and 1992 over North and South America based on analyses from European Center for Medium Range Weather Forecasts, US National Meteorological Center, and United Kingdom Meteorological Office. The calculations are performed over 20 deg by 20 deg box-shaped midlatitude domains located to the east of the Rocky Mountains in North America, and to the east of the Andes Mountains in South America. The fluxes are determined from operational center gridded analyses of wind and moisture. Differences in the monthly mean moisture flux divergence determined from these analyses are as large as 7 cm/month precipitable water equivalent over South America, and 3 cm/month over North America. Gridded analyses at higher spatial and temporal resolution exhibit better agreement in the moisture budget study. However, significant discrepancies of the moisture flux divergence computed from different gridded analyses still exist. The conclusion is more pessimistic than Rasmusson's estimate based on station data. Further analysis reveals that the most significant sources of error result from model surface elevation fields, gaps in the data archive, and uncertainties in the wind and specific humidity analyses. Uncertainties in the wind analyses are the most important problem. The low-level jets, in particular, are substantially different in the different data archives. Part of the reason for this may be due to the way the different analysis models parameterized physical processes affecting low-level jets. The results support the inference that the noise/signal ratio of the moisture budget may be improved more rapidly by providing better wind observations and analyses than by providing better moisture data.

Wang, Muyin↗

NASA NTRS2004

Atmospheric Soundings from AIRS/AMSU/HSB

AIRS was launched on EOS Aqua on May 4, 2002, together with AMSU A and HSB, to form a next generation polar orbiting infrared and microwave atmospheric sounding system. The primary products of AIRS/AMSU/HSB are twice daily global fields of atmospheric temperature-humidity profiles, ozone profiles, sea/land surface skin temperature, and cloud related parameters including OLR. The sounding goals of AIRS are to produce 1 km tropospheric layer mean temperatures with an rms error of lK, and 1 km tropospheric layer precipitable water with an rms error of 20%, in cases with up to 80% effective cloud cover. Pre-launch simulation studies indicated that these results should be achievable. Minor modifications have been made to the pre-launch retrieval algorithm as alluded to in this paper. Sample fields of parameters retrieved from AIRS/AMSU/HSB data are presented and temperature profiles are validated as a function of retrieved effective fractional cloud cover. As in simulation, the degradation of retrieval accuracy with increasing cloud cover is small. Select fields are also compared to those contained in the ECMWF analysis, done without the benefit of AIRS data, to demonstrate information that AIRS can add to that already contained in the ECMWF analysis. Assimilation of AIRS temperature soundings in up to 80% cloud cover for the month of January 2003 into the GSFC FVSSI data assimilation system resulted in improved 5 day forecasts globally, both with regard to anomaly correlation coefficients and the prediction of location and intensity of cyclones.

Susskind, Joel↗

DOE OSTI2022

Can reanalysis products outperform mesoscale numerical weather prediction models in modeling the wind resource in simple terrain?

Mesoscale numerical weather prediction (NWP) models are generally considered more accurate than reanalysis products in characterizing the wind resource at heights of interest for wind energy, given their finer spatial resolution and more comprehensive physics. However, advancements in the latest ERA-5 reanalysis product motivate an assessment on whether ERA-5 can model wind speeds as well as a state-of-the-art NWP model – the Weather Research and Forecasting (WRF) Model. We consider this research question for both simple terrain and offshore applications. Specifically, we compare wind profiles from ERA-5 and the preliminary WRF runs of the Wind Integration National Dataset (WIND) Toolkit Long-term Ensemble Dataset (WTK-LED) to those observed by lidars at a site in Oklahoma, United States, and in a United States Atlantic offshore wind energy area. We find that ERA-5 shows a significant negative bias (~-1ms-1) at both locations, with a larger bias at the land-based site. WTK-LED-predicted wind speed profiles show a limited negative bias (~-0.5ms-1) offshore and a slight positive bias (~+0.5ms-1) at the land-based site. On the other hand, we find that ERA-5 outperforms WTK-LED in terms of the centered root-mean-square error (cRMSE) and correlation coefficient, for both the land-based and offshore cases, in all atmospheric stability conditions. We find that WTK-LED's higher cRMSE is caused by its tendency to overpredict the amplitude of the wind speed diurnal cycle. At the land-based site, this is partially caused by wind plant wake effects not being accurately captured by WTK-LED.

17 WIND ENERGY↗

NASA NTRS2021

A Machine Learning Approach to Objective Identification of Dust in Satellite Imagery

Airborne dust has broad adverse effects on human activity, including aviation, human health, and agriculture. Remote sensing observations are used to detect dust and aerosols in the atmosphere using long established techniques. False color Red-Green-Blue (RGB) imagery using band differences sensitive to dust absorption (Dust RGB) is currently used operationally to assist forecasters and decision-makers in identifying dust at night, but there are still limitations, subjectivity, and nuances to image interpretation making night-time dust identification difficult even for experts. This study applies machine learning to the problem of night-time dust detection with a simple random forest (RF) model using Geostationary Operational Environmental Satellite-16 (GOES-16) Advanced Baseline Imager (ABI) infrared imagery, band differences sensitive to dust absorption, and Dust RGB color components as inputs to the model. The RF model achieves an Area-Under-Curve (AUC) of 0.97 with a standard deviation of 0.04 for dust cases. For images with dust present, the model correctly labels 85% of dust pixels and 99.96% of no-dust pixels for all dust images in the validation data set. The addition of a single null case to the training data set drastically reduces error in labeling no-dust pixels as dust from 45% to 14.5%. Application of the machine learning model to the April 13–14, 2019 dust event demonstrates the ability of the model to identify dust during night-time hours when visual dust detection is limited by the cooling ground surface characteristics.

dust↗

NASA NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗

NASA NTRS2003

Lidar Measurements of Wind, Moisture and Boundary Layer Evolution in a Dryline During IHOP2002

Variability in the convective boundary layer moisture, wind and temperature fields and their importance in the forecasting and understanding of storms have been discussed in the literature. These variations have been reported in relation to frontal zones, stationary boundaries and during horizontal convective rolls. While all three vary substantially in the convective boundary layer, moisture poses a particular challenge. Moisture or water vapor concentration (expressed as a mass mixing ratio, g/kg), is conserved in all meteorological processes except condensation and evaporation. The water vapor mixing ratio often remains distinct across an air -mass boundary even when the temperature difference is indistinct. These properties make it an ideal choice in visualizing and understanding many of the atmosphere's dynamic features. However, it also presents a unique measurement challenge because water vapor content can vary by more than three orders of magnitude in the troposphere. Characterization of the 3D-distribution of water vapor is also difficult as water vapor observations can suffer from large sampling errors and substantial variability both in the vertical and horizontal. This study presents groundbased measurements of wind, boundary layer structure and water vapor mixing ratio measurements observed by three co-located lidars. This presentation will focus on the evolution and variability of moisture and wind in the boundary layer during a dry line event that occurred on 22 May 2002. These data sets and analyses are unique in that they combine simultaneous measurements of wind, moisture and CBL structure to study the detailed thermal variability in and around clear air updrafts during a dryline event. It will quantify the variation caused by, in and around buoyant plumes and across a dryline. The data presented here were collected in the panhandle of Oklahoma as part of the International BO Project (IHOP-2002), a field experiment that took place over the Southern Great Plains (SGP) of the United States from 13 May to 30 June 2002. The chief goal of IHOP-2002 is to improve characterization of the four-dimensional (4-D) distribution of water vapor and its application to improving the understanding and prediction of convection

Demoz, Belay↗

NASA NTRS2010

Detailed Evaluation of MODIS Fire Radiative Power Measurements

Satellite remote sensing is providing us tremendous opportunities to measure the fire radiative energy (FRE) release rate or power (FRP) from open biomass burning, which affects many vegetated regions of the world on a seasonal basis. Knowledge of the biomass burning characteristics and emission source strengths of different (particulate and gaseous) smoke constituents is one of the principal ingredients upon which the assessment, modeling, and forecasting of their distribution and impacts depend. This knowledge can be gained through accurate measurement of FRP, which has been shown to have a direct relationship with the rates of biomass consumption and emissions of major smoke constituents. Over the last decade or so, FRP has been routinely measured from space by both the MODIS sensors aboard the polar orbiting Terra and Aqua satellites, and the SEVIRI sensor aboard the Meteosat Second Generation (MSG) geostationary satellite. During the last few years, FRP has been gaining recognition as an important parameter for facilitating the development of various scientific studies relating to the quantitative characterization of biomass burning and their emissions. Therefore, we are conducting a detailed analysis of the FRP products from MODIS to characterize the uncertainties associated with them, such as those due to the MODIS bow-tie effects and other factors, in order to establish their error budget for use in scientific research and applications. In this presentation, we will show preliminary results of the MODIS FRP data analysis, including comparisons with airborne measurements.

Ichoku, Charles↗

NASA NTRS2017

HMI Data Driven Magnetohydrodynamic Model Predicted Active Region Photospheric Heating Rates: Their Scale Invariant, Flare Like Power Law Distributions, and Their Possible Association With Flares

There are many flare forecasting models. For an excellent review and comparison of some of them see Barnes et al. (2016). All these models are successful to some degree, but there is a need for better models. We claim the most successful models explicitly or implicitly base their forecasts on various estimates of components of the photospheric current density J, based on observations of the photospheric magnetic field B. However, none of the models we are aware of compute the complete J. We seek to develop a better model based on computing the complete photospheric J. Initial results from this model are presented in this talk. We present a data driven, near photospheric, 3 D, non-force free magnetohydrodynamic (MHD) model that computes time series of the total J, and associated resistive heating rate in each pixel at the photosphere in the neutral line regions (NLRs) of 14 active regions (ARs). The model is driven by time series of B measured by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. Spurious Doppler periods due to SDO orbital motion are filtered out of the time series of B in every AR pixel. Errors in B due to these periods can be significant.

Goodman, Michael L.↗

DOE OSTI2022

Coupling localized Noah-MP-Crop model with the WRF model improved dynamic crop growth simulation across Northeast China

Croplands play a critical role in regulating the energy and moisture exchanges between the land surface and atmosphere. However, the interactions between cropland and climate are usually poorly represented due to a lack of detailed representation in crop types and field management. Here, we coupled the Noah-MP-Crop model with the state-of-the-art Weather Research and Forecasting (WRF) model to explore and evaluate the crop growth dynamics in response to climate variations across Northeast China. The default parameters of the crop model were not exactly suitable for the agricultural ecosystems in Northeast China. The detailed cropland distribution, and crop phenology parameters including growing degree days (GDD) and planting (harvesting) date were first created using multi-source remote sensing products and reanalysis data, and was then successfully used to simulate the growth and yield for corn and soybean and associated energy exchanges. We also optimized and calibrated other crop parameters using the time-series of the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface products. The modified crop model substantially improved the simulation of crop growth, plant physiology, and biomass accumulation for both corn and soybean. Coupling the localized dynamic crop model into the WRF led to considerable decreases in the simulated mean-absolute-errors (MAEs) and biases of the leaf area index, evapotranspiration, and gross primary production compared with the MODIS observed values. Compared with the statistical yield from each province, the modified crop model underestimated the corn yield from 11.1% to 48.6%, whereas overestimated the soybean yield from 16.5% to 162.6%.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2002

Lidar Measurements of Wind, Moisture, and Boundary Layer Evolution in a Dry Line during 1HOP 2002

Variability in the convective boundary layer moisture, wind and temperature fields and their importance in the forecasting and understanding of storms have been discussed in the literature. These . variations have been reported in relation to frontal zones, stationary boundaries and during horizontal convective rolls. While all three vary substantially in the convective boundary layer, moisture poses a particular challenge. Moisture or water vapor concentration (expressed as a mass mixing ratio, g/kg), is conserved in all meteorological processes except condensation and evaporation. The water vapor mixing ratio often remains distinct across an air-mass boundary even when the temperature difference is indistinct. These properties make it an ideal choice in visualizing and understanding many of the atmosphere's dynamic features. However, it also presents a unique measurement challenge because water vapor content can vary by more than three orders of magnitude in the troposphere. Characterization of the 3D-distribution of water vapor is also difficult as water vapor observations can suffer from large sampling errors and substantial variability both in the vertical and horizontal. This study presents ground-based measurements of wind, boundary layer structure and water vapor mixing ratio measurements observed by three co-located lidars. This presentation will focus on the evolution and variability of moisture and wind in the boundary layer during a dry line event that occurred on 22 May 2002. These data sets and analyses are unique in that they combine simultaneous measurements of wind, moisture and CBL structure to study the detailed thermal variability in and around clear air updrafts during a dryline event. It will quantify the variation caused by, in and around buoyant plumes and across a dryline. The data presented here were collected in the panhandle of Oklahoma as part of the International H2O Project (MOP-2002), a field experiment that took place over the Southern Great Plains (SGP) of the United States from 13 May to 30 June 2002. The chief goal of MOP-2002 is to improve characterization of the four-dimensional (4-D) distribution of water vapor and its application to improving the understanding and prediction of convection

Demoz, Belay↗

DOE OSTI2021

Bias Characterization, Vertical Interpolation, and Horizontal Interpolation for Distributed Wind Siting Using Mesoscale Wind Resource Estimates

Much like their counterparts in utility-scale wind energy, developers of industrial, small-scale and distributed wind turbine deployments need to understand and accurately characterize the wind resource to properly assess the power generation and financial ramifications during siting and planning. National Renewable Energy Laboratory’s WIND (Wind Integration National Dataset) Toolkit (WTK) provides a best-in-class wind resource dataset generated using the Weather Research and Forecasting (WRF) model. This dataset includes parameters such as the wind speed, wind direction, and temperature at various heights, plus atmospheric stability near the surface. This data is available at 2-km spatial resolution and five-minute temporal resolution across 7 years, from 2007 to 2013 through a publicly accessible API interface. The Tools Assessing Performance (TAP) project seeks to extend this dataset to allow long term resource estimates and leverage it to better equip distributed wind equipment manufacturers, owner-operators, and installation professionals with better tools for practical siting applications. In this report, we present the results from our investigation within the TAP project focused on characterization of bias in WTK-based wind speed estimates and evaluation of vertical and horizontal interpolation techniques. We discuss the tradeoffs between different techniques and their combinations, as well as describe the lower bounds we determine for the studied validation errors. While the specific estimates we present are specific to WTK and the validation dataset we have chosen for this investigation (NREL's Wind Resource Meteorological Database), the overall analysis and the studied techniques are general enough to be applied to a broader set of wind datasets, both simulation-based and observational.

17 WIND ENERGY↗

NASA NTRS2016

Impact of the Timing of a SAR Image Acquisition on the Calibration of a Flood Inundation Model

Synthetic Aperture Radar (SAR) data have proven to be a very useful source of information for the calibration of flood inundation models. Previous studies have focused on assigning uncertainties to SAR images in order to improve flood forecast systems (e.g. Giustarini et al. (2015) and Stephens et al. (2012)). This paper investigates whether the timing of a SAR acquisition of a flood has an important impact on the calibration of a flood inundation model. As no suitable time series of SAR data exists, we generate a sequence of consistent SAR images through the use of a synthetic framework. This framework uses two available ERS-2 SAR images of the study area, one taken during the flood event of interest, the second taken during a dry reference period. The obtained synthetic observations at different points in time during the flood event are used to calibrate the flood inundation model. The results of this study indicate that the uncertainty of the roughness parameters is lower when the model is calibrated with an image taken before rather than during or after the flood peak. The results also show that the error on the modeled extent is much lower when the model is calibrated with a pre-flood peak image than when calibrated with a near-flood peak or a post-flood peak image. It is concluded that the timing of the SAR image acquisition of the flood has a clear impact on the model calibration and consequently on the precision of the predicted flood extent.

Flood inundation↗

NASA NTRS2021

Evaluation of GEOS Precipitation Flagging for SMAP Soil Moisture Retrieval Accuracy

The precipitation flag in the Soil Moisture Active Passive (SMAP) Level 2 passive soil moisture (L2SMP) retrieval product indicates the presence or absence of heavy precipitation at the time of the SMAP overpass. The flag is based on precipitation estimates from the Goddard Earth Observing System (GEOS) Forward Processing numerical weather prediction system. An error in flagging during an active or recent precipitation event can produce either 1) an overestimation of soil moisture due to short-term surface wetting of vegetation and/or surface ponding (if soil moisture retrieval was attempted in the presence of rain) or 2) an unnecessary nonretrieval of soil moisture and loss of data (if retrieval is flagged due to an erroneous indication of rain). Satellite precipitation estimates from the Integrated Multisatellite Retrievals for GPM (IMERG), version 06, Early Run (latency of ~4 h) precipitationCal product are used here to evaluate the GEOS-based precipitation flag in the L2SMP product for both the 1800 local time (LT) ascending and 0600 LT descending SMAP overpasses over the first five years of the mission (2015–20). Consisting of blended precipitation measurements from the Global Precipitation Mission (GPM) satellite constellation, IMERG is treated as the “truth” when comparing to the GEOS model forecasts of precipitation used by SMAP. Key results include (i) IMERG measurements generally show higher spatial variability than the GEOS forecast precipitation, (ii) the IMERG product has a higher frequency of light precipitation amounts, and (iii) the effect of incorporating IMERG rainfall measurements in lieu of GEOS precipitation forecasts are minimal on the L2SMP retrieval accuracy (determined vs in situ soil moisture measurements at core validation sites). Our results indicate that L2SMP retrievals continue to meet the mission’s accuracy requirement [standard deviation of the unbiased RMSE (ubRMSE) less than 0.04 cu. m/cu. m].

IMERG-precipitationCal↗

NASA NTRS2013

The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling

The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems.

Maggioni, V.↗

NASA NTRS2019

Assimilation of Sentinel 1 and SMAP – Based Satellite Soil Moisture Retrievals into SWAT Hydrological Model: the Impact of Satellite Revisit Time and Product Spatial Resolution on Flood Simulations in Small Basins

In runoff generation process, soil moisture plays an important role as it controls the magnitude of the flood events to the rainfall inputs. In this study, we investigated the ability of a new era of satellite soil moisture retrievals to improve the Soil & Water Assessment Tool (SWAT) daily discharge simulations via soil moisture data assimilation for two small (< 500 km2) and hydrologically different catchments located in Central Italy. We ingested 1) the Soil Moisture Active and Passive (SMAP) Enhanced L3 Radiometer Global Daily 9 km EASE-Grid soil moisture, 2) the Advanced SCATterometer (ASCAT) H113 soil moisture product released within the EUMETSAT Satellite Application Facility on Support to Operational Hydrology and Water Management (H-SAF) which has a nearly daily temporal resolution and sampling of 12.5 km, and 3) a fused ASCAT/Sentinel-1 (S1) satellite soil moisture product named SCATSAR-SWI with temporal and spatial sampling of 1 day and 1 km, respectively into SWAT hydrological model via Ensemble Kalman Filter (EnKF). Different configurations were tested with the aim of exploring the effect of the hydrological regime, the land use conditions, the spatial sampling and the revisit time of the products (which controls the amount of available data to be potentially ingested). Results show a general improvement of SWAT discharge simulations for all products in terms of error and Nash Sutcliffe efficiency index. In particular, we found a relatively good behavior of both the active and the passive products in terms of low flows improvement especially for the catchment characterized by a higher baseflow component. The benefit of the higher spatial resolution of SCATSAR obtained via S1 over ASCAT was small likely due to very challenging areas for S1 retrieval. Eventually, better performances were obtained for the passive product in the more forested catchment. With the aim of exploring the benefit of having more frequent satellite soil moisture observations to be ingested, we tested the performance of the ASCAT product with a reduced temporal sampling obtained by temporally matching ASCAT observations to that of SMAP. The results show a significant reduction of the performance of ASCAT, suggesting that the correction frequency (due to the higher number of observations available) for small catchments is an important aspect for improving flood forecasting as it helps to adjust more frequently the soil moisture pre-storm conditions.

Small basins↗