Search NASA⌕ Search

Engineering topics

Randal D Koster

Publications and source records attributed to Randal D Koster.

At least 19 records

Late-Fall Satellite-Based Soil Moisture Observations Show Clear Connections to Subsequent Spring Streamflow

Because runoff production is more efficient over wetter soils, and because soil moisture has an intrinsic memory, soil moisture information can potentially contribute to the accuracy of streamflow predictions at seasonal leads. In this work, we use surface (0-5 cm) soil moisture retrievals obtained with the National Aeronautics and Space Administration’s Soil Moisture Active Passive satellite instrument in conjunction with streamflow measurements taken within 236 intermediate-scale (2,000 – 10,000 km2) unregulated river basins in the conterminous United States to show that late-fall satellite-based surface soil moisture estimates are indeed strongly connected to subsequent springtime streamflow. We thus show that the satellite-based soil moisture retrievals, all by themselves, have the potential to produce skillful seasonal streamflow predictions several months in advance. In poorly instrumented regions, they could perform better than reanalysis soil moisture products in this regard.

Randal D Koster↗

The Simulation and Subseasonal Forecasting of Hydrological Variables: Insights from a Simple Water Balance Model

Past work has shown that a land surface model’s (LSM’s) implicit (not explicitly coded) relationships between soil moisture and both evapotranspiration (ET) and runoff largely determine the LSM’s hydrological behavior. Here we estimate the relationships that appear to be operating in the real world and compare them to those of the LSM component of a state-of-the-art Earth system model (ESM). The two sets of relationships are determined by calibrating them within a simple water balance model (WBM): once using stream gauge observations from small, unregulated rivers over the eastern half of the U.S., and once using the runoffs generated by the LSM as part of a state-of-the-art atmospheric reanalysis. Hydrological simulations and subseasonal hydrological forecasts performed with the two calibrated versions of the WBM provide two key results. First, the version calibrated to the LSM-generated runoffs does successfully reproduce, to first order, the hydrological behavior of the full LSM within its ESM environment. Second, of the two WBM versions, the one calibrated to the observations reproduces more accurately a broad collection of fully independent streamflow observations as well as a similarly broad collection of in-situ soil moisture measurements. Taken together, the two results suggest that the observations-calibrated ET and runoff efficiency functions do successfully represent, at least to some degree, soil moisture controls over hydrological variability in Nature and can serve as potentially useful targets for further LSM development.

Water Balance Model↗

Skillful Seasonal Forecasts of Land Carbon Uptake in Northern Mid- and High Latitudes

Here we present a first look at the Gross Primary Production (GPP) forecast skill levels achievable with a state-of-the-art subseasonal-to-seasonal (S2S) forecast system. Using NASA’s retrospective S2S ensemble forecast in conjunction with a terrestrial biosphere model, and using an independent, remote sensing-based dataset for validation, we demonstrate an ability to accurately forecast spring-summer carbon uptake at multi-month leads. Averaged across mid-and high latitudes of the Northern Hemisphere land, the GPP forecast initialized on January 1 produces statistically significant skill through summer. The skill achieved, however, is spatially variable, with some regions appearing to extract skill from accurate forecasts of snowpack removal and others extracting skill from the initialization of carbon and nitrogen states. Our results reveal some heretofore unexplored facets of climate predictability and provide a look at what might be possible with future S2S forecast systems that are fully integrated with biogeochemical cycles.

Eunjee Lee↗

Improved Estimates of Pentad Precipitation through the Merging of Independent Precipitation Datasets

Three independent, quasi-global, gridded datasets of precipitation (a rain gauge-based dataset, the satellite-only component of the NASA Integrated Multi-satellitE Retrievals for Global Precipitation Measurement mission [IMERG] Final Run precipitation product, and precipitation estimates derived from NASA Soil Moisture Active Passive [SMAP] soil moisture retrievals), are objectively combined into a single pentad precipitation dataset at 36-km resolution using a unique approach based on extended triple collocation. The quality of each of the four datasets is then evaluated against independent observations. When a global land surface model at 36-km resolution is integrated four times, once utilizing the merged precipitation forcing and once with each of the three contributing datasets, the near-surface soil moisture variations produced with the merged forcing validate best against independent satellite-based soil moisture fields. In addition, the merged dataset is found to be more consistent, relative to each contributor, with estimates of air temperature variations across the globe. The merged dataset thus appears to draw successfully on the complementary strengths of each contributor: the particularly high quality of the rain gauge-based dataset in areas of high gauge density, the more uniform accuracy across the globe of the IMERG data, and the moderate accuracy, particularly in semi-arid regions, of the soil moisture retrieval-based data. Plain Language Summary Obtaining measurements of precipitation across the globe can be challenging. Rain gauges in some ways provide the most accurate measurements, but gauges are absent in many parts of the world, and even where they exist, they only measure precipitation at the gauge itself and therefore may not provide an accurate large-scale average. Satellite-based estimates of precipitation largely overcome these problems, but such data have their own issues, notably a “snapshot” (rather than a time-average) character of the measurements and difficulty associated with interpreting the measured radiances in the presence of complex land surfaces. In the present paper, we use a novel approach to generate a “merged” dataset, one that optimally combines the gauge precipitation information and the satellite-based precipitation information with a third set of estimates derived from soil moisture retrievals. The merged precipitation dataset and each of the three contributors (aggregated here to 5-day averages at a spatial resolution of about 36-km) are then evaluated for consistency with independent geophysical fields. The merged dataset is found to perform best, a clear indication that it takes proper advantage of the complementary strengths of each contributor and, accordingly, that the presented approach for merging the different contributors is indeed viable.

Precipitation↗

The Response of the Amazon Ecosystem to the Photosynthetically Active Radiation Fields: Integrating Impacts of Biomass Burning Aerosol and Clouds in the NASA GEOS Earth System Model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon’s clouds in tempering this impact. The study focuses on a seven-year period (2010-2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity—called here the aerosol diffuse radiation fertilization effect —is found to increase Amazonian Gross Primary Production (GPP) by 2.6% via a 3.8% increase in diffuse PAR (DFPAR) despite a 5.4% decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9% (occurring in August 2010). Consequently, the net primary production (NPP) in Amazon is increased by 1.5%, or ~92 Tg C a-1– equivalent to ~37% of the average carbon lost due to Amazon fires over the seven years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

biomass burning aerosols on ecosystem productivity↗

Asymmetry in Subseasonal Surface Air Temperature Forecast Error with Respect to Soil Moisture Initialization

Soil moisture (W) helps control evapotranspiration (ET), and ET variations can in turn have a distinct impact on 2-m air temperature (T2M), given that increases in evaporative cooling encourage reduced temperatures. Soil moisture is accordingly linked to T2M, and realistic soil moisture initialization has, in previous studies, been shown to improve the skill of subseasonal T2M forecasts. The relationship between soil moisture and evapotranspiration, however, is distinctly nonlinear, with ET tending to increase with soil moisture in drier conditions and to be insensitive to soil moisture variations in wetter conditions. Here, through an extensive analysis of subseasonal forecasts produced with a state-of-the-art seasonal forecast system, this nonlinearity is shown to imprint itself on T2M forecast error in the conterminous United States in two unique ways: (i) the T2M forecast bias (relative to independent observations) induced by a negative precipitation bias tends to be larger for dry initializations, and (ii) on average, the unbiased root-mean-square error (ubRMSE) tends to be larger for dry initializations. Such findings can aid in the identification of forecasts of opportunity; taken a step further, they suggest a pathway for improving bias correction and uncertainty estimation in subseasonal T2M forecasts by conditioning each on initial soil moisture state.

Air Temperature↗

On the Development and Demise of the Fall 2019 Southeast U. S. Flash Drought: Links to an Extreme Positive IOD

Much of the southeast United States experienced record dry conditions during September of 2019, with the area in abnormally dry to exceptional drought conditions growing from 25% at the beginning of the month to 80% by the end of the month. The drought ended just as abruptlydue to above normal rain that fell during the second half of October. In this study we employed MERRA-2 and the GEOS-5 AGCM to diagnose the underlying causes of the drought’s onset, maintenance, and demise. The basic approach involves performing a series of AGCM simulations in which the model is constrained to remain close to MERRA-2 over pre-specified areas that are external to the drought region. The start of the drought appears to have been forced by anomalous heating in the central/western tropical Pacific that resulted in low level anti-cyclonic flow and a tendency for descending motion over much of the southeast. An anomalous ridge associated with a Rossby wave train (emanating from the Indian Ocean region) is found to be the main source of the most intense temperature and precipitation anomalies that develop over the southeast during the last week of September. A second Rossby wave train (emanating from the same region) is responsible for the substantial rain that fell during the second half of October to end the drought. The links to the Indian Ocean Dipole (with record positive values) as well a sa waning El Nino allow some speculation as to the likelihood of similar events occurring in the future.

Siegfried D Schubert↗

The Contributions of Gauge-Based Precipitation and SMAP Brightness Temperature Observations to the Skill of the SMAP Level-4 Soil Moisture Product

Soil Moisture Active Passive (SMAP) mission L-band brightness temperature (Tb) observations are routinely assimilated into the Catchment land surface model to generate Level-4 Soil Moisture (L4_SM) estimates of global surface and root-zone soil moisture at 9-km, 3-hourly resolution with ~2.5-day latency. The Catchment model in the L4_SM algorithm is driven with ¼-degree, hourly surface meteorological forcing data from the Goddard Earth Observing System (GEOS). Outside of Africa and the high latitudes, GEOS precipitation is corrected using Climate Prediction Center Unified (CPCU) gauge-based, ½-degree, daily precipitation. L4_SM soil moisture was previously shown to improve over land model-only estimates that use CPCU precipitation but no Tb assimilation (CPCU_SIM). Here, we additionally examine the skill of model-only (CTRL) and Tb assimilation-only (SMAP_DA) estimates derived without CPCU precipitation. Soil moisture is assessed versus in situ measurements in well-instrumented regions and globally through the Instrumental Variable (IV) method using independent soil moisture retrievals from the Advanced Scatterometer. At the in situ locations, SMAP_DA and CPCU_SIM have comparable soil moisture skill improvements relative to CTRL for the unbiased root-mean-square error (surface and root-zone) and correlation metrics (root-zone only). In the global average, SMAP Tb assimilation increases the surface soil moisture anomaly correlation by 0.10-0.11 compared to an increase of 0.02-0.03 from the CPCU-based precipitation corrections. The contrast is particularly strong in central Australia, where CPCU is known to have errors and observation-minus-forecast Tb residuals are larger when CPCU precipitation is used. Validation versus streamflow measurements in the contiguous U.S. reveals that CPCU precipitation provides most of the skill gained in L4_SM runoff estimates over CTRL.

SMAP↗

Mechanisms Associated with Daytime and Nighttime Heat Waves over the Contiguous United States

Heat waves are extreme climate events that have the potential to cause immense stress on human health, agriculture, and energy systems, so understanding the processes leading to their onset is crucial. There is no single accepted definition for heat waves, but they are generally described as a sustained amount of time where temperature exceeds a local threshold. Multiple different temperature variables are potentially relevant, as high values of both daily maximum (Tmax) and minimum (Tmin) temperatures can be detrimental to human health. In this study, we focus explicitly on the different mechanisms associated with summertime heat waves manifested during daytime versus nighttime hours over the contiguous United States. Heat waves are examined using the National Aeronautics and Space Administration (NASA) Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2). Over 1980–2018, the increase in the number of heat wave days per summer was generally stronger for nighttime heat wave days than daytime heat wave days, with localized regions of significant positive trends. Processes linked with daytime and nighttime heat waves are identified through composite analysis of precipitation, soil moisture, clouds, humidity and fluxes of heat and moisture. Daytime heat waves are associated with dry conditions, reduced cloud cover, and increased sensible heating. Mechanisms leading to nighttime heat waves differ regionally across the US, but they are typically associated with increased clouds, humidity and/or low-level temperature advection. In the Midwest US, enhanced moisture is transported from the Gulf of Mexico during nighttime heat waves.

Extreme events↗

Investigation of the 2016 Eurasia heat wave as an event of the recent warming

This study investigates the physical mechanisms that contributed to the 2016 Eurasian heat wave during boreal summer season (July–August, JA), characterized by much higher than normal temperatures over eastern Europe, East Asia, and the Kamchatka Peninsula. It is found that the 2016 JA mean surface air temperature, upper-tropospheric height, and soil moisture anomalies are characterized by a tri-pole pattern over the Eurasia continent and a wave train-like structure not dissimilar to recent (1980–2016) trends in those quantities. A series of forecast experiments designed to isolate the impacts of the land, ocean, and sea ice conditions on the development of the heat wave is carried out with the Global Seasonal Forecast System version 5. The results suggest that the tri-pole blocking pattern over Eurasia, which appears to be instrumental in the development of the 2016 summer heat wave, can be viewed as an expression of the recent trends, amplified by record-breaking oceanic warming and internal land-atmosphere interactions.

Eunkyo Seo↗

Impact of a Regional US Drought on Land and Atmospheric Carbon

The impacts of drought on regional land and atmospheric carbon are still poorly understood. Here we quantify the impact of a regional US drought on land carbon fluxes (Gross Primary Production, or GPP, and Net Biosphere Production, or NBP) and atmospheric carbon (CO2) by imposing an idealized 3-month meteorological drought in an ensemble of coupled land-atmosphere climate simulations. The imposed drought, applied to the lower Mississippi River Valley (~500,000 km (exp 2)), leads to a 23% GPP reduction in the drought area in the month immediately following the drought’s termination. The drought also caused GPP reductions in some remote areas through drought-induced impacts on remote meteorology, particularly the areas adjacent to the imposed drought. In the remote areas, the induced precipitation changes are responsible for most of the anomalous land productivity. The impact of the drought-induced meteorological anomalies on GPP is greater than that of the CO2 anomalies by at least an order of magnitude. While their impact on GPP is secondary, the drought-induced atmospheric CO2 anomalies near the land surface can be as large as 3.57ppm. The significant CO2 anomalies cover an area up to three times of that of the imposed drought, suggesting that atmospheric transport needs to be considered in the interpretation of drought-induced CO2 anomalies in the atmosphere. The imposed drought also leads to column-averaged CO2 increases of up to 0.78ppm, which is at the edge of the uncertainty from single soundings of current greenhouse gas (GHG) observing satellites.

Eunjee Lee↗

Prediction skill of the 2012 U.S. Great Plains flash drought in Subseasonal Experiment (SubX) models

Rapid-onset droughts, known as flash droughts, can have devastating impacts on agriculture, water resources, and ecosystems. The ability to predict flash droughts in advance would greatly enhance our preparation for them and potentially mitigate their impacts. Here, we investigate the prediction skill of the extreme 2012 flash drought over the U.S. Great Plains at subseasonal lead times (3 weeks or more in advance) in global forecast systems participating in the Subseasonal Experiment (SubX). An additional comprehensive set of subseasonal hindcasts with NASA’s GEOS model, a SubX model with relatively high prediction skill, was performed to investigate the separate contributions of atmospheric and land initial conditions to flash drought prediction skill. The results show that the prediction skill of the SubX models is quite variable. While skillful predictions are restricted to within the first two forecast weeks in most models, skill is considerably better (3-4 weeks or more) for certain models and initialization dates. The enhanced prediction skill is found to originate from two robust sources: 1) accurate soil moisture initialization once dry soil conditions are established, and 2) the satisfactory representation of quasi-stationary cross-Pacific Rossby wave trains that lead to the rapid intensification of flash droughts. Evidence is provided that the importance of soil moisture initialization applies more generally to central U.S. summer flash droughts. Our results corroborate earlier findings that accurate soil moisture initialization is important for skillful subseasonal forecasts and highlight the need for additional research on the sources and predictability of drought-inducing quasi-stationary atmospheric circulation anomalies.

Anthony M DeAngelis↗

Using a Simple Water Balance Framework to Quantify the Impact of Soil Moisture Initialization on Subseasonal Evapotranspiration and Air Temperature Forecasts

Past studies have shown that accurate soil moisture initialization can contribute significant skill to near-surface air temperature (T2M) forecasts at subseasonal leads. The mechanisms by which soil moisture contributes such skill are examined here with a simple water balance-based model that captures the essence of soil moisture behavior in a state-of-the-art subseasonal-to-seasonal (S2S) forecasting system. The simple model successfully transforms initial soil moisture contents into average “forecasted” ET values at 16-30 day lead that agree well, during summer, with the values forecasted by the full NASA GEOS S2S system, indicating that soil moisture initialization dominates over forecasted meteorology in determining ET fluxes at subseasonal leads. When the simple model’s ET anomalies are interpreted in terms of T2M anomalies, a similar conclusion is reached for T2M: soil moisture initialization explains much (about 50% in the eastern half of the continental US) of the T2M anomalies produced by the full GEOS S2S system at 16-30 day lead, and the T2M forecasts produced by the simple model capture about half of the skill attained by the full system. The simple model’s framework is particularly conducive to an analysis of uncertainty in forecasts. Drier soils are generally found to induce larger uncertainty in ET (and thus T2M) forecasts, a result linked to the functional form relating ET to soil moisture in the simple model and verified by an analysis of the ensemble spreads within the forecasts produced by the full GEOS S2S system

Randal D Koster↗

Length Scales of Hydrological Variability as Inferred from SMAP Soil Moisture Retrievals

Hydrological variability at a given location is characterized in part by a horizontal length scale – a measure of how far one can travel from that location and still see similar time variations of a hydrological variable of interest. Here, using Level-2 soil moisture retrievals produced by the NASA Soil Moisture Active Passive (SMAP) mission, we compute global distributions of these length scales for the Northern Hemisphere warm and cold seasons (May-September and November-March, respectively). The length scales show significant spatial and seasonal variability, with, as expected, much larger values (e-folding scales of greater than 500 km) often seen in the cold season, when convective rainfall is less prominent. The SMAP-derived length scales are found to be largely consistent with those derived directly, where possible, from precipitation measurements. This suggests a unique value of the retrievals: outside of well-instrumented areas, satellite-based soil moisture datasets have the potential to provide otherwise unattainable estimates of the horizontal length scales of hydrological variability.

Hydrology↗

Version 4 of the SMAP Level-4 Soil Moisture Algorithm and Data Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with a mean latency of ~2.5 days. The underlying L4_SM algorithm assimilates SMAP radiometer brightness temperature (Tb) observations into the NASA Catchment land surface model using a spatially-distributed ensemble Kalman filter. Version 4 of the L4_SM modeling system includes a reduction in the upward recharge of surface soil moisture from below under non-equilibrium conditions, resulting in reduced bias and improved dynamic range of L4_SM surface soil moisture compared to earlier versions. This change and additional technical modifications to the system reduce the mean and standard deviation of the observation-minus-forecast Tb residuals and overall soil moisture analysis increments while maintaining the skill of the L4_SM soil moisture estimates versus independent in situ measurements; the average, bias-adjusted RMSE in Version 4 is 0.039 m(exp 3) m(exp -3) for surface and 0.026 m(exp 3) m(exp -3) for root-zone soil moisture. Moreover, the coverage of assimilated SMAP observations in Version 4 is near-global owing to the use of additional satellite Tb records for algorithm calibration. L4_SM soil moisture uncertainty estimates are biased low (by 0.01-0.02 m(exp 3) m(exp -3)) against actual errors (computed versus in situ measurements). L4_SM runoff estimates, an additional product of the L4_SM algorithm, are biased low (by 35 mm year (exp -1)) against streamflow measurements. Compared to Version 3, bias in Version 4 is reduced by 46% for surface soil moisture uncertainty estimates and by 33% for runoff estimates.

RMSE↗