Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

A high resolution,multiproxy Late‐glacial record of climate change and intrasystem responses in northwest England

A lacustrine carbonate sequence from Hawes Water, Lancashire, UK, has been studied using stable isotopic, lithological, pollen and mineral magnetic analysis. The data reveal four abrupt climatic oscillations in the Late‐glacial Interstadial leading up to the onset of the Loch Lomond Stadial. The data also point to climatic warming relatively early within the stadial, ca. 12 500 GRIP yr, prior to the onset of the Holocene. The oxygen isotope record is taken as a signature of climate forcing against which the response of the lake‐system can be monitored. By adopting this approach it is revealed that the response of the biological system to the rapid climatic oscillations is non‐linear and primarily a function of the antecedent conditions. A significant end‐Devensian isotopic excursion (A) is matched by only minor changes in the cold‐adapted floras and faunas. During the warmer interstadial, the response of the biological ecosystem (events B–D) is clearly influenced by thresholds: major changes in the catchment vegetation associated with relatively minor oscillations in the isotopic signature. The stratigraphical patterns reveal significant lag effects between the onset of climate deterioration and resulting changes in vegetation. Copyright © 2002 John Wiley & Sons, Ltd.     

Partitioning of polycyclic aromatic hydrocarbons between water and particles compared to bioaccumulation in mussels: a harbour case

Water and particles from three sites located in Halifax Harbour, near and further away from combined municipal sewage effluents (CMSE) were analysed for parental and alkylated polycyclic aromatic hydrocarbons (par and alkPAH). Bioavailability of PAH was compared for inter-tidal mussels collected at the same sites in April and November. The PAH fingerprint determined over 9 sampling times covering a period of 19 months differed more between phases (water and particles) and seasons, than between sites. In the spring, more alkPAH associated with diesel and gasoline were detected in the soluble phase, along with more bioaccumulation of alkPAH in inter-tidal mussels. A broader number of parPAH were detected in mussels collected in the fall. The mean sum of dissolved alkPAH concentrations was higher in water at the site closer to raw CMSE than at the other two sites and particles of that site. However, lowest bioconcentration factors (BCF) were determined in mussels of this more contaminated site. Similar biota-particle accumulation factors (BPAF) were determined for parPAH in mussels from the three sites, all lower than the BCF of alkPAH. The study indicates that sewage treatment plants will reduce the amount of parPAH especially larger than fluoranthene and pyrene from being deposited in the harbour; that alk naphthalenes and fluorenes present in water will continue to be discharged; that the disturbance of sediments can make particle-bound PAH available to mussels.     

A model of an exploding,radiating star in general relativity

In a previous paper Adams, Cary and Cohen (1994) presented a model of a supernova. In that paper the equations of General Relativity describing the evolution of a spherically symmetric, radiating star were solved analytically. The evolution of the star was determined by the application of boundary conditions at the center and at the edge. Due to lmitations in the presupernova model, only the very slow inward motion of an unstable, degenerate core could be considered. The solution was also limited by the need to exclude a runaway term, one that increased exponentially with time. Without the exclusion of the runaway, the luminosity would have increased without bound and the mass would have become negative.This paper presents a completely analytic solution to the equations of General Relativity describing the evolution of a Type II supernova. Professor S.E. Woosley kindly gave us data on the physical variables of a 12M ₀ presupernova star. In our model the core collapses within 1 s, leaving a 1.3M ₀ remnant. Shortly afterward 10.6M ₀ is ejected to infinity, and 0.17M ₀ is radiated away in the form of neutrinos. The distance of the edge from the center increases proportionally to the two-thirds power of the time. The luminosity decreases proportionally to the inverse four-thirds power.Although the runaway solution was modified by the exploding rather than a static envelope, it must still be excluded by adjusting initial conditions. Its character is changed from an exponential to a very large power (55) of time. The removal of a degree of freedom by this exclusion leads to physically non-sensical results such as negative luminosity. The inclusion of a term describing motion of the mantle due to neutrino interactions provides the additional degree of freedom necessary for physically reasonable results.     

Periodic orbits of an electric charge in a magnetic dipole field

Periodic orbits in the Stormer problem are studied using the symmetry lines of the Poincaré map introduced by De Vogelaere. Many known facts are explained by mean of these lines. The dynamics of four special symmetry lines when the Stormer parameter ₁ changes is presented, and we obtain a clear global view of the structure of the simple periodic orbits and their bifurcations, including the asymmetrical ones. New asymmetrical multiple periodic orbits are obtained.     

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.