Sediment movement under unidirectional flows: an assessment of empirical threshold curves

A selection of threshold data sets from previous investigations is used to re-examine, selectively, some of the empirical curves that are most commonly used for the prediction of sediment threshold. Simple analytical formulae are derived describing single line curves for mean threshold values and envelopes for the initial movement of discrete particles and the commencement of mass sediment transport phases of the critical condition.The complexity of factors governing the prevailing conditions at threshold, together with the randomness of events, limits deterministic solutions; however, the curves provide an approach that can be used readily in practical applications. In contrast, the envelope presentation is more of a stochastic approach to the definition of threshold. It recognises that there is no definitive threshold condition under which a particular particle size can be displaced, but rather that a range of threshold values exist.The empirical threshold relationships were derived on the basis of a relatively large amount of data. However, the sources utilised have interpreted the ‘critical condition’ in various ways, implying that the formulae must be treated with caution in applications.     

Application of environmental groundwater tracers at the Sulphur Bank Mercury Mine, California, USA

Boron, chloride, sulfate, δD, δ¹⁸O, and ³H concentrations in surface water and groundwater samples from the Sulphur Bank Mercury Mine (SBMM), California, USA were used to examine geochemical processes and provide constraints on evaporation and groundwater flow. SBMM is an abandoned sulfur and mercury mine with an underlying hydrothermal system, adjacent to Clear Lake, California. Results for non-³H tracers (i.e., boron, chloride, sulfate, δD, and δ¹⁸O) identify contributions from six water types at SBMM. Processes including evaporation, mixing, hydrothermal water input and possible isotopic exchange with hydrothermal gases are also discerned. Tritium data indicate that hydrothermal waters and other deep groundwaters are likely pre-bomb (before ~1952) in age while most other waters were recharged after ~1990. A boron-based steady-state reservoir model of the Herman Impoundment pit lake indicates that 71–79% of its input is from meteoric water with the remainder from hydrothermal contributions. Results for groundwater samples from six shallow wells over a 6–month period for δD and δ¹⁸O suggests that water from Herman Impoundment is diluted another 3% to more than 40% by infiltrating meteoric water, as it leaves the site. Results for this investigation show that environmental tracers are an effective tool to understand the SBMM hydrogeologic regime.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé Les concentrations en bore, chlorure, sulfate, δD, δ¹⁸O, et ³H d’échantillons d’eaux de surface et souterraine prélevés dans le banc de soufre de la mine de mercure (SBMM en anglais) en Californie, USA, ont été utilisées pour étudier les processus géochimiques et pour fournir des contraintes à l’évaporation et à l’écoulement des eaux souterraines. La SBMM est une mine de soufre et de mercure abandonnée, adjacente au lac Clear en Californie et sous laquelle se trouve un système hydrothermal.Les résultats des traceurs autres que le tritium (bore, chlorure, sulfate, δD, et δ¹⁸O) ont permis d’identifier des contributions de six types d’eaux à SBMM. Des processus tels que l’évaporation, le mélange, l’entrée d’eau hydrothermale et de possibles échanges isotopiques avec des gaz hydrothermaux ont également été identifiés. Les données de tritium montrent que les eaux hydrothermales et d’autres eaux profondes sont probablement d’age antérieure à la bombe (avant ~1952), alors que la plupart des autres eaux sont issues de la recharge après ~1990. Un modèle de réservoir représentant le lac situé dans la partie Herman Impoundment de l’ancienne mine, en régime permanent et basé sur le bore, montre que 71–79% de l’eau provient des précipitations, le reste provenant de contributions hydrothermales. Les résultats de δD et δ¹⁸O pour des échantillons d’eau souterraine de six puits peu profonds sur une période de 6 mois suggèrent que l’eau de Hermann Impoundment est encore diluée entre 3% jusqu’à plus de 40% lorsqu’elle quitte le site du fait de l’infiltration d’eau météorique. Les résultats de cette étude montrent que les traceurs environnementaux constituent un outil efficace pour comprendre le régime hydrogéologique de la SBMM.

---

Resumen Se han utilizado datos de concentraciones de boro, cloruros, sulfatos, δD, δ¹⁸O, y ³H en muestras de aguas superficiales y subterráneas procedentes de la Mina de Mercurio Sulphur Bank (SBMM), California, USA, para estudiar los procesos geoquímicos y caracterizar la evaporación y el flujo del agua subterránea. SBMM es una mina abandonada de azufre y mercurio con un sistema hidrotermal subyacente, cercano al Lago Clear, California.Los resultados de los trazadores que no son ³H (por ejemplo, boro, cloruros, sulfatos δD, y δ¹⁸O) identifican las contribuciones de seis tipos de agua en la SBMM. Se han identificado diversos procesos, como evaporación, mezcla, entrada de agua hidrotermal y la posibilidad de intercambios isotópicos con gases hidrotermales. Los datos de tritio indican que las aguas hidrotermales y otras aguas subterráneas profundas son probablemente pre-bombas (previas a ~1952) en edad mientras que la mayoría de las otras aguas se han recargado después de ~1990. Un modelo estacionario basado en el boro de la fosa del lago Herman Impoundment indica que el 71–79% de su entrada procede de agua meteórica con restos de contribución hidrotermal. Los resultados para las muestras de aguas subterráneas procedentes de seis pozos superficiales en un periodo de 6 meses para δD y δ¹⁸O sugieren que el agua procedente de Herman Impoundment se diluye en otro 3% hasta más del 40% debido a la infiltración de agua meteórica., que sale del sitio. Los resultados de esta investigación muestran que los trazadores ambientales son una herramienta efectiva para entender el régimen hidrogeológico de la SBMM.

     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

High heat flow anomalies on an old oceanic plate observed seaward of the Japan Trench

Thirty-three new measurements on the seaward slope and outer rise of the Japan Trench along a parallel of 38°45′N revealed the existence of high heat flow anomalies on the subducting Pacific plate, where the seafloor age is about 135 m.y.. The most prominent anomaly with the highest value of 114 mW/m² is associated with a small mound on the outer rise, which was reported to be a kind of mud volcano. On the seaward slope of the trench, heat flow is variable: high (70–90 mW/m²) at some locations and normal for the seafloor age (about 50 mW/m²) at others. The spatial variation of heat flow may be related to development of normal faults and horst/graben structures due to bending of the Pacific plate before subduction, with fluid flow along the fault zones enhancing the vertical heat transfer. Possible heat sources of the high heat flow anomalies are intra-plate volcanism in the last several million years like that discovered recently on the Pacific plate east of the Japan Trench.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Fracture analysis in the south-western Corinth rift (Greece) and implications on fault hydraulic behavior

This paper reviews the data concerning the fracture network and the hydraulic characteristics of faults in an active zone of the Gulf of Corinth. Pressure gap measured through fault planes shows that in this area the active normal faults (Aigion, Helike) act, at least temporarily and locally, as transversal seal. The analysis of the carbonate cements in the fractures on both the hangingwall and the footwall of the faults also suggests that they have acted as local seals during the whole fault zone evolution. However, the pressure and the characteristics of the water samples measured in the wells indicate that meteoric water circulates from the highest part of the relief to the coast, which means it goes through the fault zones. Field quantitative analysis and core studies from the AIG-10 well have been performed to define both regional and fault-related fracture networks. Then laboratory thin section observations have been done to recognize the different fault rocks characterizing the fault zone components. These two kinds of approach give information on the permeability characteristics of the fault zone. To synthesize the data, a schematic conceptual 3D fluid flow modeling has been performed taking into account fault zone permeability architecture, sedimentation, fluid flow, fault vertical offset and meteoric water influx, as well as compaction water flow. This modeling allows us to fit all the data with a model where the fault segments act as a seal whereas the relays between these segments allow for the regional flow from the Peloponnese topographic highs to the coast.     

Statistical evaluation of methods for the calculation of static formation temperatures in geothermal and oil wells using an extension of the error propagation theory

The feasibility of using the more sophisticated weighted least-squares (WLS) model, as opposed to the traditional ordinary least-squares (OLS), in linear regressions of BHT data to estimate the static formation temperatures (SFT) was investigated. The most commonly used analytical methods (line-source; spherical and radial heat flow; and cylindrical heat source) were applied. Error propagation equations were derived to calculate errors in the time function of each method. These errors were combined with the BHT measurement errors to compute weighting factors for applying the WLS. Intercept uncertainties were estimated for all regressions using sets of synthetic and actual borehole logs taken from geothermal and oil applications. SFT computed with the spherical and radial heat flow method were generally greater than those from the other two methods.     

Temperature–composition–time (T–X–t) data from authigenic K-feldspar: An integrated methodology for dating fluid flow events

Integration of fluid inclusion analysis with high spatial resolution Ar–Ar dating of K-feldspar cements has been used to resolve and reconstruct palaeo-fluid flow. Fluid inclusion analysis allows discrimination of distinct cement phases, thereby identifying discrete episodes of fluid flow. Ar–Ar dating of the same cements via high spatial resolution laserprobe establishes absolute age constraints on the framework previously constructed. Integration of these two datasets yields temperature–composition–time data.     

Up-temperature flow of surface-derived fluids in the mid-crust: the role of pre-orogenic burial of hydrated fault rocks

The Walter‐Outalpa shear zone in the southern Curnamona Province of NE South Australia is an example of a shear zone that has undergone intensely focused fluid flow and alteration at mid‐crustal depths. Results from this study have demonstrated that the intense deformation and ductile shear zone reactivation, at amphibolite facies conditions of 534 ± 20 °C and 500 ± 82 MPa, that overprint the Proterozoic Willyama Supergroup occurred during the Delamerian Orogeny (c. 500 Ma) (EPMA monazite ages of 501 ± 16 and 491 ± 19 Ma). This is in contrast to the general belief that the majority of basement deformation and alteration in the southern Curnamona Province occurred during the waning stages of the Olarian Orogeny (c. 1610–1580 Ma). These shear zones contain hydrous mineral assemblages that cut wall rocks that have experienced amphibolite facies metamorphism during the Olarian Orogeny. The shear zone rock volumes have much lower δ¹⁸O values (as low as 1‰) than their unsheared counterparts (7–9‰), and calculated fluid δ¹⁸O values (5–8‰) consistent with a surface‐derived fluid source. Hydrous minerals show a decrease in δD_(H2O) from ?14 to ?22‰, for minerals outside the shear zones, to ?28 to ?40‰, for minerals within the shear zones consistent with a contribution from a meteoric source. It is unclear how near‐surface fluids initially under hydrostatic pressure penetrate into the middle crust where fluid pressures approach lithostatic, and where fluid flow is expected to be dominantly upward because of pressure gradients. We propose a mechanism whereby faulting during basin formation associated with the Adelaidean Rift Complex (c. 700 Ma) created broad hydrous zones containing mineral assemblages in equilibrium with surface waters. These panels of fault rock were subsequently buried to depths where the onset of metamorphism begins to dehydrate the fault rock volumes evolving a low δ¹⁸O fluid that is channelled through shear zones related to Delamerian Orogenic activity.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.