Discrete element model for crack propagation in brittle materials

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical modeling of fracking fluid migration through fault zones and fractures in the North German Basin

Gas production from shale formations by hydraulic fracturing has raised concerns about the effects on the quality of fresh groundwater. The migration of injected fracking fluids towards the surface was investigated in the North German Basin, based on the known standard lithology. This included cases with natural preferential pathways such as permeable fault zones and fracture networks. Conservative assumptions were applied in the simulation of flow and mass transport triggered by a high pressure boundary of up to 50 MPa excess pressure. The results show no significant fluid migration for a case with undisturbed cap rocks and a maximum of 41 m vertical transport within a permeable fault zone during the pressurization. Open fractures, if present, strongly control the flow field and migration; here vertical transport of fracking fluids reaches up to 200 m during hydraulic fracturing simulation. Long-term transport of the injected water was simulated for 300 years. The fracking fluid rises vertically within the fault zone up to 485 m due to buoyancy. Progressively, it is transported horizontally into sandstone layers, following the natural groundwater flow direction. In the long-term, the injected fluids are diluted to minor concentrations. Despite the presence of permeable pathways, the injected fracking fluids in the reported model did not reach near-surface aquifers, either during the hydraulic fracturing or in the long term. Therefore, the probability of impacts on shallow groundwater by the rise of fracking fluids from a deep shale-gas formation through the geological underground to the surface is small.     

The constancy of galactic cosmic rays as recorded by cosmogenic nuclides in iron meteorites

We measured the He, Ne, and Ar isotopic concentrations and the ¹⁰Be, ²⁶Al, ³⁶Cl, and ⁴¹Ca concentrations in 56 iron meteorites of groups IIIAB, IIAB, IVA, IC, IIA, IIB, and one ungrouped. From ⁴¹Ca and ³⁶Cl data, we calculated terrestrial ages indistinguishable from zero for six samples, indicating recent falls, up to 562 ± 86 ka. Three of the studied meteorites are falls. The data for the other 47 irons confirm that terrestrial ages for iron meteorites can be as long as a few hundred thousand years even in relatively humid conditions. The ³⁶Cl‐³⁶Ar cosmic ray exposure (CRE) ages range from 4.3 ± 0.4 Ma to 652 ± 99 Ma. By including literature data, we established a consistent and reliable CRE age database for 67 iron meteorites. The high quality of the CRE ages enables us to study structures in the CRE age histogram more reliably. At first sight, the CRE age histogram shows peaks at about 400 and 630 Ma. After correction for pairing, the updated CRE age histogram comprises 41 individual samples and shows no indications of temporal periodicity, especially not if one considers each iron meteorite group separately. Our study contradicts the hypothesis of periodic GCR intensity variations (Shaviv 2002, 2003), confirming other studies indicating that there are no periodic structures in the CRE age histogram (e.g., Rahmstorf et al. 2004; Jahnke 2005). The data contradict the hypothesis that periodic GCR intensity variations might have triggered periodic Earth climate changes. The ³⁶Cl‐³⁶Ar CRE ages are on average 40% lower than the ⁴¹K‐K CRE ages (e.g., Voshage 1967). This offset can either be due to an offset in the ⁴¹K‐K dating system or due to a significantly lower GCR intensity in the time interval 195–656 Ma compared to the recent past. A 40% lower GCR intensity, however, would have increased the Earth temperature by up to 2 °C, which seems unrealistic and leaves an ill‐defined ⁴¹K‐K CRE age system the most likely explanation. Finally, we present new ²⁶Al/²¹Ne and ¹⁰Be/²¹Ne production rate ratios of 0.32 ± 0.01 and 0.44 ± 0.03, respectively.     

Development and validation of a finite element morphological model for shallow water basins

Recognising the importance of understanding sediment dynamics to evaluate the status of a coastal lagoon environment, this work has been focused on the investigation of the hydrodynamic and sediment transport processes occurring in such basins. In order to describe the lagoon system, a modelling approach combining hydrodynamics, waves and sediment dynamics has been developed. The framework of the numerical model consists of a finite element hydrodynamic model, a third generation finite element spectral wave model and a sediment transport and morphodynamic model for both cohesive and non-cohesive sediments. The model adopts the finite element technique for spatial integration, which has the advantage to describe more accurately complicated bathymetry and irregular boundaries for shallow water areas. The developed model has been applied to test cases and to a very shallow tidal lagoon, the Venice Lagoon, Italy. Numerical results show good agreement with water level, waves and turbidity measurements collected in several monitoring stations inside the Lagoon of Venice. Such a model represents an indispensable tool in analysing coastal problems and assessing morphological impacts of human interference.     

A cross-ecosystem comparison of spatial and temporal patterns of covariation in the recruitment of functionally analogous fish stocks

Temporal and spatial patterns of recruitment (R) and spawning stock biomass (S) variability were compared among functionally analogous species and similar feeding guilds from six marine ecosystems. Data were aggregated into four regions including the Gulf of Maine/Georges Bank, the Norwegian/Barents Seas, the eastern Bering Sea, and the Gulf of Alaska. Variability was characterized by calculating coefficients of variation and anomalies for three response variables: ln(R), ln(R/S), and stock–recruit model residuals. Patterns of synchrony and asynchrony in the response variables were examined among and between ecosystems, between- and within-ocean basins and among functionally analogous species groups using pair-wise correlation analysis corrected for within-time series autocorrelation, multivariate cross-correlation analyses and regime shift detectors. Time series trends in response variables showed consistent within basin similarities and consistent and coherent differences between the Atlantic and Pacific basin ecosystems. Regime shift detection algorithms identified two broad-scale regime shift time periods for the pelagic feeding guild (1972–1976 and 1999–2002) and possibly one for the benthic feeding guild (1999–2002). No spatial patterns in response variable coefficients of variation were observed. Results from multivariate cross-correlation analysis showed similar trends. The data suggest common external factors act in synchrony on stocks within ocean basins but temporal stock patterns, often of the same species or functional group, between basins change in opposition to each other. Basin-scale results (similar within but different between) suggest that the two geographically broad areas are connected by unknown mechanisms that, depending on the year, may influence the two basins in opposite ways. This work demonstrates that commonalities and synchronies in recruitment fluctuations can be found across geographically distant ecosystems but biophysical causes of the fluctuations remain difficult to identify.     

Experimental impact cratering: A summary of the major results of the MEMIN research unit

This paper reviews major findings of the Multidisciplinary Experimental and Modeling Impact Crater Research Network (MEMIN). MEMIN is a consortium, funded from 2009 till 2017 by the German Research Foundation, and is aimed at investigating impact cratering processes by experimental and modeling approaches. The vision of this network has been to comprehensively quantify impact processes by conducting a strictly controlled experimental campaign at the laboratory scale, together with a multidisciplinary analytical approach. Central to MEMIN has been the use of powerful two-stage light-gas accelerators capable of producing impact craters in the decimeter size range in solid rocks that allowed detailed spatial analyses of petrophysical, structural, and geochemical changes in target rocks and ejecta. In addition, explosive setups, membrane-driven diamond anvil cells, as well as laser irradiation and split Hopkinson pressure bar technologies have been used to study the response of minerals and rocks to shock and dynamic loading as well as high-temperature conditions. We used Seeberger sandstone, Taunus quartzite, Carrara marble, and Weibern tuff as major target rock types. In concert with the experiments we conducted mesoscale numerical simulations of shock wave propagation in heterogeneous rocks resolving the complex response of grains and pores to compressive, shear, and tensile loading and macroscale modeling of crater formation and fracturing. Major results comprise (1) projectile–target interaction, (2) various aspects of shock metamorphism with special focus on low shock pressures and effects of target porosity and water saturation, (3) crater morphologies and cratering efficiencies in various nonporous and porous lithologies, (4) in situ target damage, (5) ejecta dynamics, and (6) geophysical survey of experimental craters.     

High-precision stellar abundances of the elements: methods and applications

Efficient spectrographs at large telescopes have made it possible to obtain high-resolution spectra of stars with high signal-to-noise ratio and advances in model atmosphere analyses have enabled estimates of high-precision differential abundances of the elements from these spectra, i.e. with errors in the range 0.01–0.03 dex for F, G, and K stars. Methods to determine such high-precision abundances together with precise values of effective temperatures and surface gravities from equivalent widths of spectral lines or by spectrum synthesis techniques are outlined, and effects on abundance determinations from using a 3D non-LTE analysis instead of a classical 1D LTE analysis are considered. The determination of high-precision stellar abundances of the elements has led to the discovery of unexpected phenomena and relations with important bearings on the astrophysics of galaxies, stars, and planets, i.e. (i) Existence of discrete stellar populations within each of the main Galactic components (disk, halo, and bulge) providing new constraints on models for the formation of the Milky Way. (ii) Differences in the relation between abundances and elemental condensation temperature for the Sun and solar twins suggesting dust-cleansing effects in proto-planetary disks and/or engulfment of planets by stars; (iii) Differences in chemical composition between binary star components and between members of open or globular clusters showing that star- and cluster-formation processes are more complicated than previously thought; (iv) Tight relations between some abundance ratios and age for solar-like stars providing new constraints on nucleosynthesis and Galactic chemical evolution models as well as the composition of terrestrial exoplanets. We conclude that if stellar abundances with precisions of 0.01–0.03 dex can be achieved in studies of more distant stars and stars on the giant and supergiant branches, many more interesting future applications, of great relevance to stellar and galaxy evolution, are probable. Hence, in planning abundance surveys, it is important to carefully balance the need for large samples of stars against the spectral resolution and signal-to-noise ratio needed to obtain high-precision abundances. Furthermore, it is an advantage to work differentially on stars with similar atmospheric parameters, because then a simple 1D LTE analysis of stellar spectra may be sufficient. However, when determining high-precision absolute abundances or differential abundance between stars having more widely different parameters, e.g. metal-poor stars compared to the Sun or giants to dwarfs, then 3D non-LTE effects must be taken into account.     

On the occurrence and origin of anthropogenic radionuclides found in a fragment of the Chelyabinsk (LL5) meteorite

A piece of the 2013 Chelyabinsk meteorite was investigated for its content of anthropogenic radionuclides. In addition to traces of cesium‐137 that had been previously reported for this particular fragment, we found an unusually high amount of strontium‐90, which indicates that the source of this contamination was the Kyshtym accident (1957). A high Sr‐90/Cs‐137 activity ratio is characteristic for Kyshtym‐derived contaminations. Based on the cesium‐137 content in the soil from the finding site, it is estimated that the fragment was contaminated with soil particles in the milligram range upon impact. Investigation of the soil revealed very unusual ferromagnetic characteristics and an iron‐rich chemical composition. Mössbauer spectroscopy indicated the presence of steel components in this soil, suggesting that the investigated meteorite fragment was found in an industrial dumping site rather than natural soil.     

Effects of different types of mechanical disturbances on a charophyte dominated macrophyte community

Charophytes are unique and endangered macroalgae. Mechanical stress is known to damage charophyte communities, however, virtually nothing is known how such disturbance affect seasonal succession and recovery of charophyte communities. We examined the response of a shallow water charophyte community to different types and timing of disturbances by conducting seasonal in situ field experiment in a brackish water ecosystem. The disturbances associated with mechanical removal of vegetation or sediment had the largest impact on the charophyte community. Disturbances that only partly damaged the plants had minor effects and communities were characterized by fast recovery. The timing of disturbance had significant effects on the charophyte community with stronger effects when disturbed at the early stage compared to mid-stage of seasonal succession.     

Changes in Venice Lagoon dynamics due to construction of mobile barriers

The MoSE project (construction of mobile barrier to safeguard the Lagoon of Venice) entails changes to the structure of the lagoon's inlets. This could have consequences for the areas near the inlets and for the dynamics of the lagoon ecosystem as a whole. In order to predict the effects of the proposed alterations on the hydrodynamics of the lagoon, a well-tested hydrodynamic-dispersion model was applied. Simulations were carried out considering both idealised and realistic tide and wind scenarios.