Analysis of radial movement of an unconfined leaky aquifer due to well pumping and injection

Radial movement of an unconfined leaky aquifer was studied with respect to hydraulic forces that are induced by well recharge and discharge. New analytic solutions in the velocity and displacement fields were found and applied to describe transient movement in an unconfined leaky aquifer. Linear momentum and mass balance of saturated porous sediments, the Darcy-Gersevanov law, and the analytic solution of hydraulic drawdown for unsteady flow within the unconfined leaky aquifer were introduced to find the new solutions. Analytic results indicate that the nonlinear relation between the initial hydraulic head (h₀) and the well function has an insignificant effect on the aquifer transient movement when the drawdown s<0.02h ₀. When the well function is simplified with different assumptions and pumping conditions, the new solutions correspondingly reduce to cases that are similar to the Hantush-Jacob, Muskat, and Theis transient movement of a confined leaky aquifer. It was found that large leakance is important in slowing radial movement and reducing aquifer deformation. Flow velocity in the aquifer is more responsive to leakance than to cumulative displacement within the aquifer. The zones and boundary with tensile stress can be located using the same approach applied to a confined aquifer for risk assessment of earth fissuring.

---

Résumé Le mouvement radial dans un aquifère semi-perméable libre a été étudié en considérant les forces hydrauliques induites par la recharge et le prélèvement au niveau d’un puits. De nouvelles solutions analytiques des champs de vitesse et de déplacement ont été trouvées et appliquées pour décrire le mouvement transitoire dans un aquifère semi-perméable libre. Pour trouver les nouvelles solutions, le moment linéaire et le bilan de masse dans les sédiments poreux saturés, la loi de Darcy-Gersevanov, et la solution analytique du rabattement hydraulique pour un écoulement transitoire dans un aquifère semi-perméable libre sont pris en compte. Les résultats analytiques montrent que la relation non-linéaire entre le potentiel hydraulique initial (h₀) et la fonction de puits a un effet insignifiant sur le mouvement transitoire de l’aquifère quand le rabattements s <0.02h₀. Si la fonction de puits est simplifiée en tenant compte de différentes hypothèses et conditions de pompage, les nouvelles solutions se réduisent en conséquence à des cas similaires au mouvement transitoire dans un aquifère semi-perméable captif selon Hantush-Jacob, Muskat, et Theis. Il s’est avéré qu’une drainance élevée est importante pour le ralentissement du mouvement radial et la réduction de la déformation de l’aquifère. La vitesse d’écoulement dans l’aquifère est plus sensible à la drainance qu’au déplacement cumulé dans l’aquifère. Les zones et la frontière présentant des contraintes de traction peuvent être localisées en utilisant la même approche appliquée à un aquifère captif pour l’évaluation des risques de fissuration du sol.

---

Resumen Se ha estudiado el movimiento radial en un acuífero no confinado con fuga en relación a las fuerzas hidráulicas que son inducidas por la descarga y recarga del pozo. Se encontraron nuevas soluciones analíticas en los campos de desplazamiento y velocidad, las cuales se aplicaron para describir el movimiento transitorio en un acuífero no confinado con fuga. Se han introducido cálculos de balance de masa y momentum lineal de los sedimentos porosos saturados, la Ley Darcy-Gersevanov, y la solución analítica de descenso hidráulico para flujo sin régimen permanente dentro del acuífero no confinado con goteo para encontrar las nuevas soluciones. Los resultados analíticos indican que la relación no lineal entre la presión hidráulica inicial (h₀) y la función del pozo tiene efecto insignificante en el movimiento transitorio del acuífero cuando el descenso s < 0.02h₀. Cuando se simplifica la función del pozo en base a diferentes supuestos y condiciones de bombeo se reducen las nuevas soluciones a casos que son similares a el movimiento transitorio Hantush-Jacob, Muskat, y Theis de un acuífero confinado con goteo. Se encontró que las fugas grandes son importantes en la disminución del movimiento radial y en la reducción de la deformación del acuífero. La velocidad de flujo en el acuífero responde más a las fugas que al desplazamiento acumulativo dentro del acuífero. Pueden localizarse zonas y límites con esfuerzos tensionales usando el mismo enfoque aplicado a acuíferos confinados para evaluación de riesgo de generación de fracturas en el terreno.

---

     

Spatial and Temporal Variations in <Emphasis Type="Italic">Pyrodinium bahamense</Emphasis> Cyst Concentrations in the Sediments of Bioluminescent Mangrove Lagoon,St. Croix,USVI

Mangrove Lagoon, located on the island of St. Croix, US Virgin Islands (USVI), is one of few actively bioluminescent lagoons in a location experiencing significant anthropogenic impacts. The bioluminescence is due to an abundance of the dinoflagellate Pyrodinium bahamense in the water column. We recovered surface sediments and sediment cores from Mangrove Lagoon to analyze the spatial distribution and temporal variability of P. bahamense cysts in this system. Surface sediment P. bahamense cyst concentrations ranged from 0 to 466 cysts g^?1 dry sediment, with higher abundances associated with elevated surface water nutrient concentrations and a mixed terrestrial–marine organic matter source regime. In combination with available bioassay data, we hypothesize that phytoplankton utilize nutrients rapidly and subsequent decay of organic matter makes nutrients available for dinoflagellates at the sediment–water interface in the eastern and northern quadrants of the lagoon. However, the nutrients are rapidly exhausted during counterclockwise lagoon circulation resulting in the decline of primary productivity and dinoflagellate abundance in the western quadrants. Downcore profiles suggest that P. bahamense blooms have been occurring for decades, declining in recent years. No cysts were present in sediments predating dredging activities of the 1960s that created Mangrove Lagoon. Recent reductions in cyst abundance may be the result of limited primary productivity caused by restricted water exchange with Salt River Bay due to shallowing of a sill at the mouth of the lagoon. This research highlights the need for more comprehensive geochemical and fossil analyses to better understand long-term ecological variability and inform conservation efforts of these unique habitats.     

In situ investigation into fracture and subsidence of overburden strata for solid backfill mining

Controlling the subsidence of overburden strata is the main objective for backfill mining technology. In order to understand the controlling effect of filling body, a detection scheme with three direction boreholes was designed in the tailgate of T₃290 mining face in Tangshan Coal Mine to monitor the fracture developing process of overburden strata. In addition, a few displacement and stress sensors used for monitoring the roof subsidence and filling body stress were arranged behind the longwall face to monitor and record the force and deformation of gangue. The field investigation shows that the solid backfill mining technology can effectively reduce the fracture of overburden strata, i.e., the fracture height and the damage degree. The subsidence of overburden strata is mainly caused by the subsidence of roof which is divided into four stages: roof subsidence stage of mining face, rapid subsidence stage of gob roof, relatively stable stage, and long-term rheological stage. The gauge plays an essential role in slowing down the subsidence in each stage of the overburden subsidence, especially the relatively stable stage.     

Zircon U–Pb and Hf isotopic study of gneissic rocks from the Chinese Altai: Progressive accretionary history in the early to middle Palaeozoic

Gneissic rocks in the Chinese Altai Mountains have been interpreted as either Paleozoic metasedimentary rocks or Precambrian basement. This study reports geochemical and geochronological data for banded paragneisses and associated gneissic granitoids collected along a NE–SW traverse in the northwestern Chinese Altai. Petrological and geochemical data suggest that the protoliths of the banded gneisses were possibly immature sediments with significant volcanic input and that the gneissic granitoids were derived from I-type granites formed in a subduction environment. Three types of morphological features can be recognized in zircons from the banded gneisses and are interpreted to correlate with different sources. Zircons from five samples of banded paragneiss cluster predominantly between 466 and 528 Ma, some give Neoproterozoic ages, and a few yield discordant Paleoproterozoic to Archean ages. Zircon Hf isotopic compositions indicate that both juvenile/mantle and crust materials were involved in the generation of the source rocks from which these zircons were derived. In contrast, zircons occur ubiquitously as elongated euhedral prismatic crystals in the four samples of the gneissic granitoids, and define single populations for each sample with mean ages between 380 and 453 Ma. The general absence of Precambrian inheritance and positive zircon ?Hf values for these granitoids suggest insignificant crustal contribution to the generation of the precursor magmas. Our data can be interpreted in terms of a progressive accretionary history in early to middle Palaeozoic times, and the Chinese Altai may possibly represent a magmatic arc built on a continental margin dominated by Neoproterozoic rocks.     

Tides in a body librating about a spin–orbit resonance: generalisation of the Darwin–Kaula theory

The Darwin-Kaula theory of bodily tides is intended for celestial bodies rotating without libration. We demonstrate that this theory, in its customary form, is inapplicable to a librating body. Specifically, in the presence of libration in longitude, the actual spectrum of Fourier tidal modes differs from the conventional spectrum rendered by the Darwin–Kaula theory for a nonlibrating celestial object. This necessitates derivation of formulae for the tidal torque and the tidal heating rate, that are applicable under libration. We derive the tidal spectrum for longitudinal forced libration with one and two main frequencies, generalisation to more main frequencies being straightforward. (By main frequencies we understand those emerging due to the triaxiality of the librating body.) Separately, we consider a case of free libration at one frequency (once again, generalisation to more frequencies being straightforward). We also calculate the tidal torque. This torque provides correction to the triaxiality-caused physical libration. Our theory is not self-consistent: we assume that the tidal torque is much smaller than the permanent-triaxiality-caused torque, so the additional libration due to tides is much weaker than the main libration due to the permanent triaxiality. Finally, we calculate the tidal dissipation rate in a body experiencing forced libration at the main mode, or free libration at one frequency, or superimposed forced and free librations.     

Assessment of the Pacific decadal oscillation’s contribution to the occurrence of local torrential rainfall in north China

On 21–22 July 2012, torrential rains hit North China, with the daily precipitation record at Beijing station reaching 160.6 mm; this event is named the Beijing 7–21 case. This paper assesses the likelihood of the occurrence of local torrential rains, such as the Beijing 7–21 case, from the perspective of climate variability. In particular, the influence of the Pacific Decadal Oscillation (PDO) is assessed. There were five extreme events, with daily precipitation records equal to or larger than 160.6 mm, at Beijing station during the period 1951–2012; all of these events happened during negative phases of the PDO. The present analysis indicates that precipitation events more extreme than the Beijing 7–21 case should happen more than once per decade during negative phases of the PDO, but only about once every four decades during positive PDO phases. The negative phase of the PDO is found to be associated with a much greater probability of daily records of southerly winds in North China during summer. Strong southerly summer monsoons are deemed favorable for increasing the occurrence of local extreme rainfall over North China.     

Crossing Filaments

Solar filaments show the position of large-scale polarity-inversion lines and are used for the reconstruction of large-scale solar magnetic field structure on the basis of Hα synoptic charts for the periods that magnetographic measurements are not available. Sometimes crossing filaments are seen in Hα filtergrams. We analyze daily Hα filtergrams from the archive of Big Bear Solar Observatory for the period of 1999 – 2003 to find crossing and interacting filaments. A number of examples are presented and filament patterns are compared with photospheric magnetic field distributions. We have found that all crossing filaments reveal quadrupolar magnetic configurations of the photospheric field and presume the presence of null points in the corona.     

Mapping the thickness of sediments in the Ljubljana Moor basin (Slovenia) using microtremors

The Ljubljana Moor basin is characterized by moderate bedrock topography and thicknesses of Quaternary lacustrine and fluvial sediments ranging from 0 to 200 m. More than 65 boreholes which reached the bedrock were drilled in the area, but their distribution in the basin is very uneven and some data from the boreholes uncertain. There are also no data on S-velocity distribution within the basin, but seismic refraction measurements pointed out a rather uniform increase of P-velocity with depth, great impedance contrast with the bedrock and relatively small lateral velocity variations. The microtremor horizontal-to-vertical spectral ratio (HVSR) method was therefore applied as a complementary tool to seismic refraction survey to map the thickness of sediments. First, microtremors were measured at the locations of boreholes which reached the bedrock and the resonance frequencies determined. The inverse power relationship between the resonance frequency and the thickness of sediments was then determined from 53 data pairs. The quality of the correlation is moderate due to possible heterogeneities in sediments and possible 3D effects in some minor areas, but the obtained parameters correspond well to the values obtained in six other European basins. Secondly, a 16 km-long discontinuous seismic refraction profile was measured across the whole basin, leaving uncovered some larger segments where active seismic measurements were not possible. Microtremors were then measured at 64 locations along the same profile, using 250 m point spacing, without leaving any gaps. The frequency–thickness relationship was used to invert resonance frequencies to depths. These were first validated using the results of the seismic refraction survey, which showed good agreement, and finally used for interpolation in the segments of missing refraction data to obtain a continuous depth profile of the bedrock. The study has shown that the microtremor method can be used as a complementary tool for mapping the thickness of unconsolidated sediments also in areas characterized by moderate bedrock topography. As the input data are always to some extent uncertain, it is important to have a sufficiently large number of borehole data to establish a frequency–thickness relationship, as well as some additional independent geophysical information for its validation.     

Comparison of Seismic Dispersion and Attenuation Models

The frequency-dependent attenuation of seismic waves causes decreased resolution of seismic images with depth, and the difference in transmission losses induces amplitude variations with offset. Transmission losses may occur due to friction or fluid movement, or may result from scattering in thin-layer. Whatever the physical mechanism, they can often be conveniently described using an empirical formulation wherein the elastic moduli and propagation velocity are complex functions of frequency.We have compiled and compared algebraically and numerically eight different models involving complex velocity: the Kolsky-Futterman model, the power-law model, Kjartansson's model, Müller's model, Azimi's second and third model, the Cole-Cole model, and the standard linear-solid model.For two different parameter sets, the attenuation and phase velocity are computed in the seismic frequency band, and the plane-wave propagation of a Ricker wavelet for the other models is compared with that for the Kolsky-Futterman model. The first parameter set consists of parameters for each of the models calculated from expressions given in the appendix. These expressions make the different models behave similarly to the KF model. The second parameter set consists of model parameters that are numerically adapted to the KF model.By selecting proper parameters, all models, except the standard linear-solid model, show behavior similar to that of the Kolsky-Futterman model. The SLS model behaves differently from the other models as the frequency goes to zero or infinity. Broadband measurement data is needed to select a specific model for a given seismic experiment.     

Correlation between molecular absorption spectral slope ratios and fluorescence humification indices in characterizing CDOM

Ultraviolet–visible absorption spectral slope ratios S_R (slope in 275–295 nm divided by slope in 350–400 nm) and humification index (HIX, integrated fluorescence emission in 435–480 nm divided by that in 435–480 and 300–345 nm) were compared when characterizing chromophoric dissolved organic matter (CDOM) in three humic acids and 44 whole water samples. HIX increased with increasing pH for humic acids, while their S_R showed much more complicated dependencies on pH. There was a negative correlation between S_R and HIX. S_R increased in the order terrestrial coal/peat < terrestrial soil/river < seawater, while HIX increased in the order seawater < terrestrial soil/river < terrestrial coal/peat. The comparative study in this work indicates that terrestrially derived CDOM has higher HIX and lower S_R than marine CDOM. Investigators may potentially use these two indices to compare qualitatively the character of CDOM in different sources (e.g., terrestrial vs. marine).