Michigan basin regional ground water flow discharge to three Great Lakes

Ground water discharge to the Great Lakes around the Lower Peninsula of Michigan is primarily from recharge in riparian basins and proximal upland areas that are especially important to the northern half of the Lake Michigan shoreline. A steady-state finite-difference model was developed to simulate ground water flow in four regional aquifers in Michigan's Lower Peninsula: the Glaciofluvial, Saginaw, Parma-Bayport, and Marshall aquifers interlayered with the Till/"red beds," Saginaw, and Michigan confining units, respectively. The model domain was laterally bound by a continuous specified-head boundary, formed from lakes Michigan, Huron, St. Clair, and Erie, with the St. Clair and Detroit River connecting channels. The model was developed to quantify regional ground water flow in the aquifer systems using independently determined recharge estimates. According to the flow model, local stream stages and discharges account for 95% of the overall model water budget; only 50% enters the lakes directly from the ground water system. Direct ground water discharge to the Great Lakes' shorelines was calculated at 36 m3/sec, accounting for 5% of the overall model water budget. Lowland areas contribute far less ground water discharge to the Great Lakes than upland areas. The model indicates that Saginaw Bay receives only approximately 1.13 m3/sec ground water; the southern half of the Lake Michigan shoreline receives only approximately 2.83 m3/sec. In contrast, the northern half of the Lake Michigan shoreline receives more than 17 m3/sec from upland areas.     

Estimating deep recharge rates beneath an interlobate moraine using temperature logs

The Sandilands area of southeastern Manitoba contains an interlobate moraine that is a major ground water recharge area. Underlying the highly permeable sediments of the moraine are up to 100 m of till and the subcrop of the Winnipeg Formation, which contains a major sandstone aquifer. Ground water flow within the till is examined using high-resolution temperature profiles and solutions to the differential equation for heat flow in porous media. These analyses indicate that recharge to the sandstone aquifer is occurring at a rate of approximately 2 x 10(-8) m/sec beneath the moraine, which is in agreement with recharge rates determined by conventional ground water hydraulics (10(-7) to 10(-10)(m/sec) and another study using multiple environmental tracers (1 x 10(-9) to 6 X 10(-9) m/sec). The use of temperature to determine ground water flux is not limited by half-lives as many environmental tracers are, and this allows for cost-effective estimation of recharge and discharge rates over longer periods.     

Geochemical processes during five years of aquifer storage recovery

A key factor in the long-term viability of aquifer storage recovery (ASR) is the extent of mineral solution interaction between two dissimilar water types and consequent impact on water quality and aquifer stability. We collected geochemical and isotopic data from three observation wells located 25, 65, and 325 m from an injection well at an experimental ASR site located in a karstic, confined carbonate aquifer in South Australia. The experiment involved five major injection cycles of a total of 2.5 x 10(5) m3 of storm water (total dissolved solids [TDS] approximately 150 mg/L) into the brackish (TDS approximately 2400 mg/L) aquifer. Approximately 60% of the mixture was pumped out during the fifth year of the experiment. The major effect on water quality within a 25 m radius of the injection well following injection of storm water was carbonate dissolution (35 +/- 6 g of CaCO3 dissolved/m3 of aquifer) and sulfide mineral oxidation (50 +/- 10 g as FeS2/m3 after one injection). < 0.005% of the total aquifer carbonate matrix was dissolved during each injection event, and approximately 0.2% of the total reduced sulfur. Increasing amounts of ambient ground water was entrained into the injected mixture during each of the storage periods. High 14C(DIC) activities and slightly more negative delta13C(DIC) values measured immediately after injection events show that substantial CO2(aq) is produced by oxidation of organic matter associated with injectant. There were no detectable geochemical reactions while pumping during the recovery phase in the fifth year of the experiment.     

Geochemical imaging of flow near an artificial recharge facility, Orange County, California

Critical for the management of artificial recharge operations is detailed knowledge of ground water dynamics near spreading areas. Geochemical tracer techniques including stable isotopes of water, tritium/helium-3 (T/3He) dating, and deliberate gas tracer experiments are ideally suited for these investigations. These tracers were used to evaluate flow near an artificial recharge site in northern Orange County, California, where approximately 2.5 x 10(8) m3 (200,000 acre-feet) of water are recharged annually. T/3He ages show that most of the relatively shallow ground water within 3 km of the recharge facilities have apparent ages < 2 years; further downgradient apparent ages increase, reaching > 20 years at approximately 6 km. Gas tracer experiments using sulfur hexafluoride and xenon isotopes were conducted from the Santa Ana River and two spreading basins. These tracers were followed in the ground water for more than two years, allowing subsurface flow patterns and flow times to be quantified. Results demonstrate that mean horizontal ground water velocities range from < 1 to > 4 km/year. The leading edges of the tracer patch moved at velocities about twice as fast as the center of mass. Leading edge velocities are important when considering the potential transport of microbes and other "time sensitive" contaminants and cannot be determined easily with other methods. T/3He apparent ages and tracer travel times agreed within the analytical uncertainty at 16 of 19 narrow screened monitoring wells. By combining these techniques, ground water flow was imaged with time scales on the order of weeks to decades.     

Rapid evolution of redox processes in a petroleum hydrocarbon-contaminated aquifer

Ground water chemistry data collected over a six-year period show that the distribution of contaminants and redox processes in a shallow petroleum hydrocarbon-contaminated aquifer has changed rapidly over time. Shortly after a gasoline release occurred in 1990, high concentrations of benzene were present near the contaminant source area. In this contaminated zone, dissolved oxygen in ground water was depleted, and by 1994 Fe(III) reduction and sulfate reduction were the predominant terminal electron accepting processes. Significantly, dissolved methane was below measurable levels in 1994, indicating the absence of significant methanogenesis. By 1996, however, depletion of solid-phase Fe(III)-oxyhydrox ides in aquifer sediments and depletion of dissolved sulfate in ground water resulted in the onset of methanogenesis. Between 1996 and 2000, water-chemistry data indicated that methanogenic metabolism became increasingly prevalent. Molecular analysis of 16S-rDNA extracted from sediments shows the presence of a more diverse methanogenic community inside as opposed to outside the plume core, and is consistent with water-chemistry data indicating a shift toward methanogenesis over time. This rapid evolution of redox processes reflects several factors including the large amounts of contaminants, relatively rapid ground water flow (approximately 0.3 m/day [approximately foot/day]), and low concentrations of microbially reducible Fe(III) oxyhydroxides ( approximately 1 micromol/g) initially present in aquifer sediments. These results illustrate that, under certain hydrologic conditions, redox conditions in petroleum hydrocarbon-contaminated aquifers can change rapidly in time and space, and that the availability of solid-phase Fe(III)-oxyhydroxides affects this rate of change.     

Estimate of recharge of a rising water table in semiarid Niger from 3H and 14C modeling

A hydrodynamic survey carried out in semiarid southwest Niger revealed an increase in the unconfined ground water reserves of approximately 10% over the last 50 years due to the clearing of native vegetation. Isotopic samplings (3H, 18O, 2H for water and 14C, 13C for the dissolved inorganic carbon) were performed on about 3500 km2 of this silty aquifer to characterize recharge. Stable isotope analyses confirmed the indirect recharge process that had already been shown by hydrodynamic surveys and suggested the tracers are exclusively of atmospheric origin. An analytical model that takes into account the long-term rise in the water table was used to interpret 3H and 14C contents in ground water. The natural, preclearing median annual renewal rate (i.e., recharge as a fraction of the saturated aquifer volume) lies between 0.04% and 0.06%. For representative characteristics of the aquifer (30 m of saturated thickness, porosity between 10% and 25%), this implies a recharge of between 1 and 5 mm/year, which is much lower than the estimates of 20 to 50 mm/year for recent years, obtained using hydrological and hydrodynamic methods and the same aquifer parameters. Our study, therefore, reveals that land clearing in semiarid Niger increased ground water recharge by about one order of magnitude.     

Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA

Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of ‘interactive’ ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d⁻¹. ³H/³He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d⁻¹) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.     

Gravitational effect of water circulation in the northwest Yunnan

ＧｒａｖｉｔａｔｉｏｎａｌｅｆｆｅｃｔｏｆｗａｔｅｒｃｉｒｃｕｌａｔｉｏｎｉｎｔｈｅｎｏｒｔｈｗｅｓｔＹｕｎｎａｎＭｉｎ－ＹｕＪＩＡ；（贾民育）Ｓｈａｏ－ＡｎＳＵＮ；（孙少安）Ａｉ－ＭｉｎＸＩＡＮＧ；（项爱民）ａｎｄＤｏｎｇ－ＺｈｉＬＩＵ（刘冬至）（...     

The Rn content of groundwater in strata of the Plio-Pleistocene Osaka Group in the Sennan area, Japan

Studies on the ²²²Rn content of groundwater were conducted to obtain knowledge of its distribution in the geosphere just below ground surface and below the depth of 100 m in the Osaka Group (Plio-Pleistocene). Samples analyzed from ten bore holes within 10 m depth below land surface showed that a sandy layer contained groundwater with consistent radon concentrations except for a few locations where the variations in the water level were large due to precipitation and/or the water level was within 2 m of ground surface. The average ²²²Rn concentration in groundwater within 10 m of land surface was found to be approximately equal to $\text{1}\text{3}$ to $\text{3}\text{5}$ of the highest concentration found in groundwater from a well penetrating 100–200 m (420 pCi/l^?1).     

Inflow and outflow in a rotating circular container—A laboratory model for the Loop Current of the Gulf of Mexico

The Loop Current of the Gulf of Mexico is simulated in the laboratory. A circular tank is filled with water and is placed off-center on a rotating table and the flow field is generated by injecting and withdrawing water at two openings on the wall. The free surface becomes parabolic due to balance of gravitational and centrifugal forces, simulating the latitudinal change of the Coriolis parameter (-effect) in the ocean. The flow characteristics depend on the influx and the rate of rotation and can be classified according to non-dimensional parameters (Rossby, Ekman and Froude numbers denoted byR ₀,E andF, respectively). When the influx is small and the rotation rate is large (smallR ₀,E andF) the flow will be almost linear, and the fluid flows along the side-wall boundary layer under constraint of the -effect. For a very large influx (largeR ₀ andE) inertial forces become very large compared to the Coriolis force and the flow behaves like a potential flow. The flow studied had characteristics between these two extreme cases and hasR ₀ andF similar to the Gulf circulation, though similarity inE is ambiguous. Photographs of the flow indicate that the inflow penetrates further into the interior when the rotation rate is increased while the influx is kept constant. The numerical analysis of the non-linear vorticity equation confirms this for the parameters corresponding to the experiment. In addition, the photographs reveal eddies embedded on both sides of the main stream, particularly near the inflow region. These eddies are intensified and become uniform in size as the influx increases. It is pointed out that such eddies were actually observed near the Loop Current north of the Yucatan Straits.     

Uplift and submarine formation of some Melanesian porphyry copper deposits: Stable isotope evidence

Hydrogen and oxygen isotope analyses of sericites and kaolinites from four young porphyry copper deposits (Ok Tedi (1.2 Ma) and Yandera (6.5 Ma), Papua New Guinea; Koloula (1.5 Ma), Solomon Islands; and Waisoi (<5 Ma), Fiji) indicate that the fluids from which these minerals precipitated were of mixed magmatic and non-magmatic sources. The non-magmatic component of the fluid from the island arc deposits (Koloula, Waisoi) was ocean water.For Ok Tedi, the non-magmatic component was a meteoric water with an isotopic composition different from that of the present meteoric water in the region. The isotopic signature of the former meteoric water is consistent with a surface elevation of 200 m a.s.l. or less at the time of mineralization. The deposit was later exposed and supergene kaolinitization commenced at approximately 1200 m a.s.l. Uplift and erosion has continued to the present at which time the elevation of the exposed deposit is 1800 m a.s.l. This rate of uplift is consistent with that known from other geological evidence. If the rate of uplift were approximately constant during the last 1.2 Ma, the age of supergene enrichment can be dated at approximately 0.4 Ma B.P.Similarly, influx of meteoric water at Yandera occurred when the ground surface above the deposit was at an elevation of approximately 600 m a.s.l. The deposit's present elevation is 1600 m a.s.l. In this case a total uplift of approximately 2.2 km is indicated, with removal of 1.2 km of overburden by erosion.     

A Geostatistically Based Ground Water Monitoring Study of Nonpoint Source No3--N Concentrations

Geostatistical interpretations of ground water monitoring data are presented to define the spatial distributions of NO₃^--N in the ground water at two demonstration test sites in the Idaho Snake River Plain. Sequential Gaussian simulation was used to delineate monthly ground water NO₃^--N changes during and after implementation of a prescribed crop rotation at test site 1. Trend surface analyses were used to illustrate monthly ground water NO₃^--N changes during and after a prescribed irrigation practice was implemented at test site 2. These evaluations suggest that geostatistically based ground water monitoring can be effective in the delineation of changes in ground water quality in shallow, unconfined aquifers in agricultural areas such as those in southern Idaho. Geostatistical methods showed spatial and temporal changes in ground water NO₃^--N inferred to be a result of the agricultural practices implemented.     

Calibration of base flow separation methods with streamflow conductivity

The conductivity mass-balance (CMB) method can be used to calibrate analytical base flow separation methods. The principal CMB assumptions are base flow conductivity is equal to streamflow conductivity at lowest flows, runoff conductivity is equal to streamflow conductivity at highest flows, and base flow and runoff conductivities are assumed to be constants over the period of record. To test the CMB assumptions, fluid conductivities of ground water, surface runoff, and streamflow were measured during wet and dry conditions in a 12-km(2) stream basin. Ground water conductivities at wells varied an average of 6% from dry to wet conditions, while stream conductivities varied 58%. Shallow ground water conductivity varied significantly with distance from the stream, with lowest conductivities of 87 microS/cm near the divide, a maximum of 520 microS/cm 59 m from the stream, and 215 microS/cm 22 m from the stream. Runoff conductivities measured in three rain events remained nearly constant, with lower conductivities of 35 microS/cm near the divide and 50 microS/cm near the stream. The CMB method was applied to the records from 10 USGS stream-gauging stations in Texas, Kentucky, Georgia, and Florida to calibrate the USGS base flow separation technique, HYSEP, by varying the time parameter 2N*. There is a statistically significant relationship between basin areas and calibrated values of 2N*, expressed as N = 0.46A(0.44), with N in days and A in km(2). The widely accepted relationship N = 0.83A(0.2) is not valid for these basins. Other analytic methods can also be calibrated with the CMB method.     

Water percolation in a thick unsaturated loess layer considering the ground‐atmosphere interaction

The key objective of this paper is to advance our present understanding of how surface water infiltrates in thick unsaturated loess, which is found in arid and semiarid regions of the world, considering the ground‐atmosphere interaction. In situ data for a period of 1 year in thick loess layer at a site in the Loess Plateau of China that has groundwater table at 97.5 m depth were collected for achieving this objective. Climate factors, mainly rainfall and actual evaporation, were measured. In addition, variations of soil temperature and water content at different depths in the unsaturated zone were also measured. The data were used to interpret the water percolation characteristics by dividing the thick unsaturated zone into three zones; namely, (i) surface zone, which constitutes the top 1.0 m, (ii) unsteady zone, which is from 1.0 to 7.0 m, and (iii) steady zone, which is below 7.0 m. In the surface zone, soil temperature and water content are sensitive to climate factors. There is a variation of water content associated with the cumulative influence of infiltration and evaporation in the precipitation and nonprecipitation periods, respectively. In the unsteady zone, the water content is relatively constant; however, temperature varies in different seasons. Water percolation in this zone is both in liquid and vapour phases. In the steady zone, both soil temperature and water content are constant during the entire investigation period. The percolation velocity in this zone is approximately 1.23 × 10^?8 m/s or 0.39 m/year, which suggests that it will take approximately 230.8 years for surface water to pass through the thick unsaturated zone and recharge the groundwater.     

The effect of terrace geology on ground‐water movement and on the interaction of ground water and surface water on a mountainside near Mirror Lake,New Hampshire,USA

The west watershed of Mirror Lake in the White Mountains of New Hampshire contains several terraces that are at different altitudes and have different geologic compositions. The lowest terrace (FSE) has 5 m of sand overlying 9 m of till. The two next successively higher terraces (FS2 and FS1) consist entirely of sand and have maximum thicknesses of about 7 m. A fourth, and highest, terrace (FS3) lies in the north‐west watershed directly adjacent to the west watershed. This highest terrace has 2 m of sand overlying 8 m of till. All terraces overlie fractured crystalline bedrock. Numerical models of hypothetical settings simulating ground‐water flow in a mountainside indicated that the presence of a terrace can cause local ground‐water flow cells to develop, and that the flow patterns differ based on the geologic composition of the terrace. For example, more ground water moves from the bedrock to the glacial deposits beneath terraces consisting completely of sand than beneath terraces that have sand underlain by till. Field data from Mirror Lake watersheds corroborate the numerical experiments. The geology of the terraces also affects how the stream draining the west watershed interacts with ground water. The stream turns part way down the mountainside and passes between the two sand terraces, essentially transecting the movement of ground water down the valley side. Transects of water‐table wells were installed across the stream's riparian zone above, between, and below the sand terraces. Head data from these wells indicated that the stream gains ground water on both sides above and below the sand terraces. However, where it flows between the sand terraces the stream gains ground water on its uphill side and loses water on its downhill side. Biogeochemical processes in the riparian zone of the flow‐through reach have resulted in anoxic ground water beneath the lower sand terrace. Results of this study indicate that it is useful to understand patterns of ground‐water flow in order to fully understand the flow and chemical characteristics of both ground water and surface water in mountainous terrain. Copyright © 2007 John Wiley & Sons, Ltd.     

利用解析法有效快速估计将来GRACE Follow-On地球重力场的精度

本文首次利用解析法有效快速估计了将来GRACE（Gravity Recovery and Climate Experiment） Follow-On地球重力场的精度. 第一,基于功率谱原理分别建立了新的GRACE Follow-On卫星激光干涉星间测量系统星间速度、GPS接收机轨道位置和轨道速度以及加速度计非保守力误差影响累计大地水准面的单独和联合解析误差模型. 第二,利用提出的GRACE卫星关键载荷匹配精度指标和美国喷气推进实验室（JPL）公布的GRACE Level 1B实测精度指标的一致性,以及估计的GRACE累计大地水准面精度和德国波兹坦地学研究中心（GFZ）公布的EIGEN-GRACE02S地球重力场模型实测精度的符合性,验证了本文建立的解析误差模型是可靠的. 第三,论证了GRACE Follow-On卫星不同关键载荷匹配精度指标和轨道高度对地球重力场精度的影响. 在360阶处,利用轨道高度250 km、星间距离50 km、星间速度误差1×10^-9m/s、轨道位置误差3×10^-5m、轨道速度误差3×10^-8m/s和非保守力误差3×10^-13m/s²,基于联合解析误差模型估计累计大地水准面的精度为1.231×10^-1 m. 本文的研究不仅为当前GRACE和将来GRACE Follow-On地球重力场精度的有效快速确定提供了理论基础和计算保证,同时对国际将来GRAIL（Gravity Recovery and Interior Laboratory）月球卫星重力测量计划的成功实施具有重要的参考意义.     

A first estimate of ground water ages for the deep aquifer of the Kathmandu Basin, Nepal, using the radioisotope chlorine-36

The Kathmandu Basin in Nepal contains up to 550 m of Pliocene-Quaternary fluvio-lacustrine sediments which have formed a dual aquifer system. The unconfined sand and gravel aquifer is separated by a clay aquitard, up to 200 m thick, from the deeper, confined aquifer, comprised of Pliocene sand and gravel beds, intercalated with clay, peat, and lignite. The confined aquifer currently provides an important water supply to the central urban area but there are increasing concerns about its sustainability due to overexploitation. A limited number of determinations of the radioisotope 36Cl have been made on bore waters in the basin, allowing us to postulate on the age of ground water in the deeper, confined aquifer. Ground water evolution scenarios based on radioisotope decay, gradual dissolution of formational salts as the ground waters move downgradient, and flow velocity estimations produce comparable ground water ages for the deep waters, ranging from 200,000 to 400,000 years. From these ages, we deduce a mean ground water flow velocity of only 45 mm/year from recharge in the northeast to the main extraction region 15 km to the southwest. We thus estimate current recharge at about 5 to 15 mm/year, contributing 40,000 to 1.2 million m3/year to the ground water system. Current ground water extraction is estimated to be 20 times this amount. The low specific discharge confirms that the resource is being mined, and, based on current projections, reserves will be used up within 100 years.     

Using temperature to test models of flow near Yucca Mountain,Nevada

Ground water temperatures in the fractured volcanic aquifer near Yucca Mountain, Nevada, have previously been shown to have significant spatial variability with regions of elevated temperatures coinciding roughly with near-vertical north-south trending faults. Using insights gained from one-dimensional models, previous investigators have suggested upwelling along faults from an underlying aquifer as a likely explanation for this ground water temperature pattern. Using a three-dimensional coupled flow and heat-transport model, we show that the thermal high coinciding with the Paintbrush fault zone can be explained without significant upwelling from the underlying aquifer. Instead, the thermal anomaly is consistent with thermal conduction enhanced slightly by vertical ground water movement within the volcanic aquifer sequence. If more than approximately 400 m3/day of water enters the volcanic aquifer from below along a 10 km fault zone, the calculated temperatures at the water table are significantly greater than the measured temperatures. These results illustrate the potential limitations in using one-dimensional models to interpret ground water temperature data, and underscore the value in combining temperature data with fully coupled three-dimensional simulations.     

Estimating ground water discharge by hydrograph separation

Iron Mountain is located in the West Shasta Mining District in California. An investigation of the generation of acid rock drainage and metals loading to Boulder Creek at Iron Mountain was conducted. As part of that investigation, a hydrograph separation technique was used to determine the contribution of ground water to total flow in Boulder Creek. During high-flow storm events in the winter months, peak flow in Boulder Creek can exceed 22.7 m3/sec, and comprises surface runoff, interflow, and ground water discharge. A hydrograph separation technique was used to estimate ground water discharge into Boulder Creek during high-flow conditions. Total ground water discharge to the creek approaches 0.31 m3/sec during the high-flow season. The hydrograph separation technique combined with an extensive field data set provided reasonable estimates of ground water discharge. These estimates are useful for other investigations, such as determining a corresponding metals load from the metal-rich ground water found at Iron Mountain and thus contributing to remedial alternatives.