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.     

Diel Patterns in Chlorophyll <Emphasis Type="Italic">a</Emphasis> Specific Absorption Coefficient and Absorption Efficiency Factor of Picoplankton

Diel patterns in the chlorophyll a specific absorption coefficient of surface picoplankton, a*_pico (γ) (m²·[mg chlorophyll a]⁻¹), were studied at 7 stations under daily cycle of in situ light condition in the western subarctic Pacific and Japan Sea. All the data were normalized by dividing the anomaly with daily averaged a*_pico (γ). Opposite diel patterns were observed for the normalized a*_pico (443) and a*_pico (675) with maximum toward dawn or dusk and minimum toward midday at 4 stations under low-irradiance (LI) conditions and vice versa at 3 stations under high-irradiance (HI) conditions. The absorption efficiency factors at red absorption peak, Q _a (675), were determined by reconstruction with intracellular chlorophyll a concentration and cell diameter. The normalized Q _a (675) also showed diel pattern with maximum toward midday and minimum toward dawn or dusk under LI. The diel pattern in a*_pico (675) and Q _a (675) were primarily caused by changes in intracellular chlorophyll a concentration due to photoadaptation under LI. The diel pattern in a*_pico (443) was influenced by pigmentation, as recognized by blue to red ratio [a*_pico (443)/a*_pico (675)] under HI. This study proposed that the opposite diel pattern in a*_pico (γ) might occur for a wide range of algal species. The results presented here have important consequences for the interpretation of diel variations in optical properties observed in the open ocean.     

CO<Subscript>2</Subscript> Flux Measurements over the Sea Surface by Eddy Correlation and Aerodynamic Techniques

The measurements of the vertical transport of CO₂ were carried out over the Sea of Japan using the specially designed pier of Kyoto University on September 20 to 22, 2000. CO₂ fluxes were measured by the eddy correlation and aerodynamic techniques. Both techniques showed comparable CO₂ fluxes during sea breeze conditions: −0.001 to −0.08 mg m⁻²s⁻¹ with the mean of −0.05 mg m⁻²s⁻¹. This means that the measuring site satisfies the fetch requirement for meteorological observations under sea breeze conditions. Moreover, the eddy diffusivity coefficient used in the aerodynamic technique is found to be consistent with the coefficient used in the eddy correlation technique. The present result leads us to conclude that the aerodynamic technique may be applicable to underway CO₂ flux measurements over the ocean and may be used in place of the bulk technique. The important point is the need to maintain a measuring accuracy of CO₂ concentration difference of the order of 0.1 ppmv on the research vessels or the buoys.     

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.     

Seasonal and Spatial Variations of Total Mass Flux around Coral Reefs in the Southern Ryukyus,Japan

Total mass flux, size distribution of sediment particles and some chemical components such as total carbon (TC), total nitrogen (TN) and calcium carbonate (CaCO₃) were monitored monthly using a multi-cup sediment traps at seven coral reef sites (6 reef flat and 1 reef slope) of the Marine Protected Areas around Ishigaki, Kohama, Kuroshima and Iriomote Islands in the southern Ryukyus, Japan from September 2000 to September 2001. The size distribution of trapped sediments revealed mostly uni-modal fine sand to mud in the reef flat and gravelly to coarse sand in the reef slope. The total mass flux ranged between 0.54 to 872 gm⁻²d⁻¹, and showed a pronounced seasonality (high in summer-autumn and low in spring) at each site, which was consistent with the rainfall and typhoon regime. Exceptionally high values were observed on the reef slope (Iriomote) in February–March 2001 (1533 gm⁻²d⁻¹) owing to a large amount of bottom sediment re-suspension. On the reef flat (Todoroki South and North; Ishigaki), values obtained in July–August 2001 (872 gm⁻²d⁻¹) and August–September 2001 (800 gm^{− 2}d⁻¹) indicate the high terrestrial discharge from Todoroki River. Trapped sediment particles consist of CaCO₃ (1.2–27.1%) and a non-carbonate fraction (98.8–72.9%), which contains total carbon (4.9–26%), carbonate carbon (CO₂-C) (0.2–3.1%) and non-carbonate carbon (NC-C) (7.9–25.6%). Total nitrogen content was in the range 0.02–0.48%. TN is contained mainly in the carbonate fraction and NC-C may be contained in the non-carbonate fraction. The low TN/OC ratio of the trapped sediments suggests that they were mostly of terrestrial origin and that both fractions migrated. The high total mass flux derived from Todoroki River exceeded the threshold at which a lethal effect on coral community is caused. The results stress the importance of conducting seasonal studies of sedimentation over more than one year and at more than one location in south Japan coral reef ecosystems to gain an understanding of the processes controlling the total mass fluxes and their nutrients content, also to develop an awareness of how to prevent the damage of coral reef ecosystems and, if it does occur, to allow mitigation measures to be undertaken.     

Seismic structure and tectonics of the Shackleton Fracture Zone (Drake Passage,Scotia Sea)

The structural framework of the southern part of the Shackleton Fracture Zone has been investigated through the analysis of a 130-km-long multichannel seismic reflection profile acquired orthogonally to the fracture zone near 60° S. The Shackleton Fracture Zone is a 800-km-long, mostly rectilinear and pronounced bathymetric lineation joining the westernmost South Scotia Ridge to southern South America south of Cape Horn, separating the western Scotia Sea plate from the Antarctic plate. Conventional processing applied to the seismic data outlines the main structures of the Shackleton Fracture Zone, but only the use of enhanced techniques, such as accurate velocity analyses and pre-stack depth migration, provides a good definition of the acoustic basement and the architecture of the sedimentary sequences. In particular, a strong and mostly continuous reflector found at about 8.0 s two-way traveltime is very clear across the entire section and is interpreted as the Moho discontinuity. Data show a complex system of troughs developed along the eastern flank of the crustal ridge, containing tilted and rotated blocks, and the presence of a prominent listric normal fault developed within the oceanic crust. Positive flower structures developed within the oceanic basement indicate strike-slip tectonism and partial reactivation of pre-existing faults. Present-day tectonic activity is found mostly in correspondence to the relief, whereas fault-induced deformation is negligible across the entire trough system. This indicates that the E–W-directed stress regime present in the Drake Passage region is mainly dissipated along a narrow zone within the Shackleton Ridge axis. A reappraisal of all available magnetic anomaly identifications in the western Scotia Sea and in the former Phoenix plate, in conjunction with new magnetic profiles acquired to the east of the Shackleton Fracture Zone off the Tierra del Fuego continental margin, has allowed us to propose a simple reconstruction of Shackleton Fracture Zone development in the general context of the Drake Passage opening.