The kinetics of phytoplankton is complicated function of phytoplankton present and their differing reactions to solar radiation, temperature and balance between nutrient availabilities and phytoplankton requirements. Computers program for water quality problems helps users interpret and predict water quality responses to natural phenomena. But model parameters and initial condition needs hydrodynamic and hydrologic measures in long time interval, with combination of laboratory studies. Lake Gracanka (Serbia, Yugoslavia) is good example for description of the processes kinetics of phytoplankton. But, simulation model can not be used to guide future monitoring efforts without needing measuring and developing network of measuring station.

The paper presents one of the methods for calculation of temperature of a water course, downstream of the relase of a thermo-electric power-plant, taking into consideration the well-known temperatures and discharge values of the tributaries under different meteorological conditions. In view of the importance of temperature regime for the water eco-system, understanding the influence of the power plants on thermic capacity of the water course, is essential in design and use of the power plant. The author's intention is to help towards improving computer models for the computation of a water course thermic capacity. To this end a model developed at the Institute for Meteorology and Water Management in Poland and applied in many coutries has been presented.

Mathematical modles and adequate investigation works enables successful water quality modeling. These models includes all relevant abiotic components (concentration of dissolved oxygen, temperature regimes and temperature stratification, concentration of organic and inorganic components of physical-chemical water conditions), while the biotic components are present only with phyto- and zooplankton. In new water reservoir, ichthyofauna mostly undergo mancontrolled successions, starting with planned stocking a water reservoir with fish. The paper presents analytical approach for mathematical modeling of ichthyofauna and application of Lotka - Volterr's equations.

MPI - Metal Pollution Index is being proposed with an aim improving the freshwater pollution control, monitoring and classification based on fish metal body burden. This is a simple mathematical model which enables presentation of Cd, Pb, Cu, Zn and Al fish liver concentration as a single value. After thorough studies of bioaccumulation, biomagnification and tissue distribution patterns, liver of Carassius auratus gibelio has been chosen as sentinel organ. Normal distribution and biological variation assessment is being achieved by using normalizers - reference values (liver concentrations of selected elements in fish from reference sites Vruje and Moharac micro water reservoirs, Serbia (Yugoslavia)). Selected representative sites have been evaluated according to new classification of freshwater ecosystem based on MPI values.

Stable isotope fractionation of toluene under dynamic phase exchange was studied aiming at ascertaining the effects of gas-liquid partitioning and biodegradation of toluene stable isotope composition in liquid-air phase exchange reactors (Laper). The liquid phase consisted of a mixture of aqueous minimal media, a known amount of a mixture of deuterated (toluene-d) and non-deuterated toluene (toluene-h), and bacteria of toluene degrading strain Pseudomonas putida KT2442. During biodegradation experiments, the liquid and air-phase concentrations of both toluene isotopologues were monitored to determine the observable stable isotope fractionation in each phase. The results show a strong fractionation in both phases with apparent enrichment factors beyond −800‰. An offset was observed between enrichment factors in the liquid and the gas phase with gas-phase values showing a stronger fractionation in the gas than in the liquid phase. Numerical simulation and parameter fitting routine was used to challenge hypotheses to explain the unexpected experimental data. The numerical results showed that either a very strong, yet unlikely, fractionation of the phase exchange process or a – so far unreported – direct consumption of gas phase compounds by aqueous phase microorganisms could explain the observed fractionation effects. The observed effect can be of relevance for the analysis of volatile contaminant biodegradation using stable isotope analysis in unsaturated subsurface compartments or other environmental compartment containing a gas and a liquid phase.

The effectiveness and feasibility of the three biochar materials for remediation of arsenic (As) and lead (Pb) contaminated soil were explored in this study. Significant reduction of bioaccessibility and migration risks of both heavy metals have been explained mechanistically by incubation, column experiments and numerical simulation. Langmuir equation fitted As and Pb sorption isotherms better in the control and biochar (BC) amended soils, while Freundlich model was more suitable for iron modified biochar (Fe-BC) and sulfur/iron modified biochar (S/Fe-BC) amended soils, indicating that modified biochar promoted chemical adsorption process for As and Pb. For the three biochar materials, S/Fe-BC showed the best effects on reducing the bioavailability of As and Pb, with a decrease of 40.42%–64.21%. The reduction in bioaccessibility by metal portioning into available and non-available fractions was better for illustrating the mechanisms including adsorption, precipitation/coprecipitation and As(III) oxidation behind S/Fe-BC efficacy. Moreover, S/Fe-BC can effectively inhibit the leaching behavior of As and Pb under acid rain, which increased by 99.89% and 90.18%, respectively, compared with the control. The HYDRUS-1D modeling indicated that S/Fe-BC could continuously treat As (100 mg/L) and Pb (1000 mg/L) contaminated water for 16.22 years and 40.86 years, respectively, and ensure the groundwater quality criteria being met. Based on these insights, we believe that our study will provide meaningful information about the potentials of biochar derived materials for soil heavy metals’ remediation.

The influence of regional transport on aerosol pollution has been explored in previous studies based on numerical simulation or surface observation. Nevertheless, owing to inhomogeneous vertical distribution of air pollutants, vertical observations should be conducted for a comprehensive understanding of regional transport. Here we obtained the vertical profiles of aerosol and its precursors using ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) at the Nancheng site in suburban Beijing on the southwest transport pathway of the Beijing-Tianjin-Hebei (BTH) region, China, and then estimated the vertical profiles of transport fluxes in the southwest-northeast direction. The maximum net transport fluxes per unit cross-sectional area, calculated as pollutant concentration multiply by wind speed, of aerosol extinction coefficient (AEC), NO₂, SO₂ and HCHO were 0.98 km⁻¹ m s⁻¹, 24, 14 and 8.0 μg m⁻² s⁻¹ from southwest to northeast, which occurred in the 200–300 m, 100–200 m, 500–600 m and 500–600 m layers, respectively, due to much higher pollutant concentrations during southwest transport than during northeast transport in these layers. The average net column transport fluxes were 1200 km⁻¹ m² s⁻¹, 38, 26 and 15 mg m⁻¹ s⁻¹ from southwest to northeast for AEC, NO₂, SO₂ and HCHO, respectively, in which the fluxes in the surface layer (0–100 m) accounted for only 2.3%–4.2%. Evaluation only based on surface observation would underestimate the influence of the transport from southwest cities to Beijing. Northeast or weak southwest transports dominated in clean conditions with PM₂.₅ <75 μg m⁻³ and intense southwest transport dominated in polluted conditions with PM₂.₅ >75 μg m⁻³. Southwest transport through the middle boundary layer was a trigger factor for aerosol pollution events in urban Beijing, because it not only directly bringing air pollutants, but also induced an inverse structure of aerosols, which resulted in stronger atmospheric stability and aggravated air pollution in urban Beijing.

Oil spills are a serious environmental problem. To better support risk assessment and pollution control for oil spills, a good understanding of oil transport in the environment is required. This study focused on the numerical simulation of the nearshore oil behaviors based on computational fluid dynamics. Based on the Navier-Stokes momentum equations for an incompressible viscous fluid and volume of fluid (VOF) method, a 3D numerical model of three-phase transient flow was developed. The wave number, averaged flow velocity and oil properties would affect the oil spread extent and the oil volume fraction. The higher the averaged flow velocity and wave number, the lower the oil concentration, and the faster the horizontal movement of the oil. The spilled oil may move to contact the seafloor by increasing the averaged flow velocity at the inlet boundary. Through increasing the wave number, the oil would stay near the water surface. In the nearshore, where the wave is the main seawater motion, the oil containment boom should be set preferentially to the direction of wave transmission for oil cleaning. This study shows that by doubling the wave number and increasing the averaged flow velocity (ten times) at the same time, the maximum oil volume fraction would be reduced by around 32%. Finally, the water temperature had no significant impact on the oil migration, and the impact of evaporation should be considered in the simulation.

The growing use of coastal areas for different economic purposes is responsible for increasing pollution by hydrocarbons in marine environments. As a consequence of these activities, accidents during fuel extraction, transport, and storage can occur, causing intense environmental degradation. Numerical modeling of the trajectory of oil stains becomes an important tool with low operational costs, providing powerful support to the government agencies in charge of risk management associated with possible oil accidents, by helping to generate scenarios and strategies for containment and cleaning of affected environments. In this sense, the aim of this study is to estimate environmental vulnerability to oil at beaches located in the Itapuã State Park (PEI), a Protection Conservation Unit. This work focused on describing a methodology to estimate the vulnerability of coastal areas, with emphasis on the fact that the study was carried out in a closed environment. For that, an approach was used based on the integration of: (1) an intrinsic variable to the environment; (2) a dynamic variable determined through diesel oil surface dispersion scenarios. Four hypothetical accident scenarios with 20 m³ of diesel oil were simulated in 2018, during five days of simulations with instant dumping in the navigation channel of the local waterway near the PEI. The results suggest the forcing of the field of intensity and direction of the local winds as preponderant for the dynamics of movement and structure of the spots, with the zonal and meridional components of the fields of superficial currents acting in this process as a secondary factor. The study showed that all beaches in the park are susceptible to contact with oil throughout the simulated year, with Pombas Beach, Pedreira Beach, and Onça Beach being affected in all simulated scenarios, which classifies them as very high vulnerability and defines them as priority protection areas.