Closed Analytic Elements with Flexible Geometry

A simple and fast treatment of hydrogeologic features with irregularly shaped boundaries in two‐dimensional analytic element groundwater flow models is presented. The star domain shapes of the features are restricted to closed shapes represented as smooth and continuous single‐valued functions of distance from a focus point, . The element can be used to treat a variety of boundary and continuity conditions, including those of irregularly shaped lakes or heterogeneities in hydraulic conductivity. The new element is demonstrated via some simple illustrative test cases and shown to be efficient, accurate, and much simpler to implement than existing solutions for irregular shapes.     

Disconnected Surface Water and Groundwater: From Theory to Practice

When describing the hydraulic relationship between rivers and aquifers, the term disconnected is frequently misunderstood or used in an incorrect way. The problem is compounded by the fact that there is no definitive literature on the topic of disconnected surface water and groundwater. We aim at closing this gap and begin the discussion with a short introduction to the historical background of the terminology. Even though a conceptual illustration of a disconnected system was published by Meinzer (1923) , it is only within the last few years that the underlying physics of the disconnection process has been described. The importance of disconnected systems, however, is not widely appreciated. Although rarely explicitly stated, many approaches for predicting the impacts of groundwater development on surface water resources assume full connection. Furthermore, management policies often suggest that surface water and groundwater should only be managed jointly if they are connected. However, although lowering the water table beneath a disconnected section of a river will not change the infiltration rate at that point, it can increase the length of stream that is disconnected. Because knowing the state of connection is of fundamental importance for sustainable water management, robust field methods that allow the identification of the state of connection are required. Currently, disconnection is identified by showing that the infiltration rate from a stream to an underlying aquifer is independent of the water table position or by identifying an unsaturated zone under the stream. More field studies are required to develop better methods for the identification of disconnection and to quantify the implications of heterogeneity and clogging processes in the streambed on disconnection.     

Discontinuous Galerkin methods for modeling Hurricane storm surge

Storm surge due to hurricanes and tropical storms can result in significant loss of life, property damage, and long-term damage to coastal ecosystems and landscapes. Computer modeling of storm surge can be used for two primary purposes: forecasting of surge as storms approach land for emergency planning and evacuation of coastal populations, and hindcasting of storms for determining risk, development of mitigation strategies, coastal restoration and sustainability.Storm surge is modeled using the shallow water equations, coupled with wind forcing and in some events, models of wave energy. In this paper, we will describe a depth-averaged (2D) model of circulation in spherical coordinates. Tides, riverine forcing, atmospheric pressure, bottom friction, the Coriolis effect and wind stress are all important for characterizing the inundation due to surge. The problem is inherently multi-scale, both in space and time. To model these problems accurately requires significant investments in acquiring high-fidelity input (bathymetry, bottom friction characteristics, land cover data, river flow rates, levees, raised roads and railways, etc.), accurate discretization of the computational domain using unstructured finite element meshes, and numerical methods capable of capturing highly advective flows, wetting and drying, and multi-scale features of the solution.The discontinuous Galerkin (DG) method appears to allow for many of the features necessary to accurately capture storm surge physics. The DG method was developed for modeling shocks and advection-dominated flows on unstructured finite element meshes. It easily allows for adaptivity in both mesh (h) and polynomial order (p) for capturing multi-scale spatial events. Mass conservative wetting and drying algorithms can be formulated within the DG method.In this paper, we will describe the application of the DG method to hurricane storm surge. We discuss the general formulation, and new features which have been added to the model to better capture surge in complex coastal environments. These features include modifications to the method to handle spherical coordinates and maintain still flows, improvements in the stability post-processing (i.e. slope-limiting), and the modeling of internal barriers for capturing overtopping of levees and other structures. We will focus on applications of the model to recent events in the Gulf of Mexico, including Hurricane Ike.     

210Pb in the Pacific: the GEOSECS measurements of particulate and dissolved concentrations

We report here on particulate and dissolved²¹⁰Pb profiles at 16 stations, and on total²¹⁰Pb profiles at 3 stations, all occupied during the Pacific GEOSECS expedition. Comparison with measurements at Yale on GEOSECS library samples indicates that during separation of particulate lead from dissolved lead, our filtered water samples suffered some loss of²¹⁰Pb in the filtration system; this effect appears to have reduced the dissolved²¹⁰Pb activities by ～ 20% in stations where the water was filtered. However, for these first Pacific data on the²¹⁰Pb distribution between the two phases, this effect does not significantly interfere with our recognition of the major features of both particulate and dissolved²¹⁰Pb distributions.The dissolved²¹⁰Pb profiles in general vary geographically, following the²²⁶Ra profiles. In deep water,²²⁶Ra increases northward and eastward from the southwest Pacific, from ～ 22dpm/100kg, to over 40 dpm/100 kg in the northeast Pacific. Our dissolved²¹⁰Pb profiles show a similar increase in deep water, varying from about 10 to 20 dpm/100 kg along this line, and are commonly characterized by a mid-depth maximum. This²¹⁰Pb maximum reflects the mid-depth²²⁶Ra maximum of the Pacific Deep Water observed along the western boundary current.In surface water at low latitudes there is a significant²¹⁰Pb flux from the atmosphere, which produces a²¹⁰Pb/²²⁶Ra activity ratio generally greater than unity. This flux penetrates as deep as 600 m, as indicated by an “induced”²¹⁰Pb minimum caused by the surface maximum. The surface water²¹⁰Pb excess decreases toward high southern latitudes and vanishes in the Circumpolar region.The particulate²¹⁰Pb profiles show a general increase with depth, from ～ 0.3dpm/100kg in subsurface water to ～ 1.5dpm/100kg in bottom water, with or without a mid-depth maximum that reflects the²²⁶Ra or dissolved²¹⁰Pb maximum. The particulate²¹⁰Pb normally comprises about 2% of the total²¹⁰Pb in subsurface water, and this fraction increases to about 10% near the bottom. As the filtration loss is not taken into account, the fraction of particulate²¹⁰Pb quoted here is an upper limit. Since the particulate matter concentrations are quite uniform in the water column below a few hundred meters, the²¹⁰Pb activity of the particulate matter also increases with depth. The particulate matter has a²¹⁰Pb concentration of ～ 100dpm/g in subsurface water, but the concentration increases to ～ 500dpm/g or more toward the bottom. This indicates that there is a cumulative adsorption of Pb onto the suspended particles as they are sinking through the water column.     

Total mercury,methyl mercury and sulphide in river carron sediments

Total mercury, methyl mercury and sulphide contents of River Carron sediments (Lothian, Scotland) have been determined. Total mercury concentrations are comparable to other mercury polluted estuaries in the UK, but the methyl mercury values for low-sulphide Carron sediments are generally higher. It has been found that methyl mercury levels are initially in direct proportion to the sulphide concentrations of the sediments but beyond sulphide concentrations of 1.8 mg g^?1 the methyl mercury levels decline sharply.     

Estimating Flow Using Tracers and Hydraulics in Synthetic Heterogeneous Aquifers

Regional ground water flow is most usually estimated using Darcy's law, with hydraulic conductivities estimated from pumping tests, but can also be estimated using ground water residence times derived from radioactive tracers. The two methods agree reasonably well in relatively homogeneous aquifers but it is not clear which is likely to produce more reliable estimates of ground water flow rates in heterogeneous systems. The aim of this paper is to compare bias and uncertainty of tracer and hydraulic approaches to assess ground water flow in heterogeneous aquifers. Synthetic two-dimensional aquifers with different levels of heterogeneity (correlation lengths, variances) are used to simulate ground water flow, pumping tests, and transport of radioactive tracers. Results show that bias and uncertainty of flow rates increase with the variance of the hydraulic conductivity for both methods. The bias resulting from the nonlinearity of the concentration–time relationship can be reduced by choosing a tracer with a decay rate similar to the mean ground water residence time. The bias on flow rates estimated from pumping tests is reduced when performing long duration tests. The uncertainty on ground water flow is minimized when the sampling volume is large compared to the correlation length. For tracers, the uncertainty is related to the ratio of correlation length to the distance between sampling wells. For pumping tests, it is related to the ratio of correlation length to the pumping test's radius of influence. In regional systems, it may be easier to minimize this ratio for tracers than for pumping tests.     

Coordinate mapping of analytical contaminant transport solutions to non-uniform flow fields

Existing analytical solutions to 2D and 3D contaminant transport problems are limited by the mathematically convenient assumption of uniform flow. An approximate method is developed herein for coordinate mapping of 2D (vertically-averaged) transport solutions to non-uniform steady-state irrotational and divergence-free flow fields in single-layer aquifers. The method enables existing analytical transport solutions to be applied to aquifer systems with wells, non-uniform saturated thickness, surface water features, and (to a limited degree) heterogeneous hydraulic conductivity and recharge. This mass-conservative coordinate mapping approach is inexact in its approximation of the dispersion process but is still sufficiently accurate for many simple flow systems. The degree of model error is directly proportional to the variation of velocity magnitude within the domain. These mapped analytical solutions are compared to numerical simulation results and the coordinate mapping errors are investigated. The methods described herein may be used in the traditional capacity of analytical transport models, i.e., screening and preliminary site assessment, without sacrificing accuracy by assuming locally uniform flow conditions or applying an ad-hoc coordinate transformation. The solutions benefit from the traditional advantages of analytical methods, particularly the removal of artifacts due to spatial and temporal discretization: no time-stepping or numerical discretization is required.     

Triggered earthquakes and deep well activities

Earthquakes can be triggered by any significant perturbation of the hydrologic regime. In areas where potentially active faults are already close to failure, the increased pore pressure resulting from fluid injection, or, alternatively, the massive extraction of fluid or gas, can induce sufficient stress and/or strain changes that, with time, can lead to sudden catastrophic failure in a major earthquake. Injection-induced earthquakes typically result from the reduction in frictional strength along preexisting, nearby faults caused by the increased formation fluid pressure. Earthquakes associated with production appear to respond to more complex mechanisms of subsidence, crustal unloading, and poroelastic changes in response to applied strains induced by the massive withdrawal of subsurface material. As each of these different types of triggered events can occur up to several years after well activities have begun (or even several years after all well activities have stopped), this suggests that the actual triggering process may be a very complex combination of effects, particularly if both fluid extraction and injection have taken place locally. To date, more than thirty cases of earthquakes triggered by well activities can be documented throughout the United States and Canada. Based on these case histories, it is evident that, owing to preexisting stress conditions in the upper crust, certain areas tend to have higher probabilities of exhibiting such induced seismicity.