Integrated Surface and Subsurface Hydrological Modeling with Snowmelt and Pore Water Freeze–Thaw

For the simulation of winter hydrological processes a gap in the availability of flow models existed: one either had the choice between (1) physically-based and fully-integrated, but computationally very intensive, or (2) simplified and compartamentalized, but computationally less expensive, simulators. To bridge this gap, we here present the integration of a computationally efficient representation of winter hydrological processes (snowfall, snow accumulation, snowmelt, pore water freeze–thaw) in a fully-integrated surface water-groundwater flow model. This allows the efficient simulation of catchment-scale hydrological processes in locations significantly influenced by winter processes. Snow accumulation and snowmelt are based on the degree-day method and pore water freeze–thaw is calculated with a vertical heat conduction approach. This representation of winter hydrological processes is integrated into the fully-coupled surface water-groundwater flow model HydroGeoSphere. A benchmark for pore water freeze–thaw as well as two illustrative examples are provided.     

Homogenization of spatial patterns of hydrologic response in artificially drained agricultural catchments

Anthropogenic modifications to the landscape, with agricultural activities being a primary driver, have resulted in significant alterations to the hydrologic cycle. Artificial drainage, including surface and subsurface drainage (tile drains), is one of the most extensive manipulations in agricultural landscapes and thus is expected to provide a distinct signature of anthropogenic modification. This study adopts a data synthesis approach in an effort to characterize the signature of artificial subsurface drainage. Daily discharge data from 24 basins across the state of Iowa, which encapsulate a range of anthropogenic modifications, are assessed using a variety of flow metrics. Results indicate that the presence of artificial subsurface drainage leads to a homogenization of landscape hydrologic response. Non‐tiled watersheds exhibit a decrease in the area‐normalized peak discharge and an increase in the baseflow ratio (baseflow/streamflow) with increases in the spatial scale, while scale invariance is apparent in tiled basins. Within‐basin variability in hydrograph recession coefficients also appears to decrease with increases in the proportion of the catchment that is artificially drained. Finally, the differences between tiled and non‐tiled landscapes disappear at scales greater than approximately 2200 km², indicating that this may be a threshold scale for studying the effects of tile drainage. This decrease in within‐basin variability and the scale invariance of hydrologic metrics in artificially drained watersheds are attributed to the creation of a bypass flow hydrologic pathway that bypasses the complexity of the catchment travel paths. Spatial homogeneity in responses implies that it may be possible to develop more parsimonious hydrologic models for these regions. Copyright © 2013 John Wiley & Sons, Ltd.     

Groundwater monitoring at the watershed scale: An evaluation of recharge and nonpoint source pollutant loading in the Clear Creek Watershed,Iowa

Determining the groundwater contribution of nonpoint source pollution at a watershed scale is a challenging issue. In this study, we utilized a top‐down approach to characterize representative groundwater response units (GRUs) based on land use and landscape position (e.g., upland, sideslope, or floodplain) in the 275‐km² Clear Creek Watershed, Iowa. Groundwater monitoring wells were then established along downslope transects in representative GRUs. This unique combination of top‐down/bottom‐up approaches allowed us to estimate groundwater pollutant loads at the watershed scale with minimal monitoring. For the 2015 study period, results indicated that more groundwater recharge occurred in the floodplain (404 mm) compared to the uplands or sideslopes (281 and 165 mm, respectively), irrespective of land use. Recharge in the floodplains consisted of 37% of the annual precipitation, whereas upland wells averaged 26% and sideslopes averaged 15% of the annual precipitation. Less recharge was found to occur beneath perennial grass compared to row crop and urbanized areas. Baseflow discharge accounted for 69% of the total NO₃‐N exported from the Clear Creek Watershed, with row crop areas contributing approximately 95% of the annual load. Orthophosphorus (OP) yields were approximately 0.72 kg/ha beneath urban and suburban areas, three times higher than those in row crop or perennial areas. Urban and suburban areas accounted for 21.4% of groundwater orthophosphorus and chloride loads in the watershed compared to only 8.5% of the land area. Overall, the groundwater load allocation model for baseflow nutrient discharge to Clear Creek can be used to target future nonpoint source load reduction strategies at the watershed scale. The use of GRUs can pinpoint better areas of concern for controlling nutrient loads.     

Rapid,low-cost photogrammetry to monitor volcanic eruptions: an example from Mount St. Helens,Washington, USA

We describe a low-cost application of digital photogrammetry using commercially available photogrammetric software and oblique photographs taken with an off-the-shelf digital camera to create sequential digital elevation models (DEMs) of a lava dome that grew during the 2004–2008 eruption of Mount St. Helens (MSH) volcano. Renewed activity at MSH provided an opportunity to devise and test this method, because it could be validated against other observations of this well-monitored volcano. The datasets consist of oblique aerial photographs (snapshots) taken from a helicopter using a digital single-lens reflex camera. Twelve sets of overlapping digital images of the dome taken during 2004–2007 were used to produce DEMs and to calculate lava dome volumes and extrusion rates. Analyses of the digital images were carried out using photogrammetric software to produce three-dimensional coordinates of points identified in multiple photos. The evolving morphology of the dome was modeled by comparing successive DEMs. Results were validated by comparison to volume measurements derived from traditional vertical photogrammetric surveys by the US Geological Survey Cascades Volcano Observatory. Our technique was significantly less expensive and required less time than traditional vertical photogrammetric techniques; yet, it consistently yielded volume estimates within 5% of the traditional method. This technique provides an inexpensive, rapid assessment tool for tracking lava dome growth or other topographic changes at restless volcanoes.     

Objective rapid delineation of areas at risk from block-and-ash pyroclastic flows and surges

Assessments of pyroclastic flow (PF) hazards are commonly based on mapping of PF and surge deposits and estimations of inundation limits, and/or computer models of varying degrees of sophistication. In volcanic crises a PF hazard map may be sorely needed, but limited time, exposures, or safety aspects may preclude fieldwork, and insufficient time or baseline data may be available for reliable dynamic simulations. We have developed a statistically constrained simulation model for block-and-ash type PFs to estimate potential areas of inundation by adapting methodology from Iverson et al. (Geol Soc America Bull 110:972–984, 1998) for lahars. The predictive equations for block-and-ash PFs are calibrated with data from several volcanoes and given by A = (0.05 to 0.1)V ^2/3, B = (35 to 40)V ^2/3, where A is cross-sectional area of inundation, B is planimetric area and V is deposit volume. The proportionality coefficients were obtained from regression analyses and comparison of simulations to mapped deposits. The method embeds the predictive equations in a GIS program coupled with DEM topography, using the LAHARZ program of Schilling (1998). Although the method is objective and reproducible, any PF hazard zone so computed should be considered as an approximate guide only, due to uncertainties on the coefficients applicable to individual PFs, the authenticity of DEM details, and the volume of future collapses. The statistical uncertainty of the predictive equations, which imply a factor of two or more in predicting A or B for a specified V, is superposed on the uncertainty of forecasting V for the next PF to descend a particular valley. Multiple inundation zones, produced by simulations using a selected range of volumes, partly accommodate these uncertainties. The resulting maps show graphically that PF inundation potentials are highest nearest volcano sources and along valley thalwegs, and diminish with distance from source and lateral distance from thalweg. The model does not explicitly consider dynamic behavior, which can be important. Ash-cloud surge impact limits must be extended beyond PF hazard zones and we provide several approaches to do this. The method has been used to supply PF and surge hazard maps in two crises: Merapi 2006; and Montserrat 2006–2007.     

Field Observation of Diurnal Dissolved Oxygen Fluctuations in Shallow Groundwater

Dissolved oxygen (DO) concentrations influence many biogeochemical processes in groundwater systems but studies of temporal variability in DO are lacking. In this study, we used an optical DO probe to measure rapid changes in concentration due to plant‐groundwater interaction at an alluvial aquifer field site in Iowa. Diurnal DO concentrations were observed during mid‐ to late‐summer when soil conditions were dry, fluctuating approximately 0.2 to 0.3 mg/L on a daily basis. DO fluctuations in groundwater were out‐of‐phase with diurnal water table fluctuations, increasing during the day and decreasing at night. DO consumption at night is likely due to increased soil autotrophic and heterotrophic respiration linked with patterns of carbon supply derived from daytime photosynthetic activity, and consistent with available literature on diurnal soil respiration patterns. Although more work is needed to quantify specific processes, our results indicate the potential usefulness of the new optical DO technology to reveal insights regarding many ecohydrological processes.     

An experimental investigation of Iceland and Reykjanes Ridge tholeiites: I. Phase relations

Tholeiite basalts from 60° N to 65° N on the Mid-Atlantic Ridge were melted and recrystallized at atmospheric pressure in a CO₂-H₂ gas mixture. Seven basalts are from the Langjokull-Thingvellir volcanic zone and the Reykjanes Peninsula of Iceland and nine are from the Reykjanes Ridge. The crystallization sequence in both Iceland and Reykjanes Ridge basalts with (Total Fe as FeO)/(Total Fe as FeO+ MgO) [F/F + M] less than 0.6 is olivine, plagioclase, clinopyroxene. Chromian spinel crystallizes before plagioclase in one Iceland and one Reykjanes Ridge basalt with F/F+M less than 0.57. Chemical differences of the two groups of basalts (lower SiO₂ and higher alkalis in Iceland basalts) can not simply be a result of low pressure fractional crystallization. Liquidus temperatures of the seven Iceland basalts decreases from 1,230° C to 1,170° C as the F/F+M of the rock increases from 0.52 to 0.70. The liquidus temperatures of the Reykjanes Ridge basalts are about 10° C lower than those of the Iceland basalts for the same F/F+M value. The profile of measured liquidus temperatures from 65° N on Iceland to 60° N on the Reykjanes Ridge has a minimum value at 63.2° N on the Reykjanes Ridge just south of Iceland. Model calculations of the pressure of phenocryst crystallization indicate that olivine and plagioclase in Langjokull basalts could have equilibrated between 2.0 and 6.2 kb (200 to 620 MPa). Phenocryst assemblages in Reykjanes Ridge basalts at 60° N could have crystallized together at greater than 2 kb (200 MPa) and probably less than 8 kb (800 MPa). A minimum in the equilibrium pressure of phenocryst crystallization occurs between 62.9° and 64° N and coincides with the minimum in the experimentally determined liquidus temperatures. The more extensive fractionation at low pressure in this area could be related to the shift of the Mid-Atlantic Ridge axis along the leaky transform fault from the Reykjanes Ridge to the Thingvellir volcanic zone.     

The Hydraulic Effects of Lost Circulation and the Implications for Contaminant Migration During Drilling

The impact of lost circulation during rotary drilling near an existing monitoring well cluster was evaluated by periodic measurements of water levels and contaminant concentrations at the well cluster. Due to regulatory concerns, changes in water levels or VOC concentration in the well cluster during drilling would trigger monitoring well redevelopment. The borehole was drilled approximately 30 feet northeast of four nested monitoring wells that screen Devonian and Silurian carbonate bedrock at depths of 15, 60, 130, and 190 feet. Following complete circulation loss at depths of 177 and 1 S3 feet in the borehole, a rapid decrease in water levels was observed in the upper three monitoring wells. The water level in the well that was screened through the lost circulation zones increased slightly.
Decreasing water levels in formations located above the point of circulation loss appear to occur in response to a sudden decrease in borehole fluid pressure caused by the flow of drilling fluid into the formation. The relative contribution of contaminated formation water lo the borehole can be estimated by using the time-drawdown relationship and estimates of transmissivity. At the point of circulation loss, significant dilution of contaminant concentrations occurs from the loss of drilling fluid into the contaminated zone. Contaminated formation water entering the borehole during periods of complete lost circulation may mobilize contaminants from upper lo lower formations. Lost circulation into a formation would be signaled by a water level increase in monitoring wells. The wells would subsequently require development to remove the volume of fluid lost to the formation, including both drilling fluid and contaminated formation water. Monitoring wells exhibiting declining water levels following lost circulation would not require development since drilling water has not entered the zones screened by these wells.