Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

Computer simulations of seasonal outbreak and diurnal vertical migration of cyanobacteria

Algal blooms caused by cyanobacteria are characterized by two features with different time scales: one is seasonal outbreak and collapse of a bloom and the other is diurnal vertical migration. Our two-component mathematical model can simulate both phenomena, in which the state variables are nutrients and cyanobacteria. The model is a set of one-dimensional reaction-advection-diffusion equations, and temporal changes of these two variables are regulated by the following five factors: (1) annual variation of light intensity, (2) diurnal variation of light intensity, (3) annual variation of water temperature, (4) thermal stratification within a water column and (5) the buoyancy regulation mechanism. The seasonal change of cyanobacteria biomass is mainly controlled by factors, (1), (3) and (4), among which annual variations of light intensity and water temperature directly affect the maximum growth rate of cyanobacteria. The latter also contributes to formation of the thermocline during the summer season. Thermal stratification causes a reduction in vertical diffusion and largely prevents mixing of both nutrients and cyanobacteria between the epilimnion and the hypolimnion. Meanwhile, the other two factors, (2) and (5), play a significant role in diurnal vertical migration of cyanobacteria. A key mechanism of vertical migration is buoyancy regulation due to gas-vesicle synthesis and ballast formation, by which a quick reversal between floating and sinking becomes possible within a water column. The mechanism of bloom formation controlled by these five factors is integrated into the one-dimensional model consisting of two reaction-advection-diffusion equations.

Determining the Reaeration Coefficient and Hydrodynamic Properties of Rivers Using Inert Gas Tracers

Various contaminants which can be aerobically degraded find their way directly or indirectly into surface water bodies. The reaeration coefficient (K2) characterises the rate at which oxygen can transfer from the atmosphere across the air-water interface following oxygen depletion in a water body. Other mechanisms (like advection, dispersion and transient storage) determine how quickly the contaminants can spread in the water, affecting their spatial and temporal concentrations. Tracer methods involving injection of a gas into the water body have traditionally been used for direct (in-situ) measurement of K2 in a given reach. This paper shows how additional modelling of tracer test results can be used to quantify also hydrodynamic mechanisms (e.g. dispersion and storage exchange coefficients, etc.). Data from three tracer tests conducted in the River Lagan (Northern Ireland) using an inert gas (krypton, Kr) are re-analysed using two solute transport models (ADM, TSM) and an inverse-modelling framework (OTIS-P). Results for K2 are consistent with previously published values for this reach (K2(20)~10-40 d-1). The storage area constituted 30-60% of the main cross-section area and the storage exchange rate was between 2.5×10-3-3.2×10-3s-1. The additional hydrodynamic parameters obtained give insight into transport and dispersion mechanisms within the reach.

Dynamics of the Indonesian seas circulation. Part I – The influence of bottom topography on temperature and salinity distributions

The influence of bottom topography on the distribution of temperature and salinity in the Indonesian seas region has been studied with a high-resolution model based on the Princeton Ocean Model. One of the distinctive properties of the model is an adequate reproduction of all major topographic features in the region by the model bottom relief. The three major routes of flow of Pacific water through the region have been identified. The western route follows the flow of North Pacific Water through the Sulawesi Sea, Makassar Strait, Flores Sea, and Banda Sea. This is the main branch of the Indonesian Throughflow. The eastern routes follow the flow of South Pacific water through the eastern Indonesian seas. This water enters the region either through the Halmahera Sea or by flowing to the north around Halmahera Island into the Morotai Basin and then into the Maluku Sea. A deep southward flow of South Pacific Water fills the Seram Sea below 1200 m through the Lifamatola Passage. As it enters the Seram Sea, this overflow turns eastward at depths greater than 2000 m, then upwells in the eastern part of the Seram Sea before returning westward at ~1500-2000 m. The flow continues westward across the Seram Sea, spreading to greater depths before entering the Banda Sea at the Buru-Mangole passage. It is this water that shapes the temperature and salinity of the deep Banda Sea. Topographic elevations break the Indonesian seas region down into separate basins. The difference in the distributions of potential temperature, ?, and salinity, S, in adjacent basins is primarily due to specific properties of advection of ? and S across a topographic rise. By and large, the topographic rise blocks deep flow between basins whereas water shallower than the depth of the rise is free to flow between basins. To understand this process, the structure of simulated fields of temperature and salinity has been analyzed. To identify a range of advected ? or S, special sections over the sills with isotherms or isohalines and isotachs of normal velocity have been considered. Following this approach the impact of various topographic rises on the distribution of ? and S has been identified. There are no substantial structural changes of potential temperature and salinity distributions between seasons, though values of some parameters of temperature and salinity distributions, e.g., magnitudes of maxima and minima, can change. It is shown that the main structure of the observed distributions of temperature and salinity is satisfactorily reproduced by the model throughout the entire domain.

STEADY-STATE SOLUTIONS IN NONLINEAR DIFFUSIVE SHOCK ACCELERATION

Stationary solutions to the equations of nonlinear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean free path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cutoff in energy space-particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for the substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven nonresonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off-a diagnostic that is in principle accessible to observation.

Ecological and societal benefits of jellyfish blooms

Jellyfish are often considered as stressors on marine ecosystems or as indicators of highly perturbed systems. Far less attention is given to the potential of such species to provide beneficial ecosystem services in their own right. In an attempt to redress this imbalance we take the liberty of portraying jellyfish in a positive light and suggest that the story is not entirely one of doom and gloom. More specifically, we outline how gelatinous marine species contribute to the four categories of ecosystem services (regulating, supporting, provisioning and cultural) defined by the Millennium Ecosystem Assessment. This discussion ranges from the role of jellyfish in carbon capture and advection to the deep ocean through to the creation of micro habitat for developing fishes and the advancement of citizen science programmes. Attention is paid also to incorporation of gelatinous species into fisheries or ecosystem level models and the mechanisms by which we can improve the transfer of information between jellyfish researchers and the wider non-specialist community.

Assessing the environmental fate of selected polybrominated diphenyl ethers in the region surrounding the Zhuoshui River of Taiwan based on an Equilibrium Constant fugacity model

Polybrominated diphenyl ethers (PBDEs) are a group of flame retardants that have been in use since the 1970s. They are included in the list of hazardous substances known as persistent organic pollutants (POPs) because they are extremely hazardous to the environment and human health. PBDEs have been extensively used in industry and manufacturing in Taiwan, thus its citizens are at high risk of exposure to these chemicals.
An assessment of the environmental fate of these compounds in the Zhuoshui river and Changhua County regions of western Taiwan, and also including the adjacent area of the Taiwan Strait, was conducted for three high risk congeners, BDE-47, -99 and -209, to obtain information regarding the partitioning, advection, transfer and long range transport potential of the PBDEs in order to identify the level of risk posed by the pollutants in this region.
The results indicate that large amounts of PBDEs presently reside in all model compartments – air, soil, water, and sediment – with particularly high levels found in air and especially in sediment. The high levels found in sediment, particularly for BDE-209, are significant, since there is the threat of these pollutants entering the food chain, either directly through benthic feeding, or through resuspension and subsequent feeding in the pelagic region of the water column which is a distinct possibility in the strong currents found within the Taiwan Strait. Another important result is that a substantial portion of emissions leave the model domain directly through advection, particularly for BDE-47 (58%) and BDE-209 (75%), thus posing a risk to adjacent communities.
Model results were generally in reasonable agreement with available measured concentrations. In air, model concentrations are in reasonably good agreement with available measured values. For both BDE-47 and -99, model concentrations are a factor of 2-3 higher and BDE-209 within the range of measured values. In soil, model results are somewhat less than measured values. In sediment, model results are at the high end of measured values.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.