Modeling the 2D behavior of dry‐stone retaining walls by a fully discrete element method

The dry‐stone retaining walls (DSRW) have been tipped as a promising solution for sustainable development. However, before recently, their behavior is relatively obscure. In this study, discrete element method (DEM) approach was applied to simulate the plane strain failure of these walls. A commercial DEM package (PFC2D™) was used throughout this study. The authors used a fully discrete approach; thus, both the wall and the backfill were modeled as discrete elements. The methodology for obtaining the micromechanical parameters was discussed in detail; this includes the three mechanical sub‐systems of DSRWs: wall, backfill and interface. The models were loaded progressively until failure, and then the results were compared with the full‐scale experimental results where the walls were loaded, respectively, with hydrostatic load and backfill. Despite its complexity and its intensive calculation time, DEM model can then be used to validate a more simplified approach. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of a finite element model to the computation of tides in the Mersey Estuary and Eastern Irish Sea

A variable mesh finite element model of the Irish and Celtic Sea regions with/without the inclusion of the Mersey estuary is used to examine the influence of grid resolution and the Mersey upon the higher harmonics of the tide in the region. Comparisons are made with observations and published results from finite difference models of the area. Although including a high resolution representation of the Mersey had little effect upon computed tides in the western Irish Sea it had a significant effect upon tidal currents in the eastern Irish Sea. In addition the higher harmonics of the M₂ tide in near-shore regions of the eastern Irish Sea particularly the Solway and Mersey estuary together with Morecambe Bay showed significant small scale variability. The Mersey was used to test the sensitivity to including estuaries because high resolution accurate topography was available. The results presented here suggest that comparable detailed topographic data sets are required in all estuaries and near-shore regions. In addition comparisons clearly show the need for an unstructured grid model of the region that can include all the estuaries. Such an unstructured grid solution was developed here within a finite element approach, although other methods in particular the finite volume, or coordinate transformations/curvilinear grids and nesting could be applied.     

Oblique oceanic opening and passive margin irregularity,as inherited in the Zagros fold‐and‐thrust belt

The Zagros fold‐and‐thrust belt of SW Iran represents deformation of the former Arabian passive margin since Permian–Triassic opening of the Neo‐Tethys ocean. The Zagros belt is characterized by a present‐day structural salient‐recess setting inherited from past marginal embayment‐promontory geometry, which was involved in discontinuous ophiolite obduction and diachronous continental collision. We examine outcrop‐scale Mesozoic extensional brittle tectonics, preserved as syn‐depositional normal faults within the folded strata, in terms of stress tensor inversion. The result is then integrated with belt‐scale isopach, seismic and topographical data to delineate the geometry of a major irregularity along the passive margin originating from oblique oceanic opening. The implication of this configuration within the tectonic framework of oceanic closure is discussed.     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

The groundwater recharge regime of some slightly metamorphosed neoproterozoic sedimentary rocks: an application of natural environmental tracers

Isotope data of precipitation and groundwater in parts of the Voltaian Basin in Northern Ghana were used to explain the groundwater recharge regime in the area. Groundwater recharge is an important parameter in the development of a decision support system for the management and efficient utilization of groundwater resources in the area. It is therefore important to establish the processes and sources of groundwater recharge. δ¹⁸O and δ²H data for local precipitation suggest enrichment relative to the Global Meteoric Water Line (GMWL) and indicate that precipitation takes place at a relative humidity less than 100%. The groundwater data plot on an evaporation line with a slope of 5, suggesting a high degree of evaporative enrichment of the precipitation in the process of vertical infiltration and percolation through the unsaturated zone into the saturated zone. This finding is consistent with the observation of high evapotranspiration rates in the area and ties in with the fact that significant clay fraction in the unsaturated zone limits vertical percolation and thus exposes the percolating rainwater to the effects of high temperatures and low humidities resulting in high evapotranspiration rates. Groundwater recharge estimates from the chloride mass balance, CMB, method suggest recharge in the range of 1.8–32% of the annual average precipitation in the form of rainfall. The highest rates are associated with areas where open wells encourage significant amount of groundwater recharge from precipitation in the area. In the northern parts of the study area, groundwater recharge is lower than 12%. The recharge so computed through the application of the CMB methodology takes on a spatial distribution akin to the converse of the spatial pattern of both δ¹⁸O and δ²H in the area. As such, the locations of the highest recharge are associated with the most depleted values of the two isotopes. This observation is consistent with the assertion that low vertical hydraulic conductivities slow down vertical percolation of precipitation down to the groundwater water. The percolating precipitation water thus gets enriched in the heavier isotopes through high evapotranspiration rates. At the same time, the amount of water that finally reaches the water table is considerably reduced. Copyright © 2013 John Wiley & Sons, Ltd.     

Formation of mud clast breccias and the process of sedimentary autobrecciation in the hominin-bearing (Homo naledi) Rising Star Cave system,South Africa

Unconsolidated mud clast breccia facies in the hominin-bearing (Homo naledi) Rising Star Cave, Cradle of Humankind, South Africa, are interpreted to have formed through a process termed sedimentary autobrecciation in this study. This process, by which most of the angular mud clast breccia deposits are thought to have formed autochthonously to para-autochthonously via a combination of erosion, desiccation, diagenesis and microbial alteration of laminated mud deposits, is thought to have taken place under relatively dry (i.e. non-flooded) conditions inside the cave. Subsequently, gravitational slumping and collapse was the dominant mechanism that produced the mud clast breccia deposits, which commonly accumulate into debris aprons. The mud clast breccia is typically associated with (micro) mammal fossils and is a common facies throughout the cave system, occurring in lithified and unlithified form. This facies has not been described from other cave localities in the Cradle of Humankind. Additionally, sedimentary autobrecciation took place during the deposition of some of the fossils within the Rising Star Cave, including the abundant Homo naledi skeletal remains found in the Dinaledi Subsystem. Reworking of the mud clast breccia deposits occurs in some chambers as they slump towards floor drains, resulting in the repositioning of fossils embedded in the breccias as evidenced by cross-cutting manganese staining lines on some Homo naledi fossil remains. The formation of the unlithified mud clast breccia deposits is a slow process, with first order formation rates estimated to be ca 8 × 10⁻⁴ mm year⁻¹. The slow formation of the unlithified mud clast breccia facies sediments and lack of laminated mud facies within these deposits, indicates that conditions in the Dinaledi Chamber were probably stable and dry for at least the last ca 300 ka, meaning that this study excludes Homo naledi being actively transported by fluvial mechanisms during the time their remains entered the cave.     

Coseismic Fault Slip of the September 16, 2015 <Emphasis Type="Italic">Mw</Emphasis> 8.3 Illapel,Chile Earthquake Estimated from InSAR Data

The complete surface deformation of 2015 Mw 8.3 Illapel, Chile earthquake is obtained using SAR interferograms obtained for descending and ascending Sentinel-1 orbits. We find that the Illapel event is predominantly thrust, as expected for an earthquake on the interface between the Nazca and South America plates, with a slight right-lateral strike slip component. The maximum thrust-slip and right-lateral strike slip reach 8.3 and 1.5 m, respectively, both located at a depth of 8 km, northwest to the epicenter. The total estimated seismic moment is 3.28 × 10²¹ N.m, corresponding to a moment magnitude Mw 8.27. In our model, the rupture breaks all the way up to the sea-floor at the trench, which is consistent with the destructive tsunami following the earthquake. We also find the slip distribution correlates closely with previous estimates of interseismic locking distribution. We argue that positive coulomb stress changes caused by the Illapel earthquake may favor earthquakes on the extensional faults in this area. Finally, based on our inferred coseismic slip model and coulomb stress calculation, we envision that the subduction interface that last slipped in the 1922 Mw 8.4 Vallenar earthquake might be near the upper end of its seismic quiescence, and the earthquake potential in this region is urgent.     

A comparison of shelf observation platforms for assimilation in an eddy-resolving ocean model

An assessment of the likely benefits of assimilating in situ temperature (T) and salinity (S) observations from repeat glider transects and surface velocity observations from high-frequency radar arrays into an eddy-resolving ocean model is presented. The deployment of new shelf observation platforms around Australia is being undertaken through the Australian Integrated Marine Observing System program. In this study, various options for an observing system along the coast of New South Wales, Australia, are assessed for their benefits to an ocean forecast and reanalysis system. The forecast system considered here uses ensemble optimal interpolation (EnOI) for data assimilation. Using error estimates from the EnOI scheme, estimates of the theoretical analysis errors are calculated for different observing systems that include a range of remotely sensed and in situ observations. The results demonstrate that if HF radar observations are assimilated along with the standard components of the global ocean observing system, the analysis errors are likely to reduce by as much as 80% for velocity and 60% for T, S and sea-level in the vicinity of the observations. Owing to the relatively short along-shore decorrelation length-scales for T and S near the shelf, the glider observations are likely to provide the forecast system with a more modest gain.     

The Geochemistry of Lake Joyce, McMurdo Dry Valleys, Antarctica

Lake Joyce is one of the least studied lakes of the McMurdo Dry Valleys. Similar to other lakes in this region, Lake Joyce is a closed-basin, permanently ice-covered, meromictic lake. We present here a detailed investigation of major ions, nutrients, and dissolved trace elements for Lake Joyce. Specifically, we investigate the role of iron and manganese oxides and hydrous oxides in trace metal cycling.Lake Joyce is characterized by fresh, oxic waters overlying an anoxic brine, primarily Na–Cl. Surface waters have a maximum nitrate concentration of 26M with a molar dissolved inorganic nitrogen to phosphorus ratio of 477. The supply of nitrogen is attributed to atmospheric deposition, possibly from polar stratospheric clouds. Dissolved phosphorus is scavenged by hydrous iron oxides. The pH is highest (10.15) just beneath the 7-m thick ice cover and decreases to a minimum of 7.29 in the redox transition zone. Dissolved Al exceeds 8M in surface waters, and appears to be controlled by equilibrium with gibbsite. In contrast, concentrations of other trace elements in surface waters are quite low (e.g., 5.4nM Cu, 0.19nM Co, <20pM La). Dissolved Fe, Mn, Ni and Cd were below our detection limits of 13 nM, 1. 8 nM, 4.7 nM and 15pM (respectively) in surface waters. There was a 6-m vertical separation in the onset of Mn and Fe reduction, with dissolved Mn appearing higher in the water column than Fe. Based on thermodynamic calculations, dissolved Mn appears to be controlled by equilibrium with hausmannite (Mn3O4). Co tracks the Mn profile closely, suggesting Co(III) is bound in the lattice of Mn oxides, whereas the Ce profile is similar, yet the Ce anomaly suggests oxidative scavenging of Ce. Release of Cu, Ni, Cd and trivalent REE appears to be controlled by pH-induced desorption from Fe and Mn oxides, although Cu (and perhaps Ni) may be scavenged by organic matter in surface waters.     

Mechanisms of coal metamorphism: case studies from Paleozoic coalfields

Controls on coal metamorphism can be complex. In this paper, we examine four Paleozoic coalfields: the western Kentucky portion of the Illinois Basin, the Pennsylvania anthracite fields, the South Wales Coalfield, and the Bowen Basin. An increase in temperature with depth of burial is certainly a factor in coal metamorphism. In many coalfields, however, including the coalfields reviewed here, it has become apparent that such a simple mechanism does not explain the coal rank patterns observed. The flow of hydrothermal fluids through the coals has been proposed as a cause of coal metamorphism. Evidence includes inverted rank gradients, elevated CFL as an indicator of brine fluids, isotopic evidence for hydrothermal fluids, and vein and cleat mineral assemblages. In any case, multiple hypotheses must often be evaluated in the examination of any coalfield since the simple paradigm of coal rank increases with a simple increase in temperature with increasing depth does not fit the evidence observed in many cases.