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 qualitative and quantitative evaluation of experimental model catchment evolution

Due to the geological time scales required for observation of catchment evolution, surrogates or analogues of field data are necessary to understand long‐term processes. To investigate long‐term catchment behaviour, two experimental model catchments that developed without rigid boundaries under controlled conditions are examined and a qualitative and quantitative analysis of their evolution is presented. Qualitatively, the experimental catchments have the visual appearance of field scale data. Observation demonstrates that changes in catchment shape and network form are conservative. Quantitative analysis suggests that the catchments reach an equilibrium form while a reduction in the channel network occurs. While the catchments are laboratory scale models, the results provide insights into field scale behaviour. Copyright © 2003 John Wiley & Sons, Ltd.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

Benthic photosynthesis in submerged Wadden Sea intertidal flats

In this study we compare benthic photosynthesis during inundation in coarse sand, fine sand, and mixed sediment (sand/mud) intertidal flats in the German Wadden Sea. In situ determinations of oxygen-, DIC- and nutrient fluxes in stirred benthic chamber incubations were combined with measurements of sedimentary chlorophyll, incident light intensity at the sediment surface and scalar irradiance within the sediment. During submergence, microphytobenthos was light limited at all study sites as indicated by rapid response of gross photosynthesis to increasing incident light at the sea floor. However, depth integrated scalar irradiance was 2 to 3 times higher in the sands than in the mud. Consequently, gross photosynthesis in the net autotrophic fine sand and coarse sand flats during inundation was on average 4 and 11 times higher than in the net heterotrophic mud flat, despite higher total chlorophyll concentration in mud. Benthic photosynthesis may be enhanced in intertidal sands during inundation due to: (1) higher light availability to the microphytobenthos in the sands compared to muds, (2) more efficient transport of photosynthesis-limiting solutes to the microalgae with pore water flows in the permeable sands, and (3) more active metabolic state and different life strategies of microphytobenthos inhabiting sands.     

基于圆特征的地图要素自动综合算法研究

分析了目前地图综合的地位、作用和发展情况,介绍了“优胜劣汰”思想,分析对比出自然界“优胜劣汰”基本思想和地图制图综合在某些方面存在共同的相似之处,从而尝试把它运用到地图制图综合中来,并结合Circle技术及其基本特点,提出了基于Circle原理和“优胜劣汰”思想的地图综合新算法,在地图综合的选取、化简、合并等方面得到了较好的应用,相应的例子在文中详细地进行了列举。     

The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

Effects of crab–halophytic plant interactions on creek growth in a S.W. Atlantic salt marsh: A Cellular Automata model

The Bahía Blanca Estuary (38° 50′ S, and 62° 30′ W) presents salt marshes where interactions between the local flora (Sarcocornia perennis) and fauna (Chasmagnathus granulatus) generate some kind of salt pans that alter the normal water circulation and condition its flow and course towards tidal creeks. The crab–vegetation dynamics in the salt marsh presents variations that cannot be quantified in a reasonable period of time. The interaction between S. perennis plant and C. granulatus crab is based on simple laws, but its result is a complex biological mechanism that causes an erosive process on the salt marsh and favors the formation of tidal creeks. To study it, a Cellular Automata model is proposed, based on the laws deduced from the observation of these phenomena in the field, and then verified with measurable data within macroscale time units. Therefore, the objective of this article is to model how the interaction between C. granulatus and S. perennis modifies the landscape of the salt marsh and influences the path of tidal creeks. The model copies the basic laws that rule the problem based on purely biological factors.The Cellular Automata model proved capable of reproducing the effects of the interaction between plants and crabs in the salt marsh. A study of the water drainage of the basins showed that this interaction does indeed modify the development of tidal creeks. Model dynamics would likewise follow different laws, which would provide a different formula for the probability of patch dilation. The patch shape can be obtained changing the pattern that dilates.     

Freshening of the Intermediate Waters in the Northern South China Sea over the PastSix Decades

The properties of salinity in the South China Sea (SCS), a significant marginal sea connecting the Pacific andIndian Oceans, are greatly influenced by the transport of fresh water flux between the two oceans. However, the long-termchanges in the intermediate water in the SCS have not been thoroughly studied due to limited data, particularly in relationto its thermodynamic variations. This study utilized reanalysis data products to identify a 60-year trend of freshening in theintermediate waters of the northern South China Sea (NSCS), accompanied by an expansion of low-salinity water. Thestudy also constructed salinity budget terms, including advection and entrainment processes, and conducted an analysis ofthe salinity budget to understand the impacts of external and internal dynamic processes on the freshening trend of theintermediate water in the NSCS. The analysis revealed that the freshening in the northwest Pacific Ocean and theintensification of intrusion through the Luzon Strait at intermediate levels are the primary drivers of the salinity changes inthe NSCS. Additionally, a weakened trend in the intensity of vertical entrainment also contributes to the freshening in theNSCS. This study offers new insights into the understanding of regional deep sea changes in response to variations in boththermodynamics and oceanic dynamic processes.     

基于不同GCP数量和获取方法的数字正射影像的制作与评定

SPOT-5影像的高分辨率和高信息量的性质,使其得到了广泛的应用,特别是在全国第二次土地利用调查中得到了更加深入的应用。而数字正射影像的制作是SPOT-5遥感影像应用的关键所在。本文在介绍DOM制作的基础上,依据当前获取GCP主要的三种途径采用ERDAS遥感图像处理软件进行了实验,并对最终的结果精度进行分析,总结了DOM制作过程的相关问题。