Bed composition generation for morphodynamic modeling: case study of San Pablo Bay in California, USA

Applications of process-based morphodynamic models are often constrained by limited availability of data on bed composition, which may have a considerable impact on the modeled morphodynamic development. One may even distinguish a period of “morphodynamic spin-up??in which the model generates the bed level according to some ill-defined initial bed composition rather than describing the realistic behavior of the system. The present paper proposes a methodology to generate bed composition of multiple sand and/or mud fractions that can act as the initial condition for the process-based numerical model Delft3D. The bed composition generation (BCG) run does not include bed level changes, but does permit the redistribution of multiple sediment fractions over the modeled domain. The model applies the concept of an active layer that may differ in sediment composition above an underlayer with fixed composition. In the case of a BCG run, the bed level is kept constant, whereas the bed composition can change. The approach is applied to San Pablo Bay in California, USA. Model results show that the BCG run reallocates sand and mud fractions over the model domain. Initially, a major sediment reallocation takes place, but development rates decrease in the longer term. Runs that take the outcome of a BCG run as a starting point lead to more gradual morphodynamic development. Sensitivity analysis shows the impact of variations in the morphological factor, the active layer thickness, and wind waves. An important but difficult to characterize criterion for a successful application of a BCG run is that it should not lead to a bed composition that fixes the bed so that it dominates the “natural??morphodynamic development of the system. Future research will focus on a decadal morphodynamic hindcast and comparison with measured bathymetries in San Pablo Bay so that the proposed methodology can be tested and optimized.     

Electrochemical Treatment of Dye Solution by Oxalate Catalyzed Photoelectro‐Fenton Process Using a Carbon Nanotube‐PTFE Cathode: Optimization by Central Composite Design

Decolorization of C.I. Basic Blue 3 (BB3) by oxalate catalyzed photoelectro‐Fenton process based on carbon nanotube‐polytetrafluoroethylene (CNT‐PTFE) electrode as cathode under visible light was studied. A comparison of electro‐Fenton, photoelectro‐Fenton, and photoelectro‐Fenton/oxalate processes for decolorization of the solution containing BB3 has been performed. The results showed that color removal follows the decreasing order: photoelectro‐Fenton/oxalate > photoelectro‐Fenton > electro‐Fenton. Response surface methodology (RSM) was employed to assess individual and interactive effects of the four main independent parameters on the decolorization efficiency. A central composite design (CCD) was employed for optimization of photoelectro‐Fenton/oxalate treatment of BB3. The analysis of variance (ANOVA) showed a high coefficient of determination value (R² = 0.958) and satisfactory prediction second‐order regression. This study clearly showed that RSM was one of the suitable methods to optimize the operating conditions.     

Cr(VI) Immobilization Process in a Cr‐Spiked Soil by Zerovalent Iron Nanoparticles: Optimization Using Response Surface Methodology

The response surface methodology involving the five‐level central composite design (CCD) was employed to model and optimize the Cr(VI) immobilization process in a Cr‐spiked soil using starch‐stabilized zerovalent iron nanoparticles (ZVIn). ZVIn were synthesized via a borohydride reduction method and characterized by X‐ray diffraction (XRD) and scanning electron microscopy (SEM). All Cr(VI) immobilization experiments were conducted in a batch system. The variables for the CCD optimization were the ZVIn dosage (% w/w), reaction time (min), and initial Cr(VI) concentration in soil (mg/kg). The predicted response values by the second‐order polynomial model were found to be in good agreement with experimental values (R² = 0.968 and adj‐R² = 0.940). The optimization result showed that the Cr(VI) immobilization efficiency presented the maximal result (90.63%) at the following optimal conditions: ZVIn dosage of 1.5% w/w, reaction time of 60 min, and an initial Cr(VI) concentration of 400 mg/kg.     

Characterization of fractured reservoirs using a consistent stiffness‐permeability model: focus on the effects of fracture aperture

In this paper we propose a method for the characterization of naturally fractured reservoirs by quantitative integration of seismic and production data. The method is based on a consistent theoretical frame work to model both effective hydraulic and elastic properties of fractured porous media and a (non‐linear) Bayesian method of inversion that provides information about uncertainties as well as mean (or maximum likelihood) values. We model a fractured reservoir as a porous medium containing a single set of vertical fractures characterized by an unknown fracture density, azimuthal orientation and aperture. We then look at the problem of fracture parameter estimation as a non‐linear inverse problem and try to estimate the unknown fracture parameters by joint inversion of seismic amplitude versus angle and azimuth data and dynamic production data. Once the fracture parameters have been estimated the corresponding effective stiffness and permeability tensors can be estimated using consistent models. A synthetic example is provided to clearly explain and test the workflow. It shows that seismic and production data complement each other, in the sense that the seismic data resolve a non‐uniqueness in the fracture orientation and the production data help to recover the true fracture aperture and permeability, because production data are more sensitive to the fracture aperture than the seismic data.     

Sensitivity of the strong ground motion time histories to a finite source model: A case study for the January 12, 2010 Haiti earthquake (Mw=7.0)

We analyze the waveforms generated by the January 12, 2010 Haiti earthquake (M_w=7.0) for its source characteristics. A 60 to 25 km source model is retrieved by the Kikuchi and Kanamori finite source inversion technique that uses broadband teleseismic body wave records. The derived rupture model points out unilateral rupture propagation commenced at the eastern side of the fault plane where the major seismic moment release occurred. The rupture front propagated westward and terminated at a site where the largest aftershocks occurred. Our estimates yield a seismic moment of M_o=8.17×10¹⁹ N m released on a 60 km-long fault plane. A patch at the eastern side of the ruptured fault plane inferred as a region of maximum moment release.     

Selecting and scaling real ground motion records using harmony search algorithm

In this study, a solution model is proposed to obtain input ground motion datasets compatible with given design spectra based on meta-heuristic harmony search algorithm. The utility of the solution model is demonstrated by generating ground motion datasets matching the Eurocode-8 design spectra for different soil types out of an extensive database of recorded motions. A total of 352 records are selected from the Pacific Earthquake Engineering Center (PEER) Strong Motion Database based on magnitude, distance, and site conditions to form the original ground motion domain. Then, the proposed harmony search based solution algorithm is applied on the pre-selected 352 time-series to obtain the ground motion record sets compatible with design spectra. The results demonstrate that the proposed HS based solution model provides an efficient way to develop input ground motion record sets that are consistent with code-based design spectra.     

Pb isotopic variability in the modern-Pleistocene Indus River system measured by ion microprobe in detrital K-feldspar grains

The western Himalaya, Karakoram and Tibet are known to be heterogeneous with regard to Pb isotope compositions in K-feldspars, which allows this system to be used as a sediment provenance tool. We used secondary ion mass spectrometry to measure the isotopic character of silt and sand-sized grains from the modern Sutlej and Chenab Rivers, together with Thar Desert sands, in order to constrain their origin. The rivers show a clear Himalayan provenance, contrasting with grains from the Indus Suture Zone, but with overlap to known Karakoram compositions. The desert dunes commonly show ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁶Pb/²⁰⁴Pb values that are much higher than those seen in the rivers, most consistent with erosion from Nanga Parbat. This implies at least some origin from the trunk Indus, probably reworked by summer monsoon winds from the SW, a hypothesis supported by bulk Nd and U-Pb zircon dating. Further data collected from Holocene and Pleistocene sands shows that filled and abandoned channels on the western edge of the Thar Desert were sourced from Himalayan rivers before and at 6-8 ka, but that after that time the proportion of high isotopic ratio grains rose, indicating increased contribution from the Thar Desert dunes prior to ∼4.5 ka when flow ceased entirely. This may be linked to climatic drying, northward expansion of the Thar Desert, or changes in drainage style including regional capture, channel abandonment, or active local Thar tributaries. Our data further show a Himalayan river channel east of the present Indus, close to the delta, in the Nara River valley during the middle Holocene. While this cannot be distinguished from the Indus it is not heavily contaminated by reworking from the desert. The Pb system shows some use as a provenance tool, but is not effective at demonstrating whether these Nara sediments represent a Ghaggar-Hakra stream independent from the Indus. Our study highlights an important role for eolian reworking of floodplain sediments in arid rivers such as the Indus.     

CFD modeling of the launder of settler of an industrial copper solvent extraction plant: A case study on Sarcheshmeh copper complex,Iran

In this study, a computational fluid dynamics (CFD) modeling accompanied by experimental field measurements has been applied to study the behavior of aqueous–organic dispersion on the performance of the launder of the settler of the copper solvent extraction plant at the Sarcheshmeh copper complex, Iran. Fluid flow field has been calculated by solving the continuity and Navier Stokes equations along with the standard k-ε turbulence model. The simulation results have been compared with the field measurements data to check the accuracy of the CFD work. The effects of picket fence on the launder performance and pressure drop over the picket fences in the launder have been investigated. The model predicts that by setting the picket fences, flow pattern becomes uniform, and turbulent eddies disappear. The results also show that by installing two picket fences without dam in the launder, the phase separation is improved and the performance of the launder is optimized.     

Widespread overspill from a saline density‐current channel and its interaction with topography on the south‐west Black Sea shelf

Seaward of the Bosphorus Strait, the south‐west Black Sea shelf is dominated by the world's largest channel network maintained by a quasi‐continuous saline (ca 35 → 31 psu) underflow. Calculations indicate that >85% of the initial discharge of ca 10⁴ m³s^?1 spills overbank before the shelf edge. This paper documents interaction of the overspill with sea bed topography using multibeam bathymetry, echo‐sounder images of the water column, conductivity–temperature–depth profiles and sediment cores. Overbank spill is widespread, particularly through crevasse channels and on the middle shelf where confinement by channel banks is negligible. Towards the outer shelf, the wind‐driven Rim Current advects mud along the shelf, contributing to levée successions and deposition on stoss sides of elongate transverse ridges. Echo‐sounder profiles reveal metre‐scale eddies over megaflutes, and breaking lee waves and internal hydraulic jumps over ridges. Megaflutes reach 600 m long and 7 m deep, yet form where the underflow, outside the flute, is no thicker than ca 2 to 5 m. Two types of elongate seaward‐facing ridges are recognized. Type 1 ridges, 2 to 5 m high, consist of bivalve‐rich muddy sand in low‐angle (3·5° to 6°) cross‐sets created by the underflow. Type 2 ridges, ca 5 m high, have crests up to 2 km long and a buried wedge‐shaped foundation (the ‘ridge‐core’) comprised of facies similar to Type 1 ridges. These ridge‐cores are blanketed on the landward side by stratified muds, and are capped by obliquely oriented ribs supporting a diverse benthic community. This facies distribution is interpreted to result from stoss‐side and lee‐side velocity and turbulence fluctuations induced by internal hydraulic jumps and breaking lee waves in overspilling portions of the underflow. Experimental results published by W.H. Snyder and co‐workers effectively explain ridge evolution and flow across the ridges, and therefore can be applied with confidence to less easily studied deep‐marine settings swept by turbidity currents.     

Public perception of the risk of disasters in a developing economy: the case of Saudi Arabia

Saudi Arabia has in recent years experienced frequent disasters, including flooding, epidemics, and dust storms, while many parts of the country are subject to regular earthquake and volcanic activity. The paper examines public perception of the risk of disasters in this interesting socio-cultural and regional environment not already covered by existing literature. A wide national survey conducted between March and May 2012 resulted in 1,164 responses across the 13 regions of Saudi Arabia. The study showed that the majority of the participants have faith that God is in control of the world and that disasters may be a punishment from him. However, this does not hinder their desire to be prepared to cope with disasters. It also highlighted that direct experience with such disasters does not directly influence perception. The research findings lead to the emergence of a number of recommendations regarding raising awareness of hazards and the risk of disasters, including education, training, encouraging voluntary work, and improving public access to vital information resources. A requirement for research into resilience also emerges, to prepare communities to cope with disasters; this is the focus of the authors’ future research.