Flows through forest canopies in complex terrain

Recent progress on boundary layer flow within and above tall forest canopies in complex terrain is reviewed from the perspective of developing methods to interpret carbon dioxide fluxes from tower measurements in real terrain. Two examples of complex terrain are considered in detail: a forest edge, which exemplifies nonuniform forests, and hilly terrain, which can lead to drainage currents at night. Dynamical arguments show that, when boundary layer winds approach a forest edge, the mean wind adjusts on a length scale of approximately 3L(c), where L(c) is the canopy drag length scale, which depends inversely on the leaf area density of the forest. Over a further distance that also scales on L(c), turbulence in the flow adjusts, and the mixing and transport in the canopy approaches the homogeneous limit. Even low hills change the neutral flow within and above the forest canopy substantially. When the canopy is tall, pressure gradients drive flow up both the upwind and downwind slopes of the hill, leading to an ejection of air out of the top of the canopy just downwind of the crest. This flow at the crest can then advect scalar out of the top of the forest, leading to large variations in the flux of scalar across the hill. At night, when the air near the ground cools and becomes stably stratified, turbulence within the canopy can collapse, even when the flow above the canopy remains turbulent. This leads to a decoupling of the air motions within the canopy from those above. The air above the canopy can then continue to pass up and over the hill, as it does in the neutral case, but at the same time, air within the canopy drains down the hill slopes as drainage currents. These analyses will help us understand when flux towers are reliably measuring the net ecosystem exchange and suggest ways of correcting the flux tower data in more complex situations.     

Analysis of a two-layer cloud system with RAMS model and comparison to airborne observations

A three-dimensional numerical model (Regional Atmospheric Modeling System—-RAMS) was used to study the formation and evolution of water forms in a two-layer cloud structure observed during a field campaign over Brest (France). The model performance in regular operations, using conventional meteorological data as initial and lateral boundary conditions, was also examined. Remote sensing observations of the cloud system and in-situ aircraft data, selected during the campaign, were used to validate the model outputs. The model simulations showed that the lower cloud formation was characterized by high number concentration of pristine ice and snow, while the concentration of aggregates, graupel and hail were considerably lower. Hydrometeors in liquid phase appeared demonstrating high number concentration and water content on the top of this layer. The upper cloud layer consisted only of frozen water substances in lower amounts. The qualitative and quantitative comparison of the model-calculated meteorological and microphysical fields to the available observational data revealed that the model reflected fairly well the cloud structure (e.g., the spatio-temporal variability of the cloud parameters, the geometry of the cloud system). However, there were deviations as far as the model underestimating the ice water content (IWC) and number concentration (N_t) fields is concerned, especially at the atmospheric layer between 2.5 and 4 km of altitude. These deviations of the model simulated quantities from the measured ones may be attributed either to the performance of the model’s microphysics scheme, to instrument inaccuracies and to the local disturbance caused by the aircraft.     

Using nocturnal cold air drainage flow to monitor ecosystem processes in complex terrain.

This paper presents initial investigations of a new approach to monitor ecosystem processes in complex terrain on large scales. Metabolic processes in mountainous ecosystems are poorly represented in current ecosystem monitoring campaigns because the methods used for monitoring metabolism at the ecosystem scale (e.g., eddy covariance) require flat study sites. Our goal was to investigate the potential for using nocturnal down-valley winds (cold air drainage) for monitoring ecosystem processes in mountainous terrain from two perspectives: measurements of the isotopic composition of ecosystem-respired CO2 (delta13C(ER)) and estimates of fluxes of CO2 transported in the drainage flow. To test if this approach is plausible, we monitored the wind patterns, CO2 concentrations, and the carbon isotopic composition of the air as it exited the base of a young (approximately 40 yr-old) and an old (>450 yr-old) steeply sided Douglas-fir watershed. Nocturnal cold air drainage within these watersheds was strong, deep, and occurred on more than 80% of summer nights. The depth of cold air drainage rapidly increased to tower height or greater when the net radiation at the top of the tower approached zero. The carbon isotope composition of CO2 in the drainage system holds promise as an indicator of variation in basin-scale physiological processes. Although there was little vertical variation in CO2 concentration at any point in time, we found that the range of CO2 concentration over a single evening was sufficient to estimate delta 13C(ER) from Keeling plot analyses. The seasonal variation in delta 13C(ER) followed expected trends: during the summer dry season delta 13C(ER) became less negative (more enriched in 13C), but once rain returned in the fall, delta 13C(ER) decreased. However, we found no correlation between recent weather (e.g., vapor pressure deficit) and delta 13C(ER) either concurrently or with up to a one-week lag. Preliminary estimates suggest that the nocturnal CO2 flux advecting past the 28-m tower is a rather small fraction (<20%) of the watershed-scale respiration. This study demonstrates that monitoring the isotopic composition and CO2 concentration of cold air drainage at the base of a watershed provides a new tool for quantifying ecosystem metabolism in mountainous ecosystems on the basin scale.     

Influence of low and decreasing food levels on Daphnia-algal interactions: Numerical experiments with a new dynamic energy budget model

Based on numerical experiments with a new physiologically structured population model we demonstrate that predator physiology under low food and under starving conditions can have substantial implications for population dynamics in predator-prey interactions. We focused on Daphnia-algae interactions as model system and developed a new dynamic energy budget (DEB) model for individual daphnids. This model integrates the κ-rule approach common to net assimilation models into a net-production model, but uses a fixed allocation of net-productive energy in juveniles. The new DEB-model agrees well with the results of life history experiments with Daphnia. Compared to a pure κ-rule model the new allocation scheme leads to significant earlier maturation at low food levels and thus is in better agreement with the data. Incorporation of the new DEB-model into a physiologically structured population model using a box-car elevator technique revealed that the dynamics of Daphnia-algae interactions are highly sensitive to the assumptions on the energy allocation of juveniles under low food conditions. Additionally we show that also other energy allocation rules of our DEB-model concerning decreasing food levels and starving conditions at the individual level have strong implications for Daphnia-algae interactions at the population level. With increasing carrying capacity of algae a stable equilibrium with coexistence of Daphnia occurs and algae shifts to limit cycles. The amplitudes of the limit cycles increase with increasing percentage of sustainable weight loss. If a κ-rule energy allocation is applied to juveniles, the stable equilibrium occurs for a much narrower range of algal carrying capacities, the algal concentration at equilibrium is about 2 times larger, and the range of algae carrying capacities at which daphnids become extinct extends to higher carrying capacities than in the new DEB-model. Because predator-prey dynamics are very sensitive to predator physiology under low food and starving conditions, empirical constraints of predator physiology under these conditions are essential when comparing model results with observations in laboratory experiments or in the field.     

Simulation analysis of environmental risk accident and management of high-sulfur gas field development in complex terrain

Environmental risk of high sulfur gas field exploitation has become one of the hot spots of environmental management studies. Severe gas H₂S blowout accidents in recent years have shown that poor understanding and estimates of the poisonous gas movement could lead to dangerous evacuation delays. It is important to evaluate the real concentration of H₂S, especially in complex terrain. Traditional experiential models are not valid in the case of rough terrain, especially in low-lying areas where the gas accumulates. This study, using high sulfur content gas field of Sichuan “Pu Guang gas field” as study object and adopting objective diagnosis of wind field of land following coordinate three dimensions, applied Lagrangian Puff Model and breaking up technique of puffs to simulate the H₂S diffusion condition of blowout accidents produced in the high sulfur content gas field of complex terrain area. The results showed that the H₂S distribution did not occur mainly in low wind direction, and due to the obstruction of the mountain’s body, it accumulated in front of mountain on produced turn over, flowed around submitted jumping type distribution. The mountain waist near the hilltop and low hollow river valley site rapture points simulating contrast showed that the higher the rapture point, the better the diffusing condition of pollutant, the distribution of risk sensitive point decided piping rupture environmental risk size combining the H₂S diffusion result and residential area dispersing in the study area, synthetic judge located in the high rapture point environmental risk was smaller than the low hollow point, thus it was suggested to carryout laying of lining build of equal high line of higher terrain. According to simulation results, the environmental risk management measures aimed at putting down adverse effects were worked out.     

Intercomparison of Two Models,ETA and RAMS,with TRACT Field Campaign Data

In this work a model intercomparison between RAMS and ETA models is carried out, with the aim of evaluating the quality and accuracy of these mesoscale models in reproducing the time evolution of the meteorology in real complex terrain. This is of great importance not only for meteorological forecast but also for air quality assessment. Numerical simulations are performed to reproduce the mean variables' fields and to compare them with measurements collected during the field campaign TRACT. The domain covers the Rhine valley and surrounding mountainous region and we consider a time period of two days. Results from simulations are compared to observations relative to ground stations and radiosoundings. A qualitative analysis is joined to a quantitative estimation of some reference statistical indexes. Both RAMS and ETA models performances are satisfactory when compared to the measured data and also their relative agreement is good. The mean variable fields are reproduced with a satisfactory degree of reliability, even if the simulated profiles are not able to describe the largest fluctuations of the variables. At the surface stations, the best agreement between predictions and observations is obtained for the wind velocity, while the quality of the results is lower for temperature and humidity.     

Mass conservation above slopes in the Regional Atmospheric Modelling System (RAMS)

This paper examines the mass balance in calculations with the Regional Atmospheric Modelling System (RAMS). An error is pointed out that concerns the calculation of the surface fluxes on slopes. This error affects all the prognostic variables in RAMS when sloping terrain is involved. Here we explain how the error can be corrected. To study the impact of the error, we compared simulations with the uncorrected and corrected model. The model contains CO₂ transport, and online mass balance calculations were performed for this tracer. Without correction, effective surface CO₂-fluxes on mountain slopes were found to be enhanced under certain common conditions to several times the parameterized fluxes. Neglecting this error may cause substantial deviations in both forward and inverse model calculations. After the correction a very good closure of the mass balance is obtained. The correction also modifies the meteorological parameters, although the consequences were limited compared to the CO₂-fluxes.     

Mass conservation and atmospheric dynamics in the Regional Atmospheric Modeling System (RAMS)

This paper examines the spatio-temporal patterns of atmospheric carbon dioxide transport predicted by the Regional Atmospheric Modeling System (RAMS). Forty-eight hour simulations over northern New England incorporating a simple representation of the diurnal summertime surface carbon dioxide forcing arising from biological activity indicate that, in its native formulation, RAMS exhibits a significant degree of mass non-conservation. Domain-wide rates of non-physical mass gain and mass loss are as large as three percent per day which translates into approximately eleven parts per million per day for carbon dioxide — enough to rapidly dilute the signature of carbon dioxide fluxes arising from biological activity. Analysis shows that this is due to the approximation used by RAMS to compute the Exner function. Substitution of the exact, physically complete equation improves mass conservation by two orders of magnitude. In addition to greatly improving mass conservation, use of the complete Exner function equation has a substantial impact on the spatial pattern of carbon dioxide predicted by the model, yielding predictions differing from a conventional RAMS simulation by as much as forty parts per million. Such differences have important implications both for comparisons of modeled atmospheric carbon dioxide concentrations to observations and for carbon dioxide inversion studies, which use estimates of atmospheric transport of carbon dioxide in conjunction with measurements of atmospheric carbon dioxide concentrations to infer the spatio-temporal distribution of surface carbon dioxide fluxes. Furthermore, use of the complete Exner function equation affects the vertical velocity and water mixing ratio fields, causing significant changes in accumulated precipitation over the region.     

冀北山区森林群落草本多样性及其与地形关系研究

通过野外调查取样,应用Menhinick丰富度指数、Simpson指数、Shannon-Wiener指数和Pielou均匀度指数研究了冀北山区12种典型森林群落内的草本植物多样性规律,并采用典范对应分析（CCA）方法研究其分布与地形的关系。结果表明：不同的群落类型草本植物的组成和多样性指数不同,山杨（Populus davidiana）纯林的种数最多为34种,黑桦（Betula dahurica Pall）纯林的种数最少为19种,其Simpson指数、Shannon-Wiener指数最高,分别为0.927 3和2.879 6;在1 200～1 700 m之间,随海拔升高物种多样性指数、均匀度指数和丰富度指数呈正相关性,但没有线性上升的趋势关系;选取的地形因子对样点草本的影响程度大小为：坡向〉海拔〉坡度〉坡位〉坡形,坡向是样点草本植物空间差异的最主要地形制约因子,环境解释率为98.7%,证明排序可信。     

The performance of RAMS in representing the convective boundary layer structure in a very steep valley

Data from a comprehensive field study in the Riviera Valley of Southern Switzerland are used to investigate convective boundary layer structure in a steep valley and to evaluate wind and temperature fields, convective boundary layer height, and surface sensible heat fluxes as predicted by the mesoscale model RAMS. Current parameterizations of surface and boundary layer processes in RAMS, as well as in other mesoscale models, are based on scaling laws strictly valid only for flat topography and uniform land cover. Model evaluation is required to investigate whether this limits the applicability of RAMS in steep, inhomogeneous terrain. One clear-sky day with light synoptic winds is selected from the field study. Observed temperature structure across and along the valley is nearly homogeneous while wind structure is complex with a wind speed maximum on one side of the valley. Upvalley flows are not purely thermally driven and mechanical effects near the valley entrance also affect the wind structure. RAMS captured many of the observed boundary layer characteristics within the steep valley. The wind field, temperature structure, and convective boundary layer height in the valley are qualitatively simulated by RAMS, but the horizontal temperature structure across and along the valley is less homogeneous in the model than in the observations. The model reproduced the observed net radiation, except around sunset and sunrise when RAMS does not take into account the shadows cast by the surrounding topography. The observed sensible heat fluxes fall within the range of simulated values at grid points surrounding the measurement sites. Some of the scatter between observed and simulated turbulent sensible heat fluxes are due to sub-grid scale effects related to local topography.     

Evolution of a stream ecosystem in recently deglaciated terrain

Climate change and associated glacial recession create new stream habitat that leads to the assembly of new riverine communities through primary succession. However, there are still very few studies of the patterns and processes of community assembly during primary succession for stream ecosystems. We illustrate the rapidity with which biotic communities can colonize and establish in recently formed streams by examining Stonefly Creek in Glacier Bay, Alaska (USA), which began to emerge from a remnant glacial ice mass between 1976 and 1979. By 2002, 57 macroinvertebrate and 27 microcrustacea species had become established. Within 10 years of the stream's formation, pink salmon and Dolly Varden charr colonized, followed by other fish species, including juvenile red and silver salmon, Coast Range sculpin, and sticklebacks. Stable-isotope analyses indicate that marine-derived nitrogen from the decay of salmon carcasses was substantially assimilated within the aquatic food web by 2004. The findings from Stonefly Creek are compared with those from a long-term study of a similarly formed but older stream (12 km to the northeast) to examine possible similarities in macroinvertebrate community and biological trait composition between streams at similar stages of development. Macroinvertebrate community assembly appears to have been initially strongly deterministic owing to low water temperature associated with remnant ice masses. In contrast, microcrustacean community assembly appears to have been more stochastic. However, as stream age and water temperature increased, macroinvertebrate colonization was also more stochastic, and taxonomic similarity between Stonefly Creek and a stream at the same stage of development was <50%. However the most abundant taxa were similar, and functional diversity of the two communities was almost identical. Tolerance is suggested as the major mechanism of community assembly. The rapidity with which salmonids and invertebrate communities have become established across an entire watershed has implications for the conservation of biodiversity in freshwater habitats.     

Evaluation and application of a three-dimensional water quality model in a shallow lake with complex morphometry

Fundamental hydrodynamic and ecological processes of a lake or reservoir could be adequately depicted by one-dimensional (1D) numerical simulation models. Whereas, lakes with significant horizontal water quality and hydrodynamic gradients due to their complex morphometry, inflow or water level fluctuations require a three-dimensional (3D) hydrodynamics and ecological analyses to accurately simulate their temporal and spatial dynamics. In this study, we applied a 3D hydrodynamic model (ELCOM) coupled with an ecological model (CAEDYM) to simulate water quality parameters in three bays of the morphologically complex Lake Minnetonka. A considerable effort was made in setting up the model and a systematic parameterization approach was adopted to estimate the value of parameters based on their published values. Model calibration covered the entire length of the simulation periods from March 29 to October 20, 2000. Sensitivity analysis identified the top parameters with the largest contributions to the sensitivity of model results. The model was next verified with the same setup and parameter values for the period of April 25 to October 10, 2005 against field data. Spatial and temporal dynamics were well simulated and model output results of water temperature (T), dissolved oxygen (DO), total phosphorus (TP) and one group of algae (Cyanobacteria) represented as chlorophyll a (Chla) compared well with an extensive field data in the bays. The results show that the use of the model along with an accurate bathymetry, a systematic calibration and corroboration (verification) process will help to analyze the hydrodynamics and geochemical processes of the morphologically complex Lake Minnetonka. An example of an ecological application of the model for Lake Minnetonka is presented by examining the effect of spatial heterogeneity on coolwater fish habitat analysis in 3D and under a scenario where horizontal spatial heterogeneity was eliminated (1D). Both analyses captured seasonal fish habitat changes and the total seasonal averages differed moderately. However, the 1D analysis did not capture local and short duration variabilities and missed suitable fish habitat variations of as much as 20%. The experiment highlighted the need for a 3D analysis in depicting ecological hot spots such as unsuitable fish habitats in Lake Minnetonka.     

An isobole-based statistical model and test for synergism/antagonism in binary mixture toxicity experiments

Synergism and antagonism are often defined in relation to the model of Concentration Addition (CA). Hence, it is vital for the conclusion of mixture toxicity studies to be able to test whether an observed deviation from CA reflects a true deviation or whether it is simply due to random variation. In this paper we consider a non-linear regression model for the classical ray designs for binary mixture experiments. The model combines dose–response curves for each mixture in the experiment with an isobole model, describing possible deviations from CA. The method allows us to test whether the chosen isobole model is reasonable for the data and to test the hypothesis of CA. Furthermore, it provides us with a measure of the degree of synergism/antagonism. The method is flexible since both the dose–response relationships and the isobole model can be chosen arbitrarily. We demonstrate the use of the method on datasets where combinations of pesticides are tested on a floating plant, Lemna minor, and an algae, Pseudokirchneriella subcapitata. Furthermore, we conduct a simulation study in order to explore the power with which a specific deviation from CA can be distinguished in different test-systems.     

Linking species- and ecosystem-level impacts of climate change in lakes with a complex and a minimal model

To study the interaction between species- and ecosystem-level impacts of climate change, we focus on the question of how climate-induced shifts in key species affect the positive feedback loops that lock shallow lakes either in a transparent, macrophyte-dominated state or, alternatively, in a turbid, phytoplankton-dominated state. We hypothesize that climate warming will weaken the resilience of the macrophyte-dominated clear state. For the turbid state, we hypothesize that climate warming and climate-induced eutrophication will increase the dominance of cyanobacteria. Climate change will also affect shallow lakes through a changing hydrology and through climate change-induced eutrophication. We study these phenomena using two models, the full ecosystem model PCLake and a minimal dynamic model of lake phosphorus dynamics. Quantitative predictions with the complex model show that changes in nutrient loading, hydraulic loading and climate warming can all lead to shifts in ecosystem state. The minimal model helped in interpreting the non-linear behaviour of the complex model. The main output parameters of interest for water quality managers are the critical nutrient loading at which the system will switch from clear to turbid and the much lower critical nutrient loading – due to hysteresis – at which the system switches back. Another important output parameter is the chlorophyll-a level in the turbid state. For each of these three output parameters we performed a sensitivity analysis to further understand the dynamics of the complex model PCLake. This analysis showed that our model results are most sensitive to changes in temperature-dependence of cyanobacteria, planktivorous fish and zooplankton. We argue that by combining models at various levels of complexity and looking at multiple aspects of climate changes simultaneously we can develop an integrated view of the potential impact of climate change on freshwater ecosystems.     

漓江流域红壤侵蚀区植被演替与复合农林试验

漓江季节性的干旱和洪涝并存成为地方政府、有关专家和旅游行业密切关注的环境问题,漓江流域的环境建设尤其是上游森林植被恢复,成为生态治理重点内容。漓江上游红壤区植被类型主要包括针叶林、阔叶林、竹林、灌丛和草丛,植被的水分特征差异较大,植被的生态调节、水源涵养功能各异。在漓江上游红壤区开展4种模式的植被恢复和复合农林试验示范,3年的试验表明,乡土常绿阔叶树适应性强、有效保持水土,果林经济效益显著、但生态功能较差,木本药材林具有较好的生态功能和潜在的经济效益。通过调查和试验示范,提出亚热带红壤退化山地农林经营的植物种类和复合模式,对漓江流域水资源调控和水土流失治理提出建议和对策。     

Emergent properties modelled with the functional structural tree growth model ALMIS: Computer experiments on resource gain and use

The functional structural tree growth model ALMIS uses the individual based modelling approach and is implemented in the object-oriented programming language SIMULA. All features (state variables) and functions (processes) are specified locally, on the level of the single plant organs. Increasing numbers of “copies” (objects) of these elementary units, Internodes, Leaves, Meristems, Roots, and Root tips, form the growing tree. Various procedures (e.g. Photosynthesis, Nutrient_uptake, Transport, Storage, Mobilisation, Respiration, Growth) are employed to describe carbon and nutrient uptake, and matter fluxes between the different plant organs. Combining plant physiology and architecture, ALMIS allows studying the effect of single ecophysiological and structural processes on whole tree growth and in the tree–environment system. Some of these effects, driven by microclimate, self-shading, variable nutrient availability, variable transport dynamics, and branching patterns are exemplified. From the interactions at the organ and sub-organ levels new features emerge at higher levels of plant organisation. These so-called emergent properties are, for example, lifetime spectrum of single organs, space filling (“architecture”) and self-thinning of the crown. The most prominent emergent properties are the different growth forms of trees resulting from simulations under various conditions. Their causal interrelations are discussed in detail.     

237Pu experiments with the plaice Pleuronectes platessa

Some experiments have been made on the metabolism of ²³⁷Pu — a high specific activity, gamma-emitting isotope of plutonium — by the plaice Pleuronectes platessa L. Very little of the isotope was incorporated into fish after 2 months exposure to labelled water. The oral retention from a variety of labelled foods was also very poor, and apart from the gut itself no incorporation of the isotope could, in fact, be demonstrated. Injected ²³⁷Pu was eliminated very slowly, and some evidence was obtained for an inverse relationship between the rate of ²³⁷Pu loss and the rate of growth of the fish. Growing fish incorporated a relatively larger fraction of the ²³⁷Pu body burden into skeletal material, and this was attained at the expense of the isotopic content of the liver. Very little ²³⁷Pu was incorporated into muscle.