Variably saturated groundwater modelling for optimizing managed aquifer recharge using trench infiltration

Spreading‐basin methods have resulted in more than 130 million cubic metres of recharge to the unconfined Navajo Sandstone of southern Utah in the past decade, but infiltration rates have slowed in recent years because of reduced hydraulic gradients and (or) clogging. Trench infiltration is a promising alternative technique for increasing recharge and minimizing evaporation. This paper uses a variably saturated flow model to further investigate the relative importance of the following variables on rates of trench infiltration to unconfined aquifers: saturated hydraulic conductivity, trench spacing and dimensions, initial water‐table depth, alternate wet/dry periods, and number of parallel trenches. Modelling results showed (1) increased infiltration with higher hydraulic conductivity, deeper initial water tables, and larger spacing between parallel trenches, (2) deeper or wider trenches do not substantially increase infiltration, (3) alternating wet/dry periods result in less overall infiltration than keeping the trenches continuously full, and (4) larger numbers of parallel trenches within a fixed area increases infiltration but with a diminishing effect as trench spacing becomes tighter. An empirical equation for estimating expected trench infiltration rates as a function of hydraulic conductivity and initial water‐table depth was derived and can be used for evaluating feasibility of trench infiltration in other hydrogeologic settings. Copyright © 2014 John Wiley & Sons, Ltd.     

Selecting best mapping strategies for storm runoff modeling in a mountainous semi‐arid area

Accurate runoff and soil erosion modeling is constrained by data availability, particularly for physically based models such as OpenLISEM that are data demanding, as the processes are calculated on a cell‐by‐cell basis. The first decision when using such models is to select mapping units that best reflect the spatial variability of the soil and hydraulic properties in the catchment. In environments with limited data, available maps are usually generic, with large units that may lump together the values of the soil properties, affecting the spatial patterns of the predictions and output values in the outlet. Conversely, the output results may be equally acceptable, following the principle of equifinality. To studyhow the mapping method selected affects the model outputs, four types of input maps with different degrees of complexity were created: average values allocated to general soil map units (ASG1), average values allocated to detailed map units (ASG2), values interpolated by ordinary kriging (OK) and interpolated by kriging with external drift (KED). The study area was Ribeira Seca, a 90 km² catchment located in Santiago Island, Cape Verde (West Africa), a semi‐arid country subject to scarce but extreme rainfall during the short tropical summer monsoon. To evaluate the influence of rainfall on runoff and erosion, two storm events with different intensity and duration were considered. OK and KED inputs produced similar results, with the latter being closer to the observed hydrographs. The highest soil losses were obtained with KED (43 ton ha^{? 1} for the strongest event). To improve the results of soil loss predictions, higher accurate spatial information on the processes is needed; however, spatial information of input soil properties alone is not enough in complex landscapes. The results demonstrate the importance of selecting the appropriate mapping strategy to obtain reliable runoff and erosion estimates. Copyright © 2013 John Wiley & Sons, Ltd.     

Measurement of radon potential from soil using a special method of sampling

Soil radon gas and/or its exhalation rate are used as indicators for some applications, such as uranium exploration, indoor radon concentration, seismic activity, location of subsurface faults, etc., and also in the studies where the main interest is the field verification of radon transport models. This work proposes a versatile method for the soil radon sampling using a special manner of pumping. The soil gas is passed through a column of charcoal by using passive pumping. A plastic bottle filled with water is coupled to an activated charcoal column and the flow of water through an adjustable hole made at the bottom of bottle assures a controlled gas flow from the soil. The results obtained for the activity of activated charcoal are in the range of 20–40 kBq/m³, for a depth of approximately 0.8 m. The results obtained by this method were confirmed by simultaneous measurements using LUK 3C device for soil radon measurements. Possible applications for the estimation of radon soil potential are discussed.     

Improving a regional model using reduced complexity and parameter estimation

The availability of powerful desktop computers and graphical user interfaces for ground water flow models makes possible the construction of ever more complex models. A proposed copper-zinc sulfide mine in northern Wisconsin offers a unique case in which the same hydrologic system has been modeled using a variety of techniques covering a wide range of sophistication and complexity. Early in the permitting process, simple numerical models were used to evaluate the necessary amount of water to be pumped from the mine, reductions in streamflow, and the drawdowns in the regional aquifer. More complex models have subsequently been used in an attempt to refine the predictions. Even after so much modeling effort, questions regarding the accuracy and reliability of the predictions remain. We have performed a new analysis of the proposed mine using the two-dimensional analytic element code GFLOW coupled with the nonlinear parameter estimation code UCODE. The new model is parsimonious, containing fewer than 10 parameters, and covers a region several times larger in areal extent than any of the previous models. The model demonstrates the suitability of analytic element codes for use with parameter estimation codes. The simplified model results are similar to the more complex models; predicted mine inflows and UCODE-derived 95% confidence intervals are consistent with the previous predictions. More important, the large areal extent of the model allowed us to examine hydrological features not included in the previous models, resulting in new insights about the effects that far-field boundary conditions can have on near-field model calibration and parameterization. In this case, the addition of surface water runoff into a lake in the headwaters of a stream while holding recharge constant moved a regional ground watershed divide and resulted in some of the added water being captured by the adjoining basin. Finally, a simple analytical solution was used to clarify the GFLOW model's prediction that, for a model that is properly calibrated for heads, regional drawdowns are relatively unaffected by the choice of aquifer properties, but that mine inflows are strongly affected. Paradoxically, by reducing model complexity, we have increased the understanding gained from the modeling effort.     

The uniform,regular differential equations of the KS transformed perturbed two-body problem

The Newtonian differential equations of motion for the two-body problem can be transformed into four, linear, harmonic oscillator equations by simultaneously applying the regularizing time transformation dt/ds=r and the Kustaanheimo-Stiefel (KS) coordinate transformation. The time transformation changes the independent variable from time to a new variables, and the KS transformation transforms the position and velocity vectors from Cartesian space into a four-dimensional space. This paper presents the derivation of uniform, regular equations for the perturbed twobody problem in the four-dimensional space. The variation of parameters technique is used to develop expressions for the derivatives of ten elements (which are constants in the unperturbed motion) for the general case that includes both perturbations which can arise from a potential and perturbations which cannot be derived from a potential. These element differential equations are slightly modified by introducing two additional elements for the time to further improve long term stability of numerical integration.Originally presented at the AAS/AIAA Astrodynamics Specialists Conference, Vail, Colorado, July 1973     

Extraterrestrial hydrogeology

Subsurface water processes are common for planetary bodies in the solar system and are highly probable for exoplanets (planets outside the solar system). For many solar system objects, the subsurface water exists as ice. For Earth and Mars, subsurface saturated zones have occurred throughout their planetary histories. Earth is mostly clement with the recharge of most groundwater reservoirs from ample precipitation during transient ice- and hot-house conditions, as recorded through the geologic and fossilized records. On the other hand, Mars is mostly in an ice-house stage, which is interrupted by endogenic-driven activity. This activity catastrophically drives short-lived hydrological cycling and associated climatic perturbations. Regional aquifers in the Martian highlands that developed during past, more Earth-like conditions delivered water to the northern plains. Water was also cycled to the South Polar Region during changes in climate induced by endogenic activity and/or by changes in Mars orbital parameters. Venus very likely had a warm hydrosphere for hundreds of millions of years, before the development of its current extremely hot atmosphere and surface. Subsequently, Venus lost its hydrosphere as solar luminosity increased and a run-away moist greenhouse took effect. Subsurface oceans of water or ammonia-water composition, induced by tidal forces and radiogenic heating, probably occur on the larger satellites Europa, Ganymede, Callisto, Titan, and Triton. Tidal forces operating between some of the small bodies of the outer solar system could also promote the fusion of ice and the stability of inner liquid-water oceans.

---

Resumen Los procesos hídricos subsuperficiales son comunes en cuerpos planetarios del sistema solar y son altamente probables para exoplanetas (planetas fuera del sistema solar). Para muchos cuerpos del sistema solar, el agua subsuperficial existe como hielo. Para la Tierra y Marte han ocurrido zonas saturadas subsuperficiales a través de sus historias planetarias. La Tierra es principalmente generosa con la recarga de la mayoría de reservorios de aguas subterráneas a partir de amplia precipitación reconocida en condiciones transitorias calientes y heladas, tal y como aparece en los registros fósiles y geológicos. Por otro lado, Marte se encuentra principalmente en una etapa de cámara de hielo la cual es interrumpida por actividad de tipo endogénico. Esta actividad pone en funcionamiento catastróficamente ciclos hidrológicos de vida corta y perturbaciones climáticas asociadas. Acuíferos regionales en las montañas de Marte que se desarrollaron en el pasado en condiciones similares a la Tierra distribuyen agua a las planicies del norte. El agua ha sido transportada hacia el sur de la región polar durante cambios en el clima inducidos por actividad endogénica y/o cambios en los parámetros orbitales de Marte. Venus muy probablemente tuvo una hidrósfera caliente durante cientos de millones de años, antes de que se desarrollara su atmósfera y superficie actual extremadamente caliente. Subsecuentemente, Venus perdió su hidrósfera a medida que la luminosidad solar aumentó y un efecto de invernadero húmedo escapatorio se llevó a cabo. Océanos subsuperficiales de composición agua o amoniaco-agua, inducidos por fuerzas de marea y calentamiento radiogénico, probablemente ocurren en los satélites más grandes como Europa, Ganimeda, Callisto, Titan y Triton. Las fuerzas de marea que operan entre los cuerpos pequeños del sistema solar externo podrían también promover la fusión de hielo y la estabilidad de líquido interno-aguas de los océanos.

---

Résumé Les processus de subsurface impliquant leau sont communs pour les corps planétaires du système solaire et sont très probables sur les exoplanètes (planètes en dehors du système solaire). Pour plusieurs objets du systèmes solaire, leau de subsurface est présente sous forme de glace. Pour la Terre et Mars, les zones saturées de subsurface apparaissent à travers toute leur histoire planétaire. La Terre est particulièrement clémente avec la recharge des réservoirs, avec de amples précipitations, des conditions glaciaires et de fortes chaleurs, comme latteste les enregistrements géologiques et paléontologiques. Dun autre côté, Mars se trouve dans une phase essentiellement glaciaire, qui est interrompue par des activités contraintes par les phénomènes endogéniques. Cette activité conduit de manière catastrophique à des cycles hydrologiques et à des perturbations climatiques brutaux. Les aquifères régionaux dans les haute terres martiennes qui se sont formés dans des conditions similaires aux conditions terrestres, alimentent les plaines du Nord. Leau a également été déplacée vers le Pôle Sud martien durant des changements marqués par une forte activité endogénique et une modification des paramètres de lorbite de Mars. Venus possèdait vraisemblablement une hydrosphère chaude durant des millions dannée, avant le développement de son atmosphère et sa surface particulièrement chaude. Par après Venus a perdit son hydrosphère alors que la luminosité solaire augmentait et quune humidité liée à un effet de serre sinstallait. Les océans de subsurface deau ou deau ammoniacale, induits par les forces de marée et le chauffage radiogénique, apparaissent probablement sur les satellites les plus importants (Europa, Ganymede, Callisto, Titan, Triton). Les forces de marée entre les petits corps externes du système solaire peuvent également occasionner la fusion de glace et la stabilité des océans internes deau liquide.

     

Gas-partitioning tracer test to quantify trapped gas during recharge

Dissolved helium and bromide tracers were used to evaluate trapped gas during an infiltration pond experiment. Dissolved helium preferentially partitioned into trapped gas bubbles, or other pore air, because of its low solubility in water. This produced observed helium retardation factors of as much as 12 relative to bromide. Numerical simulations of helium breakthrough with both equilibrium and kinetically limited advection/dispersion/retardation did not match observed helium concentrations. However, better fits were obtained by including a decay term representing the diffusive loss of helium through interconnected, gas-filled pores. Calculations indicate that 7% to more than 26% of the porosity beneath the pond was filled with gas. Measurements of laboratory hydraulic properties indicate that a 10% decrease in saturation would reduce the hydraulic conductivity by at least one order of magnitude in the well-sorted sandstone, but less in the overlying soils. This is consistent with in situ measurements during the experiment, which show steeper hydraulic gradients in sandstone than in soil. Intrinsic permeability of the soil doubled during the first six months of the experiment, likely caused by a combination of dissolution and thermal contraction of trapped gas. Managers of artificial recharge basins may consider minimizing the amount of trapped gas by using wet, rather than dry, tilling to optimize infiltration rates, particularly in well-sorted porous media in which reintroduced trapped gas may cause substantial reductions in permeability. Trapped gas may also inhibit the amount of focused infiltration that occurs naturally during ephemeral flood events along washes and playas.     

Canonical orbital elements in terms of an arbitrary independent variable

Within the framework of the Hamiltonian mechanics in the extended phase space, a set of canonical elements of the Delaunay type is developed in terms of an arbitrary independent angular variable. Application to the four classical anomalies-eccentric, true, elliptic, and mean-is presented. Particular attention is given to the generalized time equation and its conjugate energy equation.     

An electron temperature distribution derived from optical and radio measurements of HII regions

Radio measurements of the electron temperature ofHii regions are obtained from the ratio of the brightness temperature of a hydrogen recombination line to that of the adjacent continuum, while optical measurements are obtained from the ratio of [Oiii] forbidden-line intensities. The radio and optical measurements made under the assumption of an isothermalHii region,T _R andT _opt respectively, are combined to derive a temperature distribution for an entire nebula. A sphericalHii region in local thermodynamic equilibrium with constant density which is optically thin in both the line and the continuum is used as a model. Assuming linear temperature gradients withT _R=6000K andT _opt=10000K, it is found thatT=12000K (1–0.74r/R), wherer is the distance from the center andR is the radius of the nebula.