Design of a Domain Specific Language for modelling processes in landscapes

The modelling of processes that occur in landscapes is often confronted to issues related to the representation of space and the difficulty of properly handling time and multiple scales. In order to investigate these issues, a flexible modelling environment is required. We propose to develop such a tool based on a Domain Specific Language (DSL) that capitalises on the service-oriented architecture (SOA) paradigm. The modelling framework around the DSL is composed of a model building environment, a code generator and compiler, and a program execution platform. The DSL introduces five language elements (entity, service, relation, scenario and datafacer) that can be combined to offer a wide range of possibilities for modelling in space and time at different scales. When developing a model, model parts are either built using the DSL or taken from libraries of previously built ones, and adapted to the specific model. The practical usage of the DSL is illustrated first with the Lotka–Volterra model, and then with a landscape modelling experiment on the spread of a mosquito-borne disease in the Sahelian region of West Africa. An interesting characteristic of this approach is the possibility of adding new elements into an existing model, and replacing others with more appropriate ones, thus allowing potentially complex models to be built from simpler parts.     

Eliciting cognitive processes underlying patterns of human-wildlife interactions for agent-based modelling

Integrating humans in our perception of ecosystems is of critical importance to adequately protect natural resources. This poses the challenge of understanding human decision making in the context of decisions potentially threatening nature's integrity. We developed a spatially explicit agent-based model that simulates commercial whale-watching vessel movements based on a representation of the captains’ decision making process when observing marine mammals in and around the Saguenay-St. Lawrence Marine Park in Québec, Canada. We focus here on the human part of the global model, the submodel of whale movements having been developed and validated independently (Lamontagne, 2009). The objective of this study is to select and validate a model of whale-watching captains’ decision making using the pattern-oriented modelling approach (POM): three models of cognitive heuristics (satisficing, tallying and Take The Best) along with a null model (random choice) were tested. These concurrent decision making models were built upon knowledge extracted from data collected during field investigations, including interviews with whale-watching captains and park wardens, onboard and shore-based observations, and analyses of a multi-year dataset of sampled whale-watching excursions. Model selection is performed by statistically comparing simulated and real patterns of boat trajectories (excursion length), spatial hotspots (kernel home range 50%), and excursion content (species observed, time allocated to different activities). The selection process revealed that the Take The Best heuristic was the best performing model. We used the distribution of the number of whale-watching boats in the vicinity (2000 m) of each vessel as a secondary pattern to validate the ability of each decision making model to reproduce real observations. Given the prevalence of the species attribute in the choice of which whale to observe, the Take The Best heuristic's ability to deal with non-compensatory information partly explains its overall best performance. Moreover, implementation of communication abilities between modelled captains led to the emergence of persistent observation sites in the park, which is a well-known collective spatiotemporal characteristic of the whale-watching industry; thus validating the fundamental assumption that cooperation is an important mechanism behind the pattern of whale-watching boat dynamics. The relatively good performance of the satisficing and tallying heuristics supports both field evidence and literature on bounded rationality in that humans likely use collections of heuristics (adaptive toolbox) to solve decision problems in different contexts. The POM strategy appears suitable to build up an informative ABM regarding the management of human activities in a natural environment so that further developments will be assessed following the same approach.     

Quantification of advective solute travel times and mass transport through hydrological catchments

This study has investigated and outlined the possible quantification and mapping of the distributions of advective solute travel times through hydrological catchments. These distributions are essential for understanding how local water flow and solute transport and attenuation processes affect the catchment-scale transport of solute, for instance with regard to biogeochemical cycling, contamination persistence and water quality. The spatial and statistical distributions of advective travel times have been quantified based on reported hydrological flow and mass-transport modeling results for two coastal Swedish catchments. The results show that the combined travel time distributions for the groundwater-stream network continuum in these catchments depend largely on the groundwater system and model representation, in particular regarding the spatial variability of groundwater hydraulic parameters (conductivity, porosity and gradient), and the possible contributions of slower/deeper groundwater flow components. Model assumptions about the spatial variability of groundwater hydraulic properties can thus greatly affect model results of catchment-scale solute spreading. The importance of advective travel time variability for the total mass delivery of naturally attenuated solute (tracer, nutrient, pollutant) from a catchment to its downstream water recipient depends on the product of catchment-average physical travel time and attenuation rate.     

Two-phase modeling of sediment clouds

A sediment cloud release in stagnant ambient fluid occurs in many engineering applications. Examples include land reclamation and disposal of dredged materials. The detailed modeling of the distinct characteristics of both the solid and fluid phases of the sediment cloud is hitherto unavailable in the literature despite their importance in practice. In this paper, the two-phase mixing characteristics of the sediment cloud are investigated both experimentally and theoretically. Experiments were carried out to measure the transient depth penetration and the lateral spread of the sediment cloud and its entrained fluid using the laser induced fluorescence technique, with a range of particle sizes frequently encountered in the field (modeled at laboratory scale). A two-phase model of the sediment cloud that provides detailed predictions of the mixing characteristics of the individual phases is also proposed. The entrained fluid characteristics are solved by an integral model accounting for the buoyancy loss (due to particle separation) in each time step. The flow field induced by the sediment cloud is approximated by a Hill’s spherical vortex centered at the centroid and with the size of the entrained fluid. The particle equation of motion under the effect of the induced flow governs each computational particle. A random walk model using the hydrodynamic diffusion coefficient is used to account for the random fluctuation of particles in the dispersive regime. Overall, the model predictions of the two-phase mixing characteristics are in good agreement with the experimental data for a wide range of release conditions.     

Comparative study of different advective schemes: Application to the MECCA model

A 3-D version of the MECCA model (Model of Estuarine and Coastal Circulation Assessment) is used to simulate the dynamics of the Eastern part of the English Channel. This area is characterized by a strong tidal turbulent regime and a frontal zone identified near the French coast by a low-salinity band. The model uses the upwind scheme to approximate the advective terms. Results show that the model overestimates the band width of the frontal zone and that this anomaly is definitely caused by the numerical diffusion introduced by the so-called upwind scheme (first-order approximation). In this paper, we study the flux-limiter schemes as an alternative to the upwind method in order to reduce this non-physical diffusion. To illustrate the improvements provided by this type of schemes, a comparison in 2-D schematic cases is made between the upwind and centered scheme with more recent higher-order schemes combined with limiter namely Minmod, Superbee, Van Leer and Monotonized Centered (MC) (called also MUSCL scheme). Respecting the CFL condition, our numerical simulations show that the flux-limiter schemes reduce the numerical diffusion and eliminate the oscillations caused by the non-limited higher-order schemes. For the schematic and realistic cases, the Superbee limiter is a good compromise between shape preservation and computational cost.     

The contribution of advective fluxes to net ecosystem exchange in a high-elevation, subalpine forest

The eddy covariance technique, which is used in the determination of net ecosystem CO2 exchange (NEE), is subject to significant errors when advection that carries CO2 in the mean flow is ignored. We measured horizontal and vertical advective CO2 fluxes at the Niwot Ridge AmeriFlux site (Colorado, USA) using a measurement approach consisting of multiple towers. We observed relatively high rates of both horizontal (F(hadv)) and vertical (F(vadv)) advective fluxes at low surface friction velocities (u(*)) which were associated with downslope katabatic flows. We observed that F(hadv) was confined to a relatively thin layer (0-6 m thick) of subcanopy air that flowed beneath the eddy covariance sensors principally at night, carrying with it respired CO2 from the soil and lower parts of the canopy. The observed F(vadv) came from above the canopy and was presumably due to the convergence of drainage flows at the tower site. The magnitudes of both F(hadv) and F(vadv) were similar, of opposite sign, and increased with decreasing u(*), meaning that they most affected estimates of the total CO2 flux on calm nights with low wind speeds. The mathematical sign, temporal variation and dependence on u(*) of both F(hadv) and F(vadv) were determined by the unique terrain of the Niwot Ridge site. Therefore, the patterns we observed may not be broadly applicable to other sites. We evaluated the influence of advection on the cumulative annual and monthly estimates of the total CO2 flux (F(c)), which is often used as an estimate of NEE, over six years using the dependence of F(hadv) and F(vadv) on u(*). When the sum of F(hadv) and F(vadv) was used to correct monthly F(c), we observed values that were different from the monthly F(c) calculated using the traditional u(*)-filter correction by--16 to 20 g C x m(-2) x mo(-1); the mean percentage difference in monthly Fc for these two methods over the six-year period was 10%. When the sum of F(hadv) and F(vadv) was used to correct annual Fc, we observed a 65% difference compared to the traditional u(*)-filter approach. Thus, the errors to the local CO2 budget, when F(hadv) and F(vadv) are ignored, can become large when compounded in cumulative fashion over long time intervals. We conclude that the "micrometeorological" (using observations of F(hadv) and F(vadv)) and "biological" (using the u(*) filter and temperature vs. F(c) relationship) corrections differ on the basis of fundamental mechanistic grounds. The micrometeorological correction is based on aerodynamic mechanisms and shows no correlation to drivers of biological activity. Conversely, the biological correction is based on climatic responses of organisms and has no physical connection to aerodynamic processes. In those cases where they impose corrections of similar magnitude on the cumulative F(c) sum, the result is due to a serendipitous similarity in scale but has no clear mechanistic explanation.     

糖蜜酒精废水两相UASB处理工艺的酸化段特征

文章考察了糖蜜酒精废水二相UASB处理工艺的酸化段特征 ,试验结果表明酸化段的容积负荷达到3 0kgCOD/m3 时 ,系统仍能正常运行 ,TOC去除率在 3 5%以上 ,酸化段的SO4 2 -去除率在 70 %。糖蜜酒精废水酸化率在 50 % ,二相厌氧生物处理系统并未将产酸微生物与产甲烷微生物截然分开 ,在产酸相中仍含有一定数量的产甲烷菌。     

Spring abundance and distribution of the ctenophore Pleurobrachia pileus in the Seine estuary: advective transport and diel vertical migration

Several sets of plankton and hyperbenthos samples were collected from the Seine estuary in late May-early June of 1991, 1992 and 1994, and the abundance, size-frequency, horizontal and vertical distribution, and short-term fluctuations of Pleurobrachia pileus (O.F. Müller, 1776) were determined. In the Seine estuary, P. pileus was more abundant than reported for any other European shallow waters. Its depth-averaged density was 33 individuals m^-3 in 1991 and 32 individuals m^-3 in 1992, with a maximum density of 808 individuals m^-3. The ctenophores displayed passive tidal advection and active diel vertical migration. In the outer estuary, high abundances were observed around low tide, and diel vertical migration was very pronounced. However, many individuals remained aggregated near the bottom at all times. In spring, the distribution of P. pileus coincided with the 33 and 150 S isohalines, with high abundances on the marine side of the estuary. It is suggested that aggregation of ctenophores near the bottom and diel vertical migration play an important role in the retention of P. pileus in the estuary. Two schemes of P. pileus life cycle in the Bay of Seine are proposed.     

Thermodynamic modelling of bio-species accommodation

An ecological model has been constructed. A chemical analogy to the biological system has been developed, and the system's driving forces have been identified.     

糖蜜酒精废水两相UASB处理有机物去除特征

应用气相色谱研究了糖蜜酒精废水两相ＵＡＳＢ处理工艺段中有机物的去除情况,结果表明,酸化反应器内,易降解有机物的去除率上为１００％,产甲烷段反应器内,因易降解化合物在酸化反应器内大部分被去除,剩下的化合物的在该段反应器内去除率偏小,仅为８０％。质谱结果证实,进水中的酚类化合物经两级处理后,被有效去除。     

Two-phase flow insights into open-channel flows with suspended particles of different densities

The effect of particle density on the turbulent open-channel flow carrying dilute particle suspensions is investigated using two specific gravities and three concentrations of solid particles. The particles, identical in size and similar in shape, were natural sand and a neutrally buoyant plastic. The particles were fully suspended, and formed no particle streaks on the channel’s bed. Accordingly, the changes in the flow are attributed to the interactions between suspended particles and flow turbulence structures. Measurements were obtained by means of image velocimetry enabling simultaneous, but distinct, measurement of liquid and particle velocities. The experimental results show that, irrespective of particle specific gravity, particle suspension influences bulk velocity of flow and the Kármán coefficient, while friction velocity essentially remains constant. The results also show that particles in suspension modify local water turbulence over the flow depth, but in ways not accurately predicted using the customary parameters for characterizing turbulence modification.     

Numerical modelling of horizontal sediment-laden jets

Sediment-laden turbulent flows are commonly encountered in natural and engineered environments. It is well known that turbulence generates fluctuations to the particle motion, resulting in modulation of the particle settling velocity. A novel stochastic particle tracking model is developed to predict the particle settling out and deposition from a sediment-laden jet. Particle velocity fluctuations in the jet flow are modelled from a Lagrangian velocity autocorrelation function that incorporates the physical mechanism leading to a reduction of settling velocity. The model is first applied to study the settling velocity modulation in a homogeneous turbulence field. Consistent with basic experiments using grid-generated turbulence and computational fluid dynamics (CFD) calculations, the model predicts that the apparent settling velocity can be reduced by as much as 30 % of the stillwater settling velocity. Using analytical solution for the jet mean flow and semi-empirical RMS turbulent velocity fluctuation and dissipation rate profiles derived from CFD predictions, model predictions of the sediment deposition and cross-sectional concentration profiles of horizontal sediment-laden jets are in excellent agreement with data. Unlike CFD calculations of sediment fall out and deposition from a jet flow, the present method does not require any a priori adjustment of particle settling velocity.     

Mathematical modelling of a fish pond ecosystem

A mathematical model is constructed for a fish breeding pond for carp, silver carp and bighead. The model is a system of ordinary differential equations describing the material transformations in the ecosystem. It allows a choice of optimal regimes of the aeration, feeding and fertilization of a pond for different climatic conditions in order to maximize the yield.     

On the integrated modelling of coastal changes

This paper introduces a possible approach to integrated modelling of coastal change, focusing on coastal land use and cover change. Some of the mostimportant open issues in the context of integrated modelling of coastal change are introduced. The paper focuses on methdological aspects. Specific reference is made to Physiographic Unit Modelling as an approach to better handle spatial variability and ‘morphogenesis’, and as a way to focus on coastal change mechanisms instead of absolute coastal dynamics for achieving an important simplification of the problem. The application is briefly discussed with reference to a ‘minimal model’. The methodological structure introduced is considered particularly suitable to represent, according to a variable degree of simplification, the integrative dynamics between resources and uses of the resources.     

3D numerical modelling of turbidity currents

During floods, the density of river water usually increases due to a subsequent increase in the concentration of the suspended sediment that the river carries, causing the river to plunge underneath the free surface of a receiving water basin and form a turbidity current that continues to flow along the bottom. The study and understanding of such complex phenomena is of great importance, as they constitute one of the major mechanisms for suspended sediment transport from rivers into oceans, lakes or reservoirs. Unlike most of the previous numerical investigations on turbidity currents, in this paper, a 3D numerical model that simulates the dynamics and flow structure of turbidity currents, through a multiphase flow approach is proposed, using the commercial CFD code FLUENT. A series of numerical simulations that reproduce particular published laboratory flows are presented. The detailed qualitative and quantitative comparison of numerical with laboratory results indicates that apart from the global flow structure, the proposed numerical approach efficiently predicts various important aspects of turbidity current flows, such as the effect of suspended sediment mixture composition in the temporal and spatial evolution of the simulated currents, the interaction of turbidity currents with loose sediment bottom layers and the formation of internal hydraulic jumps. Furthermore, various extreme cases among the numerical runs considered are further analyzed, in order to identify the importance of various controlling flow parameters.     

A review of recent developments in lake modelling

This paper reviews the lake models published the last five years, mainly in Ecological Modelling. The review shows that structurally dynamic modelling and coupling between hydrodynamic and ecological models are applied increasingly. A number of processes that have not been included in lake models before have been proposed. It has been shown that these additional processes in specific case studies are significant, for instance the competition between phytoplankton and macrophytes or cyanobacteria growth and growth of mussels. It is recommended to study these models for the development of models for case studies where these processes are relevant.