Search USGSSearch

SEARCH · Search USGS

Results for “Hydrology and Earth System Sciences”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 55 records · Page 3Linked to original sources

Changing suspended sediment in United States rivers and streams: Linking sediment trends to changes in land use/cover, hydrology and climate

Sediment is one of the leading pollutants in rivers and streams across the United States (US) and the world. Between 1992 and 2012, concentrations of annual mean suspended sediment decreased at over half of the 137 stream sites assessed across the contiguous US. Increases occurred at less than 25 % of the sites, and the direction of change was uncertain at the remaining 25 %. Sediment trends were characterized using the Weighted Regressions on Time, Discharge, and Season (WRTDS) model, and decreases in sediment ranged from −95 % to −8.5 % of the 1992 concentration. To explore potential drivers of these changes, the sediment trends were (1) parsed into two broad contributors of change, changes in land management versus changes in the streamflow regime, and (2) grouped by land use of the watershed and correlated to concurrent changes in land use or land cover (land use/cover), hydrology and climate variables and static/long-term watershed characteristics. At 83 % of the sites, changes in land management (captured by changes in the concentration–streamflow relationship over time; C–Q relationship) contributed more to the change in the sediment trend than changes in the streamflow regime alone (i.e., any systematic change in the magnitude, frequency or timing of flows). However, at >50 % of the sites, changes in the streamflow regime contributed at least a 5 % change in sediment, and at 11 sites changes in the streamflow regime contributed over half the change in sediment, indicating that at many sites changes in streamflow were not the main driver of changes in sediment but were often an important supporting factor. Correlations between sediment trends and concurrent changes in land use/cover, hydrology and climate were often stronger at sites draining watersheds with more homogenous, human-related land uses (i.e., agricultural and urban lands) compared to mixed-use or undeveloped lands. At many sites, decreases in sediment occurred despite small-to-moderate increases in the amount of urban or agricultural land in the watershed, suggesting conservation efforts and best-management practices (BMPs) used to reduce sediment runoff to streams may be successful, up to a point, as lands are converted to urban and agricultural uses.

Hydrology and Earth System Sciences

Implications of model selection: A comparison of publicly available, conterminous US-extent hydrologic component estimates

Spatiotemporally continuous estimates of the hydrologic cycle are often generated through hydrologic modeling, reanalysis, or remote sensing (RS) methods and are commonly applied as a supplement to, or a substitute for, in situ measurements when observational data are sparse or unavailable. This study compares estimates of precipitation ( P ), actual evapotranspiration (ET), runoff ( R ), snow water equivalent (SWE), and soil moisture (SM) from 87 unique data sets generated by 47 hydrologic models, reanalysis data sets, and remote sensing products across the conterminous United States (CONUS). Uncertainty between hydrologic component estimates was shown to be high in the western CONUS, with median uncertainty (measured as the coefficient of variation) ranging from 11 % to 21 % for P , 14 % to 26 % for ET, 28 % to 82 % for R , 76 % to 84 % for SWE, and 36 % to 96 % for SM. Uncertainty between estimates was lower in the eastern CONUS, with medians ranging from 5 % to 14 % for P, 13 % to 22 % for ET, 28 % to 82 % for R , 53 % to 63 % for SWE, and 42 % to 83 % for SM. Interannual trends in estimates from 1982 to 2010 show common disagreement in R, SWE, and SM. Correlating fluxes and stores against remote-sensing-derived products show poor overall correlation in the western CONUS for ET and SM estimates. Study results show that disagreement between estimates can be substantial, sometimes exceeding the magnitude of the measurements themselves. The authors conclude that multimodel ensembles are not only useful but are in fact a necessity for accurately representing uncertainty in research results. Spatial biases of model disagreement values in the western United States show that targeted research efforts in arid and semiarid water-limited regions are warranted, with the greatest emphasis on storage and runoff components, to better describe complexities of the terrestrial hydrologic system and reconcile model disagreement.

Hydrology and Earth System Sciences

Determination of vadose zone and saturated zone nitrate lag times using long-term groundwater monitoring data and statistical machine learning

In this study, we explored the use of statistical machine learning and long-term groundwater nitrate monitoring data to estimate vadose zone and saturated zone lag times in an irrigated alluvial agricultural setting. Unlike most previous statistical machine learning studies that sought to predict groundwater nitrate concentrations within aquifers, the focus of this study was to leverage available groundwater nitrate concentrations and other environmental variables to determine mean regional vertical velocities (transport rates) of water and solutes in the vadose zone and saturated zone (3.50 and 3.75 m yr −1 , respectively). The statistical machine learning results are consistent with two primary recharge processes in this western Nebraska aquifer, namely ( 1 ) diffuse recharge from irrigation and precipitation across the landscape and ( 2 ) focused recharge from leaking irrigation conveyance canals. The vadose zone mean velocity yielded a mean recharge rate (0.46 m yr −1 ) consistent with previous estimates from groundwater age dating in shallow wells (0.38 m yr −1 ). The saturated zone mean velocity yielded a recharge rate (1.31 m yr −1 ) that was more consistent with focused recharge from leaky irrigation canals, as indicated by previous results of groundwater age dating in intermediate-depth wells (1.22 m yr −1 ). Collectively, the statistical machine learning model results are consistent with previous observations of relatively high water fluxes and short transit times for water and nitrate in the primarily oxic aquifer. Partial dependence plots from the model indicate a sharp threshold in which high groundwater nitrate concentrations are mostly associated with total travel times of 7 years or less, possibly reflecting some combination of recent management practices and a tendency for nitrate concentrations to be higher in diffuse infiltration recharge than in canal leakage water. Limitations to the machine learning approach include the non-uniqueness of different transport rate combinations when comparing model performance and highlight the need to corroborate statistical model results with a robust conceptual model and complementary information such as groundwater age.

Nebraska

Bedrock depth influences spatial patterns of summer baseflow, temperature and flow disconnection for mountainous headwater streams

In mountain headwater streams, the quality and resilience of summer cold-water habitat is generally regulated by stream discharge, longitudinal stream channel connectivity and groundwater exchange. These critical hydrologic processes are thought to be influenced by the stream corridor bedrock contact depth (sediment thickness), a parameter often inferred from sparse hillslope borehole information, piezometer refusal and remotely sensed data. To investigate how local bedrock depth might control summer stream temperature and channel disconnection (dewatering) patterns, we measured stream corridor bedrock depth by collecting and interpreting 191 passive seismic datasets along eight headwater streams in Shenandoah National Park (Virginia, USA). In addition, we used multi-year stream temperature and streamflow records to calculate several baseflow-related metrics along and among the study streams. Finally, comprehensive visual surveys of stream channel dewatering were conducted in 2016, 2019 and 2021 during summer low flow conditions (124 total km of stream length). We found that measured bedrock depths along the study streams were not well-characterized by soils maps or an existing global-scale geologic dataset where the latter overpredicted measured depths by 12.2 m (mean) or approximately four times the average bedrock depth of 2.9 m. Half of the eight study stream corridors had an average bedrock depth of less than 2 m. Of the eight study streams, Staunton River had the deepest average bedrock depth (3.4 m), the coldest summer temperature profiles and substantially higher summer baseflow indices compared to the other study steams. Staunton River also exhibited paired air and water annual temperature signals suggesting deeper groundwater influence, and the stream channel did not dewater in lower sections during any baseflow survey. In contrast, Paine Run and Piney River did show pronounced, patchy channel dewatering, with Paine Run having dozens of discrete dry channel sections ranging from 1 to greater than 300 m in length. Stream dewatering patterns were apparently influenced by a combination of discrete deep bedrock (20+ m) features and more subtle sediment thickness variation (1–4 m) depending on local stream valley hydrogeology. In combination, these unique datasets show the first large-scale empirical support for existing conceptual models of headwater stream disconnection based on spatially variable underflow capacity and shallow groundwater supply.

Virginia

Benchmarking high-resolution hydrologic model performance of long-term retrospective streamflow simulations in the contiguous United States

Because use of high-resolution hydrologic models is becoming more widespread and estimates are made over large domains, there is a pressing need for systematic evaluation of their performance. Most evaluation efforts to date have focused on smaller basins that have been relatively undisturbed by human activity, but there is also a need to benchmark model performance more comprehensively, including basins impacted by human activities. This study benchmarks the long-term performance of two process-oriented, high-resolution, continental-scale hydrologic models that have been developed to assess water availability and risks in the United States (US): the National Water Model v2.1 application of WRF-Hydro (NWMv2.1) and the National Hydrologic Model v1.0 application of the Precipitation–Runoff Modeling System (NHMv1.0). The evaluation is performed on 5390 streamflow gages from 1983 to 2016 ( ∼ 33 years) at a daily time step, including both natural and human-impacted catchments, representing one of the most comprehensive evaluations over the contiguous US. Using the Kling–Gupta efficiency as the main evaluation metric, the models are compared against a climatological benchmark that accounts for seasonality. Overall, the model applications show similar performance, with better performance in minimally disturbed basins than in those impacted by human activities. Relative regional differences are also similar: the best performance is found in the Northeast, followed by the Southeast, and generally worse performance is found in the Central and West areas. For both models, about 80 % of the sites exceed the seasonal climatological benchmark. Basins that do not exceed the climatological benchmark are further scrutinized to provide model diagnostics for each application. Using the underperforming subset, both models tend to overestimate streamflow volumes in the West, which could be attributed to not accounting for human activities, such as active management. Both models underestimate flow variability, especially the highest flows; this was more pronounced for NHMv1.0. Low flows tended to be overestimated by NWMv2.1, whereas there were both over and underestimations for NHMv1.0, but they were less severe. Although this study focused on model diagnostics for underperforming sites based on the seasonal climatological benchmark, metrics for all sites for both model applications are openly available online.

Hydrology and Earth System Sciences

Increasing seasonal variation in the extent of rivers and lakes from 1984 to 2022

Knowledge of the spatial and temporal distribution of surface water is important for water resource management, flood risk assessment, monitoring ecosystem health, constraining estimates of biogeochemical cycles and understanding our climate. While global-scale spatiotemporal change detection of surface water has significantly improved in recent years due to planetary-scale remote sensing and computing, it has remained challenging to distinguish the changing characteristics of rivers and lakes. Here we analyze the spatial extent of permanent and seasonal rivers and lakes globally over the past 38 years based on new data of river system extents and surface water trends. Results show that while the total permanent surface area of both rivers and lakes has remained relatively constant, the areas with intermittent seasonal coverage have increased by 12 % and 27 % for rivers and lakes, respectively. The increase is statistically significant in over 84 % of global water catchments based on Spearman's rank correlations (rho) above 0.05 and p values less than 0.05. The seasonal river extent is nearly 32 % larger than the previously observed annual mean river extent, suggesting large seasonal variations that impact not only ecosystem health but also estimations of terrestrial biogeochemical cycles of carbon. The outcomes of our analysis are shared as the Surface Area of Rivers and Lakes (SARL) database, serving as a valuable resource for monitoring and research of hydrological cycles, ecosystem accounting, and water management.

Hydrology and Earth System Sciences

Isotopic evaluation of the National Water Model reveals missing agricultural irrigation contributions to streamflow across the western United States

The National Water Model (NWM) provides critical analyses and projections of streamflow that support water management decisions. However, the NWM performs poorly in lower-elevation rivers of the western United States (US). The accuracy of the NWM depends on the fidelity of the model inputs and the representation and calibration of model processes and water sources. To evaluate the NWM performance in the western US, we compared observations of river water isotope ratios ( 18 O / "> 16 O and 2 H / "> 1 H expressed in δ notation) to NWM-flux-estimated (model) river reach isotope ratios. The modeled estimates were calculated from long-term (2000–2019) mean summer (June, July, and August) NWM hydrologic fluxes and gridded isotope ratios using a mass balance approach. The observational dataset comprised 4503 in-stream water isotope observations in 877 reaches across 5 basins. A simple regression between observed and modeled isotope ratios explained 57.9 % ( δ 18 O) and 67.1 % ( δ 2 H) of variance, although observations were 0.5 ‰ ( δ 18 O) and 4.8 ‰ ( δ 2 H) higher, on average, than mass balance estimates. The unexplained variance suggest that the NWM does not include all relevant water fluxes to rivers. To infer possible missing water fluxes, we evaluated patterns in observation–model differences using δ 18 O diff ( δ 18 O obs − δ 18 O mod ) and d diff ( δ 2 H diff - 8 ⋅ δ 18 O diff "> ). We detected evidence of evaporation in observations but not model estimates (negative d diff and positive δ 18 O diff ) at lower-elevation, higher-stream-order, arid sites. The catchment actual-evaporation-to-precipitation ratio, the fraction of streamflow estimated to be derived from agricultural irrigation, and whether a site was reservoir-affected were all significant predictors of d diff in a linear mixed-effects model, with up to 15.2 % of variance explained by fixed effects. This finding is supported by seasonal patterns, groundwater levels, and isotope ratios, and it suggests the importance of including irrigation return flows to rivers, especially in lower-elevation, higher-stream-order, arid rivers of the western US.

Hydrology and Earth Systems Science

Pluvial and potential compound flooding in a coupled coastal modeling framework: New York City during post-tropical Cyclone Ida (2021)

Many coastal urban areas are prone to extreme pluvial flooding due to limitations in stormwater system capacity, with the additional potential for flooding compounded by storm surge, tides, and waves. Understanding and simulating these processes can improve prediction and flood risk management. Here, we adapt the Coupled Ocean–Atmosphere–Wave–Sediment Transport modeling framework (COAWST) to simulate pluvial flooding from post-tropical Cyclone Ida (2021) in the Jamaica Bay watershed of New York City (NYC). We modify the model to capture the volumetric effects of rainfall and parameterize soil infiltration and a stormwater conveyance system as the drainage rate. We generate a spatially continuous flood map of Ida with a root-mean-square error (RMSE) of 20 cm when compared to high-water marks, useful for understanding Ida's impacts and subsequent mitigation planning. Results show that over 23 km 2 and 4621 buildings were flooded deeper than 0.3 m during Ida. Sensitivity analyses are used to study the broader risk from events like Ida (pluvial flooding) as well as potential compound (pluvial–coastal) flooding. Spatial shifting of the storm track within a typical 12 h forecast uncertainty reveals a worst-case scenario that increases this flooded area to 62 km 2 (5907 buildings). Shifting Ida's rainfall to coincide with high tide increases this flooded area by 1 km 2 , a relatively small change due to the lack of significant storm surge. The application of COAWST to this storm event addresses a broader goal of developing the capability to model compound pluvial–coastal flooding by simultaneously representing coastal storm processes such as rain, tide, waves, erosion, and atmosphere–wave–ocean interactions. The sensitivity analysis results underscore the need for detailed flood risk assessments, showing that Ida, already NYC's worst rain event, could have been even more devastating with slight shifts in the storm track.

New York

Interrogating process deficiencies in large-scale hydrologic models with interpretable machine learning

Large-scale hydrologic models are increasingly being developed for operational use in the forecasting and planning of water resources. However, the predictive strength of such models depends on how well they resolve various functions of catchment hydrology, which are influenced by gradients in climate, topography, soils, and land use. Most assessments of hydrologic model uncertainty have been limited to traditional statistical methods. Here, we present a proof-of-concept approach that uses interpretable machine learning techniques to provide post hoc assessment of model sensitivity and process deficiency in hydrologic models. We train a random forest model to predict the Kling–Gupta efficiency (KGE) of National Water Model (NWM) and National Hydrologic Model (NHM) streamflow predictions for 4383 stream gauges in the conterminous United States. Thereafter, we explain the local and global controls that 48 catchment attributes exert on KGE prediction using interpretable Shapley values. Overall, we find that soil water content is the most impactful feature controlling successful model performance, suggesting that soil water storage is difficult for hydrologic models to resolve, particularly for arid locations. We identify nonlinear thresholds beyond which predictive performance decreases for NWM and NHM. For example, soil water content less than 210 mm, precipitation less than 900 mm yr −1 , road density greater than 5 km km −2 , and lake area percent greater than 10 % contributed to lower KGE values. These results suggest that improvements in how these influential processes are represented could result in the largest increases in NWM and NHM predictive performance. This study demonstrates the utility of interrogating process-based models using data-driven techniques, which has broad applicability and potential for improving the next generation of large-scale hydrologic models.

conterminous United States

Technical note: A low-cost approach to monitoring relative streamflow dynamics in small headwater streams using time lapse imagery and a deep learning model

Despite their ubiquity and importance as freshwater habitat, small headwater streams are under-monitored by existing stream gage networks. To address this gap, we describe a low-cost, non-contact, and low-effort method that enables organizations to monitor relative streamflow dynamics in small headwater streams. The method uses a camera to capture repeat images of the stream from a fixed position. A person then annotates pairs of images, in each case indicating which image has more apparent streamflow or indicating equal flow if no difference is discernible. A deep learning modeling framework called streamflow rank estimation (SRE) is then trained on the annotated image pairs and applied to rank all images from highest to lowest apparent streamflow. From this result a relative hydrograph can be derived. We found that our modeled relative hydrograph dynamics matched the observed hydrograph dynamics well for 11 cameras at 8 streamflow sites in western Massachusetts. Higher performance was observed during the annotation period (median Kendall's Tau rank correlation of 0.75, with a range of 0.6–0.83) than after it (median Kendall's Tau of 0.59, with range 0.34–0.74). We found that annotation performance was generally consistent across the 11 camera sites and 2 individual annotators and was positively correlated with streamflow variability at a site. A scaling simulation determined that model performance improvements were limited after 1000 annotation pairs. Our model's estimates of relative flow, while not equivalent to absolute flow, may still be useful for many applications, such as ecological modeling and calculating event-based hydrological statistics (e.g., the number of out-of-bank floods). We anticipate that this method will be a valuable tool to extend existing stream monitoring networks and provide new insights on dynamic headwater systems.

Massachusetts

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States

Antecedent flow conditions and nitrate concentrations in the Mississippi River basin

The relationship between antecedent flow conditions and nitrate concentrations was explored at eight sites in the 2.9 million square kilometers (km 2 ) Mississippi River basin, USA. Antecedent flow conditions were quantified as the ratio between the mean daily flow of the previous year and the mean daily flow from the period of record (Qratio), and the Qratio was statistically related to nitrate anomalies (the unexplained variability in nitrate concentration after filtering out season, long-term trend, and contemporaneous flow effects) at each site. Nitrate anomaly and Qratio were negatively related at three of the four major tributary sites and upstream in the Mississippi River, indicating that when mean daily streamflow during the previous year was lower than average, nitrate concentrations were higher than expected. The strength of these relationships increased when data were subdivided by contemporaneous flow conditions. Five of the eight sites had significant negative relationships ( p ≤ 0.05) at high or moderately high contemporaneous flows, suggesting nitrate that accumulates in these basins during a drought is flushed during subsequent high flows. At half of the sites, when mean daily flow during the previous year was 50 percent lower than average, nitrate concentration can be from 9 to 27 percent higher than nitrate concentrations that follow a year with average mean daily flow. Conversely, nitrate concentration can be from 8 to 21 percent lower than expected when flow during the previous year was 50 percent higher than average. Previously documented for small, relatively homogenous basins, our results suggest that relationships between antecedent flows and nitrate concentrations are also observable at a regional scale. Relationships were not observed (using all contemporaneous flow data together) for basins larger than 1 million km 2 , suggesting that above this limit the overall size and diversity within these basins may necessitate the use of more complicated statistical approaches or that there may be no discernible basin-wide relationship with antecedent flow. The relationships between nitrate concentration and Qratio identified in this study serve as the basis for future studies that can better define specific hydrologic processes occurring during and after a drought (or high flow period) which influence nitrate concentration, such as the duration or magnitude of low flows, and the timing of low and high flows.

Mississippi River Basin

Assessing winter cover crop nutrient uptake efficiency using a water quality simulation model

Winter cover crops are an effective conservation management practice with potential to improve water quality. Throughout the Chesapeake Bay Watershed (CBW), which is located in the Mid-Atlantic US, winter cover crop use has been emphasized and federal and state cost-share programs are available to farmers to subsidize the cost of winter cover crop establishment. The objective of this study was to assess the long-term effect of planting winter cover crops at the watershed scale and to identify critical source areas of high nitrate export. A physically-based watershed simulation model, Soil and Water Assessment Tool (SWAT), was calibrated and validated using water quality monitoring data and satellite-based estimates of winter cover crop species performance to simulate hydrological processes and nutrient cycling over the period of 1991–2000. Multiple scenarios were developed to obtain baseline information on nitrate loading without winter cover crops planted and to investigate how nitrate loading could change with different winter cover crop planting scenarios, including different species, planting times, and implementation areas. The results indicate that winter cover crops had a negligible impact on water budget, but significantly reduced nitrate leaching to groundwater and delivery to the waterways. Without winter cover crops, annual nitrate loading was approximately 14 kg ha −1 , but it decreased to 4.6–10.1 kg ha −1 with winter cover crops resulting in a reduction rate of 27–67% at the watershed scale. Rye was most effective, with a potential to reduce nitrate leaching by up to 93% with early planting at the field scale. Early planting of winter cover crops (~30 days of additional growing days) was crucial, as it lowered nitrate export by an additional ~2 kg ha −1 when compared to late planting scenarios. The effectiveness of cover cropping increased with increasing extent of winter cover crop implementation. Agricultural fields with well-drained soils and those that were more frequently used to grow corn had a higher potential for nitrate leaching and export to the waterways. This study supports the effective implement of winter cover crop programs, in part by helping to target critical pollution source areas for winter cover crop implementation.

Chesapeake Bay Watershed

Predicting East African spring droughts using Pacific and Indian Ocean sea surface temperature indices

In southern Ethiopia, Eastern Kenya, and southern Somalia poor boreal spring rains in 1999, 2000, 2004, 2007, 2008, 2009 and 2011 contributed to severe food insecurity and high levels of malnutrition. Predicting rainfall deficits in this region on seasonal and decadal time frames can help decision makers support disaster risk reduction while guiding climate-smart adaptation and agricultural development. Building on recent research that links more frequent droughts to a stronger Walker Circulation, warming in the Indo-Pacific warm pool, and an increased western Pacific sea surface temperature (SST) gradient, we explore the dominant modes of East African rainfall variability, links between these modes and sea surface temperatures, and a simple index-based monitoring-prediction system suitable for drought early warning.

Hydrology and Earth System Sciences

Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model

Soil moisture is a crucial control on surface water and energy fluxes, vegetation, and soil carbon cycling. Earth-system models (ESMs) generally represent an areal-average soil-moisture state in gridcells at scales of 50–200 km and as a result are not able to capture the nonlinear effects of topographically-controlled subgrid heterogeneity in soil moisture, in particular where wetlands are present. We addressed this deficiency by building a subgrid representation of hillslope-scale topographic gradients, TiHy (Tiled-hillslope Hydrology), into the Geophysical Fluid Dynamics Laboratory (GFDL) land model (LM3). LM3-TiHy models one or more representative hillslope geometries for each gridcell by discretizing them into land model tiles hydrologically coupled along an upland-to-lowland gradient. Each tile has its own surface fluxes, vegetation, and vertically-resolved state variables for soil physics and biogeochemistry. LM3-TiHy simulates a gradient in soil moisture and water-table depth between uplands and lowlands in each gridcell. Three hillslope hydrological regimes appear in non-permafrost regions in the model: wet and poorly-drained, wet and well-drained, and dry; with large, small, and zero wetland area predicted, respectively. Compared to the untiled LM3 in stand-alone experiments, LM3-TiHy simulates similar surface energy and water fluxes in the gridcell-mean. However, in marginally wet regions around the globe, LM3-TiHy simulates shallow groundwater in lowlands, leading to higher evapotranspiration, lower surface temperature, and higher leaf area compared to uplands in the same gridcells. Moreover, more than four-fold larger soil carbon concentrations are simulated globally in lowlands as compared with uplands. We compared water-table depths to those simulated by a recent global model-observational synthesis, and we compared wetland and inundated areas diagnosed from the model to observational datasets. The comparisons demonstrate that LM3-TiHy has the capability to represent some of the controls of these hydrological variables, but also that improvement in parameterization and input datasets are needed for more realistic simulations. We found large sensitivity in model-diagnosed wetland and inundated area to the depth of conductive soil and the parameterization of macroporosity. With improved parameterization and inclusion of peatland biogeochemical processes, the model could provide a new approach to investigating the vulnerability of Boreal peatland carbon to climate change in ESMs.

Hydrology and Earth System Sciences

Using multi-source satellite data for lake level modelling in ungauged basins: A case study for Lake Turkana, East Africa

Managing limited surface water resources is a great challenge in areas where ground-based data are either limited or unavailable. Direct or indirect measurements of surface water resources through remote sensing offer several advantages of monitoring in ungauged basins. A physical based hydrologic technique to monitor lake water levels in ungauged basins using multi-source satellite data such as satellite-based rainfall estimates, modelled runoff, evapotranspiration, a digital elevation model, and other data is presented. This approach is applied to model Lake Turkana water levels from 1998 to 2009. Modelling results showed that the model can reasonably capture all the patterns and seasonal variations of the lake water level fluctuations. A composite lake level product of TOPEX/Poseidon, Jason-1, and ENVISAT satellite altimetry data is used for model calibration (1998-2000) and model validation (2001-2009). Validation results showed that model-based lake levels are in good agreement with observed satellite altimetry data. Compared to satellite altimetry data, the Pearson's correlation coefficient was found to be 0.81 during the validation period. The model efficiency estimated using NSCE is found to be 0.93, 0.55 and 0.66 for calibration, validation and combined periods, respectively. Further, the model-based estimates showed a root mean square error of 0.62 m and mean absolute error of 0.46 m with a positive mean bias error of 0.36 m for the validation period (2001-2009). These error estimates were found to be less than 15 % of the natural variability of the lake, thus giving high confidence on the modelled lake level estimates. The approach presented in this paper can be used to (a) simulate patterns of lake water level variations in data scarce regions, (b) operationally monitor lake water levels in ungauged basins, (c) derive historical lake level information using satellite rainfall and evapotranspiration data, and (d) augment the information provided by the satellite altimetry systems on changes in lake water levels. ?? Author(s) 2011.

Hydrology and Earth System Sciences Discussions

Assessing the potential hydrological impact of the Gibe III Dam on Lake Turkana water level using multi-source satellite data

Lake Turkana, the largest desert lake in the world, is fed by ungauged or poorly gauged river systems. To meet the demand of electricity in the East African region, Ethiopia is currently building the Gibe III hydroelectric dam on the Omo River, which supplies more than 80% of the inflows to Lake Turkana. On completion, the Gibe III dam will be the tallest dam in Africa with a height of 241 m. However, the nature of interactions and potential impacts of regulated inflows to Lake Turkana are not well understood due to its remote location and unavailability of reliable in-situ datasets. In this study, we used 12 years (1998–2009) of existing multi-source satellite and model-assimilated global weather data. We use calibrated multi-source satellite data-driven water balance model for Lake Turkana that takes into account model routed runoff, lake/reservoir evapotranspiration, direct rain on lakes/reservoirs and releases from the dam to compute lake water levels. The model evaluates the impact of Gibe III dam using three different approaches such as (a historical approach, a knowledge-based approach, and a nonparametric bootstrap resampling approach) to generate rainfall-runoff scenarios. All the approaches provided comparable and consistent results. Model results indicated that the hydrological impact of the dam on Lake Turkana would vary with the magnitude and distribution of rainfall post-dam commencement. On average, the reservoir would take up to 8–10 months, after commencement, to reach a minimum operation level of 201 m depth of water. During the dam filling period, the lake level would drop up to 2 m (95% confidence) compared to the lake level modelled without the dam. The lake level variability caused by regulated inflows after the dam commissioning were found to be within the natural variability of the lake of 4.8 m. Moreover, modelling results indicated that the hydrological impact of the Gibe III dam would depend on the initial lake level at the time of dam commencement. Areas along the Lake Turkana shoreline that are vulnerable to fluctuations in lake levels were also identified. This study demonstrates the effectiveness of using existing multi-source satellite data in a basic modeling framework to assess the potential hydrological impact of an upstream dam on a terminal downstream lake. The results obtained from this study could also be used to evaluate alternate dam-filling scenarios and assess the potential impact of the dam on Lake Turkana under different operational strategies.

Hydrology and Earth System Sciences Discussions

Linking soil moisture balance and source-responsive models to estimate diffuse and preferential components of groundwater recharge

Results are presented of a detailed study into the vadose zone and shallow water table hydrodynamics of a field site in Shropshire, UK. A conceptual model is developed and tested using a range of numerical models, including a modified soil moisture balance model (SMBM) for estimating groundwater recharge in the presence of both diffuse and preferential flow components. Tensiometry reveals that the loamy sand topsoil wets up via macropore flow and subsequent redistribution of moisture into the soil matrix. Recharge does not occur until near-positive pressures are achieved at the top of the sandy glaciofluvial outwash material that underlies the topsoil, about 1 m above the water table. Once this occurs, very rapid water table rises follow. This threshold behaviour is attributed to the vertical discontinuity in the macropore system due to seasonal ploughing of the topsoil, and a lower permeability plough/iron pan restricting matrix flow between the topsoil and the lower outwash deposits. Although the wetting process in the topsoil is complex, a SMBM is shown to be effective in predicting the initiation of preferential flow from the base of the topsoil into the lower outwash horizon. The rapidity of the response at the water table and a water table rise during the summer period while flow gradients in the unsaturated profile were upward suggest that preferential flow is also occurring within the outwash deposits below the topsoil. A variation of the source-responsive model proposed by Nimmo (2010) is shown to reproduce the observed water table dynamics well in the lower outwash horizon when linked to a SMBM that quantifies the potential recharge from the topsoil. The results reveal new insights into preferential flow processes in cultivated soils and provide a useful and practical approach to accounting for preferential flow in studies of groundwater recharge estimation.

Shropshire County