Linking water-resource models to ecosystem-response models to guide water-resource planning -- an example from the Murray--Darling Basin, Australia

Objectively assessing ecological benefits of competing watering strategies is difficult. We present a framework of coupled models to compare scenarios, using the Coorong, the estuary for the MurrayDarling River system in South Australia, as a case study. The framework links outputs from recent modelling of the effects of climate change on water availability across the MurrayDarling Basin to a hydrodynamic model for the Coorong, and then an ecosystem-response model. The approach has significant advantages, including the following: (1) evaluating management actions is straightforward because of relatively tight coupling between impacts on hydrology and ecology; (2) scenarios of 111 years reveal the impacts of realistic climatic and flow variability on Coorong ecology; and (3) ecological impact is represented in the model by a series of ecosystem states, integrating across many organisms, not just iconic species. We applied the approach to four flow scenarios, comparing conditions without development, current water-use levels, and two predicted future climate scenarios. Simulation produced a range of hydrodynamic conditions and consequent distributions of ecosystem states, allowing managers to compare scenarios. This approach could be used with many climates and/or management actions for optimisation of flow delivery to environmental assets.

Current-oriented swimming by jellyfish and its role in bloom maintenance

Cross-flows (winds or currents) affect animal movements [1-3]. Animals can temporarily be carried off course or permanently carried away from their preferred habitat by drift depending on their own traveling speed in relation to that of the flow [1]. Animals able to only weakly fly or swim will be the most impacted (e.g., [4]). To circumvent this problem, animals must be able to detect the effects of flow on their movements and respond to it [1, 2]. Here, we show that a weakly swimming organism, the jellyfish Rhizostoma octopus, can orientate its movements with respect to currents and that this behavior is key to the maintenance of blooms and essential to reduce the probability of stranding. We combined insitu observations with first-time deployment of accelerometers on free-ranging jellyfish and simulated the behavior observed in wild jellyfish within a high-resolution hydrodynamic model. Our results show that jellyfish can actively swim countercurrent in response to current drift, leading to significant life-history benefits, i.e., increased chance of survival and facilitated bloom formation. Current-oriented swimming may be achieved by jellyfish either directly detecting current shear across their body surface [5] or indirectly assessing drift direction using other cues (e.g., magnetic, infrasound). Our coupled behavioral-hydrodynamic model provides new evidence that current-oriented swimming contributes to jellyfish being able to form aggregations of hundreds to millions of individuals for up to several months, which may have substantial ecosystem and socioeconomic consequences [6, 7]. It also contributes to improve predictions of jellyfish blooms' magnitude and movements in coastal waters. Current drift can have major and potentially negative effects on the lives of weakly swimming species in particular. Fossette etal. show that jellyfish modulate their swimming behavior in relation to current. Such oriented swimming has significant life-history benefits, such as increased bloom formation and a reduction of probability of stranding.

Calibration Strategies Designed For Long Term Sediment Transport Modelling

The presented work deals with the calibration of a 2D numerical model for the simulation of long term bed load transport. A settled basin along an alpine stream was used as a case study. The focus is to parameterise the used multi fractional transport model such that a dynamically balanced behavior regarding erosion and deposition is reached. The used 2D hydrodynamic model utilizes a multi-fraction multi-layer approach to simulate morphological changes and bed load transport. The mass balancing is performed between three layers: a top mixing layer, an intermediate subsurface layer and a bottom layer. Using this approach bears computational limitations in calibration. Due to the high computational demands, the type of calibration strategy is not only crucial for the result, but as well for the time required for calibration. Brute force methods such as Monte Carlo type methods may require a too large number of model runs. All here tested calibration strategies used multiple model runs utilising the parameterization and/or results from previous run. One concept was to reset to initial bed elevations after each run, allowing the resorting process to convert to stable conditions. As an alternative or in combination, the roughness was adapted, based on resulting nodal grading curves, from the previous run. Since the adaptations are a spatial process, the whole model domain is subdivided in homogeneous sections regarding hydraulics and morphological behaviour. For a faster optimization, the adaptation of the parameters is made section wise. Additionally, a systematic variation was done, considering results from previous runs and the interaction between sections. The used approach can be considered as similar to evolutionary type calibration approaches, but using analytical links instead of random parameter changes.

Evolução da biodegradabilidade da matéria orgânica em um sistema de lagoas de estabilização

Waste stabilization ponds (WSP) have been widely used for sewage treatment in hot climate regions because they are economic and environmentally sustainable. In the present study a WSP complex comprising a primary facultative pond (PFP) followed by two maturation ponds (MP-1 and MP-2) was studied, in the city of Natal-RN. The main objective was to study the bio-degradability of organic matter through the determination of the kinetic constant k throughout the system. The work was carried out in two phases. In the first, the variability in BOD, COD and TOC concentrations and an analysis of the relations between these parameters, in the influent raw sewage, pond effluents and in specific areas inside the ponds was studied. In the second stage, the decay rate for organic matter (k) was determined throughout the system based on BOD tests on the influent sewage, pond effluents and water column samples taken from fixed locations within the ponds, using the mathematical methods of Least Squares and the Thomas equation. Subsequently k was estimated as a function of a hydrodynamic model determined from the dispersion number (d), using empirical methods and a Partial Hydrodynamic Evaluation (PHE), obtained from tracer studies in a section of the primary facultative pond corresponding to 10% of its total length. The concentrations of biodegradable organic matter, measured as BOD and COD, gradually reduced through the series of ponds, giving overall removal efficiencies of 71.95% for BOD and of 52.45% for COD. Determining the values for k, in the influent and effluent samples of the ponds using the mathematical method of Least Squares, gave the following values respectively: primary facultative pond (0,23 day-1 and 0,09 day-1), maturation 1 (0,04 day-1 and 0,03 day-1) and maturation 2 (0,03 day-1 and 0,08 day-1). When using the Thomas method, the values of k in the influents and effluents of the ponds were: primary facultative pond (0,17 day-1 and 0,07 day-1), maturation 1 (0,02 day-1 and 0,01 day-1) and maturation 2 (0,01 day-1 and 0,02 day-1). From the Partial Hydrodynamic Evaluation, in the first section of the facultative pond corresponding to 10% of its total length, it can be concluded from the dispersion number obtained of d = 0.04, that the hydraulic regime is one of dispersed flow with a kinetic constant value of 0.20 day-1

Influência dos modelos digitais de terreno na modelagem de circulação hidrodinâmica 2DH: um estudo de caso

Digital Elevation Models (DEM) are numerical representations of a portion of the earth surface. Among several factors which affect the quality of a DEM, it should be emphasized the attention on the input data and the choice of the interpolating algorithm. On the other hand, several numerical models are used nowadays to characterize nearshore hydrodynamics and morphological changes in coastal areas, whose validation is based on field data collection. Independent on the complexity of the physical processes which are modeled, little attention has been given to the intrinsic bathymetric interpolation built within the numerical models of the specific application. Therefore, this study aims to investigate and to quantify the influence of the bathymetry, as obtained by a DEM, on the hydrodynamic circulation model at a coastal stretch, off the coast of the State of Rio Grande do Norte, Northeast Brazil. This coastal region is characterized by strong hydrodynamic and littoral processes, resulting in a very dynamic morphology with shallow coastal bathymetry. Important economic activities, such as oil exploitation and production, fisheries, salt ponds, shrimp farms and tourism, also bring impacts upon the local ecosystems and influence themselves the local hydrodynamics. This fact makes the region one of the most important for the development of the State, but also enhances the possibility of serious environmental accidents. As a hydrodynamic model, SisBaHiA® - Environmental Hydrodynamics System ( Sistema Básico de Hidrodinâmica Ambiental ) was chosen, for it has been successfully employed at several locations along the Brazilian coast. This model was developed at the Coastal and Oceanographical Engineering Group of the Ocean Engineering Program at the Federal University of Rio de Janeiro. Several interpolating methods were tested for the construction of the DEM, namely Natural Neighbor, Kriging, Triangulation with Linear Interpolation, Inverse Distance to a Power, Nearest Neighbor, and Minimum Curvature, all implemented within the software Surfer®. The bathymetry which was used as reference for the DEM was obtained from nautical charts provided by the Brazilian Hydrographic Service of the Brazilian Navy and from a field survey conducted in 2005. Changes in flow velocity and free surface elevation were evaluated under three aspects: a spatial vision along three profiles perpendicular to the coast and one profile longitudinal to the coast as shown; a temporal vision from three central nodes of the grid during 30 days; a hodograph analysis of components of speed in U and V, by different tidal cycles. Small, but negligible, variations in sea surface elevation were identified. However, the differences in flow and direction of velocities were significant, depending on the DEM

Mean-field description of a dynamical collapse of a fermionic condensate in a trapped boson-fermion mixture

We suggest a time-dependent dynamical mean-field-hydrodynamic model for the collapse of a trapped boson-fermion condensate and perform numerical simulation based on it to understand some aspects of the experiment by Modugno et al. [Science 297, 2240 (2002)] on the collapse of the fermionic condensate in the K-40-Rb-87 mixture. We show that the mean-field model explains the formation of a stationary boson-fermion condensate at zero temperature with relative sizes compatible with experiment. This model is also found to yield a faithful representation of the collapse dynamics in qualitative agreement with experiment. In particular we consider the collapse of the fermionic condensate associated with (a) an increase of the number of bosonic atoms as in the experiment and (b) an increase of the attractive boson-fermion interaction using a Feshbach resonance. Suggestion for experiments of fermionic collapse using a Feshbach resonance is made.

Black soliton in a quasi-one-dimensional trapped fermion-fermion mixture

Employing a time dependent mean-field-hydrodynamic model we study the generation of black solitons in a degenerate fermion-fermion mixture in a cigar-shaped geometry using variational and numerical solutions. The black soliton is found to be the first stationary vibrational excitation of the system and is considered to be a nonlinear continuation of the vibrational excitation of the harmonic oscillator state. We illustrate the stationary nature of the black soliton, by studying different perturbations on it after its formation.

Formation of bright solitons and soliton trains in a fermion-fermion mixture by modulational instability

We employ a time- dependent mean- field- hydrodynamic model to study the generation of bright solitons in a degenerate fermion - fermion mixture in a cigar- shaped geometry using variational and numerical methods. Due to a strong Pauli- blocking repulsion among identical spin- polarized fermions at short distances there cannot be bright solitons for repulsive interspecies interactions. Employing a linear stability analysis we demonstrate the formation of stable solitons due to modulational instability of a constant-amplitude solution of the model equations for a sufficiently attractive interspecies interaction. We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton trains by jumping the effective interspecies interaction from repulsive to attractive. These fermionic solitons can be formed and studied in laboratory with present technology.

Fermionic bright soliton in a boson-fermion mixture

We use a time-dependent dynamical mean-field-hydrodynamic model to study the formation of fermionic bright solitons in a trapped degenerate Fermi gas mixed with a Bose-Einstein condensate in a quasi-one-dimensional cigar-shaped geometry. Due to a strong Pauli-blocking repulsion among spin-polarized fermions at short distances there cannot be bright fermionic solitons in the case of repulsive boson-fermion interactions. However, we demonstrate that stable bright fermionic solitons can be formed for a sufficiently attractive boson-fermion interaction in a boson-fermion mixture. We also consider the formation of fermionic solitons in the presence of a periodic axial optical-lattice potential. These solitons can be formed and studied in the laboratory with present technology.

Tightly bound gap solitons in a Fermi gas

Within the framework of the mean-field hydrodynamic model of a degenerate Fermi gas ( DFG), we study, by means of numerical methods and variational approximation ( VA), the formation of fundamental gap solitons ( FGSs) in a DFG ( or in a BCS superfluid generated by weak interaction between spin- up and spin- down fermions), which is trapped in a periodic optical- lattice ( OL) potential. An effectively one- dimensional ( 1D) con. guration is considered, assuming strong transverse confinement; in parallel, a proper 1D model of the DFG ( which amounts to the known quintic equation for the Tonks- Girardeau gas in the OL) is considered too. The FGSs found in the first two bandgaps of the OL- induced spectrum ( unless they are very close to edges of the gaps) feature a ( tightly bound) shape, being essentially confined to a single cell of the OL. In the second bandgap, we also find antisymmetric tightly bound subfundamental solitons ( SFSs), with zero at the midpoint. The SFSs are also confined to a single cell of the OL, but, unlike the FGSs, they are unstable. The predicted solitons, consisting of similar to 10(4) - 10(5) atoms, can be created by available experimental techniques in the DFG of Li-6 atoms.

Bright solitons and soliton trains in a fermion-fermion mixture

We use a time-dependent dynamical mean-field-hydrodynamic model to predict and study bright solitons in a degenerate fermion-fermion mixture in a quasi-one-dimensional cigar-shaped geometry using variational and numerical methods. Due to a strong Pauli-blocking repulsion among identical spin-polarized fermions at short distances there cannot be bright solitons for repulsive interspecies fermion-fermion interactions. However, stable bright solitons can be formed for a sufficiently attractive interspecies interaction. We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton trains. These fermionic solitons can be formed and studied in laboratory with present technology.

Free expansion of fermionic dark solitons in a boson-fermion mixture

We use a time-dependent dynamical mean-field-hydrodynamic model to study the formation of fermionic dark solitons in a trapped degenerate Fermi gas mixed with a Bose-Einstein condensate in a harmonic as well as a periodic optical-lattice potential. The dark soliton with a 'notch' in the probability density with a zero at the minimum is simulated numerically as a nonlinear continuation of the first vibrational excitation of the linear mean-field-hydrodynamic equations, as suggested recently for pure bosons. We study the free expansion of these dark solitons as well as the consequent increase in the size of their central notch and discuss the possibility of experimental observation of the notch after free expansion.

Mixing-demixing transition and collapse of a vortex state in a quasi-two-dimensional boson-fermion mixture

We investigate the mixing-demixing transition and the collapse in a quasi-two-dimensional degenerate boson-fermion mixture (DBFM) with a bosonic vortex. We solve numerically a quantum-hydrodynamic model based on a new density functional which accurately takes into account the dimensional crossover. It is demonstrated that with the increase of interspecies repulsion, a mixed state of DBFM could turn into a demixed state. The system collapses for interspecies attraction above a critical value which depends on the vortex quantum number. For interspecies attraction just below this critical limit there is almost complete mixing of boson and fermion components. Such mixed and demixed states of a DBFM could be experimentally realized by varying an external magnetic field near a boson-fermion Feshbach resonance, which will result in a continuous variation of interspecies interaction.

Miscibility in a degenerate fermionic mixture induced by linear coupling

We consider a one-dimensional mean-field-hydrodynamic model of a two-component degenerate Fermi gas in an external trap, each component representing a spin state of the same atom. We demonstrate that the interconversion between them (linear coupling), imposed by a resonant electromagnetic wave, transforms the immiscible binary gas into a miscible state, if the coupling constant, kappa, exceeds a critical value, kappa(cr). The effect is predicted in a variational approximation, and confirmed by numerical solutions. Unlike the recently studied model of a binary Bose-Einsten condensate with the linear coupling, the components in the immiscible phase of the binary fermion mixture never fill two separated domains with a wall between them, but rather form antilocked (pi-phase-shifted) density waves. Another difference from the bosonic mixture is spontaneous breaking of symmetry between the two components in terms of the numbers of atoms in them, N(1) and N(2). The latter effect is characterized by the parameter nu equivalent to(N(1)-N(2))/(N(1)+N(2)) (only N(1)+N(2) is a conserved quantity), the onset of miscibility at kappa >=kappa(cr) meaning a transition to nu equivalent to 0. At kappa model, with a chemical-potential difference between the two components. The relation between the imbalance in the spin population, induced by the linear coupling, and the developing spatial structure resembles the known Larkin-Ovchinnikov-Fulde-Ferrell states in the Fermi mixture. Dynamical states, when kappa is suddenly switched from zero to a value exceeding kappa(cr), are considered too. In the latter case, the system features oscillatory relaxation to the mixed state.

Study of a degenerate dipolar Fermi gas of Dy-161 atoms

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)