Lattice Boltzmann Models for Flow and Transport in Saturated Karst

Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions.     

Functioning of intertidal flats inferred from temporal and spatial dynamics of O<Subscript>2</Subscript>, H<Subscript>2</Subscript>S and pH in their surface sediment

In this article, we describe the dynamics of pH, O₂ and H₂S in the top 5–10 cm of an intertidal flat consisting of permeable sand. These dynamics were measured at the low water line and higher up the flat and during several seasons. Together with pore water nutrient data, the dynamics confirm that two types of transport act as driving forces for the cycling of elements (Billerbeck et al. 2006b): Fast surface dynamics of pore water chemistry occur only during inundation. Thus, they must be driven by hydraulics (tidal and wave action) and are highly dependent on weather conditions. This was demonstrated clearly by quick variation in oxygen penetration depth: Seeps are active at low tide only, indicating that the pore water flow in them is driven by a pressure head developing at low tide. The seeps are fed by slow transport of pore water over long distances in the deeper sediment. In the seeps, high concentrations of degradation products such as nutrients and sulphide were found, showing them to be the outlets of deep-seated degradation processes. The degradation products appear toxic for bioturbating/bioirrigating organisms, as a consequence of which, these were absent in the wider seep areas. These two mechanisms driving advection determine oxygen dynamics in these flats, whereas bioirrigation plays a minor role. The deep circulation causes a characteristic distribution of strongly reduced pore water near the low water line and rather more oxidised sediments in the centre of the flats. The two combined transport phenomena determine the fluxes of solutes and gases from the sediment to the surface water and in this way create specific niches for various types of microorganisms.     

On Constraining Pilot Point Calibration with Regularization in PEST

Ground water model calibration has made great advances in recent years with practical tools such as PEST being instrumental for making the latest techniques available to practitioners. As models and calibration tools get more sophisticated, however, the power of these tools can be misapplied, resulting in poor parameter estimates and/or nonoptimally calibrated models that do not suit their intended purpose. Here, we focus on an increasingly common technique for calibrating highly parameterized numerical models—pilot point parameterization with Tikhonov regularization. Pilot points are a popular method for spatially parameterizing complex hydrogeologic systems; however, additional flexibility offered by pilot points can become problematic if not constrained by Tikhonov regularization. The objective of this work is to explain and illustrate the specific roles played by control variables in the PEST software for Tikhonov regularization applied to pilot points. A recent study encountered difficulties implementing this approach, but through examination of that analysis, insight into underlying sources of potential misapplication can be gained and some guidelines for overcoming them developed.     

A heat perturbation flow meter for application in soft sediments

A prototype flow meter has been developed, based upon the heat perturbation principle, to monitor groundwater specific discharge in soft sediments. The device is designed for use in spatially intensive, long-term monitoring campaigns in remote or inconvenient locations, and is cheap, robust and capable of being logged automatically. The results of the laboratory tests indicate that the heat perturbation principle is suitable for determining the magnitude of specific discharge to a degree of accuracy that would be useful in practical applications in dynamic groundwater systems with rapidly changing flows of approximately 1 md⁻¹ or more and that the groundwater flow direction can generally be determined to a high level of precision. The accuracy and reliability of the estimates of specific discharge have been shown to depend strongly upon the geometrical precision of manufacture and the quality of the temperature monitoring system. These factors become most significant in the estimation of lower flows and further investigation is required to determine the detection limit of the device. Specific discharge estimates have been shown to be insensitive to dispersivity values appropriate to the scale of the device. Unlike the majority of heat perturbation devices, calibration is unnecessary.     

Interactions between freshwater input, light, and phytoplankton dynamics on the Louisiana continental shelf

We examined the effects of freshwater flow and light availability on phytoplankton biomass and production along the Louisiana continental shelf in the region characterized by persistent spring–summer stratification and widespread summer hypoxia. Data were collected on 7 cruises from 2005 to 2007, and spatially-averaged estimates of phytoplankton and light variables were calculated for the study area using Voronoi polygon normalization. Shelf-wide phytoplankton production ranged from 0.47 to 1.75 mg C m⁻² d⁻¹ across the 7 cruises. Shelf-wide average light attenuation (k_d) ranged from 0.19–1.01 m⁻¹ and strongly covaried with freshwater discharge from the Mississippi and Atchafalaya Rivers (R²=0.67). Interestingly, we observed that the euphotic zone (as defined by the 1% light depth) extended well below the pycnocline and to the bottom across much of the shelf. Shelf-wide average chlorophyll a (chl a) concentrations ranged from 1.4 to 5.9 mg m⁻³ and, similar to k_d, covaried with river discharge (R²=0.83). Also, chl a concentrations were significantly higher in plume versus non-plume regions of the shelf. When integrated through the water-column, shelf-wide average chl a ranged from 26.3 to 47.6 mg m⁻², but did not covary with river discharge, nor were plume versus non-plume averages statistically different. The high integrated chl a in the non-plume waters resulted from frequent sub-pycnocline chl a maxima. Phytoplankton production rates were highest in the vicinity of the Mississippi River bird's foot delta, but as with integrated chl a were not statistically different in plume versus non-plume waters across the rest of the shelf. Based on the vertical distribution of light and chl a, a substantial fraction of phytoplankton production occurred below the pycnocline, averaging from 25% to 50% among cruises. These results suggest that freshwater and nutrient inputs regulate shelf-wide k_d and, consequently, the vertical distribution of primary production. The substantial below-pycnocline primary production we observed has not been previously quantified for this region, but has important implications about the formation and persistence of hypoxia on the Louisiana continental shelf.     

Merging eruption datasets: building an integrated Holocene eruptive record for Mt Taranaki,New Zealand

Acquiring detailed eruption frequency datasets for a volcano system is essential for realistic eruption forecasts. However, accurate datasets are inherently difficult to compile, even if one or more well-dated eruption records are available. A single record typically under-represents the eruption frequency, while combining two or more records may result in an overrepresentation. Although glass compositions have proven to be successful in tephrochronological studies of dominantly rhyolitic tephras; microlitic growth and thin glass shards inhibit their application to andesitic tephras. A method consisting of a combination of two techniques for correlating syn-eruptive deposits is demonstrated on data from the typical andesitic stratovolcano of Mt. Taranaki, New Zealand. Firstly, tentative matches are identified using the radiocarbon age and associated error of each event. Secondly, the compositions of titanomagnetite micro-phenocrysts are used as an independent check, and shown to be a useful correlation tool where age data is available. Using two lake-core records containing tephra layers in an overlapping time-frame, the radiocarbon age-correlation procedure suggested 31 tephra matches. Geochemistry data were available for 15 of these pairs. In three of these cases, the titanomagnetite compositions did not match. Hence, these “paired” tephras were from compositionally distinct magmas and therefore likely represent separate events. An additional three matches were reassigned within the temporal uncertainty limits of the dating procedure, based on better geochemical pairing. The final combined dataset suggests that there have been at least 138 separate ash fall-producing eruptions between 96 and 10 150 years B.P. from Taranaki. Using the combined dataset the mixture of Weibulls renewal model forecasts a probability of 0.52 for an eruption occurring in the next 50 years at this volcano. The present annual eruption probability is estimated at 1.6%. This likelihood is almost double that obtained when relying on a single stratigraphic record.     

Preliminary attempt to distinguish the domesticated pigs from wild boars by the methods of carbon and nitrogen stable isotope analysis

Despite great achievements in the origins of domestic pigs made by the methods of zooarchaeology and molecular biology, how to scientifically distinguish the domesticated pigs from wild boars during the early stage of pig domestication is still poorly understood. Compared to wild boar’s diets which come from the natural environment, the diets of domestic pigs are more easily influenced by human feeding activities. Therefore, in principle, exploration of the dietary differences among pigs and under-standing the impact on pig diets fed by humans can have great potential to differentiate between wild boars and domesticated pigs. To reveal dietary differences among pigs and distinguish the domesticated pigs from wild boars based on comparison with the diets of humans and other animals, we analyzed the carbon and nitrogen stable isotopes of human bones from Xiaojingshan Site and animal bones from Yuezhuang Site, both of which belong to Houli Culture in Shandong Province and date to about 8500-7500 years ago. The mean δ¹³C value ((−17.8 ± 0.3)‰) and δ¹⁵N value ((9.0±0.6)‰) in human collagen indicate that although millet agriculture began it was not the main subsistence strategy as millets are typical of C₄ plants and that humans made a living mainly by gathering, hunting or raising some domesticated animals. The δ¹³C value (−16.1‰) and δ¹⁵N value (6.9‰) in the bovine suggest that C₃ plants were dominant in its diet with some C₄ plants complemented. The fish has lower δ¹³C value (−24.9‰) and higher δ¹⁵N value (8.8‰) than the bovine, which is the characteristic of the isotopic values from Eurasian freshwater fish. Based on the differences in carbon and nitrogen isotope values, the pigs can be divided into three groups. A group, composed of two pigs, has low δ¹³C values (−18.1‰, −20.0‰) and low δ¹⁵N values (4.7‰, 6.0‰). B group, only one pig, has the highest δ¹³C value (−10.6‰) and mediate δ¹⁵N value (6.4‰). As for the C group, also only one pig, low δ¹³C value (−19.0‰) and the highest δ¹⁵N value (9.1‰) are observed. Previous studies on the stable isotopes from modern or ancient wild boars’ bones have suggested that C₃ plants are predominated in their diets and that their δ¹⁵N values are close to those in herbivores and far from those in carnivores. Based on the comparison with the isotope values from humans, the wild boars and the domestic pigs from Xipo Site in Henan 6000-5500 years ago and Kangjia Site in Shaanxi 4500-4000 years ago, we conclude that A pig group belongs to wild boars while B and C groups can be attributed to domesticated pigs. Supported by Max-Planck Society and Chinese Academy of Sciences, Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KJCX3.SYW.N12), National Natural Science Foundation of China (Grant No. 40702003) and President Funding of Graduate University of Chinese Academy of Sciences     

Reconstruction of paleocoastlines for the northwestern South China Sea since the Last Glacial Maximum

The range of relative sea level rise in the northwestern South China Sea since the Last Glacial Maximum was over 100 m. As a result, lowland regions including the Northeast Vietnam coast, Beibu Gulf, and South China coast experienced an evolution from land to sea. Based on the principle of reconstructing paleogeography and using recent digital elevation model, relative sea level curves, and sediment accumulation data, this paper presents a series of paleogeographic scenarios back to 20 cal. ka BP for the northwestern South China Sea. The scenarios demonstrate the entire process of coastline changes for the area of interest. During the late glacial period from 20 to 15 cal. ka BP, coastline slowly retreated, causing a land loss of only 1×10⁴ km², and thus the land-sea distribution remained nearly unchanged. Later in 15–10 cal. ka BP coastline rapidly retreated and area of land loss was up to 24×10⁴km², causing lowlands around Northeast Vietnam and South China soon to be underwater. Coastline retreat continued quite rapidly during the early Holocene. From 10 to 6 cal. ka BP land area had decreased by 9×10⁴km², and during that process the Qiongzhou Strait completely opened up. Since the mid Holocene, main controls on coastline change are from vertical crustal movements and sedimentation. Transgression was surpassed by regression, resulting in a land accretion of about 10×10⁴km². Supported by Key Laboratory of Marginal Sea Geology, Chinese Academy of Sciences (Grant No. MSGL0711), the Guangdong Natural Science Foundation (Grant No. 04001309) and Open Fund of the Key Laboratory of Marine Geology and Environment, Chinese Academy of Sciences (Grant No. MGE2007KG04)