Avalanche and landslide simulation using the material point method: flow dynamics and force interaction with structures

In this paper, the material point method (MPM) is presented as a tool for simulating large deformation, gravity-driven landslides. The primary goal is to assess the interaction of these flow-like events with the built environment. This includes an evaluation of earthen mounds when energy dissipating devices are placed in the path of a snow avalanche. The effectiveness of the embankments is characterized using displacement, velocity, mass, and energy measures. A second example quantifies the force interaction between a landslide and a square rigid column. Multiple slide approach angles are considered, and various aspects of the impact force are discussed.     

Buried pockmarks on the Top Chalk surface of the Danish North Sea and their potential significance for interpreting palaeocirculation patterns

Three-dimensional (3D) seismic data from the southern Danish North Sea were used to analyse the morphology and spatial distribution of depressions in the Danian Chalk deposit. Previously, these depressions were either interpreted as karst structures or pockmarks. The observed depressions occur in an interval from 25 ms below to 12 ms above the Top Chalk surface. Three types of depressions were differentiated based on their plan-view geometry and their degree of symmetry: Type 1, comprising sub-circular and symmetrical depressions, is the dominant group (ca. 70 %). Type 2, elongated and symmetrical depressions, represents only a small fraction (ca. 5 %). The elongated and asymmetrical depressions of Type 3 compose ca. 25 %. In cross section, each depression type can be either characterised by a V- or a U-shape. The maximum size of the depressions ranges from 50 to 580 m, with an average internal depth of 10 m. We interpret the depressions as pockmarks formed by the expulsion of biogenic or thermogenic fluids at the Danian seafloor. Likely, the initial form of the pockmarks has been circular (Type 1) and was subsequently modified for Types 2 and 3 to an elongated form by currents. The long axis of the pockmarks is interpreted to represent the effective current direction. The inferred direction is sub-parallel to the palaeobathymetric contours. The 3D seismic interpretation of pockmarks presented in this paper contributes to the understanding of fluid migration and palaeocirculation patterns during the sedimentation of the terminal Chalk Group in the southern Danish North Sea.     

Experimental study,numerical modeling of and axial prediction approach to base grouted drilled shafts in cohesionless soils

The pressure grouting of drilled shaft tips has become popular worldwide due to its effectiveness in mobilizing a larger portion of the available tip resistance under service displacements. This paper presents experimental and numerical studies on the load transfer mechanism and factors controlling the axial response of base grouted drilled shafts in cohesionless soils. The study found that the increased axial capacity of grout-tipped drilled shafts under service loads and displacements depended mainly on preloading effects and the increased tip area provided by the grouting process. A simple prediction approach for estimating the tip capacity of grouted shafts utilizing cone penetration resistance was suggested based on the results of the study. The validity of the proposed approach was verified by the analysis of full-scale case studies of grouted shafts reported in the literature.     

Stochastic Fracture Propagation Modelling for Enhanced Geothermal Systems

Fractures and fracture networks are the fundamental components of enhanced geothermal systems and determine their technical and economic viability. A realistic fracture model that can adequately describe a fracture-stimulated reservoir is critical for subsequent flow and heat transfer analyses of the system. Fractures in these systems are essentially the product of hydraulic stimulations of the reservoir that, together with ground conditions and the local stress regime, determine how fractures are formed and propagated. This paper describes three methods for generating realistic fracture models for enhanced geothermal systems; two of them incorporate the fracture propagation process in the modelling and hence provide a stochastic fracture propagation model. The methods are: a Bayesian framework in the form of Markov Chain Monte Carlo simulation, an extended Random Sampling Consensus model and a Point and Surface Association Consensus model. The conditioning data used in these methods are seismic events recorded during fracture stimulation. Geodynamics’ Habanero reservoir in the Cooper Basin of South Australia is used as a case study to test these methods.     

Responses of microorganisms and enzymes to soil contamination with metazachlor

A greenhouse pot experiment was conducted at the University of Warmia and Mazury in Olsztyn, with the aim of describing the influence of metazachlor on counts and biodiversity of soil microorganisms, soil enzymatic activity, physicochemical properties of soil and yield of spring oilseed rape. The first experimental factor was soil contamination with increasing rates of metazachlor: 0 (soil without the herbicide), 0.333 (recommended by the manufacturer), 6.666, 13.332, 26.665, 53.328, 106.656 and 213.312 mg kg^?1 dm of soil, while the second factor comprised two dates of determination: 30 and 60 days after starting the experiment. The tested herbicide had an adverse effect on reproduction of all analyzed microorganisms: oligotrophic bacteria and their endospore forms, Azotobacter spp. bacteria, organotrophic bacteria, actinomycetes and fungi. The values of the colony development and eco-physiological diversity indices decreased under the influence of excessive doses of the herbicide. Also, significant fluctuation in the enzymatic activity of soil was observed in response to the higher doses of metazachlor and depending on the date of determinations. The herbicide had an invariably negative influence on the activity of soil enzymes, causing the inhibition of dehydrogenases, catalase, urease, acid phosphatase, alkaline phosphatase, arylsulfatase and β-glucosidase. The physicochemical status of soil depended significantly on the degree of soil contamination with the herbicide, same as the yield of spring oilseed rape.     

Streamflow response to regional climate model output in the mountainous watershed: a case study from the Swiss Alps

Regional climate model (RCM) outputs are often used in hydrological modeling, in particular for streamflow forecasting. The heterogeneity of the meteorological variables such as precipitation, temperature, wind speed and solar radiation often limits the ability of the hydrological model performance. This paper assessed the sensitivity of RCM outputs from the PRUDENCE project and their performance in reproducing the streamflow. The soil and water assessment tool was used to simulate the streamflow of the Rhone River watershed located in the southwestern part of Switzerland, with the climate variables obtained from four RCMs. We analyzed the difference in magnitude of precipitation, maximum and minimum air temperature, and wind speed with respect to the observed values from the meteorological stations. In addition, we also focused on the impact of the grid resolution on model performance, by analyzing grids with resolutions of 50 × 50 and 25 × 25 km². The variability of the meteorological inputs from various RCMs is quite severe in the studied watershed. Among the four different RCMs, the Danish Meteorological Institute provided the best performance when simulating runoff. We found that temperature lapse rate is significantly important in the mountainous snow and glacier dominated watershed as compared to other variables like precipitation, and wind speed for hydrological performance. Therefore, emphasis should be given to minimum and maximum temperature in the bias correction studies for downscaling climatic data for impact modeling in the mountainous snow and glacier dominated complex watersheds.     

Behaviour of radiogenic Pb in zircon during ultrahigh-temperature metamorphism: an ion imaging and ion tomography case study from the Kerala Khondalite Belt,southern India

Zircon crystals from a locally charnockitized Paleoproterozoic high-K metagranite from the Kerala Khondalite Belt (KKB) of southern India have been investigated by high-spatial resolution secondary ion mass spectrometry analysis of U–Th–Pb and rare earth elements (REE), together with scanning ion imaging and scanning ion tomography (depth-profiled ion imaging). The spot analyses constrain the magmatic crystallization age of the metagranite to ca. 1,850 Ma, with ultrahigh-temperature (UHT) metamorphism occurring at ca. 570 Ma and superimposed charnockite formation at ca. 520–510 Ma, while the ion imaging reveals a patchy distribution of radiogenic Pb throughout the zircon cores. Middle- to heavy-REE depletion in ca. 570 Ma zircon rims suggests that these grew in equilibrium with garnet and therefore date the UHT metamorphism in the KKB. The maximum apparent ²⁰⁷Pb/²⁰⁶Pb age obtained from the unsupported radiogenic Pb concentrations is also consistent with formation of the Pb patches during this event. The superimposed charnockitization event appears to have caused additional Pb-loss in the cores and recrystallization of the rims. The results of depth-profiling of the scanning ion tomography image stack show that the Pb-rich domains range in size from <5 nm to several 10 nm (diameter if assumed to be spherical). The occurrence of such patchy Pb has previously been documented only from UHT metamorphic zircon, where it likely results from annealing of radiation-damaged zircon. The formation of a discrete, heterogeneously distributed and subsequently immobile Pb phase effectively arrests the normal Pb-loss process seen at lower grades of metamorphism.     

Petrogenesis of a voluminous Quaternary adakitic volcano: the case of Baru volcano

The origin of adakite magmas remains controversial because initially the term adakite had petrogenetic significance implying an origin by direct melting of the eclogitized subducting oceanic crust. Many models have been produced for their origin, and until now there has not been a straightforward method to discriminate between these models in a given adakite suite. Here, we use detailed chronological and geochemical studies of selected adakitic edifices that allows for the determination of the magmatic output rate parameter (Qe), which has been correlated with the rates of magma generation deep within subduction zones. By providing temporal and eruption rate estimates, we provide constraints on the possible petrogenetic processes involved in the generation of adakite-like signatures. Adakite magmas derived from the melting of the subducting slab should be volumetrically insignificant when compared to the adakite-like magmas produced by typical arc magma generation processes. In this study, we use this observation and the extraordinary stratigraphic exposure from Miocene to present in an adakitic volcano in Panama and to study the temporal and chemical variation in erupted magmas to estimate rates of magma generation. Detailed chemical and geochronological analyses of Baru volcano indicate that the volcanic edifice was constructed in its entirety during the Quaternary and magmas display adakite-like features such as steep rare earth elements patterns, pronounced depletions in the heavy rare earth elements, low Y, high Sr, and high Sr/Y. The magmatic output rates (Qe) that we have calculated show that compared to other typical adakitic volcanoes, most of the volcanic edifice of Baru volcano was constructed extremely rapidly (<~213 k.a.) and in time frames that are similar to typical arc volcanoes. The observed chemical and mineralogical variation, coupled with the high magma production rates, indicate that Baru volcano is more representative of a typical arc volcano than a small-volume melt of the subducting oceanic crust. The technique we outline may have broader application in determining the petrogenetic conditions of other adakite suites.     

Origin and significance of poikilitic and mosaic peridotite xenoliths in the western Pannonian Basin: geochemical and petrological evidences

Peridotite xenoliths erupted by alkali basaltic volcanoes in the western Pannonian Basin can be divided into two fundamentally contrasting groups. Geochemical characteristics of the abundant protogranular, porphyroclastic and equigranular nodules suggest that these samples originate from an old consolidated and moderately depleted lithospheric mantle domain. In contrast, the geochemical features of the worldwide rare, but in the Pannonian Basin relatively abundant, poikilitic xenoliths attest to a more complex evolution. It has been argued that the origin of the peculiar texture and chemistry may be intimately linked to melt/rock reactions at successively decreasing liquid volumes in a porous melt flow system. The most likely site where such reactions can take place is the asthenosphere–lithosphere boundary. In this context, poikilitic xenoliths may provide petrological and geochemical evidence for reactions between magmatic liquids issued from the uprising asthenosphere and the solid mantle rocks of the lithosphere. These reactions are important agents of the thermal erosion of the lithosphere; thus, they could have considerably contributed to the thinning of the lithosphere in the Pannonian region. We suggest that in the Pannonian Basin, there could be a strong relation between the unusual abundance of poikilitic mantle xenoliths and the strongly eroded lithosphere.     

Overstepping the garnet isograd: a comparison of QuiG barometry and thermodynamic modeling

The consequences of overstepping the garnet isograd reaction have been investigated by comparing the composition of garnet formed at overstepped P–T conditions (the overstep or “OS” model) with the P–T conditions that would be inferred by assuming garnet nucleated in equilibrium with the matrix assemblage at the isograd (the equilibrium or “EQ” model). The garnet nucleus composition formed at overstepped conditions is calculated as the composition that produces the maximum decrease in Gibbs free energy from the equilibrated, garnet-absent, matrix assemblage for the bulk composition under study. Isopleths were then calculated for this garnet nucleus composition assuming equilibrium with the matrix assemblage (the EQ model). Comparison of the actual P–T conditions of nucleation (the OS model) with those inferred from the EQ model reveals considerable discrepancy between the two. In general, the inferred garnet nucleation P–T conditions (the EQ model) are at a lower temperature and higher or lower pressure (depending on the coexisting calcic phase(s)) than the actual (OS model) nucleation conditions. Moreover, the degree of discrepancy increases with the degree of overstepping. Independent estimates of the pressure of nucleation of garnet were made using the Raman shift of quartz inclusions in garnet (quartz-in-garnet or QuiG barometry). To test the validity of this method, an experimental synthesis of garnet containing quartz inclusions was made at 800 °C, 20 kbar, and the measured Raman shift reproduced the synthesis conditions to within 120 bars. Raman band shifts from three natural samples were then used to calculate an isochore along which garnet was presumed to have nucleated. Model calculations were made at several temperatures along this isochore (the OS model), and these P–T conditions were compared to those computed assuming equilibrium nucleation (the EQ model) to estimate the degree of overstepping displayed by these samples. A sample from the garnet isograd in eastern Vermont is consistent with overstepping of around 10 degrees and 0.6 kbar (affinities of around 2 kJ/mole garnet). A sample from the staurolite–kyanite zone in the same terrane requires overstepping of around 50 °C and 2–5 kbar (affinities of around 10–18 kJ/mole garnet). A similar amount of overstepping was inferred for a blueschist sample from Sifnos, Greece. These results indicate that overstepping of garnet nucleation reactions may be common and pronounced in regionally metamorphosed terranes, and that the P–T conditions and paths inferred from garnet zoning studies may be egregiously in error.