A comparison of two physics-based numerical models for simulating surface water–groundwater interactions

Problems in hydrology and water management that involve both surface water and groundwater are best addressed with simulation models that can represent the interactions between these two flow regimes. In the current generation of coupled models, a variety of approaches is used to resolve surface–subsurface interactions and other key processes such as surface flow propagation. In this study we compare two physics-based numerical models that use a 3D Richards equation representation of subsurface flow. In one model, surface flow is represented by a fully 2D kinematic approximation to the Saint–Venant equations with a sheet flow conceptualization. In the second model, surface routing is performed via a quasi-2D diffusive formulation and surface runoff follows a rill flow conceptualization. The coupling between the land surface and the subsurface is handled via an explicit exchange term resolved by continuity principles in the first model (a fully-coupled approach) and by special treatment of atmospheric boundary conditions in the second (a sequential approach). Despite the significant differences in formulation between the two models, we found them to be in good agreement for the simulation experiments conducted. In these numerical tests, on a sloping plane and a tilted V-catchment, we examined saturation excess and infiltration excess runoff production under homogeneous and heterogeneous conditions, the dynamics of the return flow process, the differences in hydrologic response under rill flow and sheet flow parameterizations, and the effects of factors such as grid discretization, time step size, and slope angle. Low sensitivity to vertical discretization and time step size was found for the two models under saturation excess and homogeneous conditions. Larger sensitivity and differences in response were observed under infiltration excess and heterogeneous conditions, due to the different coupling approaches and spatial discretization schemes used in the two models. For these cases, the sensitivity to vertical and temporal resolution was greatest for processes such as reinfiltration and ponding, although the differences between the hydrographs of the two models decreased as mesh and step size were progressively refined. In return flow behavior, the models are in general agreement, with the largest discrepancies, during the recession phase, attributable to the different parameterizations of diffusion in the surface water propagation schemes. Our results also show that under equivalent parameterizations, the rill and sheet flow conceptualizations used in the two models produce very similar responses in terms of hydrograph shape and flow depth distribution.     

Production of oxalates In Vitro by Microbes Isolated from Rock Surfaces with prehistoric paints in the Lower Pecos Region,Texas

Calcium oxalate‐rich rock coatings are ubiquitous on limestone inside dry rock shelters and under bluff overhangs along canyon walls in southwestern Texas. Prehistoric pictographs occur in more than 250 such sites, and the ancient paints are encapsulated within the natural rock coating. Previous studies suggest lichens were the source of the oxalate; however, we report here that microbes cultured and isolated from samples of the coating produce oxalate in vitro. Twenty different bacteria species have been identified in samples from eight different sites, with Bacillus the most common genus, represented by five species. HPLC analyses of inoculated R2B medium after eight months of bacterial growth revealed the presence of oxalate ions in the solid phase of the growth medium. © 2008 Wiley Periodicals, Inc.     

An isobaric–isenthalpic magma mixing model for the Hasan Dagi volcano,Central Anatolia,Turkey

Calc-alkaline rocks of the Hasan Dagi volcano (Central Anatolia, Turkey) are products of arc volcanism triggered by continental collision. Volcanic rocks of the Hasan Dagi range in composition from basalt to rhyolite but are dominated by andesite and dacite. Considering only the mass transfer part of the process leads to an incomplete picture of magma chamber processes. The exclusion of simultaneous calculations of heat and mass transfer between mixing magmas, however, has prevented petrologists from gaining new insights into the magma mixing process. Thus, we report our experimental results in conjunction with modeling with MELTS to test the ideas concerning the petrogenesis of Hasan Dagi volcanic rocks and quantitatively model the relevant petrogenetic processes. Our results demonstrate that the chemical diversity of Hasan Dagi volcano is inconsistent with the closed-system crystallization and differentiation. Thus, (1) our experimental modeling, (2) the agreement between the liquid line of descent defined by the natural rock data and the MELTS calculations, and (3) the agreement between the mineralogy of the rocks and calculated mineralogy corroborate the conclusion that the isobaric–isenthalpic magma mixing of basalt and rhyolite is the major controlling process in the petrogenesis of the Hasan Dagi magmas.     

Quantifying the effects of three-dimensional subsurface heterogeneity on Hortonian runoff processes using a coupled numerical,stochastic approach

The impact of three-dimensional subsurface heterogeneity in the saturated hydraulic conductivity on hillslope runoff generated by excess infiltration (so-called Hortonian runoff) is examined. A fully coupled, parallel subsurface–overland flow model is used to simulate runoff from an idealized hillslope. Ensembles of correlated, Gaussian random fields of saturated hydraulic conductivity are used to create uncertainty in spatial structure. A large number of cases are simulated in a parametric manner with the variance of the hydraulic conductivity varied over orders of magnitude. These cases include rainfall rates above, equal and below the geometric mean of the hydraulic conductivity distribution. These cases are also compared to theoretical representations of runoff production based on simple assumptions regarding (1) the rainfall rate and the value of hydraulic conductivity in the surface cell using a spatially-indiscriminant approach; and (2) a percolation-theory type approach to incorporate so-called runon. Simulations to test the ergodicity of hydraulic conductivity on hillslope runoff are also performed. Results show that three-dimensional stochastic representations of the subsurface hydraulic conductivity can create shallow perching, which has an important effect on runoff behavior that is different than previous two-dimensional analyses. The simple theories are shown to be very poor predictors of the fraction of saturated area that might runoff due to excess infiltration. It is also shown that ergodicity is reached only for a large number of integral scales (∼30) and not achieved for cases where the rainfall rate is less than the geometric mean of the saturated hydraulic conductivity.     

Quantifying the effects of subsurface heterogeneity on hillslope runoff using a stochastic approach

The role of heterogeneity and uncertainty in hydraulic conductivity on hillslope runoff production was evaluated using the fully integrated hydrologic model ParFlow. Simulations were generated using idealized high-resolution hillslopes configured both with a deep water table and a water table equal to the outlet to isolate surface and subsurface flow, respectively. Heterogeneous, correlated random fields were used to create spatial variability in the hydraulic conductivity. Ensembles, generated by multiple realizations of hydraulic conductivity, were used to evaluate how this uncertainty propagates to runoff. Ensemble averages were used to determine the effective runoff for a given hillslope as a function of rainfall rate and degree of subsurface heterogeneity. Cases where the water table is initialized at the outlet show runoff behavior with little sensitivity to variance in hydraulic conductivity. A technique is presented that explicitly interrogates individual realizations at every simulation timestep to partition overland and subsurface flow contributions. This hydrograph separation technique shows that the degree of heterogeneity can play a role in determining proportions of surface and subsurface flow, even when effective hillslope outflow is seen. This method is also used to evaluate current hydrograph separation techniques and demonstrates that recursive filters can accurately proportion overland and base-flow for certain cases.     

An Inter-Laboratory Assessment of the Thorium Isotopic Composition of Synthetic and Rock Reference Materials

We present a concerted international effort to cross-calibrate five synthetic Th isotope reference materials (UCSC Th "A", OU Th "U", WUN, IRMM-35 and IRMM-36), and six rock reference materials (UCSC TML, Icelandic ATHO, USGS BCR-2, USGS W-2, USGS BHVO-2, LV18) using multi-collector inductively coupled plasma-mass spectrometry (MC-ICP-MS). We then compare our new values with a compilation of literature mass spectrometric data for these reference materials and derive recommended "consensus"²³⁰Th/²³²Th values for each. We also present isotope dilution U and Th concentration data for four rock reference materials (UCSC TML, Icelandic ATHO, USGS BCR-2, USGS W-2).     

Patterns of sediment deposition in subsiding coastal salt marshes,Terrebonne Bay,Louisiana: The role of winter storms

High rates of wetland loss in the Mississippi deltaic plain have been attributed to a combination of insufficient marsh sedimentation and relative sea-level rise rates of over 1.2 cm yr^?1. This study examines contemporary patterns of sediment delivery to the marsh surface by evaluating the contribution of individual marsh flooding events. Strong meteorological effects on water level in Terrebonne Bay often mask the usual microtidal fluctuations in water level and cause flood events to be of unpredictable frequency and duration. Sediment deposited on the marsh surface was collected weekly at two sites. Preliminary results allow the relative contributions of tidal and storm inundations to be calculated. Maximum sedimentation is associated with strong southerly winds both causing increased flooding and mobilizing sediment from open bay areas. Sediment deposition is limited by the availability of suspended sediment and the opportunity for its transport onto the marsh surface.