Direct and Adjoint Oil Spill Estimates

Propagation of the oil spilling from a damaged oil tanker is considered in a limited sea area. The accident consequences are evaluated by means of direct and adjoint oil concentration estimates in ecologically sensitive zones. While the direct estimates are preferable to get a comprehensive idea of the oil spill impact on the whole area, the adjoint ones are useful and economical in studying the sensitivity of the oil concentration in some zones to variations in the accident site and oil spill rate from the tanker. Thanks to special boundary conditions set at the inflow and outflow parts of the open boundary, the main and adjoint oil transport problems are both well-posed according to Hadamard (1923). The estimates obtained in Skiba (1996a) are generalized to the three dimensions. Balanced, absolutely stable 2nd-order finite-difference schemes based on the splitting method are constructed for the two- and three-dimensional cases, both. The main and adjoint schemes are compatible in the sense that at every fractional step of the splitting algorithm, the one-dimensional split operators of both the schemes satisfy a discrete form of the Lagrange identity (Marchuk, 1995). In the special unforced and non-dissipative case, the schemes have two conservation laws each. Each split problem is solved by the factorization method.     

Application of Adjoint Approach to Oil Spill Problems

Adjoint method is applied to various oil spill problems. A three-dimensional model for describing the dispersion of a quasi-passive substance (a pollutant or a nutrient) and its adjoint model are considered in a limited sea region. Direct and adjoint estimates are used to get dual (equivalent) estimates of the mean concentration of the substance in important zones of the region. The role of dual estimates is illustrated with a few examples. They include such oil spill problems as the search of the most dangerous point of the oil tanker route, the oil dispersion with a climatic velocity, and the dependence of the oil concentration estimates on the oil spill rate. One more example is the application of optimal bioremediation strategy for cleaning a few zones polluted by oil. In this case, instead of oil, the model describes the dispersion of a nutrient released to marine environment. Balanced, unconditionally stable second-order finite-difference schemes based on the splitting method for the solution of the dispersion model and its adjoint are suggested. The main and adjoint difference schemes are compatible in the sense that at every fractional step of the splitting algorithm, the one-dimensional split operators of both schemes satisfy a discrete form of Lagrange identity. In the special unforced and non-dissipative case, each scheme has two conservation laws. Every split one-dimensional problem is solved by Thomas’ factorization method.     

A Linear-Programming-Based Strategy for Bioremediation of Oil-Polluted Marine Environments

A linear programming problem is considered with the aim to determine the optimal discharge point and the optimal discharge rate of a nutrient to be released to a marine environment polluted with oil. The objective is to minimize the total discharge of nutrient into the system provided that the concentrations of nutrient will reach critical values sufficient to eliminate oil residuals in certain affected zones through bioremediation. An initial boundary-value 3D problem for the advection–diffusion equation and its adjoint problems are considered to model, estimate, and control the dispersion of nutrient in a limited region. It is shown that the advection–diffusion problem is well posed, and its solution satisfies the mass balance equation. In each oil-polluted zone, the mean concentration of nutrient is determined by means of an integral formula in which the adjoint model solution serves as a weight function. Critical values of these mean concentrations are used as the constraints of linear programming problem. Some additional constraints are posed in order to limit not only the local discharge of the nutrient, but also the mean concentration of this substance in the whole region. Both constraints serve for environmental protection. The ability of the new method is demonstrated by numerical experiments on the remediation in oil-polluted channel using three control zones. The experiments show that the optimal discharge rate can always be got with a simple combination of step functions.     

Industrial pollution transport. Part 1. Formulation of the problem and air pollution estimates

A pollution transport problem is formulated in a limited area. As the pollution sources we take emissions from industrial plants. Physically and mathematically suitable conditions are prescribed on the open boundaries. We show that the problem (as well as its adjoint) is well posed in the sense that a weak solution exists, is unique and depends continuously on its data. Direct and adjoint estimates of the average pollution concentration in an ecologically important zone are given, and the sensitivity of these estimates to perturbations in model parameters is analyzed.     

Environmental genotoxicity and cytotoxicity studies in mussels before and after an oil spill at the marine oil terminal in the Baltic Sea

Environmental genotoxicity and cytotoxicity effects in the gills of mussels Mytilus edulis, from the Baltic Sea areas close to the Būtingė oil terminal (Lithuania) before and after accidental oil spill in 31 January 2008 were studied. Mussels from the oil spillage zones were collected in 12 days, in 3 and 6 months after the spill to determine the effects of the spill. Mussels sampled in 2006–2007 were used for the assessment of the background levels of genotoxicity and cytotoxicity in the Būtingė oil terminal area. Comparison of the responses in M. edulis before and after the oil spill revealed significant elevation of frequencies of micronuclei (MN), nuclear buds (NB) and fragmented-apoptotic (FA) cells. Environmental genotoxicity and cytotoxicity levels in mussels from the Palanga site before the accident (in June 2007) served as a reference. Six months after the accident, in July 2008, 5.6-fold increase of MN, 2.9-fold elevation of NB, and 8.8-fold elevation of FA cells were observed in mussels from the same site.     

Air Pollution Estimates in Guadalajara City

A method for estimating the impact of industrial emissions is suggested and applied to the Guadalajara City Metropolitan Area (GCMA). The method is based on solutions to the pollution transport model and its adjoint. Two equivalent direct and adjoint mean pollution concentration estimates are considered for ecologically important zones of the GCMA. The dependence of these estimates on the number, positions and emission rates of industrial plants, as well as on the wind and initial pollution distribution in the GCMA is qualitatively and quantitatively examined. It is shown that the adjoint model solutions serve as the influence functions providing valuable information on the role of each of the industrial plants in polluting different zones within the GCMA. These solutions have been calculated with a balanced and absolutely stable second-order finite-difference scheme based on the splitting method. A method for an optimal allocation of a new industrial plant is considered.     

Mathematical Modeling and Numerical Algorithms for Simulation of Oil Pollution

This paper deals with the mathematical modeling and algorithms for the problem of oil pollution. For solving this task, we derive the adjoint problem for the advection–diffusion equation describing the propagation of oil slick after an accident, which we call the main problem. We prove a fundamental equality between the solutions of the main and the adjoint problems. Based on this equality, we propose a novel method for the identification of the pollution source location and the accident time of oil emission. This approach is illustrated on an example for an accident in the offshore of the central part of the Vietnamese coast. Numerical simulations demonstrate the effectiveness of the proposed method. Besides, the method is verified for 1D model of substance propagation.     

An assessment of acute biomarker responses in the demersal catfish Cathorops spixii after the Vicuña Oil Spill in a harbour estuarine area in Southern Brazil

The Vicuña oil tanker exploded in Paranaguá Bay (South of Brazil), during methanol unloading operations in front of Paranaguá Harbour, on November 15th, 2004, releasing a large amount of bunker oil and methanol. Two weeks after the accident, the acute effects of the Vicuña Oil Spill (VOS) were evaluated in the demersal catfish Cathorops spixii, comparing a contaminated (at the spill site) and a reference site inside the Bay. Data were compared to previous measurements, taken before the accident, in the same species, from the same sites. The physiological biomarkers were the ones that best reflected acute effects of the spill: plasma osmolality, chloride, calcium, magnesium, and potassium. Morphological (liver and gill histopathology) and genetic (piscine micronucleus and DNA strand breaks) biomarkers revealed that damage was already present in fishes from both reference and contaminated sites inside the Bay. Thus, the reference site is not devoid of contamination, as water circulation tends to spread the contaminants released into other areas of the Bay. Acute field surveys of oil spill effects in harbour areas with a long history of contamination should thus be viewed with caution, and whenever possible previous evaluations should be considered for proper appraisal of biomarker sensitivity, especially in mobile bioindicators such as fish.     

Effects of an Oil Spill on Soil Physico-Chemical Properties of a Spill Site in a Typic Udipsamment of the Niger Delta Basin of Nigeria

Physico-chemical analysis of soil samples at an oil spill site in a Typic Udipsamment of the Niger Delta Basin of Nigeria showed that the total hydrocarbon content of top soil layers ranged from 18.6 to 23.6 ppm in the heavy impact zone and the oil had penetrated to a depth op 8.4 m. The concentration of hydrocarbons in the medium impact zone ranged from 10.04 to 10.38 ppm while hydrocarbons were not detected in 85% of samples from the unimpacted reference zone. Heavy metal concentration measurements in the soil revealed a significant build-up (P < 0.05) of lead, copper and zinc in the heavy impact zone. Other quality parameters including electrical conductity, exchangeable cations, total nitrogen and available phosphorus in impacted soils were relatively low, while the total organic carbon was high compared with the reference site. Textural class of soil from the different depths showed a predominantly brown sand at the topsoil, loamy sand and grey coarse sand at medium depths, and grey coarse sand and greyish sandy clay at greater depths.     

Effects of an oil spill on soil physico-chemical properties of a spill site in the Niger Delta Area of Nigeria

Physico-chemical analysis of soil samples at an oil spill site in the Niger Delta Area of Nigeria showed that the total hydrocarbon content of top soil layers ranged from 0.8 to 12.4 ppm in the heavy impact zone and the oil had penetrated to a depth of 7.2 m. Hydrocarbon concentration in the medium impact zone ranged from 0.02 to 0.40 ppm while hydrocarbons were not detected in 75% of samples from the unimpacted reference zone. Measurement of heavy metal concentrations in the soils revealed a significant build-up (p<0.05) of lead, iron and zinc in the heavy impact zone. Other parameters including electrical conductivity, exchangeable cations, available phosphorus and total nitrogen in impacted soils were comparatively low while the total organic carbon was high, compared with the reference site. Textural class of soil from the different depths showed a predominantly brown sand at the top soil, loamy sand at medium depths and grey coarse sand at greater depths.     

EFFECTS OF AN OIL SPILL ON SOIL PHYSICO-CHEMICAL PROPERTIES OF A SPILL SITE IN A TYPIC PALEUDULT OF MIDWESTERN NIGERIA

Physico-chemical analysis of soil samples at an oil spill sitein a Typic Paleudult of midwestern Nigeria showed that thetotal hydrocarbon content of top soil layers ranged from 0.6to 12.6 ppm in the heavy impact zone and the oil hadpenetrated to a depth of 9.2 m. Hydrocarbon concentration inthe medium impact zone ranged from 0.36 ppm while hydrocarbonswere not detected in 80% of samples from the unimpactedreference zone. Measurement of heavy metal concentrations inthe soil revealed a significant build-up (p< 0.05) of lead,iron and zinc in the heavy impact zone. Other parametersincluding electrical conductivity, exchangeable cations,available phosphorus and total nitrogen in impacted soils werecomparatively low while the total organic carbon was high,compared with the reference site. Textural class of soil fromthe different depths showed a predominantly brown sand at thetopsoil, loamy sand and grey fine sand at medium depths andred sandy clay at greater depths.     

Air Quality Assessment And Control Of Emission Rates

Mathematical methods based on the adjoint model approach are given for the air-pollution estimation and control in an urban region. A simple advection–diffusion-reaction model and its adjoint are used to illustrate the application of the methods. Dual pollution concentration estimates in ecologically important zones are derived and used to develop two non-optimal strategies and one optimal strategy for controlling the emission rates of enterprises. A linear convex combination of these strategies represents a new sufficient strategy. A method for detecting the enterprises, which violate the emission rates prescribed by a control, is given. A method for determining an optimal position for a new enterprise in the region is also described.     

Protecting Sensitive Coastal Areas with Exclusion Booms during Oil Spill Events

Oil spills at sea remain a serious threat to coastal settlements and sensitive ecosystems. Although the impacts of spills are contingent upon a variety of environmental factors and the chemical composition of the oil itself, spill effects can be long lasting in the pelagic zone with broad impacts on sensitive bacterial, microbial, plant, and animal communities. Efforts to contain, deflect, protect, and mitigate the effects of oil are increasingly important, given the massive social, economic, and environmental fallout connected to large spills. The purpose of this paper is to provide geographic perspective for protecting coastal areas with exclusion booms during oil spill events. Specifically, we introduce a generalized, extendable, spatial optimization model that simultaneously minimizes spill effects on vulnerable shorelines and the total costs associated with dispatching booms. The multiobjective model is solved with a weighting method to produce a Pareto optimal curve that reveals how the costs and protection operations change under different priorities. A simulated tanker spill near Mobile Bay, AL, USA, is used as an illustrative example.

     

Industrial pollution transport. Part 2. Control of industrial emissions

Solutions of the pollution transport problem and its adjoint are used to monitor mean pollution concentration in an ecologically important zone. Four strategies of control over pollutants released into the atmosphere by industrial plants are suggested. They differ by the restrictions imposed on the emission rate of each plant. All the strategies use solutions of the adjoint transport problem and assure the fulfillment of the sanitary norm in the zone. A linear interpolation of these strategies also brings pollution level in the zone down to the sanitary norm. A method of detecting the plants violating the prescribed emission rates is also given. A simple example is given to illustrate the strategies suggested.     

用于环境风险评价的大连湾溢油预测模型

基于POM建立了M2分潮的大连湾海域三维潮流场动力模式。利用此动力模式,以"油粒子"模型来模拟溢油在大连湾中的漂移扩散过程。通过欧拉-拉格朗日法,将油膜中心的轨迹和各个时刻的油膜分布情况可视化。溢油预测模型综合考虑了三维流场、湍流、风导输移和波浪破等环境因素对海上溢油行为的影响,考虑了蒸发和乳化过程,预测了溢油的归宿。本模型旨在为大连湾海上溢油风险评价提供信息。     

Groundwater contamination evolution in the Guadiamar and Agrio aquifers after the Aznalcóllar spill: assessment and environmental implications

In 1998, the Agrio and Guadiamar rivers underwent an enormous environmental disaster caused by the rupture of the Aznalcóllar tailings dam and the release of 6 hm³ of pyrite sludge and acidic water. Both rivers run over recent alluvial materials which form a small-sized aquifer which is however important because underground water feeds the flow of the rivers. This work analyzes the state of groundwater 10 years after the spill. Before the dam failure, this aquifer was already contaminated in the zone nearest to the mine, to which the impact of the spill was added. Contamination levels in the alluvial aquifer of the Agrio River have decreased remarkably. However, they are still important, with acidic pH values and high concentrations of toxic elements (maximum values of 16 mg/L of Zn and 15 mg/L of Al). There are also important levels of contamination in the Guadiamar alluvial area closest to the mine, as well as in specific zones located further south. The concentration of toxic elements is mainly controlled by pH. The evolution of contaminant levels show a sharp decrease after the first years following the spill, followed by a subsequent stabilization. It is necessary to take measures for the recovery of the aquifer because, otherwise, groundwater will continue contributing contaminants into the Agrio and Guadiamar rivers.     

Spatial and temporal trends of polycyclic aromatic hydrocarbons in wild mussels from the Cantabrian coast (N Spain) after the Prestige oil spill

Polycyclic aromatic hydrocarbon (PAH) levels were determined in tissues of wild mussels (Mytilus galloprovincialis) collected at 17 stations along the Cantabrian coast (N Spain), from Navia (Asturias) to Fuenterrabía (Basque Country), in order to assess the extent of the environmental impact caused by the Prestige oil spill (November 13, 2002). Six sampling campaigns were carried out in April, June and November in 2003 and 2004. The comparison of PAH data with those obtained earlier in 2000 showed a widespread pyrolytic and petrogenic contamination and allowed an estimation, for the first time, of the background pollution in the region and identification of the chronic hotspots. The spatial distribution found in the first samples after the oil spill revealed the eastern area as the most affected due to the continuous arrival of fuel slicks since early summer 2003. Several stations in this area showed increased total PAH concentrations of up to 15 times the pre-spill levels, which did not recover until April 2004, more than one year after the accident. Molecular parameters within the aliphatic and aromatic fractions were determined to assess the presence of Prestige oil in these samples.     

Simulation of impact of oil spill in the ocean – a case study of Arabian Gulf

To meet the growing energy demand worldwide, oil and gas exploration and production activities have increased rapidly both in onshore and offshore areas. The produced oil from the ocean bed is transported onshore either by ship or pipeline. This has increased the risk of oil spill in the coastal area. In order to prepare an emergency preparedness plan and to assess the magnitude of risk involved in transporting and offloading oil, oil spill simulation studies play an important role. This paper describes a simulation of oil spill in coastal bay of Arabian Gulf where new developments are taking place using MIKE 21 model. The developments include a diesel based thermal power plant near Sir Baniyas Island, which is an ecological fragile area. Based on the project activity, two probable scenarios, one for diesel leak (250 m3/h) for 1 h and the other for instantaneous spill (500 m3) are considered. The MIKE 21 model was calibrated for hydrodynamics using measured field data followed by diesel-spill simulation to track its movement in the Arabian Gulf. The results for both leak and instantaneous spill indicate that spilled diesel will not move towards the Sir Banyas Island and more than 45% of the diesel will be evaporated within 48 h of oil spill. Based on the results, a clean up and contingency plan is proposed to mitigate the adverse impacts arising due to diesel spill in the study area.     

Recovery of a freshwater wetland from chemical contamination after an oil spill

In March 2009, a cargo ship spilled 250 tons of heavy fuel oil off the Queensland coast of Australia. The pristine National Park Moreton Island, seven nautical miles to the east of the spill site, was most affected by the oil slick. Contamination of the island's shoreline was widespread, with freshwater wetlands particularly slow to recover as clean-up needed to be carefully managed to avoid damage to this sensitive ecosystem. During the clean-up process on Moreton Island a monitoring program was initiated using traditional chemical analysis in combination with bioanalytical techniques to assess the extent and variability in contamination at sites on the shoreline and freshwater wetlands. Water accommodated fractions (WAF) of oil residues from samples taken directly after the spill on the shoreline showed the same level of toxic potency as samples from the wetland while baseline-toxicity equivalent concentrations (baseline-TEQ) and 2,3,7,8-tetrachlorodibenzodioxin equivalent concentrations (TCDDEQ) were much lower in oil collected from the sandy beach. The umuC assay for genotoxicity and the E-SCREEN assay for estrogenic effects indicated the extracts were not genotoxic or estrogenic. PAH concentrations and toxicity in grab water samples were below detectable levels, however, extracts from time integrated silicone passive samplers deployed for several weeks at the contaminated sites gave measurable responses in the bioassays with TCDDEQ levels increased relative to the control site. The low levels of baseline-TEQ and TCDDEQ present after 8 months had further decreased 6 months later indicating satisfactory recovery of this pristine ecosystem after an oil spill.     

新疆典型城-郊-乡梯度带土壤微量元素空间分布与来源

从新疆某地典型城-郊-乡梯度带采集了77个表层(0~20 cm)土壤样品,基于GIS技术与多元统计分析方法,研究各梯度带表层土壤中Hg、As、V、Co、Ni、Cu、Zn、Cd、Pb和Sb等10种微量元素的空间分布特征与主要来源。结果表明:Hg元素在城区、郊区和乡村表层土壤中的平均含量均超出研究区土壤背景值,As元素在城区和郊区表层土壤中的平均含量超出背景值,Zn和Pb元素在城区表层土壤中的平均含量超出背景值,其余元素在3个梯度带表层土壤中的平均含量均低于相应的背景值。研究区表层土壤中,V、Co、Ni、Zn、Cd、Pb和Sb等7种元素的空间分布格局基本相似,均呈现沿城市化梯度带分布的地带性格局;As、Cu和Hg等3种元素的空间分布呈现岛状格局。来源分析结果表明,各梯度带表层土壤中的微量元素的来源各不相同。