Mathematical modeling and analysis of the depletion of dissolved oxygen in water bodies

In this paper, a nonlinear mathematical model is proposed to study the depletion of dissolved oxygen in a water body caused by industrial and household discharges of organic matters (pollutants). The problem is formulated as a food chain model by considering various interaction processes (biodegradation and biochemical) involving organic pollutants, bacteria, protozoa, an aquatic population and dissolved oxygen. Using stability theory, it is shown that as the rate of introduction of organic pollutants in a water body increases, the concentration of dissolved oxygen decreases due to various interaction processes. It is found that if the organic pollutants are continuously discharged into water body, the concentration of dissolved oxygen may become negligibly small, thus threatening the survival of aquatic populations. However, by using some effort to control the cumulative discharge of these pollutants into the water body, the concentration of dissolved oxygen can be maintained at a desired level.     

Mathematical modeling and analysis of the depletion of dissolved oxygen in eutrophied water bodies affected by organic pollutants

In this paper, an ecological type nonlinear mathematical model is proposed to study the simultaneous effect of water pollution and eutrophication on the concentration of dissolved oxygen (DO) in a water body. It is assumed that the organic pollutants and the nutrients are discharged into water body from outside with constant rates. The system is modeled by considering the variables such as cumulative concentration of organic pollutants, the densities of bacteria, nutrients, algae, detritus and the concentration of DO. The analysis of the model shows that the decrease in the concentration of DO due to simultaneous effect of water pollution and eutrophication is much more than when only single effect is present in the water body, thus leading to more uncertainty about the survival of DO-dependent species.     

Modeling the depletion of dissolved oxygen in a water body located near a city

Water bodies located nearby cities are much prone to pollution, especially in the developing countries, where effluents treatment facilities are generally lacking. The main reason for this phenomenon is the increasing population in the cities, and the large number of industries located near them. This leads to generation of huge amounts of domestic and industrial sewage that is discharged into the water bodies, increasing their organic pollutant load and resulting in the depletion of dissolved oxygen. In this paper, we propose a mathematical model for this situation, focusing especially on the resulting quality of the water, determined by the level of dissolved oxygen. The model also accounts for resources needed for the population survival and for the industrial operations. In addition, we describe also the decomposition of organic pollutants by bacteria in the aquatic medium. Feasibility conditions and stability criteria of the system's equilibria are determined analytically. The results show that human population and industries are relevant influential factors responsible for the increase in organic pollutants and the decrease in dissolved oxygen in the water body, in the sense that they may exert a destabilizing effect on the system. The numerical simulations confirm the analytical results. Copyright © 2016 John Wiley & Sons, Ltd.     

Modeling and analysis of the removal of an organic pollutant from a water body using fungi

In this paper, a non linear mathematical model for removing an organic pollutant such as a dye from a water body is proposed and analyzed. In the modeling process four variables are considered, namely, (i) the concentration of the dye, (ii) the density of fungus population, (iii) the concentration of a nutrient and (iv) the concentration of dissolved oxygen (DO). It is assumed that an organic pollutant is present in water with given concentration or discharged with a constant rate in water. It is assumed further that fungus population is kept alive and active due to supply of a nutrient. It is considered that nutrient and oxygen are supplied to the water body from outside with constant rates. The model is analyzed by using the stability theory of differential equations. The model analysis shows that organic pollutant can be removed from the water body by fungus population and the level of degradation depends upon the concentration of organic pollutant, the density of fungal population and the interaction processes involved.The simulation analysis of the proposed model confirms the analytical results. It is also found that these results are qualitatively in line with the experimental observations of one of the authors (Sanghi).     

Modeling the depletion of dissolved oxygen due to algal bloom in a lake by taking Holling type-III interaction

In this paper a non-linear mathematical model for depletion of dissolved oxygen due to algal bloom in a lake is proposed and analyzed. The model is formulated by considering four variables namely, cumulative concentration of nutrients, density of algal population, density of detritus and concentration of dissolved oxygen. In the modeling process it is assumed that nutrients are continuously coming with a constant rate to the lake through water runoff from agricultural fields and domestic drainage. The Holling type-III interaction between nutrients and algal population is considered. Equilibrium values have been obtained and their stability analysis has also been performed. Numerical simulations are carried out to explain the mathematical results.     

A nonlinear mathematical model to study the interactions of hot gases with cloud droplets and raindrops

In this paper, a nonlinear mathematical model is proposed and analyzed to study the interactions of hot gases with cloud droplets as well as with raindrops and their removal by rain from the stable atmosphere. The atmosphere, during rain, is assumed to consist of five nonlinearly interacting phases i.e. the vapour phase, the phase of cloud droplets, the phase of raindrops, the phase of hot gaseous pollutants and the absorbed phase of hot gases in the raindrops (if it exists). It is further assumed that these phases undergo ecological type growth and nonlinear interactions. The proposed model is analyzed using stability theory of differential equations and by numerical simulation. It is shown that the cumulative concentration of gaseous pollutants decreases due to rain and its equilibrium level depends upon the density of cloud droplets, the rate of formation of raindrops, emission rate of pollutants, the rate of falling absorbed phase on the ground, etc. It is noted here that if gases are very hot, cloud droplets are not formed and rain may not take place. In such a case gaseous pollutants may not be removed from the atmosphere due to non-occurrence of rain.     

A solution to the distribution problems arising in the studies of a two-sex population process

Given a population of two sexes, the birth rate of one sex of which depends upon the population size of the other, it is very difficult to find an explicit expression for the probability distribution of the former. In this paper we have explicitly found the probability generating function of the joint distribution from which individual probability distributions and, in particular, moments of all orders in each case can be obtained in principle. As an example, using this probability generating function we have worked out explicitly the first and second order moments of the male and female populations and the explicit expression for the distribution of the male population in a particular case. This method can be successfully applied for the same purpose in the studies of chemical and biological processes where the synthesis or production of one species depends upon the concentration of another species.     

Velocity Vector Field Optimization in Bioventing

Bioventing is a technology used to remove some kinds of pollutants from the subsoil and it is based on the capability of some bacteria species to biodegrade contaminants. The biochemical reaction requires, among other things, oxygen and, therefore, oxygen is inflated into the subsoil by wells. The mathematical model describes the movement of the different fluids which are present in the subsoil – air, water, pollutants, oxygen and so on – and the bacteria population dynamics. The presence of source reactive terms in the continuity equations allows the contaminant biodegradation to be described. The design of a subsoil decontamination intervention concerns bioavailability problems and, in particular, the oxygen concentration. Therefore, in order to enhance the biodegradation phenomenon, the optimization of the subsoil oxygen velocity field in the polluted area is required, by an appropriate choice of the well positions and of the well air inflating rates. In mathematical terms, the goal is to obtain the decontamination of the subsoil with an optimal value of an objective function by acting on some control variables which, in this case, are the well positions and the inflating rates. In this paper several kind of objective function are proposed. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling the removal of gaseous pollutants and particulate matters from the atmosphere of a city

In this paper, a nonlinear mathematical model is proposed and analyzed to study the removal of gaseous pollutants and particulate matters from the atmosphere of a city by precipitation. The atmosphere consists of four interacting phases i.e. the raindrops phase, the gaseous pollutants phase, the phase of gaseous pollutants absorbed (dissolved) in rain drops and the phase of particulate matters. The dynamics of these phases is assumed to be governed by ordinary differential equations with source, interaction, removal and recycle terms. The proposed model is analyzed by using stability theory of differential equations. It is shown that the pollutants can be removed from the atmosphere and their removal rates depend mainly upon the rates of emission of the pollutants, rate of rain drops formation and the rate of falling rain drops on the ground. If the rate of precipitation is very high, the pollutants may be removed completely from the atmosphere.     

Asymptotic Analysis of Reacting Materials with Saturated Explosion. II. High-Frequency Waves

The interaction of small-scale material inhomogeneities with high-frequency acoustic waves is known to have a prominent role in accelerating the heat-release rate in liquid and solid explosive materials. In the present paper, simplified asymptotic equations are studied which incorporate the above interaction, and which include reactant depletion at leading order. Because fuel may be completely exhausted, singularities do not always form in the model equations; it is conjectured that when a singularity does form, the material has initiated. The detailed mechanisms by which shock formation and resonant wave interaction can either enhance or retard reaction are explored. In a realistic model for inhomogeneous condensed-phase reaction, with pressure-dependent reaction rate and nonconstant initial fuel concentration, initiation of the material depends on correct placement of the fuel relative to the acoustic waves.     

Random signal analysis in an environmental sciences problem

Two major difficulties are encountered in the identification of wastewater treatment plant and river water quality dynamics: process behaviour can neither be easily observed, nor easily experimented upon, and the underlying biological nature of the processes involved is only partially understood. This paper describes the derivation of a model for the interaction between dissolved oxygen (DO), biochemical oxygen demand (BOD), and algae in a freshwater river. Noisy measurements from a stretch of the River Cam downstream of Cambridge are analysed using various techniques of identification, parameter estimation and filtering. An important feature of the paper is the interpretation of system identification as a hypothesis testing/decision making procedure.     

UNSTEADY STATE DISPERSION OF AIR POLLUTANTS UNDER THE EFFECTS OF DELAYED AND NONDELAYED REMOVAL MECHANISMS

Abstract In this paper, we present a two‐dimensional time‐dependent mathematical model for studying the unsteady state dispersion of air pollutants emitted from an elevated line source in the atmosphere under the simultaneous effects of delayed (slow) and nondelayed (instantaneous) removal mechanisms. The wind speed and coefficient of diffusion are taken as functions of the vertical height above the ground. The deposition of pollutants on the absorptive ground and leakage into the atmosphere at the inversion layer are also included in the model by applying appropriate boundary conditions. The model is solved numerically by the fractional step method. The Lagrangian approach is used to solve the advection part, whereas the Eulerian finite difference scheme is applied to solve the part with the diffusion and removal processes. The solutions are analyzed to observe the effects of coexisting delayed and nondelayed removal mechanisms on overall dispersion. Comparison of delayed and nondelayed removal processes of equal capacity shows that the latter (nondelayed) process is more effective than the former (delayed removal) in the removal of pollutants from the atmosphere.     

Large time behaviour of solutions to a class of non-autonomous, degenerate parabolic equations

We consider a class of non-autonomous, degenerate parabolic equations and we study the asymptotic behaviour of the solutions. Even if the equation depends explicitly upon the time, we prove that several asymptotic properties, valid for the autonomous case, are preserved in this more general situation. To our knowledge, it is the first time that the asymptotic behaviour of solutions to non-autonomous equations is studied.     

MODELING EFFECTS OF TWO INTERACTING POPULATIONS ON A RENEWABLE RESOURCE

Abstract In this paper, a general mathematical model is presented to study the effect of two populations on a resource biomass. The interaction between two populations is assumed to be competition, predation, or cooperation. These two populations may depend on the resource biomass partially, wholly, or they may predate on the resource. In each case, criteria for local stability, instability, and global stability are obtained. Numerical simulation is presented to illustrate the theoretical results obtained in each case. It is shown that the depletion of resource biomass is maximum in the case of cooperation and is minimum in the case of competition.     

Dynamic modelling of dissolved oxygen in the creeks of Sagar island,Hooghly–Matla estuarine system,West Bengal,India

Hooghly–Matla estuarine ecosystem is one of the largest estuarine ecosystems of the world. Sagar island is the largest delta in this estuarine complex. This island is criss-crossed by small and large creeks with mangrove vegetation and all are connected to the principal estuarine water. Decomposition of mangrove litter in soil is major source of inorganic nutrient to phytoplankton of the adjacent creeks. Deforestation of mangrove affects the primary production, which in turn reduces the availability of dissolved oxygen for the organisms residing in the estuary. Considering the importance of dissolved oxygen in various aspects of aquatic life, a dynamic model of dissolved oxygen at Sagar island of Hooghly–Matla estuarine complex with the help of single dimension differential equation is proposed in the present paper. Different physical, chemical and biological factors such as solar irradiance, temperature, salinity of water, particulate organic matter, re-aeration, wind velocity, phytoplankton and zooplankton, which control the fluctuation of dissolved oxygen, are included in the present model. Most of the parameter values are collected directly from the field surveys. The parameter values which are not able to collect from the field, obtained from literatures are calibrated. To make the model realistic it is properly validated with observed data and to know the statistical significance, chi square goodness fit test is performed. Field surveys are performed over two years. During calibration and validation, two sets of data (first year and second year data) are used. Chi-square values are 5.97 and 6.17 for first and second sets of data respectively (p < 0.05). Sensitivity analysis reveals that optimal light intensity is the most sensitive parameter for dissolved oxygen dynamics. Results also show that wind velocity, solar irradiation, salinity of water and temperature are important factors for controlling the dynamics of dissolved oxygen. Macrophytes have very little contribution to oxygen production in the creeks of Sagar island. Model reveals that low dissolved oxygen in the creek water is one of the causes of decline in fish population of the estuary.     

MODELING THE EFFECTS OF WOOD AND NON‐WOOD BASED INDUSTRIES ON FORESTRY RESOURCES

Pollution due to the wood and non‐wood based industries is one of the main reasons for the depletion of forestry resources. On the onset of industrial era, the number of industries was less but now it has increased much due to rapid pace of development. Industrialization affects forestry resources in two ways: (i) using raw materials by wood based industries and cut down by non‐wood based industries and (ii) emitting pollutants by wood and non‐wood based industries. In view of this, in this paper, we have proposed and analyzed a mathematical model to assess the effects of wood and non‐wood based industries on the depletion of forestry resources. In the modeling process, it is assumed that wood and non‐wood based industries deplete the forestry resources directly, whereas pollutants emitted by both types of industries decrease the growth rate of forestry resources indirectly. The equilibria are obtained and their stability is discussed. The model analysis reveals that forestry resources decrease due to the growth in wood based industries but the overgrowth in non‐wood based industries and pollutants emitted by them in the environment adversely affect the forestry resources. Numerical simulation is provided to support analytical findings.     

MODELS FOR THE EFFECT OF HIGH SPEED WIND ON THE DEPLETION OF FERTILE TOPSOIL

ABSTRACT. In this paper a mathematical model is proposed and analyzed to study the effect of high speed wind on the depletion of fertile topsoil depth and crop yield. The system is modeled by considering a Cobb‐Douglas production function which depends on depreciating capital stock, a labor force and depth of fertile topsoil. A model to conserve the depth of fertile topsoil and to control the high speed wind is also presented.     

Management strategies for the control of eutrophication processes in Fogliano lagoon (Italy): a long-term analysis using a mathematical model

Models are developed and used to analyse and test different management strategies aimed at limiting eutrophication processes in Fogliano Lagoon: modification of lagoon hydrodynamics by tidal flow regulation, harvest of algae biomass, reclaim of sediments. Mathematical models, which have been constructed and proposed, simulate, on a multiyear time scale, the main ecological processes responsible for the most important effects of eutrophication: vegetal blooms, summer anoxia. For different management strategies, hydrodynamic fields produced by wind and tide, and three-dimensional concentration fields of significant species in the ecological phenomena, in water and into sediments, are quantified and compared. The species simulated are: in the water column dissolved oxygen, phytoplanktonic biomass, macrophytic biomass, orthophosphate, dissolved organic carbon, particulate organic carbon and hydrogen sulphide; in sediments dissolved oxygen, dissolved organic carbon, particulate organic carbon, orthophosphate, adsorbed phosphorous and hydrogen sulphide. On the basis of the results of the simulations carried out, the best management strategy limiting eutrophication processes in Fogliano lagoon has been pointed out.     

STOPPING RULES AND THE CONTROL OF STOCK POLLUTANTS

A stock pollutant is defined as a residual waste that might accumulate over time. This paper examines some of the important distinctions between degradable and nondegradable stock pollutants and between nondegradable stock pollutants with known versus uncertain environmental cost. The latter case is examined using the more recent literature on stochastic control with Brownian motion. The presence of irreversibility and uncertainty is known to lead to more conservative investment rules and places a value on the preservation of options. In the case of a nondegradable stock pollutant with Brownian environmental cost, options are preserved by stopping accumulation at a lower level than in the corresponding certainty-equivalent problem. The model presented in this paper permits the derivation of closed-form stopping rules. For a simple numerical problem, the optimal nondegradable stock with Brownian environmental cost was 20 to 45 percent lower than the optimal level with known environmental cost. The empirical study of an actual nondegradable stock pollutant will require time series data on private and social cost in order to estimate drift and variance parameters which will influence the actual extent to which the optimal stock is less than the certainty-equivalent stock.     

On the motion of a heavy body with a circular base on a horizontal plane and riddles of curling

In this paper we discuss the dynamics of an axisymmetric rigid body whose circular area moves upon a horizontal rough surface. We investigate the interaction between the character of the law of friction and the curvature of the body’s trajectory. For the case of a curling stone it is shown that the observed effects can only be explained using the dependence of the friction coefficient on the Gümbel number. The procedure for constructing the law of friction based on experimental data is developed. It is shown that the available data can only be substantiated by means of anisotropic friction. The simplest model of such friction is constructed which provides quantitative coincidence with the experiment.