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 (2000-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.

Assimilation of SMAP Brightness Temperature Observations in the GEOS Land-Atmosphere Data Assimilation System.

Errors in soil moisture adversely impact the modeling of land-atmosphere water and energy fluxes and, consequently, near-surface atmospheric conditions in atmospheric data assimilation systems (ADAS). To mitigate such errors, a land surface analysis is included in many such systems, although not yet in the currently operational NASA Goddard Earth Observing System (GEOS) ADAS. This article investigates the assimilation of L-band brightness temperature (Tb) observations from the Soil Moisture Active Passive (SMAP) mission in the GEOS weakly coupled land-atmosphere data assimilation system (LADAS) during boreal summer 2017. The SMAP Tb analysis improves the correlation of LADAS surface and root-zone soil moisture versus in situ measurements by ~0.1-0.26 over that of ADAS estimates; the unbiased root-mean-square error of LADAS soil moisture is reduced by 0.002-0.008 m3/m3 from that of ADAS. Furthermore, the global land average RMSE versus in situ measurements of screen-level air specific humidity (q2m) and daily maximum temperature (T2mmax) is reduced by 0.05 g/kg and 0.04 K, respectively, for LADAS compared to ADAS estimates. Regionally, the RMSE of LADAS q2m and T2mmax is improved by up to 0.4 g/kg and 0.3 K, respectively. Improvement in LADAS specific humidity extends into the lower troposphere (below ~700 mb), with relative improvements in bias of 15-25%, although LADAS air temperature bias slightly increases relative to that of ADAS. Finally, the root mean square of the LADAS Tb observation-minus-forecast residuals is smaller by up to ~0.1 K than in a land-only assimilation system, corroborating the positive impact of the Tb analysis on the modeled land-atmosphere coupling.

DroughtCast: A Machine Learning Forecast of the United States Drought Monitor.

Drought is one of the most ecologically and economically devastating natural phenomena affecting the United States, causing the U.S. economy billions of dollars in damage, and driving widespread degradation of ecosystem health. Many drought indices are implemented to monitor the current extent and status of drought so stakeholders such as farmers and local governments can appropriately respond. Methods to forecast drought conditions weeks to months in advance are less common but would provide a more effective early warning system to enhance drought response, mitigation, and adaptation planning. To resolve this issue, we introduce DroughtCast, a machine learning framework for forecasting the United States Drought Monitor (USDM). DroughtCast operates on the knowledge that recent anomalies in hydrology and meteorology drive future changes in drought conditions. We use simulated meteorology and satellite observed soil moisture as inputs into a recurrent neural network to accurately forecast the USDM between 1 and 12 weeks into the future. Our analysis shows that precipitation, soil moisture, and temperature are the most important input variables when forecasting future drought conditions. Additionally, a case study of the 2017 Northern Plains Flash Drought shows that DroughtCast was able to forecast a very extreme drought event up to 12 weeks before its onset. Given the favorable forecasting skill of the model, DroughtCast may provide a promising tool for land managers and local governments in preparing for and mitigating the effects of drought.

Increased high-latitude photosynthetic carbon gain offset by respiration carbon loss during an anomalous warm winter to spring transition.

Arctic and boreal ecosystems play an important role in the global carbon (C) budget, and whether they act as a future net C sink or source depends on climate and environmental change. Here, we used complementary in situ measurements, model simulations, and satellite observations to investigate the net carbon dioxide (CO2 ) seasonal cycle and its climatic and environmental controls across Alaska and northwestern Canada during the anomalously warm winter to spring conditions of 2015 and 2016 (relative to 2010-2014). In the warm spring, we found that photosynthesis was enhanced more than respiration, leading to greater CO2 uptake. However, photosynthetic enhancement from spring warming was partially offset by greater ecosystem respiration during the preceding anomalously warm winter, resulting in nearly neutral effects on the annual net CO2 balance. Eddy covariance CO2 flux measurements showed that air temperature has a primary influence on net CO2 exchange in winter and spring, while soil moisture has a primary control on net CO2 exchange in the fall. The net CO2 exchange was generally more moisture limited in the boreal region than in the Arctic tundra. Our analysis indicates complex seasonal interactions of underlying C cycle processes in response to changing climate and hydrology that may not manifest in changes in net annual CO2 exchange. Therefore, a better understanding of the seasonal response of C cycle processes may provide important insights for predicting future carbon-climate feedbacks and their consequences on atmospheric CO2 dynamics in the northern high latitudes.

Snow depth variability in the Northern Hemisphere mountains observed from space.

Accurate snow depth observations are critical to assess water resources. More than a billion people rely on water from snow, most of which originates in the Northern Hemisphere mountain ranges. Yet, remote sensing observations of mountain snow depth are still lacking at the large scale. Here, we show the ability of Sentinel-1 to map snow depth in the Northern Hemisphere mountains at 1 km² resolution using an empirical change detection approach. An evaluation with measurements from ~4000 sites and reanalysis data demonstrates that the Sentinel-1 retrievals capture the spatial variability between and within mountain ranges, as well as their inter-annual differences. This is showcased with the contrasting snow depths between 2017 and 2018 in the US Sierra Nevada and European Alps. With Sentinel-1 continuity ensured until 2030 and likely beyond, these findings lay a foundation for quantifying the long-term vulnerability of mountain snow-water resources to climate change.

Uncertainty in Soil Moisture Retrievals: an Ensemble Approach using SMOS L-Band Microwave Data.

The uncertainty of surface soil moisture (SM) retrievals from satellite brightness temperature (TB) observations depends primarily on the choice of radiative transfer model (RTM) parameters, prior SM information and TB inputs. This paper studies the sensitivity of several established and experimental SM retrieval products from the Soil Moisture Ocean Salinity (SMOS) mission to these choices at 11 reference sites, located in 7 watersheds across the United States (US). Different RTM parameter sets cause large biases between retrievals. Whereas typical RTM parameter sets are calibrated for SM retrievals, it is shown that a parameter set carefully optimized for TB forward modeling can also be used for retrieving SM. It is also shown that the inclusion of dynamic prior SM estimates in a Bayesian retrieval scheme can strongly improve SM retrievals, regardless of the choice of RTM parameters. The second part of this paper evaluates the ensemble uncertainty metrics for SM retrievals obtained by propagating a wide range of RTM parameters through the RTM, and the relationship with time series metrics obtained by comparing SM retrievals with in situ data. As expected for bounded variables, the total spread in the ensemble SM retrievals is smallest for wet and dry SM values and highest for intermediate SM values. After removal of the strong long-term SM bias associated with the RTM parameter values for individual ensemble members, the remaining anomaly ensemble SM spread shows higher values when SM deviates further from its long-term mean SM. This reveals higher-order biases (e.g. differences in variances) in the retrieval error, which should be considered when characterizing retrieval error. The time-average anomaly ensemble SM spread of 0.037 m3/m3 approximates the actual time series unbiased root-mean-square-difference of 0.042 m3/m3 between ensemble mean retrievals and in situ data across the reference sites.

Global Satellite Retrievals of the Near-Surface Atmospheric Vapor Pressure Deficit from AMSR-E and AMSR2.

Near-surface atmospheric Vapor Pressure Deficit (VPD) is a key environmental variable affecting vegetation water stress, evapotranspiration, and atmospheric moisture demand. Although VPD is readily derived from in situ standard weather station measurements, more spatially continuous global observations for regional monitoring of VPD are lacking. Here, we document a new method to estimate daily (both a.m. and p.m.) global land surface VPD at a 25-km resolution using a satellite passive microwave remotely sensed Land Parameter Data Record (LPDR) derived from the Advanced Microwave Scanning Radiometer (AMSR) sensors. The AMSR-derived VPD record shows strong correspondence (correlation coefficient ≥ 0.80, p-value < 0.001) and overall good performance (0.48 kPa ≤ Root Mean Square Error ≤ 0.69 kPa) against independent VPD observations from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data. The estimated AMSR VPD retrieval uncertainties vary with land cover type, satellite observation time, and underlying LPDR data quality. These results provide new satellite capabilities for global mapping and monitoring of land surface VPD dynamics from ongoing AMSR2 operations. Overall good accuracy and similar observations from both AMSR2 and AMSR-E allow for the development of climate data records documenting recent (from 2002) VPD trends and potential impacts on vegetation, land surface evaporation, and energy budgets.

Estimating snow mass in North America through assimilation of AMSR-E brightness temperature observations using the Catchment land surface model and support vector machines.

To estimate snow mass across North America, brightness temperature observations collected by the Advanced Microwave Scanning Radiometer from 2002 to 2011 were assimilated into the Catchment model using a support vector machine (SVM) as the observation operator and a one-dimensional ensemble Kalman filter. The performance of the assimilation system is evaluated through comparisons against ground-based measurements and reference snow products. In general, there are no statistically significant skill differences between the domain-averaged, model-only ("open loop", or OL) snow estimates and assimilation estimates. The assessment of improvements (or degradations) in snow estimates is difficult because of limitations in the measurements (or products) used for evaluation. It is found that assimilation estimates agree slightly better in terms of root-mean-square error (RMSE) and Nash-Sutcliffe model efficiency with ground-based snow depth measurements than OL estimates in 82% (56 out of 62) of pixels that are colocated with at least two ground-based stations. Assimilation estimates tend to agree slightly better in terms of mean difference with reference snow products over tun-dra snow, alpine snow, maritime snow, and sparsely-vegetated, snow covered pixels. Changes in snow mass via assimilation translate into improvements (e.g.,by 22% on average in terms of RMSE, relative to OL) in cumulative runoff estimates when compared against discharge measurements in 11 out of 13 snow-dominated basins in Alaska. These results suggest that a SVM can potentially serve as an effective observation operator for snow mass estimation within a radiance assimilation system, but a better observational baseline is required to document a statistically significant improvement.

Assimilation of MODIS Snow Cover Fraction Observations into the NASA Catchment Land Surface Model.

The NASA Catchment land surface model (CLSM) is the land model component used for the Modern-Era Retrospective Analysis for Research and Applications (MERRA). Here, the CLSM versions of MERRA and MERRA-Land are evaluated using snow cover fraction (SCF) observations from the Moderate Resolution Imaging Spectroradiometer (MODIS). Moreover, a computationally-efficient empirical scheme is designed to improve CLSM estimates of SCF, snow depth, and snow water equivalent (SWE) through the assimilation of MODIS SCF observations. Results show that data assimilation (DA) improved SCF estimates compared to the open-loop model without assimilation (OL), especially in areas with ephemeral snow cover and mountainous regions. A comparison of the SCF estimates from DA against snow cover estimates from the NOAA Interactive Multisensor Snow and Ice Mapping System showed an improvement in the probability of detection of up to 28% and a reduction in false alarms by up to 6% (relative to OL). A comparison of the model snow depth estimates against Canadian Meteorological Centre analyses showed that DA successfully improved the model seasonal bias from -0.017 m for OL to -0.007 m for DA, although there was no significant change in root-mean-square differences (RMSD) (0.095 m for OL, 0.093 m for DA). The time-average of the spatial correlation coefficient also improved from 0.61 for OL to 0.63 for DA. A comparison against in situ SWE measurements also showed improvements from assimilation. The correlation increased from 0.44 for OL to 0.49 for DA, the bias improved from -0.111 m for OL to -0.100 m for DA, and the RMSD decreased from 0.186 m for OL to 0.180 m for DA.

Estimating Basin-Scale Water Budgets with SMAP Soil Moisture Data.

Soil Moisture Active Passive (SMAP) Level-2 soil moisture retrievals collected during 2015-2017 are used in isolation to estimate 10-day warm-season precipitation and streamflow totals within 145 medium-sized (2,000-10,000 km2) unregulated watersheds in the conterminous United States. The precipitation estimation algorithm, derived from a well documented approach, includes a locally-calibrated loss function component that significantly improves its performance. For the basin-scale water budget analysis, the precipitation and streamflow algorithms are calibrated with two years of SMAP retrievals in conjunction with observed precipitation and streamflow data and are then applied to SMAP retrievals alone during a third year. While estimation accuracy (as measured by the square of the correlation coefficient, r2, between estimates and observations) varies by basin, the average r2 for the basins is 0.53 for precipitation and 0.22 for streamflow. For the subset of 22 basins that calibrate particularly well, the r2 increases to 0.63 for precipitation and to 0.51 for streamflow. The magnitudes of the estimated variables are also accurate, with sample pairs generally clustered about the 1:1 line. The chief limitation to the estimation involves large biases induced during periods of high rainfall; the accuracy of the estimates (in terms of r2 and RMSE) increases significantly when periods of higher rainfall are not considered. The potential for transferability is also demonstrated by calibrating the streamflow estimation equation in one basin and then applying the equation in another. Overall, the study demonstrates that SMAP retrievals contain, all by themselves, information that can be used to estimate large-scale water budgets.

Improved Hydrological Simulation Using SMAP Data: Relative Impacts of Model Calibration and Data Assimilation.

The assimilation of remotely sensed soil moisture information into a land surface model has been shown in past studies to contribute accuracy to the simulated hydrological variables. Remotely sensed data, however, can also be used to improve the model itself through the calibration of the model's parameters, and this can also increase the accuracy of model products. Here, data provided by the Soil Moisture Active/Passive (SMAP) satellite mission are applied to the land surface component of the NASA GEOS Earth system model using both data assimilation and model calibration in order to quantify the relative degrees to which each strategy improves the estimation of near-surface soil moisture and streamflow. The two approaches show significant complementarity in their ability to extract useful information from the SMAP data record. Data assimilation reduces the ubRMSE (the RMSE after removing the long-term bias) of soil moisture estimates and improves the timing of streamflow variations, whereas model calibration reduces the model biases in both soil moisture and streamflow. While both approaches lead to an improved timing of simulated soil moisture, these contributions are largely independent; joint use of both approaches provides the highest soil moisture simulation accuracy.

Spatial and temporal variability of root-zone soil moisture acquired from hydrologic modeling and AirMOSS P-band radar.

The accurate estimation of grid-scale fluxes of water, energy, and carbon requires consideration of sub-grid spatial variability in root-zone soil moisture (RZSM). The NASA Airborne Microwave Observatory of Subcanopy and Subsurface (AirMOSS) mission represents the first systematic attempt to repeatedly map high-resolution RZSM fields using airborne remote sensing across a range of biomes. Here we compare 3-arc-sec (~100-m) spatial resolution AirMOSS RZSM retrievals from P-band radar acquisitions over 9 separate North American study sites with analogous RZSM estimates generated by the Flux-Penn State Hydrology Model (Flux-PIHM). The two products demonstrate comparable levels of accuracy when evaluated against ground-based soil moisture products and a significant level of temporal cross-correlation. However, relative to the AirMOSS RZSM retrievals, Flux-PIHM RZSM estimates generally demonstrate much lower levels of spatial and temporal variability, and the spatial patterns captured by both products are poorly correlated. Nevertheless, based on a discussion of likely error sources affecting both products, it is argued that the spatial analysis of AirMOSS and Flux-PIHM RZSM fields provide meaningful upper and lower bounds on the potential range of RZSM spatial variability encountered across a range of natural biomes.

Characterizing permafrost active layer dynamics and sensitivity to landscape spatial heterogeneity in Alaska.

An important feature of the Arctic is large spatial heterogeneity in active layer conditions, which is generally poorly represented by global models and can lead to large uncertainties in predicting regional ecosystem responses and climate feedbacks. In this study, we developed a spatially integrated modelling and analysis framework combining field observations, local scale (~ 50 m resolution) active layer thickness (ALT) and soil moisture maps derived from airborne low frequency (L+P-band) radar measurements, and global satellite environmental observations to investigate the ALT sensitivity to recent climate trends and landscape heterogeneity in Alaska. Modelled ALT results show good correspondence with in situ measurements in higher permafrost probability (PP ≥ 70%) areas (n = 33, R = 0.60, mean bias = 1.58 cm, RMSE = 20.32 cm), but with larger uncertainty in sporadic and discontinuous permafrost areas. The model results also reveal widespread ALT deepening since 2001, with smaller ALT increases in northern Alaska (mean trend = 0.32 ± 1.18 cm yr-1) and much larger increases (> 3 cm yr-1) across interior and southern Alaska. The positive ALT trend coincides with regional warming and a longer snow-free season (R = 0.60 ± 0.32). A spatially integrated analysis of the radar retrievals and model sensitivity simulations demonstrated that uncertainty in the spatial and vertical distribution of soil organic carbon (SOC) was the largest factor affecting modeled ALT accuracy, while soil moisture played a secondary role. Potential improvements in characterizing SOC heterogeneity, including better spatial sampling of soil conditions and advances in remote sensing of SOC and soil moisture, will enable more accurate predictions of active layer conditions and refinement of the modelling framework across a larger domain.

Global Assessment of the SMAP Level-4 Surface and Root-Zone Soil Moisture Product Using Assimilation Diagnostics.

The Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture and related land surface variables from 31 March 2015 to present with ~2.5day latency. The ensemble-based L4_SM algorithm assimilates SMAP brightness temperature (Tb) observations into the Catchment land surface model. This study describes the spatially distributed L4_SM analysis and assesses the observation-minus-forecast (O-F) Tb residuals and the soil moisture and temperature analysis increments. Owing to the climatological rescaling of the Tb observations prior to assimilation, the analysis is essentially unbiased, with global mean values of ~0.37 K for the O-F Tb residuals and practically zero for the soil moisture and temperature increments. There are, however, modest regional (absolute) biases in the O-F residuals (under ~3 K), the soil moisture increments (under ~0.01 m3 m-3), and the surface soil temperature increments (under ~1 K). Typical instantaneous values are ~6 K for O-F residuals, ~0.01 (~0.003) m3 m-3 for surface (root-zone) soil moisture increments, and ~0.6 K for surface soil temperature increments. The O-F diagnostics indicate that the actual errors in the system are overestimated in deserts and densely vegetated regions and underestimated in agricultural regions and transition zones between dry and wet climates. The O-F auto-correlations suggest that the SMAP observations are used efficiently in western North America, the Sahel, and Australia, but not in many forested regions and the high northern latitudes. A case study in Australia demonstrates that assimilating SMAP observations successfully corrects short-term errors in the L4_SM rainfall forcing.

A Data-Driven Approach for Daily Real-Time Estimates and Forecasts of Near-Surface Soil Moisture.

NASA's Soil Moisture Active Passive (SMAP) mission provides global surface soil moisture retrievals with a revisit time of 2-3 days and a latency of 24 hours. Here, to enhance the utility of the SMAP data, we present an approach for improving real-time soil moisture estimates ("nowcasts") and for forecasting soil moisture several days into the future. The approach, which involves using an estimate of loss processes (evaporation and drainage) and precipitation to evolve the most recent SMAP retrieval forward in time, is evaluated against subsequent SMAP retrievals themselves. The nowcast accuracy over the continental United States (CONUS) is shown to be markedly higher than that achieved with the simple yet common persistence approach. The accuracy of soil moisture forecasts, which rely on precipitation forecasts rather than on precipitation measurements, is reduced relative to nowcast accuracy but is still significantly higher than that obtained through persistence.

Benefits and Pitfalls of GRACE Data Assimilation: a Case Study of Terrestrial Water Storage Depletion in India.

This study investigates some of the benefits and drawbacks of assimilating Terrestrial Water Storage (TWS) observations from the Gravity Recovery and Climate Experiment (GRACE) into a land surface model over India. GRACE observes TWS depletion associated with anthropogenic groundwater extraction in northwest India. The model, however, does not represent anthropogenic groundwater withdrawals and is not skillful in reproducing the interannual variability of groundwater. Assimilation of GRACE TWS introduces long-term trends and improves the interannual variability in groundwater. But the assimilation also introduces a negative trend in simulated evapotranspiration whereas in reality evapotranspiration is likely enhanced by irrigation, which is also unmodeled. Moreover, in situ measurements of shallow groundwater show no trend, suggesting that the trends are erroneously introduced by the assimilation into the modeled shallow groundwater, when in reality the groundwater is depleted in deeper aquifers. The results emphasize the importance of representing anthropogenic processes in land surface modeling and data assimilation systems.

Data Assimilation to extract Soil Moisture Information from SMAP Observations.

This study compares different methods to extract soil moisture information through the assimilation of Soil Moisture Active Passive (SMAP) observations. Neural Network (NN) and physically-based SMAP soil moisture retrievals were assimilated into the NASA Catchment model over the contiguous United States for April 2015 to March 2017. By construction, the NN retrievals are consistent with the global climatology of the Catchment model soil moisture. Assimilating the NN retrievals without further bias correction improved the surface and root zone correlations against in situ measurements from 14 SMAP core validation sites (CVS) by 0.12 and 0.16, respectively, over the model-only skill and reduced the surface and root zone ubRMSE by 0.005 m3 m-3 and 0.001 m3 m-3, respectively. The assimilation reduced the average absolute surface bias against the CVS measurements by 0.009 m3 m-3, but increased the root zone bias by 0.014 m3 m-3. Assimilating the NN retrievals after a localized bias correction yielded slightly lower surface correlation and ubRMSE improvements, but generally the skill differences were small. The assimilation of the physically-based SMAP Level-2 passive soil moisture retrievals using a global bias correction yielded similar skill improvements, as did the direct assimilation of locally bias-corrected SMAP brightness temperatures within the SMAP Level-4 soil moisture algorithm. The results show that global bias correction methods may be able to extract more independent information from SMAP observations compared to local bias correction methods, but without accurate quality control and observation error characterization they are also more vulnerable to adverse effects from retrieval errors related to uncertainties in the retrieval inputs and algorithm. Furthermore, the results show that using global bias correction approaches without a simultaneous re-calibration of the land model processes can lead to a skill degradation in other land surface variables.

Evaluation and enhancement of permafrost modeling with the NASA Catchment Land Surface Model.

Besides soil hydrology and snow processes, the NASA Catchment Land Surface Model (CLSM) simulates soil temperature in six layers from the surface down to 13m depth. In this study, to examine CLSM's treatment of subsurface thermodynamics, a baseline simulation produced subsurface temperatures for 1980-2014 across Alaska at 9-km resolution. The results were evaluated using in situ observations from permafrost sites across Alaska. The baseline simulation was found to capture the broad features of inter- and intra-annual variations in soil temperature. Additional model experiments revealed that: (i) the representativeness of local meteorological forcing limits the model's ability to accurately reproduce soil temperature, and (ii) vegetation heterogeneity has a profound influence on subsurface thermodynamics via impacts on the snow physics and energy exchange at surface. Specifically, the profile-average RMSE for soil temperature was reduced from 2.96°C to 2.10°C at one site and from 2.38°C to 2.25°C at another by using local forcing and land cover, respectively. Moreover, accounting for the influence of soil organic carbon on the soil thermal properties in CLSM leads to further improvements in profile-average soil temperature RMSE, with reductions of 16% to 56% across the different study sites. The mean bias of climatological ALT is reduced by 36% to 89%, and the RMSE is reduced by 11% to 47%. Finally, results reveal that at some sites it may be essential to include a purely organic soil layer to obtain, in conjunction with vegetation and snow effects, a realistic "buffer zone" between the atmospheric forcing and soil thermal processes.