Vegetation competition model for water and light limitation. II: Spatial dynamics of groundwater and vegetation

In temperate climates groundwater can have a profound effect on vegetation, because it strongly influences the spatio-temporal distribution of soil moisture in the rootzone and therefore the occurrence of water and oxygen stress of vegetation. This article focuses on vegetation and groundwater dynamics along a hill slope by developing and evaluating a fully coupled hydrological-vegetation model for a temperate forest ecosystem. The vegetation model is described in part 1 of this series of two papers. To simulate the hydrology an extended version of the saturated-unsaturated hydrological model STARWARS has been used. The coupled model is used to investigate both the short and long-term dynamics for a system of two species. Both compete for light and water where one is adapted to wet conditions and the other to dry conditions. The daily dynamics show that the influence of groundwater is particularly strong in spring when waterlogging occurs due to decreased evapotranspiration in winter. Long simulation runs of 1000 years were performed to study the equilibrium state for the two species. Comparison of simulation results with observations of groundwater depth and vegetation types along a dry-wet gradient in a natural forest shows that a reductionist approach is able to capture these patterns well. Sensitivity analysis shows that the border between wet- and dry-adapted species moves upslope with increased rainfall, decreased slope angle and decreased aquifer thickness. These results are similar to previous findings which were based on global maximization of ecosystem evaporation or minimizing ecosystem stress. Comparison of runs with a fixed and a dynamic groundwater table shows that a dynamic groundwater table facilitates a wider transition zone between vegetation types along the hill slope. In this transition the biomass of vegetation is higher in the case of a dynamic groundwater than in case of a static groundwater table. This underlines the importance of incorporating spatial groundwater dynamics in models of groundwater influenced ecosystems.     

Mechanisms of moisture stress in a mid-latitude temperate forest: Implications for feedforward and feedback controls from an irrigation experiment

While it is well established that stomata close during moisture stress, strong correlations among environmental (e.g., vapor pressure deficit, soil moisture, air temperature, radiation) and internal (e.g., leaf water potential, sap flow, root-shoot signaling) variables obscure the identification of causal mechanisms from field experiments. Models of stomatal control fitted to field data therefore suffer from ambiguous parameter identification, with multiple acceptable (i.e., nearly optimal) model structures emphasizing different moisture status indicators and different processes. In an effort to minimize these correlations and improve parameter and process identification, we conducted an irrigation experiment on red maples (Acer rubrum L.) at Harvard Forest (summers of 2005 and 2006). Control and irrigated trees experienced similar radiative and boundary layer forcings, but different soil moisture status, and thus presumably different diurnal cycles of internal leaf water potential. Measured soil moisture and atmospheric forcing were used to drive a transient tree hydraulic model that incorporated a Jarvis-type leaf conductance in a Penman–Monteith framework with a Cowan-type (resistance and capacitance) tree hydraulic representation. The leaf conductance model included dependence on both leaf matric potential, Ψ_L (so-called feedback control) and on vapor pressure deficit, D (so-called feedforward control). Model parameters were estimated by minimizing the error between predicted and measured sap flow. The whole-tree irrigation treatment had the effect of elevating measured transpiration during summer dry-downs, demonstrating the limiting effect that subsurface resistance may have on transpiration during these times of moisture stress. From the best fitted model, we infer that during dry downs, moisture stress manifests itself in an increase of soil resistance with a resulting decrease in Ψ_L, leading to both feedforward and feedback controls in the control trees, but only feedforward control for the irrigated set. Increases in the sum-of-squares error when individual model components were disabled allow us to reject the following three null hypotheses: (1) the f(D) stress is statistically insignificant (p = 0.01); (2) the f(Ψ_L) stress is statistically insignificant (p = 0.07); and (3) plant storage capacitance is independent of moisture status (p = 0.07).     

Modeling seasonal course of carbon fluxes and evapotranspiration in response to low temperature and moisture in a boreal Scots pine ecosystem

Environmental conditions act above and below ground, and regulate carbon fluxes and evapotranspiration. The productivity of boreal forest ecosystems is strongly governed by low temperature and moisture conditions, but the understanding of various feedbacks between vegetation and environmental conditions is still unclear. In order to quantify the seasonal responses of vegetation to environmental factors, the seasonality of carbon and heat fluxes and the corresponding responses for temperature and moisture in air and soil were simulated by merging a process-based model (CoupModel) with detailed measurements representing various components of a forest ecosystem in Hyytiälä, southern Finland. The uncertainties in parameters, model assumptions, and measurements were identified by generalized likelihood uncertainty estimation (GLUE). Seasonal and diurnal courses of sensible and latent heat fluxes and net ecosystem exchange (NEE) of CO₂ were successfully simulated for two contrasting years. Moreover, systematic increases in efficiency of photosynthesis, water uptake, and decomposition occurred from spring to summer, demonstrating the strong coupling between processes. Evapotranspiration and NEE flux both showed a strong response to soil temperature conditions via different direct and indirect ecosystem mechanisms. The rate of photosynthesis was strongly correlated with the corresponding water uptake response and the light use efficiency. With the present data and model assumptions, it was not possible to precisely distinguish the various regulating ecosystem mechanisms. Our approach proved robust for modeling the seasonal course of carbon fluxes and evapotranspiration by combining different independent measurements. It will be highly interesting to continue using long-term series data and to make additional tests of optional stomatal conductance models in order to improve our understanding of the boreal forest ecosystem in response to climate variability and environmental conditions.     

Water budget and the consequent duration of canopy carbon gain in a teak plantation in a dry tropical region: Analysis using a soil–plant–air continuum multilayer model

A soil–plant–air continuum multilayer model was used to numerically simulate canopy net assimilation (A_n), evapotranspiration (ET), and soil moisture in a deciduous teak plantation in a dry tropical climate of northern Thailand to examine the influence of soil drought on A_n. The timings of leaf flush and the end of the canopy duration period (CDP) were also investigated from the perspective of the temporal positive carbon gain. Two numerical experiments with different seasonal patterns of leaf area index (LAI) were carried out using above-canopy hydrometeorological data as input data. The first experiment involved seasonally varying LAI estimated based on time-series of radiative transmittance through the canopy, and the second experiment applied an annually constant LAI. The first simulation captured the measured seasonal changes in soil surface moisture; the simulated transpiration agreed with seasonal changes in heat pulse velocity, corresponding to the water use of individual trees, and the simulated A_n became slightly negative. However, in the second simulation, A_n became negative in the dry season because the decline in stomatal conductance due to severe soil drought limited the assimilation, and the simultaneous increase in leaf temperature increased dark respiration. Thus, these experiments revealed that the leaflessness in the dry season is reasonable for carbon gain and emphasized the unfavorable soil water status for carbon gain in the dry season. Examining the duration of positive A_n (DPA) in the second simulation showed that the start of the longest DPA (LDPA) in a year approached the timing of leaf flush in the teak plantation after the spring equinox. On the other hand, the end appeared earlier than that of all CDPs. This result is consistent with the sap flow stopping earlier than the complete leaf fall, implying that the carbon assimilation period ends before the completion of defoliation. The model sensitivity analysis in the second simulation suggests that a smaller LAI and slower maximum rate of carboxylation likely extend the LDPA because soil water from the surface to rooting depth is maintained longer at levels adequate for carbon gain by decreased canopy transpiration. The experiments also suggest that lower soil hydraulic conductivity and deeper rooting depth can postpone the end of the LDPA by increasing soil water retention and the soil water capacity, respectively.     

Fluoride distribution in the environment along the gradient of a phosphate-fertilizer production emission (southern Brazil)

Airborne fluoride was determined in the rainwater, surface soil and groundwater along a gradient of emission of a phosphate fertilizer factory in Rio Grande, southern Brazil. Concentrations of fluoride in rainwater and groundwater achieved 3 mg l⁻¹ and 5 mg l⁻¹, respectively, and were dependent on pH. The fluoride deposited from emissions accumulated in a superficial horizon of soil in quantities comparable to those in the manufactured end-products—up to 23,000 mg kg⁻¹. Fluoride distribution in the environment is controlled by physical–chemical parameters of emission, rain intensity and soil properties. The highest fluoride concentrations were registered at a close distance of up to 2 km from the factory. The distribution of fluoride in groundwater resembled the same distribution in rainwater due to the high permeability of the local soils. Fluoride penetration to the groundwater also depended on the type of vegetation cover. The groundwater in woodland areas was less affected by contamination of fluoride than in the grassland areas, most probably because of the influence of eucalyptus throughfall, which increases the pH of wet precipitates.     

Raising groundwater differentially affects mineralization and plant species abundance in dune slacks.

The experience with restoring high water levels (i.e., rewetting) within restoration ecology is limited, and information on changes in soil nutrient supply is scarce. A reduction in nutrient supply is needed to restore the desired oligotrophic vegetation. We determined the effects of restoration of high water levels on decomposition and net carbon (C), nitrogen (N), and phosphorus (P) mineralization rates in wet dune slacks and its consequences for the relative abundance of eutrophic vs. oligotrophic species in the vegetation. This was done by analyzing these variables for valleys that experienced a large groundwater rise vs. valleys that had a small groundwater rise but the same current water level. In addition, the influences of underlying factors (waterlogging, vegetation dieback, and soil dynamics prior to groundwater rise) were separated in a transplantation experiment. Short-term effects of large groundwater rise were a massive dieback of vegetation, increased thickness of the fermentation layer, increased microbial decomposition activity, increased C mineralization, and decreased net N mineralization. Net P mineralization was not affected. The relative abundance of oligotrophic vs. eutrophic species was greater at large groundwater rise. Changes in decomposition and mineralization by large groundwater rise were, however, not caused by the vegetation dieback, but due to previous soil conditions. Soils experiencing waterlogged conditions for 3-4 years or more prior to large groundwater rise had lower C and higher net N mineralization rates at waterlogged conditions than soils that had experienced aerobic conditions, presumably due to differences in labile soil C contents. Contrary to expectations induced by previously determined nutrient pulses and measured vegetation dieback, large groundwater rise resulted in lower soil nutrient supply rates and more oligotrophic vegetation. If these trends continue on the longer term, restoration of high water levels may be effective in restoration ecology to establish oligotrophic, wet vegetation in dune slacks.     

Estimating the potential impact of vegetation on the water cycle requires accurate soil water parameter estimation

It is well known that vegetation dynamics at the catchment scale depends on the prevailing weather and soil moisture conditions. Soil moisture, however, is not equally distributed in space due to differences in topography, weather patterns, soil properties and the type and amount of vegetation cover. To elucidate the complex interaction between vegetation and soil moisture, the dynamic vegetation model LPJ-GUESS (Smith et al., 2001), which provides estimations of vegetation dynamics, but does not consider lateral water fluxes was coupled with the hydrological TOPMODEL (cf. Beven, 2001) in order to be able to evaluate the importance of these lateral fluxes. The new model LG-TM was calibrated and validated in two climatically different mountain catchments. The estimations of runoff were good, when monthly and weekly time scales were considered, although the low flow periods at winter time were somewhat underestimated. The uncertainty in the climate induced change vegetation carbon storage caused by the uncertainty in soil parameters was up to 3-5 kg C m⁻² (depending on elevation and catchment), compared to the total change in vegetation carbon storage of 5-9 kg C m⁻². Therefore accurate estimates of the parameters influencing the water holding capacity of the soil, for example depth and porosity, are necessary when estimating future changes in vegetation carbon storage. Similarly, changes in plant transpiration due to climatic changes could be almost double as high (88 mm m⁻²) in the not calibrated model compared to the new model version (ca 50 mm m⁻² transpiration change). The uncertainties in these soil properties were found to be more important than the lateral water exchange between grid cells, even in steep topography at least for the temporal and spatial resolution used here.     

盐渍土系统土壤水-地下水转化规律研究

通过野外定位观测和室内分析,探讨了苏打盐渍土微域特征(32m长的横截面)及其水分和地下水之间的转化规律。结果表明:土壤特征以及微地貌格局是制约土壤表层水分迁移的主要因素,并由此形成了大汽降水-土壤水-地下水转化的特殊模式。盐化草甸土分布在相对低洼的部位,苏打碱土分布在微坡地和高平地。苏打碱土质地为粉砂质亚粘土-粘土-粉砂-砂砾垂直结构,盐化草甸土以粉砂和砂砾结构为主。地表径流为水分迁移的关键环节,盐化草甸土成为地表水和地下水转化的通道。承接微坡地和高平地的径流,盐化草甸土在雨季含水率较多(甚至大到饱和)且迅速补给地下水,使地下水位逐渐上升。坡地土壤(苏打盐渍土)整个土层的含水率变化不显著。运用Vensim对盐渍土系统的土壤水-地下水的转化过程进行了模拟,结果与观测数据基本一致。在现代盐渍化过程研究中,以坡面系统为研究单元并考虑坡面径流,才能认识到盐渍化地区土壤水-地下水之间转化的实质。     

Modeling the vegetation-atmosphere carbon dioxide and water vapor interactions along a controlled CO2 gradient

Ecosystem functioning is intimately linked to its physical environment by complex two-way interactions. These two-way interactions arise because vegetation both responds to the external environment and actively regulates its micro-environment. By altering stomatal aperture, and therefore the transpiration rate, plants modify soil moisture and atmospheric humidity and these same physical variables, in return, modify stomatal conductance. Relationships between biotic and abiotic components are particularly strong in closed, managed environments such as greenhouses and growth chambers, which are used extensively to investigate ecosystem responses to climatic drivers. Model-assisted designs that account for the physiological dynamics governing two-way interactions between biotic and abiotic components are absent from many ecological studies. Here, a general model of the vegetation-atmosphere system in closed environments is proposed. The model accounts for the linked carbon-water physiology, the turbulent transport processes, and the energy and radiative transfer within the vegetation. Leaf gas exchange is modeled using a carbon gain optimization approach that is coupled to leaf energy balance. The turbulent transport within the canopy is modeled in two-dimensions using first-order closure principles. The model is applied to the Lysimeter CO₂ Gradient (LYCOG) facility, wherein a continuous gradient of atmospheric CO₂ is maintained on grassland assemblages using an elongated chamber where the micro-climate is regulated by variation in air flow rates. The model is employed to investigate how species composition, climatic conditions, and the imposed air flow rate affect the CO₂ concentration gradient within the LYCOG and the canopy micro-climate. The sensitivity of the model to key physiological and climatic parameters allows it to be used not only to manage current experiments, but also to formulate novel ecological hypotheses (e.g., by modeling climatic regimes not currently employed in LYCOG) and suggest alternative experimental designs and operational strategies for such facilities.     

Are savannas patch-dynamic systems? A landscape model

Savannas are ecosystems characterized by the coexistence of woody species (trees and bushes) and grasses. Given that savanna characteristics are mainly formed from competition, herbivory, fire, woodcutting, and patchy soil and precipitation characteristics, we propose a spatially explicit model to examine the effects of the above-mentioned parameters on savanna vegetation dynamics in space and time. Furthermore, we investigate the effects of the above-mentioned parameters on tree–bush–grass ratios, as well as the degrees of aggregation of tree–bush–grass biomass. We parameterized our model for an arid savanna with shallow soil depth as well as a mesic one with generally deeper and more variable soil depths. Our model was able to reproduce savanna vegetation characteristics for periods of time over 2000 years with daily updated time steps. According to our results, tree biomass was higher than bush biomass in the arid savanna but bush biomass exceeded tree and grass biomass in the simulated mesic savanna. Woody biomass increased in our simulations when the soil's porosity values were increased (mesic savanna), in combination with higher precipitation. Savanna vegetation varied from open savanna to woodland and back to open savanna again. Vegetation cycles varied over ∼300-year cycles in the arid and ∼220-year cycles in the mesic-simulated savanna. Autocorrelation values indicated that there are both temporal and spatial vegetation cycles. Our model indicated cycling savanna vegetation at the landscape scale, cycles in cells, and patchiness, i.e. patch dynamics.     

Coupling soil water and shoot dynamics in three grass species: A spatial stochastic model on water competition in Neotropical savanna

Savannas are ecosystems known for their high environmental and economic value. They cover at least 20% of the global land surface and, in some cases, can act as a boundary between tropical rainforest and deserts. Water is an important determinant of savanna ecosystems.In this paper, we present a theoretical stochastic model of root competition for water, which couples, soil water availability, phenology, and root and shoot architecture applied to three Neotropical savanna grasses. Soil moisture was simulated using a daily balance, as proposed by Rodriguez-Iturbe et al. [Rodriguez-Iturbe, I., Porporato, A., Ridolfi, L., Isham, V., Cox, D.R., 1999. Probabilistic modelling of water balance at a point: the role of climate, soil and vegetation. Proc. R. Soc. London, Ser. A 455, 3789–3805.]. To simulate rainfall stochasticity, we used daily precipitation data from the airport weather station in the State of Barinas, Venezuela, for the period 1991–2007. Competition among neighbouring plants took into account the spatial distribution of the individuals. As a final step, the model allowed us to calculate the shoot dynamic of the species as a function of soil water availability.Using these data, we compared the behaviour of isolated plants, pairs and trios, and we found below-ground competition to be a fundamental component of global (shoot + root) competition. Finally, our model suggests various circumstances that allow poor competitor plants to coexist in competition for water with more successful competitors. Apparently, this is not only due to transpiration rates, but also to differences in shoot emergence and shoot growth.     

徐淮平原农田防护林系统土壤水分动态

对小麦不同物候期农田防护林系统土壤含水量动态变化的测定结果表明:小麦不同物候期田间平均含水量的变化规律为幼苗期>返青期>拔节期>灌浆期,幼苗期和返青期杨树防护林对小麦的胁地影响不大,灌浆期是农田防护林与小麦争夺水分最为激烈的时期;幼龄防护林带能有效增加农田土壤含水量,对小麦生长有利,成熟林与农作物争夺水分的能力比幼龄林强,且对农作物的水分竞争主要集中在0.2倍树高范围内,造成此范围内土壤含水量明显不足。     

我国中亚热带人工林水热通量对干旱差异响应的模拟研究

季节性干旱现象在我国中亚热带地区时有发生,为了研究该区域大气-生态系统之间的相互作用关系及其碳水收支状况,2002年起在江西省千烟洲（26.7°N,115.1°E）人工林生态系统建立了通量观测塔。2003年7月该人工林生态系统遭遇了历史上少有的高温少雨天气,本研究应用基于生理生态学过程的EALCO（Ecological Assimilation of Land and Climate Observation）模型及2003和2004年通量观测数据对该生态系统的水热通量进行了模拟,同时分析了干旱胁迫对它们产生的影响。结果显示,模型能够很好的模拟该生态系统的能量通量的日变化,净辐射、显热和潜热通量模拟值与实测值相关系数的平方（R2）及标准差分别为0.99和8.05 W.m-2;0.81和41.02 W.m-2;0.90和31.49 W.m-2,模型可以解释87%的日蒸散量的变化。从模拟结果看,2003年7月下旬（发生较严重干旱胁迫）较2004年同期（干旱程度轻）相比,冠层及土壤水势下降约2倍,植物蒸腾的日变化形式改变,根系吸水滞后冠层蒸腾的时间缩短约半小时,冠层导度下降40%～60%。模拟与观测结果均表明,2003年7月下旬每天正午的波文比大都介于1～2.2,而2004年同期正午的波文比则介于0.2～0.6。EALCO模型通过Ball模型将植物碳水过程耦合在一起,从而可以很好的模拟植物的气孔行为,进而准确的模拟植物水热过程对干旱的响应。土壤水分匮乏对冠层导度的限制是2003年干旱期间冠层潜热通量模拟值下降的根本原因。     

Combined effect of fire and water scarcity on vegetation patterns in arid lands

In many arid zones around the word, the vegetation spontaneously forms regular patterns to optimize the use of the scarce water resources. The patterns act as early warning signal that fragile ecosystems may suddenly undergo irreversible shifts, thus, interpreting the structural shape of vegetation patterns is crucial to deciphering the ecosystem history and its expected further development. The sudden and irreversible shift of delicate ecosystems as a consequence of minor variation of the climatic forcing has been studied extensively in the past. The attitude of the ecosystem to recover after a catastrophic event, such as fire, did not receive as much attention so far. Here we modelled fire, as a sudden shift of the ecosystem state variables and functionality and evaluated post-fire scenarios under the hypothesis that two major feedbacks shaped the vegetation patterns: a positive feedback between preferential infiltration and plant growth, and a second feedback between infiltration and vegetation burning. A simple model solving a set of partial differential equations for soil moisture, plant biomass, surface water and dead biomass balance predicted significantly diverse post-fire vegetation patterns depending on the fire severity and on the degree of soil water repellency induced by the vegetation burning.     

SEIB–DGVM: A new Dynamic Global Vegetation Model using a spatially explicit individual-based approach

We report the development of a new spatially explicit individual-based Dynamic Global Vegetation Model (SEIB–DGVM), the first DGVM that can simulate the local interactions among individual trees within a spatially explicit virtual forest. In the model, a sample plot is placed at each grid box, and then the growth, competition, and decay of each individual tree within each plot is calculated by considering the environmental conditions for that tree as it relates to the trees that surround it. Based on these parameters only, the model simulated time lags between climate change and vegetation change. This time lags elongated when original biome was forest, because existing trees prevent newly establish trees from receiving enough sunlight and space to quickly replace the original vegetation. This time lags also elongated when horizontal heterogeneity of sunlight distribution was ignored, indicating the potential importance of horizontal heterogeneity for predicting transitional behavior of vegetation under changing climate. On a local scale, the model reproduced climate zone-specific patterns of succession, carbon dynamics, and water flux, although on a global scale, simulations were not always in agreement with observations. Because the SEIB–DGVM was formulated to the scale at which field biologists work, the measurements of relevant parameters and data comparisons are relatively straightforward, and the model should enable more robust modeling of terrestrial ecosystems.     

A mechanistic model of tree competition and facilitation for Mediterranean forests: Scaling from leaf physiology to stand dynamics

Mechanistic theories of plant competition developed to explain changes in community structure and dynamics along resource availability gradients have been mostly applied to temperate forests and grasslands where light and nutrients are the two main limiting resources. In contrast, the mechanisms underlying the structure and dynamics of water-limited plant communities have been little explored. Also previous mechanistic models rely either on complex simulators, which are difficult to interpret or on simple conceptual models, which ignore too many critical details. In this study, we develop a model of stand dynamics for light and water-limited forests of intermediate complexity and we provide an analytical framework for its analyses. The model is an individual-based simulator that describes the feedback between transpiration, stomatal function and soil water dynamics with asymmetrical competition for light and water. Trees allocate carbon to three main compartments: shoot, stem and roots. We use the model to explore general patterns that may emerge across levels of biological organization from the leaf to the stand. Model predictions are consistent with a number of features of Mediterranean forests structure and dynamics. At the plant-level the leaf-based tradeoff between carbon gain and water loss expresses as a tradeoff between mortality and growth. This tradeoff explains plant morphological changes in above-ground biomass and root to shoot allocation along a water availability gradient. At the community-level, tradeoffs among carbon acquisition and water loss govern the sign of plant interactions along the gradient. Coexistence among morphological types was not observed for the range of parameters and environmental conditions explored. Overall the model provides an unifying explanation for the observed changes in the sign of plant to plant interactions along environmental gradients as well as a process-based formulation that can be linked to empirical studies.     

Modelling and mapping topographic variations in forest soils at high resolution: A case study

This article examines the utility of a digitally derived cartographic depth-to-water (DTW) index to model and map variations in drainage, vegetation and soil type and select soil properties within a forested area (40 ha) of the Swan Hills, Alberta, Canada. This index was derived from a LiDAR (Light Detection and Ranging) derived digital elevation model (DEM), with at least 1 ground return per m². The resulting DTW pattern was set to be zero along all DEM-derived flow channels, each with a 4 ha flow-initiation threshold. Soil type (luvisol, gleysol, mesisol), drainage type (very poor to excessive), vegetation type (hydric to xeric) and forest floor depth were determined along hillslope transects. These determinations conformed more closely to the DEM-derived log₁₀(DTW) variations (R² > 60%) than to the corresponding variations of the widely used topographic wetness index (TWI) (R² < 25%). Setting log₁₀(DTW) thresholds to represent the wet to moist to dry transitions between vegetation, drainage and soil type enabled a high-resolution mapping of these types across the study area. Also determined were soil moisture content, coarse fragment and soil particle composition (sand, silt, clay), pH, total C, N, S, P, Ca, Mg, K, Fe, Al, Mn, Zn, and available Ca, Mg, K, P and NH₄, by soil layer type and depth. Most of these variables were also more correlated with log₁₀(DTW) than with TWI, with and without soil layer depth as an additional regression variable. These variables are, therefore, subject to topographic controls to at least some extent, and can be modelled and mapped accordingly, as illustrated for soil moisture, forest floor depth and pH across the study area, from ridge tops to depressions. Further examinations revealed that the DEM-produced DTW and TWI patterns complemented one another, with DTW delineating soils in relation to local water-table influences, and with TWI delineating where the water would flow and accumulate.     

Distribution,enrichment mechanism and risk assessment for fluoride in groundwater: a case study of Mihe-Weihe River Basin,China

● High fluorine is mainly HCO₃·Cl-Na and HCO₃-Na type. ● F⁻ decreases with the increase of depth to water table. ● High fluoride is mainly affected by fluorine-containing minerals and weak alkaline. ● Fluorine pollution is mainly in the north near Laizhou Bay (wet season > dry season). ● Groundwater samples have a high F⁻ health risk (children > adults). Due to the unclear distribution characteristics and causes of fluoride in groundwater of Mihe-Weihe River Basin (China), there is a higher risk for the future development and utilization of groundwater. Therefore, based on the systematic sampling and analysis, the distribution features and enrichment mechanism for fluoride in groundwater were studied by the graphic method, hydrogeochemical modeling, the proportionality factor between conventional ions and factor analysis. The results show that the fluorine content in groundwater is generally on the high side, with a large area of medium-fluorine water (0.5–1.0 mg/L), and high-fluorine water is chiefly in the interfluvial lowlands and alluvial-marine plain, which mainly contains HCO₃·Cl-Na- and HCO₃-Na-type water. The vertical zonation characteristics of the fluorine content decrease with increasing depth to the water table. The high flouride groundwater during the wet season is chiefly controlled by the weathering and dissolution of fluorine-containing minerals, as well as the influence of rock weathering, evaporation and concentration. The weak alkaline environment that is rich in sodium and poor in calcium during the dry season is the main reason for the enrichment of fluorine. Finally, an integrated assessment model is established using rough set theory and an improved matter element extension model, and the level of groundwater pollution caused by fluoride in the Mihe-Weihe River Basin during the wet and dry seasons in the Shandong Peninsula is defined to show the necessity for local management measures to reduce the potential risks caused by groundwater quality.     

Factors affecting vegetation establishment and development in a sand dune chronosequence at Newborough Warren,North Wales

Newborough Warren is a large calcareous west coast UK dune system, which has experienced rapid vegetation spread in the last 70 years. Information from two high resolution chronosequences for dry and wet dune habitats, 0–145 years, was used to answer the following questions: Does climate influence colonisation of vegetation on bare sand? What are the timescales and sequences of successional change in the vegetation? Analysis of aerial photographs showed that stabilisation of the dune system since 1945 has occurred in three main phases. The onset of stabilisation predated myxomatosis by 10 years; while stabilisation virtually halted during the period 1964–1978. Periods of rapid stabilisation were coincident with higher values of Talbot’s Mobility index (M)?>?0.3. Successional development was apparent in both dry and wet habitats. Fixed dune grassland started to replace earlier successional communities at around 40 years, and could persist to 145 years. Linear succession in dune slacks was less apparent, but a separation between communities typically regarded as ‘younger’ and ‘older’ occurred at around 40 years. Species richness in dry dune habitats increased with age to a maximum on soils around 60 years old, then declined again. Species richness was unrelated to age or soil development in wet dune slacks. The influence of climate suggests that conservation managers can only operate within the constraints imposed by natural climatic conditions. Vegetation growth and soil development are closely linked and maintaining some open areas is key to preventing soil development and over-stabilisation.     

土壤水分变化对外来入侵植物飞机草生长的影响

以外来入侵植物飞机草为研究对象,设计3个水分梯度(15mm,10mm,5.83mm),研究飞机草的叶绿素含量、生物量、光合生理特征等指标对土壤含水量变化的响应。结果表明,随着土壤水分含量的降低,飞机草的单位叶面积干质量、叶绿素含量、净光合速率、气孔导度、蒸腾速率、胞间CO2摩尔分数、光饱和点等都呈下降趋势。同时,根冠比和水分利用率随着水分胁迫的加重而有所提高,表明飞机草对土壤水分胁迫有一定的适应能力,但总体而言,土壤水分胁迫不利于飞机草的生长,限制了其向降雨量较少的地区入侵。