Multilayer shallow water flow using lattice Boltzmann method with high performance computing

A multilayer lattice Boltzmann (LB) model is introduced to solve three-dimensional wind-driven shallow water flow problems. The multilayer LB model avoids the expensive Navier–Stokes equations and obtains stratified horizontal flow velocities as vertical velocities are relatively small and the flow is still within the shallow water regime. A single relaxation time BGK method is used to solve each layer coupled by the vertical viscosity forcing term. To increase solution stability, an implicit step is suggested to obtain flow velocities. The main advantage of using the LBM is that after selecting appropriate equilibrium distribution functions, the LB algorithm is only slightly modified for each layer and retains all the simplicities of the LBM within the high performance computing (HPC) environment. The performance of the parallel LB model for the multilayer shallow water equations is investigated on CPU-based HPC environments using OpenMP. We found that the explicit loop control with cache optimization in LBM gives better performance on execution time, speedup and efficiency than the implicit loop control as the number of processors increases. Numerical examples are presented to verify the multilayer LB model against analytical solutions. We demonstrate the model’s capability of calculating lateral and vertical distributions of velocities for wind-driven circulation over non-uniform bathymetry.     

Three‐dimensional flow structure and patterns of bed shear stress in an evolving compound meander bend

Compound meander bends with multiple lobes of maximum curvature are common in actively evolving lowland rivers. Interaction among spatial patterns of mean flow, turbulence, bed morphology, bank failures and channel migration in compound bends is poorly understood. In this paper, acoustic Doppler current profiler (ADCP) measurements of the three‐dimensional (3D) flow velocities in a compound bend are examined to evaluate the influence of channel curvature and hydrologic variability on the structure of flow within the bend. Flow structure at various flow stages is related to changes in bed morphology over the study timeframe. Increases in local curvature within the upstream lobe of the bend reduce outer bank velocities at morphologically significant flows, creating a region that protects the bank from high momentum flow and high bed shear stresses. The dimensionless radius of curvature in the upstream lobe is one‐third less than that of the downstream lobe, with average bank erosion rates less than half of the erosion rates for the downstream lobe. Higher bank erosion rates within the downstream lobe correspond to the shift in a core of high velocity and bed shear stresses toward the outer bank as flow moves through the two lobes. These erosion patterns provide a mechanism for continued migration of the downstream lobe in the near future. Bed material size distributions within the bend correspond to spatial patterns of bed shear stress magnitudes, indicating that bed material sorting within the bend is governed by bed shear stress. Results suggest that patterns of flow, sediment entrainment, and planform evolution in compound meander bends are more complex than in simple meander bends. Moreover, interactions among local influences on the flow, such as woody debris, local topographic steering, and locally high curvature, tend to cause compound bends to evolve toward increasing planform complexity over time rather than stable configurations. Copyright © 2016 John Wiley & Sons, Ltd.     

Fourier spectral- and duration models for the generation of response spectra adjustable to different source-, propagation-, and site conditions

One of the major challenges related with the current practice in seismic hazard studies is the adjustment of empirical ground motion prediction equations (GMPEs) to different seismological environments. We believe that the key to accommodating differences in regional seismological attributes of a ground motion model lies in the Fourier spectrum. In the present study, we attempt to explore a new approach for the development of response spectral GMPEs, which is fully consistent with linear system theory when it comes to adjustment issues. This approach consists of developing empirical prediction equations for Fourier spectra and for a particular duration estimate of ground motion which is tuned to optimize the fit between response spectra obtained through the random vibration theory framework and the classical way. The presented analysis for the development of GMPEs is performed on the recently compiled reference database for seismic ground motion in Europe (RESORCE-2012). Although, the main motivation for the presented approach is the adjustability and the use of the corresponding model to generate data driven host-to-target conversions, even as a standalone response spectral model it compares reasonably well with the GMPEs of Ambraseys et al. (Bull Earthq Eng 3:1–53, 2005), Akkar and Bommer (Seismol Res Lett 81(2):195–206, 2010) and Akkar and Cagnan (Bull Seismol Soc Am 100(6):2978–2995, 2010).     

Impact of the configuration of stretching and ocean–atmosphere coupling on tropical cyclone activity in the variable-resolution GCM ARPEGE

This study starts by investigating the impact of the configuration of the variable-resolution atmospheric grid on tropical cyclone (TC) activity. The French atmospheric general circulation model ARPEGE, the grid of which is rotated and stretched over the North Atlantic basin, was used with prescribed sea surface temperatures. The study clearly shows that changing the position of the stretching pole strongly modifies the representation of TC activity over the North Atlantic basin. A pole in the centre of the North Atlantic basin provides the best representation of the TC activity for this region. In a second part, the variable-resolution climate model ARPEGE is coupled with the European oceanic global climate model NEMO in order to study the impact of ocean–atmosphere coupling on TC activity over the North Atlantic basin. Two pre-industrial runs, a coupled simulation and a simulation forced by the sea surface temperatures from the coupled one, are compared. The results show that the coupled simulation is more active in the Caribbean Sea and the Gulf of Mexico while the forced simulation is more active over eastern Florida and the eastern Atlantic. The difference in the distribution of TC activity is certainly linked with the location of TC genesis. In the forced simulation, tropical cyclogenesis is closer to the west African coast than in the coupled simulation. Moreover, the difference in TC activity over the eastern Atlantic seems to be related to two different mechanisms: the difference in African easterly wave activity over the west of Africa and the cooling produced, in the coupled simulation, by African easterly waves over the eastern Atlantic. Finally, the last part studies the impact of changing the frequency of ocean–atmosphere coupling on Atlantic TC activity. Increasing the frequency of coupling decreases the density of TC activity over the North Atlantic basin. However, it does not modify the spatial distribution of the TC activity. TC rainfalls are decreased by 8?% in the high frequency coupled run.     

Correction to: Development and Evaluation of Geostatistical Methods for Non-Euclidean-Based Spatial Covariance Matrices

Mathematical Geosciences - The original version of this article unfortunately contained a mistake in equation 9.     

Silicon isotopic fractionation of CAI-like vacuum evaporation residues

Calcium-, aluminum-rich inclusions (CAIs) are often enriched in the heavy isotopes of magnesium and silicon relative to bulk solar system materials. It is likely that these isotopic enrichments resulted from evaporative mass loss of magnesium and silicon from early solar system condensates while they were molten during one or more high-temperature reheating events. Quantitative interpretation of these enrichments requires laboratory determinations of the evaporation kinetics and associated isotopic fractionation effects for these elements. The experimental data for the kinetics of evaporation of magnesium and silicon and the evaporative isotopic fractionation of magnesium is reasonably complete for Type B CAI liquids (Richter F. M., Davis A. M., Ebel D. S., and Hashimoto A. (2002) Elemental and isotopic fractionation of Type B CAIs: experiments, theoretical considerations, and constraints on their thermal evolution. Geochim. Cosmochim. Acta66, 521-540; Richter F. M., Janney P. E., Mendybaev R. A., Davis A. M., and Wadhwa M. (2007a) Elemental and isotopic fractionation of Type B CAI-like liquids by evaporation. Geochim. Cosmochim. Acta71, 5544-5564.). However, the isotopic fractionation factor for silicon evaporating from such liquids has not been as extensively studied. Here we report new ion microprobe silicon isotopic measurements of residual glass from partial evaporation of Type B CAI liquids into vacuum. The silicon isotopic fractionation is reported as a kinetic fractionation factor, α_Si, corresponding to the ratio of the silicon isotopic composition of the evaporation flux to that of the residual silicate liquid. For CAI-like melts, we find that α_Si = 0.98985 ± 0.00044 (2σ) for ²⁹Si/²⁸Si with no resolvable variation with temperature over the temperature range of the experiments, 1600-1900 °C. This value is different from what has been reported for evaporation of liquid Mg₂SiO₄ (Davis A. M., Hashimoto A., Clayton R. N., and Mayeda T. K. (1990) Isotope mass fractionation during evaporation of Mg₂SiO₄. Nature347, 655-658.) and of a melt with CI chondritic proportions of the major elements (Wang J., Davis A. M., Clayton R. N., Mayeda T. K., and Hashimoto A. (2001) Chemical and isotopic fractionation during the evaporation of the FeO-MgO-SiO₂-CaO-Al₂O₃-TiO₂-REE melt system. Geochim. Cosmochim. Acta65, 479-494.). There appears to be some compositional control on α_Si, whereas no compositional effects have been reported for α_Mg. We use the values of α_Si and α_Mg, to calculate the chemical compositions of the unevaporated precursors of a number of isotopically fractionated CAIs from CV chondrites whose chemical compositions and magnesium and silicon isotopic compositions have been previously measured.     

Topographic expression of active faults in the foothills of the Northern Apennines

Active faults that rupture the earth's surface leave an imprint on the topography that is recognized using a combination of geomorphic and geologic metrics including triangular facets, the shape of mountain fronts, the drainage network, and incised river valleys with inset terraces. We document the presence of a network of active, high-angle extensional faults, collectively embedded in the actively shortening mountain front of the Northern Apennines, that possess unique geomorphic expressions. We measure the strain rate for these structures and find that they have a constant throw-to-length ratio. We demonstrate the necessary and sufficient conditions for triangular facet development in the footwalls of these faults and argue that rock-type exerts the strongest control. The slip rates of these faults range from 0.1 to 0.3 mm/yr, which is similar to the average rate of river incision and mountain front unroofing determined by corollary studies. The faults are a near-surface manifestation of deeper crustal processes that are actively uplifting rocks and growing topography at a rate commensurate with surface processes that are eroding the mountain front to base level.     

A microbial pathway for the formation of gold-anomalous calcrete

The formation of pedogenic carbonate (calcrete) in terrestrial environments is commonly mediated by microorganisms. In Australia, Au-anomalous calcrete is an important sampling medium for geochemical exploration, but current models describing its formation do not include a confirmed microbial component. This study demonstrates that bacterial communities in calcareous sands from dunes overlying the Barns Gold Deposit in semi-arid South Australia, are capable of mediating the biomineralisation of Au-anomalous carbonates. Bacterial enrichment cultures obtained from calcareous sands at three depths (0.1, 0.64 and 2.1 m, plus abiotic control) were incubated in urea and Ca²⁺-containing growth media (pH 8), unamended and amended with Au (100 parts-per-billion, ppb) as Au–aspartic-acid complex. During the incubation of the enrichment cultures urea was turned over to NH₄⁺ within 96 h to 220 h. The solution pH increased concurrently by approximately 1.2 units, and Au-anomalous Ca-carbonate crystallites were precipitated on cells, which functioned as nucleation sites; no carbonate precipitation was observed in abiotic controls. Compared to the medium, Au was strongly enriched in these carbonates and appeared to be uniformly dispersed in the individual crystallites, as shown using LA-ICP-MS; a similar distribution is present in naturally occurring Au-anomalous calcrete. Phylogenetic 16S rRNA PCR DGGE analyses, shotgun cloning and functional microbial analyses (BioLog, ureC quantitative PCR) demonstrated that naturally occurring and culture-enriched bacterial communities were dominated by alkaliphylic, halotolerant Bacillus spp. The indigenous bacterial communities were capable of utilising amino acids (including l-aspartic acid) and urea, which appears to lead to the destabilisation of the Au–amino acid complexes and concomitant co-precipitation of Au in the Ca-carbonates. In conclusion, a model combining geomicrobial– with evapotranspiration– and plant-based components is likely to best describe the formation of (Au-anomalous) calcrete in semi-arid and arid zones.