Next Article in Journal
Transferring Deep Convolutional Neural Networks for the Scene Classification of High-Resolution Remote Sensing Imagery
Next Article in Special Issue
Detection and Delineation of Localized Flooding from WorldView-2 Multispectral Data
Previous Article in Journal
Retrieval of Soil Water Content in Saline Soils from Emitted Thermal Infrared Spectra Using Partial Linear Squares Regression
Previous Article in Special Issue
Flood Hazard Mapping Combining Hydrodynamic Modeling and Multi Annual Remote Sensing data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assimilation of GRACE Terrestrial Water Storage Observations into a Land Surface Model for the Assessment of Regional Flood Potential

1
Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91125, USA
2
Department of Earth System Science, University of California, Irvine, CA 92697, USA
3
Centro de Estudios Avanzados en Zonas Áridas, La Serena Región de Coquimbo, Chile
4
NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
5
Earth System Science Interdisciplinary Center, University of Maryland, College Park, MD 20742, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2015, 7(11), 14663-14679; https://doi.org/10.3390/rs71114663
Submission received: 21 July 2015 / Revised: 21 July 2015 / Accepted: 27 October 2015 / Published: 5 November 2015
(This article belongs to the Special Issue Remote Sensing in Flood Monitoring and Management)

Abstract

:
We evaluate performance of the Catchment Land Surface Model (CLSM) under flood conditions after the assimilation of observations of the terrestrial water storage anomaly (TWSA) from NASA’s Gravity Recovery and Climate Experiment (GRACE). Assimilation offers three key benefits for the viability of GRACE observations to operational applications: (1) near-real time analysis; (2) a downscaling of GRACE’s coarse spatial resolution; and (3) state disaggregation of the vertically-integrated TWSA. We select the 2011 flood event in the Missouri river basin as a case study, and find that assimilation generally made the model wetter in the months preceding flood. We compare model outputs with observations from 14 USGS groundwater wells to assess improvements after assimilation. Finally, we examine disaggregated water storage information to improve the mechanistic understanding of event generation. Validation establishes that assimilation improved the model skill substantially, increasing regional groundwater anomaly correlation from 0.58 to 0.86. For the 2011 flood event in the Missouri river basin, results show that groundwater and snow water equivalent were contributors to pre-event flood potential, providing spatially-distributed early warning information.

Graphical Abstract

1. Introduction

Hydrologic predictability depends on knowledge of initial hydrologic conditions. While current precipitation may be the most visible hydrological driver of streamflow, antecedent states of groundwater and soil moisture are less visible, but of critical importance. These conditions serve as major factors in runoff generation and are fundamental obstacles in the accurate prediction of river discharge [1,2,3,4,5]. Traditionally, the analysis and prediction of hydrologic extremes have involved statistical assessments of long-term (30+ years) climatic and hydrological datasets to determine flood thresholds, severity, and recurrence intervals. Contemporary extremes analysis and monitoring has not deviated much from its roots, however past events are not always indicative of future conditions, especially under current climatic trends [6]. Efforts to improve hydrological prediction and reduce uncertainty can now leverage a wealth of new tools for measurement and estimation of land surface states, namely satellite remote sensing products and advanced numerical models [7].
The joint NASA/DLR Gravity Recovery and Climate Experiment (GRACE) satellite mission [8] has been an important step forward in observing terrestrial water storage with global coverage. GRACE provides monthly data on changes in the Earth’s gravity field, which we assume to be correlated with the movement of water through the Earth’s system at specific temporal and spatial scales. GRACE has been providing monthly gravity field solutions since April 2002, and has proved an effective tool to observe changes in the mass of ice sheets [9], snow mass [10], groundwater storage [11,12,13], surface water storage [14], and to hydrologic drought characterization [15,16,17].
GRACE has also contributed to flood analysis [18,19]. Generally, the land surface modulates the connection between precipitation and runoff generation through two basic mechanisms: infiltration limitation and saturation excess. During floods, one of these two mechanisms is typically identified as a driver based on the pre-flood conditions and the nature of the precipitation event. However it is difficult to ascertain the relative contribution of each. For example, preceding flood occurrence, water tables can gradually rise, saturating soils at depth and creating a dangerous precondition for sudden intense runoff generation. Though it has been a primary research topic in hydrology for years, three recent research papers have explored non-linearities in the relationship between river discharge and basin-averaged water storage using GRACE observations [19,20,21]. The idea that extreme values in terrestrial water storage have the ability to augment runoff generation during future flood events at several months lead-time has been previously referred to as “flood potential” [18]. This is distinctly different from river flood prediction, as streamflow forecasting is not employed. However, the monitoring of flood potential could provide valuable anticipatory capacity for water/natural resource managers.
GRACE terrestrial water storage anomaly (TWSA) data can contribute to an improved understanding of soil moisture and groundwater contributions to future runoff generation, but the utility of GRACE data for operational analysis or actionable prediction is challenging. Three prevailing limitations are: (1) coarse spatial resolution (one arc-degree for the level-three data product [22]); (2) the aggregated observation of multiple water storage components (i.e., snow, soil moisture, and groundwater) as a single integrated value for each grid cell; and (3) latency in data product processing and release (2–4 months lag) [23]. These obstacles make many applications difficult, as water resource management tends to occur at watershed scales and on daily to weekly timescales.
Computer modeling of land surface processes offers the potential for a level of temporal and spatial coverage beyond the capacities of most in situ observation networks. However, land surface model simulation accuracy suffers from the assumptions and simplifications that accompany an empirical representation of reality [24]. Even though many land surface models (LSMs) include detailed descriptions of vegetation and soil moisture parameters, the interactions between deeper, subsurface (i.e., groundwater) and surface water, and the accurate generation of lateral surface and subsurface flows are typically disregarded. This can cause models to have predictive inaccuracies in situations where subsurface contributions to streamflow generation are important [25].
Data assimilation is the process of merging measurements with model predictions to maximize spatial and temporal coverage, consistency, resolution, and accuracy. It addresses weaknesses in the representation of land-surface model states at the grid scale level, and can create more realistic predictions of hydrological components like runoff, river flow, and groundwater [26,27]. In hydrology, data assimilation is often assumed to be an improvement over raw observations or a model without assimilation, such that the assimilation adds unique information from both sources in order to achieve a better combined product.
The assimilation of GRACE terrestrial water storage observations has been shown to improve model outputs in several cases, decreasing RMS errors, and increasing the correlation between simulated and observed variables. Zaitchik et al. [28] discusses several potential advantages of GRACE data assimilation. For example, the pairing of GRACE estimates with a watershed-defined Catchment Land Surface Model (CLSM) domain allows for area-accurate assimilation within hydrologically-defined basins. Unlike many land surface models, CLSM provides a representation of surface and subsurface hydrologic states, including groundwater, which can be applied directly as an equivalent state variable to the integrated GRACE observations (TWSA). Additionally, the spatial resolution of modeled segments provide a physically-based mechanism for distributing GRACE observations at spatial scales smaller than the native resolution. Houborg et al. [15] suggest that data assimilation may be the key to realizing the full potential of intrinsically-coarse GRACE TWSA observations. The outcomes of these studies included moderate, but statistically significant, improvements to the hydrological modeling skill of the CLSM across major parts of the United States. Li et al. [16] also recognized the potential for GRACE data to improve the simulation of land surface processes for drought monitoring applications.
In May and June of 2011, a large flood crest moved down the Mississippi from a combination of snowmelt and heavy spring rains in the Missouri river basin. Streamflow in the Missouri were characterized as a 500-year event. Red River Landing, near the mouth of the Mississippi, went above flood stage on 19 March and did not fall below flood stage until 25 June [29].
We build upon previous work [19] that demonstrated record-length maxima in GRACE TWSA over the Missouri river basin during the months of June and July 2011 coincident with peak stream flow. This basin contributed intense snowmelt and groundwater discharge to the high river levels preceding the storm events of March–July 2011. We apply the National Climate Assessment “Northern Plains” domain to isolate the study region of the Missouri River basin (Figure 1). This entire region drains through a single point in the Missouri river and, as such, is suitable to constrain the water balance there. Following the river gage selection procedure described in Reager et al. [19], we use USGS stream gage #6818000 at St. Joseph, Missouri (Figure 1).
Figure 1. Study domain, stream gage near the mouth of the Missouri river, and monthly flow time series.
Figure 1. Study domain, stream gage near the mouth of the Missouri river, and monthly flow time series.
Remotesensing 07 14663 g001
Here we examine disaggregated water storage outputs from the GRACE-assimilated CLSM (CLSM-DA) to explore the utility of GRACE assimilation for the characterization of long lead-time river basin flood potential for an exceptional event in the Mississippi river basin. We determine the effects of assimilation by exploring disparities between CLSM-DA outputs, GRACE observations, and CLSM Open-Loop runs (CLSM-OL; obtained by running the CLSM without any data assimilation). We evaluate the improvement in the model state prediction by validation with groundwater well observations.
To provide a brief overview of the organization of the paper: In Section 2, we present results on the model validation, on the disaggregation of the water storage components, on the spatial distribution of water storage components before and during flooding, and on the contributions of water storage components to event generation. In Section 2.5 we discuss the implications of the results. Section 3 is the experimental section, and contains detailed information on the experiment setup, on the CLSM model structure, and on the assimilation methodology.

2. Results

2.1. Validation and Performance Assessment

USGS well observations were utilized for 16 monitoring well locations through the study domain (Figure 2). Well locations were chosen based on length and completeness of record, and for a geographical representation across the entire study area. Well data were processed by removing their 2003–2014 record-length mean water table depth (from 12 to 65 feet in shallow unconfined aquifers), removing the record-length seasonal climatology and then normalizing each well time series by its record-length standard deviation. This results in 14 time series of well variability that can be averaged to generate a single regional time series. USGS well stations used are: 480034105195404; 483127103373102; 472203111112602; 46540100222101; 463417099271002; 450524112380701; 473051099093601; 430027102311806; 420757104024701; 434329096521201; 420757104024701; 423730098560001; 423148098300601; 400155101521302; 392045093302401; and 383644091124901.
Figure 2. USGS well observation locations and normalized time series (grey lines), as well as the average of the normalized time series (red line).
Figure 2. USGS well observation locations and normalized time series (grey lines), as well as the average of the normalized time series (red line).
Remotesensing 07 14663 g002
The aggregated well time-series was then compared to the model groundwater estimates for the domain under the same processing. While this does not allow for evaluation of a realistic water table, the model representation of groundwater storage is empirical and is best evaluated for its relative variability and timing as opposed to absolute accuracy.
Figure 3 shows maps of GRACE observations and model outputs for March 2011, the month of highest snow melt and high water storage loading preceding flood initiation in May 2011 [19]. The bottom right panel shows the difference in terrestrial water storage before and after assimilation in the model for that month. In the assimilated model, there are generally wetter antecedent conditions in the study watershed several months prior to the spring floods. In this month, GRACE increases the model domain water storage by roughly 40 km3, or an average of 0.5 cm equivalent.
Figure 3. Pre-flood terrestrial water storage anomaly in March 2011: (a) GRACE observations; (b) CLSM open-loop simulation; (c) CLSM-DA assimilation simulation; (d) the difference between the two model simulations, CLSM-DA minus CLSM.
Figure 3. Pre-flood terrestrial water storage anomaly in March 2011: (a) GRACE observations; (b) CLSM open-loop simulation; (c) CLSM-DA assimilation simulation; (d) the difference between the two model simulations, CLSM-DA minus CLSM.
Remotesensing 07 14663 g003
Figure 4a shows the comparison of the region-averaged GRACE time-series with two model outputs after normalization. We see that the CLSM-DA run matches the GRACE observations more closely than the open-loop (CLSM). Figure 4b shows the comparison of model groundwater storage anomaly (i.e., not total water storage) with the USGS well observations after removal of the record-length climatology and normalization. The monthly climatology was removed from each of the time series to demonstrate the inter-annual persistence of above- and below-normal storage anomalies. Table 1 contains the summary statistics for the two comparisons.
Figure 4. (a) The effect of assimilation on Terrestrial Water Storage in the Northern Plains domain, showing the region-averaged GRACE TWSA, the open-loop model (CLSM), and the GRACE-assimilated model (CLSMDA); (b) well validation for the Northern Plains region, showing the groundwater component only from the open-loop model (CLSM), the GRACE-assimilated model (CLMSDA), and independent well observations averaged over the region (climatology removed in each).
Figure 4. (a) The effect of assimilation on Terrestrial Water Storage in the Northern Plains domain, showing the region-averaged GRACE TWSA, the open-loop model (CLSM), and the GRACE-assimilated model (CLSMDA); (b) well validation for the Northern Plains region, showing the groundwater component only from the open-loop model (CLSM), the GRACE-assimilated model (CLMSDA), and independent well observations averaged over the region (climatology removed in each).
Remotesensing 07 14663 g004
Table 1. Goodness of fit statistics for assimilation effects and validation of assimilation improvement of domain-averaged groundwater state estimation with well observations. All statistics are significant at the 95% confidence level.
Table 1. Goodness of fit statistics for assimilation effects and validation of assimilation improvement of domain-averaged groundwater state estimation with well observations. All statistics are significant at the 95% confidence level.
ModelMetricGRACE TWSAUSGS Groundwater Anomaly
CLSM Open-Loopr0.880.58
RMSE33 cm2.32 (normalized)
CLSM DAr0.950.86
RMSE25 cm1.81 (normalized)

2.2. Vertical Disaggregation

Figure 5 shows the vertically disaggregated components of the residual CLSM-DA time series (TWS, groundwater, RZMC, SFMC, and SWE) and GRACE residual TWSA (gold squares), for the Northern Plains region. The seasonal climatology of each time series has been removed in order to isolate inter-annual variability. Overall, CLSM-DA TWSA deficits are comparable to those from GRACE, though several events are underestimated in the study area.
The plots reveal timing differences between changes in water stores. For example, there are occasions where RZMC residuals lead groundwater changes by one to three months. RZMC also shows the most monthly variation between positive and negative values out of the storage components (i.e., as a small reservoir, it dries and wets quickly). In the Northern Plains, RZMC recovers for one to two months on several occasions during extended deficit periods while the other stores remained negative. Peak SWE residuals were often matched or immediately followed by an instance of surplus RZMC, as snow melt filled the root-zone.
For the 2011 flood event, CLSM-DA storage beneath the soil layer (i.e., “groundwater”) was a large portion of the signal before flood occurrence. Post event flood reports [19] showed that in the Missouri River basin (represented here as the Northern Plains study area), the groundwater storage climbed to record levels, accompanied by record high snow water equivalent a few months before. It is likely that record snowmelt and a high water table contributed to flooding in this event [19], and that mechanism seems consistent with the assimilation results.
Figure 5. CLSM-DA disaggregated, residual terrestrial water storage time series (i.e., climatology removed) for the Northern Plains region. Time period is from January 2003–April 2014. Negative values designate deficits and positive surplus. Variables on the left axis are: CLSM-DA total water storage (green shading), Below RZMC or “groundwater” (black), RZMC (green), and monthly GRACE TWSA (gold squares). Variables on the right axis are: CLSM-DA SFMC (red dashed) and SWE (blue dashed). Units are centimeters of equivalent water storage.
Figure 5. CLSM-DA disaggregated, residual terrestrial water storage time series (i.e., climatology removed) for the Northern Plains region. Time period is from January 2003–April 2014. Negative values designate deficits and positive surplus. Variables on the left axis are: CLSM-DA total water storage (green shading), Below RZMC or “groundwater” (black), RZMC (green), and monthly GRACE TWSA (gold squares). Variables on the right axis are: CLSM-DA SFMC (red dashed) and SWE (blue dashed). Units are centimeters of equivalent water storage.
Remotesensing 07 14663 g005
The structure of the CLSM model is such that a terrestrial water storage variable is estimated and assumed to contain snow water equivalent, soil moisture, and some residual value of subsurface water (see Section 3.2). This residual value is then assumed to be “groundwater” though, realistically, for regions with large or dynamic surface water, the residual term would also include those components, especially in the case that they are highly correlated, temporally.

2.3. Spatially Distributed View of Event Progression

Figure 6 shows the event progression for the Northern Plains floods in a series of maps for each state variable. The colors depict the number of standard deviations each grid cell is away from its record-length mean value. As could be determined from the time series, SWE achieves a record maximum in the months preceding the 2011 floods in the Missouri basin region. There are SWE values in the domain of more than 4-sigma for the month of February. This indicates the non-Gaussian nature of the time series, and that solid water storage was considerably higher than normal. Snow in the Missouri basin was reported as a record high in the months preceding the flood [29].
Groundwater and soil moisture both increase substantially and rapidly during the months after snowmelt preceding flood occurrence (April 2011). It is these two variables that indicate regions of enhanced flood potential. Visually, these two variables are highly correlated in space and time, which may be due to model structural limitations as much as it is to any physical mechanism.
Figure 6. Time progression of the 2011 flood event, showing (left) snow water equivalent; (middle) root-zone soil moisture; and (right) groundwater from the CLSMDA simulation for October 2010 through June 2011.
Figure 6. Time progression of the 2011 flood event, showing (left) snow water equivalent; (middle) root-zone soil moisture; and (right) groundwater from the CLSMDA simulation for October 2010 through June 2011.
Remotesensing 07 14663 g006

2.4. Contributions to Event Generation

Figure 7 shows the hysteresis between CLSMDA state variables and streamflow at the USGS gage. Following Sproles et al. [21], we use these relationships to understand the dynamics of the basin and how individual state variables affect a discharge response. In a normal (average) year (Figure 7a), there is nearly an equivalent curvature between soil-streamflow hysteresis and groundwater-streamflow hysteresis.
Figure 7. Hystereses between stream flow at the mouth and snow water equivalent (SWE), root-zone soil moisture (RZMC), or groundwater (GW) in the CLSMDA during (a) the climatological average of all years; and (b) the 2011 flood year. Raw GRACE data are shown for comparison. Hystereses run from October to September.
Figure 7. Hystereses between stream flow at the mouth and snow water equivalent (SWE), root-zone soil moisture (RZMC), or groundwater (GW) in the CLSMDA during (a) the climatological average of all years; and (b) the 2011 flood year. Raw GRACE data are shown for comparison. Hystereses run from October to September.
Remotesensing 07 14663 g007
During the flood year, all hysteresis loops expand along both axis due to increased storage in all state variables (horizontal) and increased discharge (vertical). However, the groundwater and snow hysteresis loops separates horizontally, demonstrating exceptional storage in the months preceding flood occurrence. This heavy skewing in the positive-storage direction shows the extent of the positive anomaly in groundwater preceding flood (nearly four standard deviations), and supports the hypothesis of groundwater as a major contributor to event generation.

2.5. Discussion

A comparison of hydrologic variables from direct observation, CLSM data-assimilated outputs, and CLSM Open-Loop runs was conducted to evaluate the accuracy and utilization of CLSM-DA for flood potential analysis. Data assimilation provides spatial resolutions currently unattainable with GRACE alone. Correlation results show that data assimilation improved CLSM TWSA simulations in the study region. CLSM-DA is a viable dataset for hydrologic analyses based on several potential benefits. This analysis suggests a useful tool for flood potential assessment—one that identifies and quantifies terrestrial water storage extremes with higher resolution (compared with GRACE alone), provides a separation of terrestrial storage components, and provides a physically-based means to obtain terrestrial water storage information beyond the latest GRACE data release.
Data assimilation improved CLSM groundwater anomaly correlations significantly with aggregated well data. The ability to represent groundwater variability in the open-loop model is limited (r = 0.58), as the model does not possess enough information from forcing and physics alone. In this case, the assimilation helps to improve model accuracy in groundwater simulation (r = 0.86). The difference in the before- and after-assimilation correlation coefficients is statistically significant at the 99.8th percentile.
Generally, storage observed by GRACE versus the CLSM includes more terrestrial features, such as surface reservoirs, stream flow, human impacts on the hydrologic system (e.g., dams, river divergences), and subsurface variations below the model’s maximum depth (i.e., below two meters). Assimilation corrects the CLSM towards GRACE but still does not address shortcomings in the model structure that could lead to a degradation of simulation accuracy without a constant constraint on accuracy (e.g., in a hydrologic forecasting application). Still, the model provides a means to generate physically-based information beyond the last assimilation time step that is presumably (based on these results) better than the open-loop model alone.
Coarse resolution makes it challenging to apply GRACE observations towards local water resource management. CLSM-DA maps improve upon this, allowing us to identify key areas being affected by flood or drought conditions. Additionally, maps of proceeding months or seasons can help identify antecedent conditions that can lead to event generation in subsequent periods. For the 2011 event, we see very high snow water equivalent (3–4 months preceding peak discharge) and soil moisture (2–3 months preceding peak discharge) that indicate the ability of this methodology to provide event early warning. Assuming assimilation were performed and then the model was run past the assimilation time step, this would certainly provide a means to generate near real-time early warning information for regional flood potential.
With vertical disaggregation from the CLSM, we were able to assess relationships between SWE, SFMC, RZMC, and groundwater. There is evidence of one-to-two month lags between storage surplus appearance in SWE, surface soils, and the groundwater layer, as water takes time to infiltrate into the subsurface. This result highlights the impact that snowmelt can have on the water table, leading to dangerous pre-conditioning for large flood event generation.
The large basin storage-discharge relationship is highly non-linear in the Missouri, and is similar to results found in other studies [19,20,21]. The hysteresis relationship is dynamic, with discharge responding to changes in basin storage at monthly time steps. This behavior is different than the traditional assumption in hydrology of a “linear reservoir”, or a near-linear proportional relationship between storage changes and discharge changes. The results shown in Figure 7 suggest that this may be especially true for the groundwater-discharge relationship: groundwater reaches extreme (3-sigma) values several months before flood occurrence, and the assimilated model seems to capture this behavior.
There may be several reasons for the existence of such strong hysteresis behavior. Riegger and Tourian [20] propose that the hysteresis effect is singularly attributable to temperature (i.e., the freeze-thaw process), which acts to load a basin with snow and frozen soils in the winter months, loading basin storage without an equivalent increase in discharge. Then, when the spring melt arrives, discharge accelerates exponentially until losses overcome inputs and storage begins to draw down. However, this explanation can be controlled for in some basins (as those authors demonstrate), resulting in a hypothetical linear relationship again between storage and discharge.
We would offer, based on the present evidence, that controlling for snow accumulation should only affect the rising limb of the hysteresis and does not really offer a full description of what is happening in basins like the Missouri, especially with respect to hysteresis in subsurface water storage. Snow and frozen soils melt fairly quickly in the Missouri, especially in the low elevation prairie (within ~3 months, Figure 6), and freeze-thaw may not provide a full explanation for the temporal displacement seen between storage and discharge. These results instead suggest a mechanism by which the storage and discharge are decoupled from the typical baseflow-driven linear reservoir process. In that case, a non-linearity in storage and discharge variability may be attributable to variable water retention in deeper soils. These are interesting phenomena for future research.

3. Experimental Section

3.1. GRACE Terrestrial Water Storage Anomaly

The NASA/DLR Gravity Recovery And Climate Experiment (GRACE) mission comprises a pair of Earth-orbiting satellites flying in tandem at over 400 km altitude, spaced about 220 kilometers apart. Distance variations between the two satellites are measured using a K-band Ranging System (KBR), which provides precise (within 10 μm) measurements of the distance changes between the two satellites caused by temporal and spatial fluctuations in the Earth’s gravity field [8]. The measurements are used to determine changes in the Earth’s mass distribution at horizontal resolutions greater than 150,000 km2, with higher measurement accuracy at larger spatial scales [30,31]. The monthly to decadal temporal changes in the gravity field are attributed to mass redistribution in the atmosphere, ocean, continents, and solid earth. After isolation and correction of “unwanted” signals, these measurements are assumed to represent the movement of water mass over time, referred to as terrestrial water storage anomalies (TWSA).
The GRACE dataset used for this project is the monthly, global, one-degree gridded, scaled GRACE land data product, processed by the University of Texas Center for Space Research (CSR; version CSR-RL05), and available for download at GRACE Tellus. [32]. The time period of data for this project is from April 2002 to December 2014. Swenson and Wahr [31] and Wahr et al. [33] provide general post-processing details and Landerer and Swenson [22] provides details on signal restoration, scaling, and regional error calculation. Several corrections to the data have been applied, including for postglacial rebound.
For temporal analyses, GRACE data were averaged in the United States National Climate Assessment (NCA) Northern Plains region. An area mask was used to spatially average monthly storage components. Time series and maps of water storage differences from the climatology were calculated following methodology in Thomas et al. [17], where surplus and deficit storage is determined by removing the monthly climatology from the original storage time series. In the Missouri basin, surface water storage variability could be an important component of the total water storage anomaly seen by GRACE (i.e., on the order of 5%–10% [34]). The surface water representation is in the level-three GRACE product based on the CESM RTM, and does not currently represent active lakes and represents no man-made reservoirs. While these surface water signals are contained within the raw GRACE observations (as a relatively small signal distributed across the resolution), they are essentially not contained in an accurate fashion in the scaled product. As such, we assume that the GRACE observations represent snow, soil moisture, and groundwater variability (i.e., the three primary reservoirs) and are suited to adjust those fields in the model.

3.2. CLSM Physical Processes

The Catchment Land Surface Model (CLSM) was developed at the NASA Goddard Space Flight Center and is a physically-based, apportioned land surface model [35]. The sub-grid horizontal structure of a rectangular atmospheric grid is divided into topographically-defined catchments with an average area of approximately 4000 km2. Water is redistributed spatially and vertically based on the topography of each basin or watershed. The model’s hydrologic processes are based on each catchment’s topographical statistics.
The strength of CLSM lies in its ability to simulate unconfined groundwater storage variations, typically 2–3 m below the terrestrial surface. This groundwater layer is an essential variable to be included in a land model if we are to generate terrestrial water storage variations that are analogous to those measured by GRACE satellites. The TOPMODEL framework of Beven and Kirkby [36] is utilized to establish root zone soil moisture distributions from the morphology of the catchment and from bulk soil moisture prognostic variables. TOPMODEL formulations allow for the calculation of both the saturated fraction and the baseflow of a watershed in addition to translating the mean watershed water table depth and a probability density function of combined topographic and soil properties [36,37]. At the core of TOPMODEL are three essential assumptions: (1) the water table is nearly parallel to the soil surface so that the local hydraulic gradient is close to the tangent of the local slope angle; (2) the saturated hydraulic conductivity declines exponentially with depth; and (3) the water table is recharged at a spatially uniform and steady rate with respect to the response timescale of the watershed. As such, recharge and baseflow are balanced in a series of steady states.
CLSM follows the founding work of Famiglietti and Wood [38], which was the first to include the TOPMODEL approach in a land surface model framework to relate water table distribution to basin topography. Sub-catchment heterogeneity of soil moisture is modeled by dividing the catchment into dynamic fractions of saturated, unsaturated, and wilting areas, each ruled by equations appropriate for its soil moisture status [28]. The distribution of soil moisture in the root zone allows the catchment to be portioned into these distinct regimes [35,39]. These surface fractions exchange heat separately with the underlying ground layer, which also exchange heat with lower ground layers. Figure 8 shows a depiction of model components [40]. Please refer to Ducharne et al. [39] and Koster et al. [35] for more detailed discussion on the philosophy, configuration, and implementation of the CLSM. Note that a transient data gap in the model forcing at the US-Canada border produces an artifact in the model storage fields which is visible in spatial maps for certain months (Figure 6).
Figure 8. CLSM model structure and state variables including: (1) soil moisture; (2) root zone; (3) surface; (46) snow; and (7) canopy. From Foreman et al. [40].
Figure 8. CLSM model structure and state variables including: (1) soil moisture; (2) root zone; (3) surface; (46) snow; and (7) canopy. From Foreman et al. [40].
Remotesensing 07 14663 g008
CLSM’s primary prognostic variable is the catchment deficit, defined as the average depth of water that would need to be added to bring the catchment to saturation [39]. The equilibrium vertical distribution of soil moisture is then diagnosed on the basis of the catchment deficit and soil parameters. This distribution includes an implicit water table, located at the depth of equilibrium saturation. Additional prognostics include reservoirs of surface excess moisture and root zone excess moisture that yield an approximate representation of non-equilibrium vertical conditions such as infiltration fronts. The surface excess moisture reservoir, comprised of the first two centimeters of the soil layer, is small relative to both root zone excess (the soil layer from 0–100 cm) and the catchment deficit.

3.3. CLSM-GRACE Assimilation

For the assimilation, the model-generated terrestrial water storage moisture components are corrected toward the GRACE observational estimate with the degree of correction determined by the levels of error associated with each using an Ensemble Kalman Smoothing Filter method (EnKS). Zaitchik et al. [28] developed the EnKS scheme to specifically assimilate GRACE observations into the CLSM in the Mississippi basin. Monthly GRACE anomaly fields are converted to absolute values by adding the time-mean total water storage field from the CLSM output. Model ensembles are generated by perturbing the model forcing (precipitation and radiative) as well as model prognostic variables [28]. Assimilation increments are calculated based on the relative uncertainty in the model and the observations where a two-step smoother is applied to handle GRACE’s monthly temporal resolution. These increments are applied directly to the column-integrated prognostic variable (the catchment deficit) and the primary non-equilibrium prognostic (the root zone excess moisture), without need for arbitrary vertical disaggregation. The CLSM-DA time series is from January 2003 to December 2014. Outputs are posted onto 0.25-degree grid cells for the intercontinental United States (domain: −126.875 23.875, −66.125 51.625), including portions of Canada and Northern Mexico. This gridded analysis is an interpolation of catchment tiles to a model grid.

3.4. Output Processing

Original CLSM outputs were averaged from millimeters of daily water storage to monthly storage anomalies in centimeters. CLSM output variables used in this analysis are: total grid cell water storage (TWS), snow water equivalent (SWE), surface moisture content (SFMC, 0–20 cm), and root zone moisture (RZMC, 0–100 cm). We also calculate storage below the root zone by subtracting SWE and RZMC from the TWS time series. The term “groundwater” is used to represent storage below the RZMC. TWS is the only variable used from the CLSM open-loop (i.e., non-assimilated) run.

Acknowledgments

This work was performed at the Jet Propulsion Laboratory (JPL), California Institute of Technology under contract with NASA. Funding sources include the NASA GRACE Science Team program and the JPL Research and Technology Development Water Initiative.

Author Contributions

John T. Reager, Alys C. Thomas and Matthew Rodell conceived of this study. Bailing Li and Hiroko K. Beaudoing conducted the model simulations and performed the assimilation. Alys C. Thomas and John T. Reager performed the analysis. All authors contributed to writing the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bierkens, M.F.P.; van Beek, L.P.H. Seasonal predictability of European discharge: NAO and hydrological response time. J. Hydromet. 2009, 10, 953–968. [Google Scholar] [CrossRef]
  2. Koster, R.D.; Reichle, R.H. Skill in streamflow forecasts derived from large-scale estimates of soil moisture and snow. Nat. Geosci. 2010, 3, 613–616. [Google Scholar] [CrossRef]
  3. Koster, D.R.; Mahanama, S.P.; Yamada, T.J.; Balsamo, G.; Berg, A.A.; Boisserie, M.; Dirmeyer, P.A.; Doblas-Reyes, F.J.; Drewitt, G.; Gordon, C.T.; et al. Contribution of land surface initialization to subseasonal forecast skill: First results from a multi-model experiment. Geophys. Res. Lett. 2010. [Google Scholar] [CrossRef]
  4. Shukla, S.; Lettenmaier, P.D. Seasonal hydrologic prediction in the United States: Understanding the role of initial hydrologic conditions and seasonal climate forecast skill. Hydrol. Earth Syst. Sci. 2011, 15, 3529–3538. [Google Scholar] [CrossRef]
  5. Kumar, S.; Dirmeyer, P.A.; Lawrence, D.M.; DelSole, T.; Altshuler, E.L.; Cash, B.A.; Fennessy, J.M.; Guo, Z.; Kinter, J.L., III; Straus, D.M. Effects of realistic land surface initializations on subseasonal to seasonal soil moisture and temperature predictability in North America and in changing climate simulated by CCSM4. J. Geophys. Res. Atmos. 2014, 119, 13250–13270. [Google Scholar] [CrossRef]
  6. Milly, P.C.D. Stationarity is dead: Whither water management? Science 2008, 319, 573–574. [Google Scholar] [CrossRef] [PubMed]
  7. Famiglietti, J.S.; Cazenave, A.; Eicker, A.; Reager, J.T.; Rodell, M.; Velicogna, I. Satellites Provide the “Big Picture” for Global Hydrology. Science 2015, 349, 684–685. [Google Scholar] [CrossRef] [PubMed]
  8. Tapley, B.D.; Bettadpur, S.; Ries, J.C.; Thompson, P.F.; Watkins, M.M. GRACE measurements of mass variability in the Earth system. Science 2004, 305, 503–505. [Google Scholar] [CrossRef] [PubMed]
  9. Velicogna, I.; Wahr, J. Measurements of time-variable gravity show mass loss in Antarctica. Science 2006, 311, 1754–1756. [Google Scholar] [CrossRef] [PubMed]
  10. Niu, G.-Y.; Seo, K.-W.; Yang, Z.-L.; Wilson, C.; Su, H.; Chen, J.; Rodell, M. Retrieving snow mass from GRACE terrestrial water storage change with a land surface model. Geophys. Res. Lett. 2007. [Google Scholar] [CrossRef]
  11. Rodell, M.; Velicogna, I.; Famiglietti, J.S. Satellite-based estimates of groundwater depletion in India. Nature 2009, 460, 999–1002. [Google Scholar] [CrossRef] [PubMed]
  12. Famiglietti, J.S.; Lo, M.; Ho, S.L.; Bethune, J.; Anderson, K.J.; Syed, T.H.; Swenson, S.C.; de Linage, C.R.; Rodell, M. Satellites measure recent rates of groundwater depletion in California’s Central Valley. Geophys. Res. Lett. 2011. [Google Scholar] [CrossRef]
  13. Richey, A.S.; Thomas, B.F.; Lo, M.-H.; Reager, J.T.; Famiglietti, J.S.; Voss, K.; Swenson, S.; Rodell, M. Quantifying renewable groundwater stress with GRACE. Water Resour. Res. 2015. [Google Scholar] [CrossRef]
  14. Kim, H.; Yeh, P.J.F.; Oki, T.; Kanae, S. Role of rivers in the seasonal variations of terrestrial water storage over global basins. Geophys. Res. Lett. 2009. [Google Scholar] [CrossRef]
  15. Houborg, R.; Rodell, M.; Li, B.; Reichle, R.; Zaitchik, B.F. Drought indicators based on model-assimilated Gravity Recovery and Climate Experiment (GRACE) terrestrial water storage observations. Water Resour. Res. 2012. [Google Scholar] [CrossRef]
  16. Li, B.; Rodell, M.; Zaitchik, B.; Reichle, R.; Koster, R.; van Dam, T. Assimilation of GRACE terrestrial water storage into a land surface model: Evaluation and potential value for drought monitoring in western and central Europe. J. Hydrol. 2012. [Google Scholar] [CrossRef]
  17. Thomas, A.C.; Reager, J.T.; Famiglietti, J.S.; Rodell, M. A GRACE-based water storage deficit approach for hydrological drought characterization. Geophys. Res. Lett. 2014, 41, 1537–1545. [Google Scholar] [CrossRef]
  18. Reager, J.T.; Famiglietti, J.S. Global terrestrial water storage capacity and flood potential using GRACE. Geophys. Res. Lett. 2009. [Google Scholar] [CrossRef]
  19. Reager, J.T.; Thomas, B.F.; Famiglietti, J.S. River basin flood potential inferred using GRACE gravity observations at several months lead time. Nat. Geosci. 2014, 7, 588–592. [Google Scholar] [CrossRef]
  20. Riegger, J.; Tourian, M.J. Characterization of runoff-storage relationships by satellite gravimetry and remote sensing. Water Resour. Res. 2014, 50, 3444–3466. [Google Scholar] [CrossRef]
  21. Sproles, E.A.; Leibowitz, S.G.; Reager, J.T.; Wigington, P.J., Jr.; Famiglietti, J.S.; Patil, S.D. GRACE storage-streamflow hystereses reveal the dynamics of regional watersheds. Hydrol. Earth Syst. Sci. 2015, 11, 12027–12062. [Google Scholar] [CrossRef]
  22. Landerer, F.W.; Swenson, S.C. Accuracy of scaled GRACE terrestrial water storage estimates. Water Resour. Res. 2012. [Google Scholar] [CrossRef]
  23. Famiglietti, J.S.; Rodell, M. Water in the Balance. Science 2013, 340, 1300–1301. [Google Scholar] [CrossRef] [PubMed]
  24. Wood, E.F.; Roundy, J.K.; Troy, T.J.; van Beek, L.P.H.; Bierkens, M.F.P.; Blyth, E.; de Roo, A.; Doll, P.; Ek, M.; Famiglietti, J.; et al. Whitehead. Hyperresolution global land surface modeling: Meeting a grand challenge for monitoring Earth’s terrestrial water. Water Resour. Res. 2011. [Google Scholar] [CrossRef]
  25. Niu, G.Y.; Yang, Z.L.; Dickinson, R.E.; Gulden, L.E. A simple TOPMODEL-based runoff parameterization (SIMTOP) for use in global climate models. J. Geophys. Res. 2005. [Google Scholar] [CrossRef]
  26. Reichle, R.H.; Walker, J.P.; Koster, R.D.; Houser, P.R. Extended versus ensemble Kalman filtering for land data assimilation. J. Hydrometeorol. 2002, 3, 728–740. [Google Scholar] [CrossRef]
  27. Gupta, H.V.; Clark, M.P.; Vrugt, J.A.; Abramowitz, G.; Ye, M. Towards a comprehensive assessment of model structural adequacy. Water Resour. Res. 2012. [Google Scholar] [CrossRef]
  28. Zaitchik, B.F.; Rodell, M.; Reichle, R.H. Assimilation of GRACE terrestrial water storage data into a land surface model: Results for the Mississippi River basin. J. Hydrometeor. 2008, 9, 535–548. [Google Scholar] [CrossRef]
  29. US Army Corps of Engineers. Post 2011 Flood Event Analysis of Missouri River. USACE: Omaha, Nebraska. Available online: http://www.nwd-mr.usace.army.mil/rcc/reports/pdfs/Post2011FloodEventAnalysisofMainstemFloodControlStorage.pdf (accessed on 30 October 2015).
  30. Wahr, J.; Swenson, S.C.; Zlotnicki, V.; Velicogna, I. Time-Variable Gravity from GRACE: First Results. Geophys. Res. Lett. 2004. [Google Scholar] [CrossRef]
  31. Swenson, S.C.; Wahr, J. Methods for inferring regional surface mass anomalies from GRACE measurements of time-variable gravity. J. Geophys. Res. 2002. [Google Scholar] [CrossRef]
  32. GRACE Tellus. Available online: grace.jpl.nasa.gov (accessed on 30 October 2015).
  33. Wahr, J.; Molenaar, M.; Bryan, F. Time-variability of the Earth’s gravity field: Hydrological and oceanic effects and their possible detection using GRACE. J. Geophys. Res. 1998, 103, 205–230. [Google Scholar] [CrossRef]
  34. Rodell, M.; Chen, J.; Kato, H.; Famiglietti, J.; Nigro, J.; Wilson, C. Estimating ground water storage changes in the Mississippi River basin (USA) using GRACE. Hydrogeol. J. 2007, 15, 159–166. [Google Scholar] [CrossRef]
  35. Koster, R.D.; Suarez, M.J.; Ducharne, A.; Stieglitz, M.; Kumar, P. A catchment-based approach to modeling land surface processes in a general circulation model 1. Model structure. J. Geophys. Res. 2000, 105, 809–822. [Google Scholar] [CrossRef]
  36. Beven, K.J.; Kirkby, M.J. A physically-based variable contributing area model of basin hydrology. Hydrol. Sci. J. 1979, 24, 43–69. [Google Scholar]
  37. Gascoin, S.; Ducharne, A.; Ribstein, P.; Carli, M.; Habets, F. Adaptation of a catchment-based land surface model to the hydrogeological setting of the Somme River basin (France). J. Hydrol. 2009, 368, 105–116. [Google Scholar] [CrossRef]
  38. Famiglietti, J.S.; Wood, E.F. Multiscale modeling of spatially variable water and energy balance processes. Water Resour. Res. 1994, 30, 3061–3078. [Google Scholar] [CrossRef]
  39. Ducharne, A.; Koster, R.D.; Suarez, M.J.; Stieglitz, M.; Kumar, P. A catchment-based approach to modeling land surface processes in a general circulation model: 2. Parameter estimation and model de-monstration. J. Geophys. Res. 2000, 105, 24823–24838. [Google Scholar] [CrossRef]
  40. Forman, B.A.; Reichle, R.H.; Rodell, M. Assimilation of terrestrial water storage from GRACE in a snow-dominated basin. Water Resour. Res. 2012. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Reager, J.T.; Thomas, A.C.; Sproles, E.A.; Rodell, M.; Beaudoing, H.K.; Li, B.; Famiglietti, J.S. Assimilation of GRACE Terrestrial Water Storage Observations into a Land Surface Model for the Assessment of Regional Flood Potential. Remote Sens. 2015, 7, 14663-14679. https://doi.org/10.3390/rs71114663

AMA Style

Reager JT, Thomas AC, Sproles EA, Rodell M, Beaudoing HK, Li B, Famiglietti JS. Assimilation of GRACE Terrestrial Water Storage Observations into a Land Surface Model for the Assessment of Regional Flood Potential. Remote Sensing. 2015; 7(11):14663-14679. https://doi.org/10.3390/rs71114663

Chicago/Turabian Style

Reager, John T., Alys C. Thomas, Eric A. Sproles, Matthew Rodell, Hiroko K. Beaudoing, Bailing Li, and James S. Famiglietti. 2015. "Assimilation of GRACE Terrestrial Water Storage Observations into a Land Surface Model for the Assessment of Regional Flood Potential" Remote Sensing 7, no. 11: 14663-14679. https://doi.org/10.3390/rs71114663

Article Metrics

Back to TopTop