A model of fish preference and mortality under hypoxic water in the coastal environment

Eutrophication and associated occurrence of hypoxic condition could cause significant damage to marine ecosystems, resulting in considerable economic losses to fisheries and aquaculture and is a major source of stress that fish often have to contend with in order to survive. This problem is likely to be exacerbated in the coming years, since the wastewater treatment facilities is unlikely to catch up with increasing human activities. Moreover, large-scale reclamation projects in coastal areas have recently been increased, and these activities certainly have adverse impacts on water quality and fisheries resources. Coastal construction has a significant role in the development of hypoxic water by changing the current and mixing pattern of water. Changes in species composition and decreases in species richness and diversity have been well documented in hypoxic systems. Hypoxia could cause endocrine disruption in fish and eliminate populations of sensitive species. Shallow coastal areas are of great importance for the special nursery of fish and shellfish and land reclamation in these areas cause strong damage to fisheries. Although the tolerance of aquatic life to hypoxia is known, there is no information about the mortality of fish caused by hypoxia because fish can swim around it and no modeling study has yet been carried out. Criteria that influence the movement of fish are: amounts of food, water temperature and depth, dissolved oxygen concentration and nature of seabed. However, among these, water temperature and dissolved oxygen are the most crucial parameters that affect survival, movement and growth of fish. In this paper, a model of fish preference and mortality for environmental conditions was developed and applied to the Hakata Bay where hypoxic water occurs every summer. For the purpose of this study, a field survey of fish behavior under hypoxic water was conducted by releasing marbled sale (Pleuronectes yokohamae) in the inner bay. Moreover, a series of preference tests for DO, salinity and temperature in the laboratory were conducted in order to decide preference parameters of fish. Using the results of both field and laboratory studies, a sub-module of fish preference and mortality was coded within an integrated hydrothermal and eutrophication model (CHEM) to predict the behavior and mortality of marbled sale when hypoxia would occur in the bay. The model could reasonably simulate the behavior of the fish under hypoxia. An assessment of the impact of the ongoing land reclamation project of about 401 ha in the Wajiro tidal flat zones at the head of the bay on the fisheries resources was also conducted using the model. The results showed that the artificial land lowered the mortality rate of fish under hypoxic condition in the bay during the summer period.     

Pits,rifts and slumps: the summit structure of Piton de la Fournaise

A clear model of structures and associated stress fields of a volcano can provide a framework in which to study and monitor activity. We propose a volcano-tectonic model for the dynamics of the summit of Piton de la Fournaise (La Reunion Island, Indian Ocean). The summit contains two main pit crater structures (Dolomieu and Bory), two active rift zones, and a slumping eastern sector, all of which contribute to the actual fracture system. Dolomieu has developed over 100 years by sudden large collapse events and subsequent smaller drops that include terrace formation. Small intra-pit collapse scars and eruptive fissures are located along the southern floor of Dolomieu. The western pit wall of Dolomieu has a superficial inward dipping normal fault boundary connected to a deeper ring fault system. Outside Dolomieu, an oval extension zone containing sub-parallel pit-related fractures extends to a maximum distance of 225 m from the pit. At the summit the main trend for eruptive fissures is N80°, normal to the north–south rift zone. The terraced structure of Dolomieu has been reproduced by analogue models with a roof to width ratio of approximately 1, suggesting an original magma chamber depth of about 1 km. Such a chamber may continue to act as a storage location today. The east flank has a convex–concave profile and is bounded by strike-slip fractures that define a gravity slump. This zone is bound to the north by strike-slip fractures that may delineate a shear zone. The southern reciprocal shear zone is probably marked by an alignment of large scoria cones and is hidden by recent aa lavas. The slump head intersects Dolomieu pit and may slide on a hydrothermally altered layer known to be located at a depth of around 300 m. Our model has the summit activity controlled by the pit crater collapse structure, not the rifts. The rifts become important on the mid-flanks of the cone, away from pit-related fractures. On the east flank the superficial structures are controlled by the slump. We suggest that during pit subsidence intra-pit eruptions may occur. During tumescence, however, the pit system may become blocked and a flank eruption is more likely. Intrusions along the rift may cause deformation that subsequently increases the slump’s potential to deform. Conversely, slumping may influence the east flank stress distribution and locally control intrusion direction. These predictions can be tested with monitoring data to validate the model and, eventually, improve monitoring.     

Impact of adsorbed petroleum hydrocarbons on marine organisms

Adsorptions of south Louisiana crude oil from seawater by clays, non-clay minerals and sediments were conducted in the laboratory. Effect of sediment-adsorbed and water dispersed crude oil on adult oysters were investigated in aquaria. Hydrocarbons in oyster tissues and surrounding water were identified by gas chromatography. Field specimens from an area of a new oil spill and an area five months after an oil spill were also analysed. Evidence of secondary chemical effect of aged oil in sediments on oyster mortality is presented. Oil coated asbestos surfaces severely reduced recruitment of sedentary larval organisms.     

Testing the suitability of frictional behaviour for pyroclastic flow simulation by comparison with a well-constrained eruption at Tungurahua volcano (Ecuador)

We use a well-monitored eruption of Tungurahua volcano to test the validity of the frictional behaviour, also called Mohr–Coulomb, which is generally used in geophysical flow modelling. We show that the frictional law is not appropriate for the simulation of pyroclastic flows at Tungurahua. With this law, the longitudinal shape of the simulated flows is a thin wedge of material progressively passing, over several hundreds of metres, from an unrealistic thickness at the front (<<1 mm) to some tens of centimetres. Simulated deposits form piles which accumulate at the foot of the volcano and are more similar to sand piles than natural pyroclastic deposits. Finally, flows simulated with a frictional rheology are not channelised by the drainage system, but affect all the flanks of the volcano. In addition, their velocity can exceed 150 m s⁻¹, allowing pyroclastic flows to cross interfluves at bends in the valley, affecting areas that would not have been affected in reality and leaving clear downstream areas that would be covered in reality. Instead, a simple empirical law, a constant retarding stress (i.e. a yield strength), involving only one free parameter, appears to be much better adapted for modelling pyroclastic flows. A similar conclusion was drawn for the Socompa debris avalanche simulation (Kelfoun and Druitt, J Geophys Res 110:B12202, 2005).     

Groundwater recharge in irrigated semi-arid areas: quantitative hydrological modelling and sensitivity analysis

For semi-arid regions, methods of assessing aquifer recharge usually consider the potential evapotranspiration. Actual evapotranspiration rates can be below potential rates for long periods of time, even in irrigated systems. Accurate estimations of aquifer recharge in semi-arid areas under irrigated agriculture are essential for sustainable water-resources management. A method to estimate aquifer recharge from irrigated farmland has been tested. The water-balance-modelling approach was based on VisualBALAN v. 2.0, a computer code that simulates water balance in the soil, vadose zone and aquifer. The study was carried out in the Campo de Cartagena (SE Spain) in the period 1999–2008 for three different groups of crops: annual row crops (lettuce and melon), perennial vegetables (artichoke) and fruit trees (citrus). Computed mean-annual-recharge values (from irrigation+precipitation) during the study period were 397 mm for annual row crops, 201 mm for perennial vegetables and 194 mm for fruit trees: 31.4, 20.7 and 20.5% of the total applied water, respectively. The effects of rainfall events on the final recharge were clearly observed, due to the continuously high water content in soil which facilitated the infiltration process. A sensitivity analysis to assess the reliability and uncertainty of recharge estimations was carried out.     

Biostratigraphy of oligocene succession in the High Folded Zone,Sulaimani, Kurdistan region,Northeastern Iraq

Oligocene and Lower Miocene sediments from High Folded Zone of Iraqi Zagros have been studied paleontologically at south of Sulaimaniyah, Kurdistan Region, NE Iraq. The identified fauna are consisted of (25) genera and species of benthonic and (16) species of planktonic foraminifera. The fauna comprises relatively abundant foraminiferal assemblages of moderate diversity. Based on the stratigraphic distribution of these species, two biozones have been recognized which are Nummulites—Rotalia and Globoquadrina dehiscens zones. These biozones indicate that the studied sections of Basara and Khewata are of Late Oligocene–Early Miocene age. Based on the microfossils, it has been found that the age of sediments is equivalent to or represents Anah and Serikagni Formations. Some previous studies described Oligocene rocks (Kirkuk Group) as interior sag basin. In the present study, the occurrence of the group inside High Folded Zone and its rich fauna content are used for the discussion of the sag basin versus normal marine water.     

Holocene rockfalls in the southern Negev Desert,Israel and their relation to Dead Sea fault earthquakes

Rockfall ages in tectonically active regions provide information regarding frequency and magnitude of earthquakes. In the hyper-arid environment of the Dead Sea fault (DSF), southern Israel, rockfalls are most probably triggered by earthquakes. We dated rockfalls along the western margin of the DSF using terrestrial cosmogenic nuclides (TCN). At each rockfall site, samples were collected from simultaneously exposed conjugate boulders and cliff surfaces. Such conjugate samples initially had identical pre-fall (“inherited”) TCN concentrations. After boulder detachment, these surfaces were dosed by different production rates due to differences in post-fall shielding and geometry. However, in our study area, pre-rockfall inheritance and post-rockfall production rates of TCN cannot be evaluated. Therefore, we developed a numerical approach and demonstrated a way to overcome the above-mentioned problems. This approach can be applied in other settings where rockfalls cannot be dated by simple exposure dating. Results suggest rockfall ages between 3.6 ± 0.8 and 4.7 ± 0.7 ka. OSL ages of sediment accumulated behind the boulders range between 0.6 ± 0.1 and 3.4 ± 1.4 ka and support the TCN results. Our ages agree with dated earthquakes determined in paleoseismic studies along the entire length of the DSF and support the observation of intensive earthquake activity around 4–5 ka.     

Pile Head Cyclic Lateral Loading of Single Pile

This paper presents an elastic continuum model using an extended nonlinear Davies and Budhu equations, which enables the nonlinear behavior of the soil around the long elastic pile to be modeled using a simple expression of pile-head stiffness method. The calculated results were validated with the measured full-scale dynamic field tests data conducted in Auckland residual clay. An idealized soil profile and soil stiffness under small strain (i.e. shear modulus, G _s and shear wave velocity, V _s of the soil) determined from in situ testing was used to model the single pile tests results. The predictions of these extended equations are also confirmed by using the three-dimensional finite-element OpenSeesPL (Lu et al. in OpenSeesPL 3D lateral pile-ground interaction: user manual, University of California, San Diego, 2010). A soil stiffness reduction factor, G _s /G _s,max of 0.36 was introduced to the proposed method and model. It was found to give a reasonable prediction for a single pile subjected to dynamic lateral loading. The reduction in soil stiffness found from the experiment arises from the cumulative effects of pile–soil separation as well as a change in the soil properties subjected to cyclic load. In summary, if the proposed method and model are accurately verified and properly used, then they are capable of producing realistic predictions. Both models provide good modelling tools to replicate the full-scale dynamic test results.     

Minimizing dredging disposal via sediment management in New York Harbor

The Port of New York/New Jersey is naturally shallow, and therefore dredging is required to maintain depths necessary for navigation. About six million cubic yards of material must be dredged annually to maintain navigation channels and berthing areas. Opportunities for disposal of dredged materials in the metropolitan region are limited. The existing ocean disposal site that has, until recently, received the majority of dredged materials is nearing capacity. Under the new, more stringent, guidance for ocean disposal, a large percentage of the dredged material is not considered acceptable for ocean disposal because of contamination. This paper explores nondredging alternatives to reduce the volume of materials to be dredged, thus reducing the disposal volume. These alternatives include short-term options (e.g., reprofiling operations, in which sediments from high spots in berths are dragged to depressions in lower spots) and long-term sedimentation minimization options (e.g., subsurface berms or air bubblers deployed around berths and interpier areas). These methods require an understanding of sedimentation/circulation patterns of the harbor as a whole and the local area in which the specific methodology may be employed.     

Mapping landslide susceptibility with frequency ratio,statistical index,and weights of evidence models: a case study in northern Iran

The main objective of this study is to investigate potential application of frequency ratio (FR), weights of evidence (WoE), and statistical index (SI) models for landslide susceptibility mapping in a part of Mazandaran Province, Iran. First, a landslide inventory map was constructed from various sources. The landslide inventory map was then randomly divided in a ratio of 70/30 for training and validation of the models, respectively. Second, 13 landslide conditioning factors including slope degree, slope aspect, altitude, plan curvature, stream power index, topographic wetness index, sediment transport index, topographic roughness index, lithology, distance from streams, faults, roads, and land use type were prepared, and the relationships between these factors and the landslide inventory map were extracted by using the mentioned models. Subsequently, the multi-class weighted factors were used to generate landslide susceptibility maps. Finally, the susceptibility maps were verified and compared using several methods including receiver operating characteristic curve with the areas under the curve (AUC), landslide density, and spatially agreed area analyses. The success rate curve showed that the AUC for FR, WoE, and SI models was 81.51, 79.43, and 81.27, respectively. The prediction rate curve demonstrated that the AUC achieved by the three models was 80.44, 77.94, and 79.55, respectively. Although the sensitivity analysis using the FR model revealed that the modeling process was sensitive to input factors, the accuracy results suggest that the three models used in this study can be effective approaches for landslide susceptibility mapping in Mazandaran Province, and the resultant susceptibility maps are trustworthy for hazard mitigation strategies.