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.     

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.     

Identifying China’s leading world city: a network approach

This paper reports our research on China’s world cities. Formal network analysis of air passenger linkages for recent years among China’s most populous cities and among many of the world’s largest cities allows us to identify the country’s leading world city from among the leading Mainland candidates, Beijing, Shanghai and Guangzhou. We theorize our findings about China’s world cities in relation to both global forces (and China’s increasing entanglement with them) and the policies and actions of the national state. We examine the national and global urban network through a longitudinal, two-level analysis of airline passenger travel for four time points between about 1990 and 2005. We show that Beijing was China’s leading world city at the beginning of the time period, a status it lost nationally in as early as 1995, and then globally 10 years later. On the other hand Shanghai became China’s leading world city, and it acquired this status first nationally in 2000, and then globally in 2005. The changing status of the Chinese capital corresponds to the country’s increasing involvement with the capitalist world economy. Shanghai’s ascendance as the leading world city in China may indicate that global forces have come to play an increasingly important role relative to that of the developmental state.

基于网络药理学探究补阳还五汤治疗颈椎病的作用机制

目的：基于网络药理学与分子对接技术分析补阳还五汤治疗颈椎病的作用机制。方法：使用TCMSP、化学专业数据库获取补阳还五汤的活性成分,并对潜在靶点进行预测及规范。分别从PharmGKB、DisGeNET、OMIM、GeneCards数据库中得到颈椎病疾病靶点,利用韦恩图获取补阳还五汤与颈椎病的交集靶点。通过CytoScape软件构建中药-活性成分-疾病靶点和蛋白质-蛋白质相互作用（PPI）网络,获得核心有效成分与关键靶点,利用David数据库对潜在靶点进行富集分析。最后运用AutoDock Vina软件对补阳还五汤核心有效成分与关键靶点进行分子对接验证。结果：共获得补阳还五汤治疗颈椎病的有效活性成分97个,包括槲皮素、山柰酚、黄芩素、木犀草素等;交集靶点64个,关键靶点有白细胞介素-6（IL-6）、肿瘤坏死因子（TNF）、丝氨酸/苏氨酸蛋白激酶1（AKT1）、白细胞介素-1B（IL-1B）等;主要涉及肿瘤坏死因子（TNF）、白细胞介素-17（IL-17）、磷脂酰肌醇3-激酶（PI3K）/蛋白激酶B（Akt）等信号通路,分子对接结果显示核心有效成分与关键靶点之间结合紧密,为补阳还五汤治疗颈椎病提供相应条件。结论：该研究在总体上预测了补阳还五汤治疗颈椎病的活性成分、靶点和信号通路,作用途径广泛,为下一步的临床应用提供参考及思路。     

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.     

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.     

Action mechanism of geotextile-reinforced cushion under breakwater on soft ground

Various measures including material experiments, centrifuge modeling tests and FEM numerical analyses were performed to study systematically the action mechanism of the geotextile-reinforced cushion under breakwater on soft ground and the effects of the strata characterization and the reinforcement condition on the stability of the breakwater-ground system. In the aspect of controlling the deformation, the geotextile-reinforced cushion under breakwater constrains the lateral displacement of both the embankment and the ground. From the viewpoint of stress, the reinforcement suppressed the range of high stress level in the system. In general, the weaker the ground is and the greater the modulus of the geotextile is, the more effective the reinforcement is. The tensile force in the geotextile is greater in the range of the main part of the embankment.     

Cross-shelf distribution of copepods and the role of event-scale winds in a northern California upwelling zone

Cross-shelf distribution and abundance of copepod nauplii and copepodids were measured during three summer upwelling seasons (2000–2002) in a coastal upwelling zone off northern California. These 3 years varied considerably in the intensity of winds, abundance of chlorophyll, and water temperature. The cruises in 2000 were characterized by relaxation conditions, with generally high levels of chlorophyll and high water temperature. The cruises in 2001 and 2002 were dominated by strong and persistent upwelling events, leading to lower chlorophyll and water temperatures. The copepod assemblage was dominated by Oithona spp., Acartia spp. and Pseudocalanus spp., with Metridia pacifica (lucens), Microsetella rosea, Oncaea spp. and Tortanus discaudatus also common during all 3 years. The cross-shelf distribution of copepods was generally shifted offshore during upwelling and onshore during relaxation events, although some variability between species occurred. Abundance of all life stages generally exhibited a negative correlation with cross-shelf transport averaged over at least 1–4 days and lagged by 0–3 days, indicating lower abundances during and immediately after active upwelling. However, copepod nauplii seemed to respond positively to wind events lasting 1–5 days followed by a period of relaxation lasting 6 or 7 days. These rapid rates of change in abundance are probably too great to be due to in situ growth and reproduction alone; physical processes must also play a role. These results suggest a highly dynamic relationship between copepods and upwelling events off northern California, with species-specific responses to upwelling to be expected.     

哀牢山-红河断裂带中段应变分析

对哀牢山—红河断裂带中段和平—水塘剖面、墨江—元江剖面和其它地段的岩石应变及磁组构进行了分析，表明应变强度的校正磁各向异性度PJ从断裂带向西至三叠系明显降低，变形强度向西迅速减弱。磁化率椭球体主轴展布反映出剪切带内、外变形方式的改变。剪切带内，以水平走滑运动为主；向西则以水平缩短为主。在应变分析中，对断裂带内的S—C组构两组面理夹角、杏仁体和石榴石应变标志体进行了测量，结果显示，剪切带内的磁组构和岩石组构间关系较差。此外，还对这两条剖面中的三叠系进行了应变和磁组构分析，找出了它们与磁化率椭球体对应轴率间的相关性。