Citizen participation for integrated assessment: new pathways in complex systems

Citizen or community participation is a guiding principle for environmental decision-making, but in practice it is often compromised and contradictory. Two themes in participation can be defined. One is in the 'product' or outcome of decision-making, and aims for more effective and appropriate evaluations and decisions. The other is participation in the 'process', where devolved decision-making itself can serve wider social and economic goals. Although both strands are crucial for the practice of 'integrated assessment' in environmental governance, a clearer view is needed of the scope and limits of such participation. This paper explores such scope and limits by applying a conceptual 'integrated assessment' framework. This is drawn from a major case study on sustainable development for a city-region in the north of the UK. Such a framework helps to identify the many roles of citizen participation from a complex systems perspective. Selected examples of current practice highlight the problems and opportunities for citizen participation in environmental governance.     

Simultaneous estimation of effective and apparent diffusion coefficients in compacted bentonite

Effective diffusion coefficients (D(e)) are usually measured by means of "through-diffusion" experiments in which steady state is reached, and the "time-lag" methods are used to estimate the apparent diffusion coefficient (D(a)). For sorbing radionuclides (as caesium), the time needed to reach steady-state conditions is very large, and the precision in D(a) determinations is not satisfactory. In this paper, a method that allows determining simultaneously effective and apparent diffusion coefficients in compacted bentonite without reaching steady-state conditions is described. Basically, this method consists of an "in-diffusion" experiment in which the concentration profile in the bentonite sample is used to estimate D(a), and the temporal evolution of the solute concentration in the reservoir is used to estimate D(e). This method has several advantages over the typical "through-diffusion" experiments, in particular: (a) experiment duration is significantly shorter, (b) D(a) values are measured with greater precision and (c) it is not necessary to maintain a constant solute concentration in the reservoir. This new method has been used to estimate the effective and apparent diffusion coefficients for caesium in FEBEX bentonite and in order to validate it, the results have been compared with results previously obtained with standard methods. Experimental results have been satisfactorily modelled using a simple model of diffusion in porewater and the measured value of D(e)(Cs) is very similar to D(e)(HTO) in the same bentonite. There is no evidence of "surface diffusion" in FEBEX bentonite for caesium.     

Eulerian derivation of the fractional advection-dispersion equation

A fractional advection-dispersion equation (ADE) is a generalization of the classical ADE in which the second-order derivative is replaced with a fractional-order derivative. In contrast to the classical ADE, the fractional ADE has solutions that resemble the highly skewed and heavy-tailed breakthrough curves observed in field and laboratory studies. These solutions, known as alpha-stable distributions, are the result of a generalized central limit theorem which describes the behavior of sums of finite or infinite-variance random variables. We use this limit theorem in a model which sums the length of particle jumps during their random walk through a heterogeneous porous medium. If the length of solute particle jumps is not constrained to a representative elementary volume (REV), dispersive flux is proportional to a fractional derivative. The nature of fractional derivatives is readily visualized and their parameters are based on physical properties that are measurable. When a fractional Fick's law replaces the classical Fick's law in an Eulerian evaluation of solute transport in a porous medium, the result is a fractional ADE. Fractional ADEs are ergodic equations since they occur when a generalized central limit theorem is employed.     

Modelling radionuclide transport for time varying flow in a channel network

Water flowrates and flow directions may change over time in the subsurface for a number of reasons. In fractured rocks flow takes place in channels within fractures. Solutes are carried by the advective flow. In addition, solutes may diffuse in and out of stagnant waters in the rock matrix and other stagnant water regions. Sorbing species may sorb on fracture surfaces and on the micropore surfaces in the rock matrix. We present a method by which solute particles can be traced in flowing water undergoing changes in flowrate and direction in a complex channel network where the solutes can also interact with the rock by diffusion in the rock matrix. The novelty of this paper is handling of diffusion in the rock matrix under transient flow conditions. The diffusive processes are stochastic and it is not possible to follow a particle deterministically. The method therefore utilises the properties of a probability distribution function for a tracer moving in a fracture where matrix diffusion is active. The method is incorporated in a three dimensional channel network model. Particle tracking is used to trace out a multitude of flowpaths, each of which consists of a large number of channels within fractures. Along each channel the aperture and velocity as well as the matrix sorption properties can vary. An efficient method is presented whereby a particle can be followed along the variable property flowpath. For stationary flow conditions and a network of channels with advective flow and matrix diffusion, a simple analytical solution for the residence time distribution along each pathway can be used. Only two parameter groups need to be integrated along each path. For transient flow conditions, a time stepping procedure that incorporates a stochastic Monte-Carlo like method to follow the particles along the paths when flow conditions change is used. The method is fast and an example is used for illustrative purposes. It is exemplified by a case where land rises due to glacial rebound. It is shown that the effects of changing flowrates and directions can be considerable and that the diffusive migration in the matrix can have a dominating effect on the results.     

Colloid migration in symmetrical non-uniform fractures: particle tracking in three dimensions

A three-dimensional model for the migration of colloids in a saturated fracture is presented, which considers the motion of colloids as a result of advection and diffusion, as well as colloid-surface interactions at the fracture walls. This model is successfully incorporated into a three-dimensional particle tracking algorithm that tracks particles within a continuum and allows consideration of the migration of colloids in symmetrical, three-dimensional, non-uniform fractures. The framework is general enough to incorporate non-local interactions that provide colloid motion relative to the fluid. The algorithm is verified against classical Taylor dispersion, and its generalization to a sorbing phase, in a uniform fracture and shows excellent agreement with theory. A simple, non-uniform fracture that has an analytically tractable velocity field is also considered, and both quantitative and qualitative comparisons are made with the uniform fracture case. The modelling of more complex fracture geometries is also discussed and a particular case is implemented within the particle tracking framework.     

Evaluation of a simplified top-down model for the spatial assessment of hot traffic emissions in mid-sized cities

Traffic emission estimation in developing countries is a key-issue for air pollution management. In most cases, comprehensive bottom-up methodologies cannot be applied in mid-sized cities because of the resource cost related to their application. In this paper, a simplified emission estimation model (SEEM) is evaluated. The model is based on a top-down approach and gives annual global hot emission. Particular attention is paid to the quality of the input traffic data. The quality of results is assessed by application of the SEEM model in the Chilean Gran Concepción urban area and by comparison with a bottom-up approach that has been led for the year 2000. The SEEM model estimates emissions with an accuracy of about $20\%$ and is related to important resource savings. The results of the SEEM model are then distributed in space with a disaggregation approach and using GIS techniques. The relevancy of the disaggregation approach is evaluated among several possibilities through statistical methods. A spatial disaggregation using principal roads density gives the best results in terms of emissions repartition and gives a globally accurate image of the distribution of hot emissions in a mid-sized city.     

Formulation of a soil-pesticide transport model based on a compartmental approach

A semianalytical soil-pesticide transport model is formulated based on a compartmental approach to determine spatial and temporal variations of pesticide residues across a soil profile. The compartmental model is implemented by drawing an analogy between a series of continuous-flow stirred tank reactors and a soil horizon that consists of multiple perfectly mixed compartments. The analogy is strengthened by exploiting a relation between the compartment series and the conventional convective-dispersive equation (CDE) for vertical transport in the soil. Consequently, the number of compartments in the model formulation is not free, but dictated as a function of transport parameters. The model formulation allows consideration of arbitrary boundary value specifications and also, for some cases, spatially varying initial concentration profiles. Sorption kinetics is represented via a two-site model that involves a linear sorption isotherm and a first-order irreversible sorption or a radial diffusive penetrating model. For these three cases, analysis of the compartmental model allows the resultant concentration profiles to be expressed in terms of the Poisson distribution. When a nonlinear kinetic sorption model is used to simulate the sorption processes, an analytical solution is not found and a numerical approach is required.     

Potential protective effects of Spirulina platensis on liver,kidney, and brain acrylamide toxicity in rats

Environmental Science and Pollution Research - Acrylamide (AA) is a hazardous chemical that is widely used in industrial practices. Spirulina platensis (SP) is a blue green alga that is rich in...     

Economic insight into municipal solid waste management in Mumbai: a critical analysis

Mumbai generates 6256 tonnes of waste every day, of which 17.20% is recyclable, but only a fraction of this is retrieved by rag-pickers. The economic value of the retrieved material is not considered by the Municipal Corporation of Greater Mumbai in valuing the waste management system as there is no retrieval mechanism except the informal rag-picking activity. Moreover, the cost of land used for the dumping of waste is also not accounted for. In the present paper, a comprehensive cost-benefit analysis for the present system of municipal solid waste management in Mumbai is carried out, with due consideration for implicit or hidden costs and benefits. Accounting for the implicit costs and benefits showed a difference of $6 per every tonne of waste disposal. This could show a considerable difference in policy development at the municipality level. Demand supply analysis proved that the present system of waste management would not yield a feasible market solution without private sector participation. With the increasing demand for improved waste management, private sector participation is essential and a Pigouvian tax is a necessary tool to make the private sector participation in solid waste management a success.     

Modeling of flow and contaminant transport in coupled stream-aquifer systems

This paper describes an integrated approach for modeling flow and contaminant transport in hydraulically connected stream-aquifer systems. The code, FTSTREAM, extended the capabilities of the ground-water model, FTWORK, to incorporate chemical fate and transport in streams. Flow in the stream network is modeled as an unsteady, spatially varying flow, while transport modeling is based on a one-dimensional advection-dispersion equation. In addition to sorption and decay during transport in ground water, the model incorporates volatilization, settling and decay during transport in surface water. The interaction between surface water and ground water is accommodated by a leakage term and is implemented in the model using an iterative Picard-type procedure to ensure mass conservation across the interface between the two systems. The modeling approach is used to simulate contaminant transport in the Mad River, Ohio, which is hydraulically connected to a buried valley aquifer of sand and gravel outwash. The river is a receiving stream in the upstream part of the modeled area. Downstream, heavy pumping from a municipal well field causes the river to become a loosing stream. Induced infiltration from the river is responsible for a considerable portion of the well yield. The flow and transport model, developed for this study, simulates coupling between flow in the aquifer and the river. Hypothetical sources of contamination are introduced at selected locations in the upstream portion of the aquifer. The model is then used to simulate the expected transport in both the aquifer and the stream. A series of simulations elucidates the role of the river in facilitating the transport of the hypothetical contaminants in ground water and surface water. Effect of sorption, retardation and volatilization on contaminant transport is also examined for the case of the volatile organic compounds.     

Scavenging of soluble trace gases by falling rain droplets in inhomogeneous atmosphere

We analyze non-isothermal absorption of trace gases by the rain droplets with internal circulation which is caused by interfacial shear stresses. It is assumed that the concentration of soluble trace gases and temperature in the atmosphere varies in a vertical direction. The rate of scavenging of soluble trace gases by falling rain droplets is determined by solving heat and mass transfer equations. In the analysis we accounted for the accumulation of the absorbate in the bulk of the falling rain droplet. The problem is solved in the approximation of a thin concentration and temperature boundary layers in the droplet and in the surrounding air. We assumed that the bulk of a droplet, beyond the diffusion boundary layer, is completely mixed and concentration of the absorbate and temperature are homogeneous and time-dependent in the bulk. By combining the generalized similarity transformation method with Duhamel's theorem, the system of transient conjugate equations of convective diffusion and energy conservation for absorbate transport in liquid and gaseous phases with time-dependent boundary conditions is reduced to a system of linear convolution Volterra integral equations of the second kind which is solved numerically. Calculations are performed using available experimental data on concentration and temperature profiles in the atmosphere.It is shown than if concentration of a trace gas in the atmosphere is homogeneous and temperature in the atmosphere decreases with height, beginning from some altitude gas absorption is replaced by gas desorption. Neglecting temperature inhomogenity in the atmosphere described by adiabatic lapse rate leads to essential overestimation of the trace gas concentration in a droplet on the ground.     

The Annual Average: An Alternative to the Second Highest Value as a Measure of Air Quality

It is common practice to use the second highest value both in determining compliance with the once-per-year air quality standards (AQS) and as a measure of air quality in pollutant trend studies and rollback calculations. A study of the variation in the second highest 8 hr carbon monoxide concentration observed at the CAMP stations 1962-72 is presented. It is shown that, for a given annual average, the second highest value can differ by a factor of 2 due to random variation. The annual average is linearly related to the observed average of the second highest value and is shown to be a good predictor of the percent of time the carbon monoxide AQS is exceeded. It is concluded that the annual average, which is less variable and not as greatly influenced by erroneous measurements, is the preferred measure of air quality for trend studies and air quality projections.     

Periphyton as a bioindicator of mercury pollution in a temperate torrential river ecosystem

Mercury presents a potential risk to the environment and humans, especially in its methylated form. It is among the highest priority environmental pollutants. River Idrijca (Slovenia) is highly contaminated with mercury due to past mercury mining. The aim of this work was to investigate whether the periphyton community in rivers such as Idrijca is a suitable indicator of Hg pollution and of changes in mercury methylation and could serve as an early warning system of increased input of MeHg in the food chain. Periphyton is the only site of primary production in temperate torrential rivers such as Idrijca and is therefore an important link in the food chain. It is also a potential site of Hg accumulation and its introduction to higher trophic levels. Our aim was to assess the response of the periphyton to seasonal and spatial variations in mercury levels and to evaluate its potential as an early warning system of changes in mercury reactivity and mobilization The results indicate that periphyton in a torrential river is too complex and unpredictable to be used as a sole indicator of mercury concentrations and changes in the river. Nevertheless, it can complement environmental measurements due to its importance in the riverine food web.     

Multi-criteria analysis and risk assessment

ELECTRE or Concordance Analysis is a French system of multicriteria analysis which is used in economics.The application of ELECTRE for hazard evaluation and risk assessment of chemical substances is suggested, and a first attempt in adapting it for this purpose is presented.     

Noncatalytic Auto Exhaust Combustor

Present methods for the determination of carbon monoxide are discussed including indicator tubes, the iodine pentoxide reaction and measurement by gas chromatography. In the gas chromatographic method an air sample is separated on a gas-solid chromatogra-phic column and the separated CO is converted to methane by hydrogenation at elevated temperature. The separated CO, in the form of methane, is passed into a hydrogen flame detector and measured. The conversion from CO to methane allows the use of a sensitive ionization detector in place of the thermal conductivity cell which is insufficiently sensitive for the measurement of trace amounts of CO. The optimum operating conditions are discussed. It is possible to determine one ppm CO in air. The iodine pentoxide reaction with CO has been combined with electrometric measurement. The iodine liberated is passed into a Ditte cell and the current generated is measured by an electrometerrecorder combination. This method is continuously direct reading with a permanent record. It is suitable for the continuous routine analysis of one ppm CO.     

Thermal and air flow characteristics in a deep pedestrian canyon under hot weather conditions

The present paper reports the results of a specific experiment carried out in a deep pedestrian canyon in Athens, Greece. The air flow inside and outside the canyon as well as the air and surface temperatures inside the canyon are measured for a period of seven continuous days during the whole day and night period. The air flow characteristics inside the canyon are analyzed. For perpendicular wind speeds, the air flow inside the canyon is characterized by either a circulatory vortex driven by the ambient air flow or a double vortex flow where the upper one is driven by the ambient air flow and the lower one is driven in the opposite direction by the circulation above. When the ambient air blows parallel to the canyon, the measured wind data inside the canyon indicate a flow in the same direction. When the air flows at a certain angle to the canyon axis, the measured data indicate that a spiral vortex is induced along the canyon length. The spatial and temporal variation of the surface and air temperature distribution inside the canyon is analyzed. The heat island intensity in the canyon is estimated based on measurements from a reference station.     

Long-term effects of soil heterogeneity on cadmium behaviour in soil

A deterministic model for long-term behaviour of contaminants in the rootzone is developed that includes sorption, leaching, and plant uptake. The model is applied to cadmium accumulation in a sandy soil and uptake of cadmium by barley. Sensitivity analysis showed that the sensitivity of the leaching rate to changes in soil chemical and soil physical parameters decreases as a function of time, and becomes zero when steady state is reached. In contrast, accumulation of cadmium in soil and the plant uptake rate of barley are increasingly sensitive to soil chemical and soil physical parameters as time preceeds. To analyse cadmium behaviour in a field that is heterogeneous with respect to soil physical properties, the interstitial flow velocity was assumed to be a random, lognormally distributed variate. Using Monte Carlo simulation, the average plant uptake rate appeared to be much higher in the stochastic analysis than in the deterministic approach. Steady state is reached after a very long period of time. For a lognormally distributed proton activity, causing heterogeneity with respect to the sorption capacity of the soil, the model predicted similar deviations from the deterministic approach. It is concluded that reference values for groundwater and crop quality are exceeded earlier in a heterogeneous field than in a homogenous soil profile. Moreover, when average values suggest an acceptable situation, variability of the leaching rate and the plant uptake rate can still cause exceedance of reference values in part of the field. Therefore, it is reasoned that environmental quality standards should take soil heterogeneity into account.     

GIS环境下流域降雨侵蚀动态模拟研究——以PCRaster系统和LISEM模型为例

本文结合ＧＩＳ和流域土壤侵蚀模型,对流域单次降雨侵蚀过程进行动态模拟研究。首先概述土壤侵蚀模型与ＧＩＳ结合的三种方式和特点;其次介绍地理信息系统ＰＣＲａｓｔｅｒ 的空间动态模型语言及其表达方法、命令函数和动态模型的逻辑组成,以及使用该语言实现的与ＧＩＳ完全结合的流域水文和土壤侵蚀动态模型－ ＬＩＳＥＭ 模型;最后分析用ＰＣＲａｓｔｅｒ 系统和ＬＩＳＥＭ 模型对黄土高原小流域次降雨过程进行动态模拟的实验结果,认为模型计算的总径流量和总土壤流失量具有较高的精度,模拟的径流过程线符合黄土高原小流域超渗产流的实际情况     

Unsteady absorption of sulfur dioxide by an atmospheric water droplet with internal circulation

Unsteady absorption characteristics of sulfur dioxide by an atmospheric water droplet in motion are predicted numerically and analyzed theoretically to recognize the physical mass transport processes inside an aerosol droplet, which is frequently encountered in the atmosphere. Considering the absorption of sulfur dioxide by a droplet in cloud or fog with various velocities, three different Reynolds numbers, viz., Re_g=0.643, 1.287, and 12.87 are studied and compared with each other. The results indicate that for the Reynolds number of 0.643, sulfur dioxide always penetrates toward the droplet centerline throughout the entire absorption period. This is due to the mass transfer dominated by diffusion along the radial direction. In contrast, when the Reynolds number is 12.87, the strength of the vortex motion inside the droplet is strong enough. It results in that, most of the time the concentration contours parallel the streamlines and the lowest SO₂ concentration is located at the vortex center. As a consequence, the diffusion distance is reduced by a factor of three and the absorption time for the droplet reaching the saturated state is shortened in a significant way. With regard to an intermediate Reynolds number such as 1.287, a two-stage mass transfer process can be clearly identified. In the first stage, it is dominated by one-dimensional diffusion, in which over 50% sulfur dioxide is absorbed before the saturated state is reached. In the second stage, the vortex motion mainly controls the mass transfer. However, the contour core is inconsistent with the vortex center. This is because the characteristic time of mass diffusion is in a comparable state with that of droplet internal circulation. The present study elucidates that the strength of a droplet's internal motion plays a vital role in determining SO₂ absorption process.     

An efficient Lagrangian stochastic model of vertical dispersion in the convective boundary layer

We consider the one-dimensional case of vertical dispersion in the convective boundary layer (CBL) assuming that the turbulence field is stationary and horizontally homogeneous. The dispersion process is simulated by following Lagrangian trajectories of many independent tracer particles in the turbulent flow field, leading to a prediction of the mean concentration. The particle acceleration is determined using a stochastic differential equation, assuming that the joint evolution of the particle velocity and position is a Markov process. The equation consists of a deterministic term and a random term. While the formulation is standard, attention has been focused in recent years on various ways of calculating the deterministic term using the well-mixed condition incorporating the Fokker–Planck equation. Here we propose a simple parameterisation for the deterministic acceleration term by approximating it as a quadratic function of velocity. Such a function is shown to represent well the acceleration under moderate velocity skewness conditions observed in the CBL. The coefficients in the quadratic form are determined in terms of given turbulence statistics by directly integrating the Fokker–Planck equation. An advantage of this approach is that, unlike in existing Lagrangian stochastic models for the CBL, the use of the turbulence statistics up to the fourth order can be made without assuming any predefined form for the probability distribution function (PDF) of the velocity. The main strength of the model, however, lies in its simplicity and computational efficiency. The dispersion results obtained from the new model are compared with existing laboratory data as well as with those obtained from a more complex Lagrangian model in which the deterministic acceleration term is based on a bi-Gaussian velocity PDF. The comparison shows that the new model performs well.