Adaptive integration of constitutive rate equations

In finite element calculations the constitutive model plays a key role. The evaluation of the stress response of the constitutive relation for a given strain increment, which is a time integration in the case of models of the rate type, is a typical sub task in such calculations. Adaptive behaviour of the time integration is essential to assure numerical stability and to control the accuracy of the solution. An adaptive second order semi-implicit method is developed in this paper. Its numerical behaviour is compared with an adaptive second order explicit scheme. The two proposed methods control the local error and guarantee numerical stability of the time integration. We include several numerical geotechnical element tests using hypoplasticity with intergranular strain. The element tests simulate the behaviour of a finite element method based on the displacement formulation.     

Polyphase scheelite and stanniferous silicates in a W-(Sn) skarn close to Felbertal tungsten mine,Eastern Alps

Mineralogy and Petrology - The scheelite exploration target Messelingscharte (Eastern Tyrol, Austria) is located in vicinity of the world-class Felbertal tungsten deposit. W-(Sn) mineralisation...     

Early burial mud diapirism and its impact on stratigraphic architecture in the Carboniferous of the Shannon Basin,County Clare,Ireland

Delta fronts are often characterized by high rates of sediment supply that result in unstable slopes and a wide variety of soft‐sediment deformation, including the formation of overpressured and mobile muds that may flow plastically during early burial, potentially forming mud diapirs. The coastal cliffs of County Clare, western Ireland, expose Pennsylvanian (Namurian) delta‐front deposits of the Shannon Basin at large scale and in three dimensions. These deposits include decametre‐scale, internally chaotic mudstone masses that clearly impact the surrounding sedimentary strata. Evidence indicates that these were true mud (unlithified sediment) diapirs that pierced overlying strata. This study documents a well‐exposed ca 20 m tall mud diapir and its impact on the surrounding mouth‐bar deposits of the Tullig Cyclothem. A synsedimentary fault and associated rollover dome, evident from stratal thicknesses and the dip of the beds, define one edge of the diapir. These features are interpreted as recording the reactive rise of the mud diapir in response to extensional faulting along its margin. Above the diapir, heterolithic sandstones and siltstones contain evidence for the creation of localized accommodation, suggesting synsedimentary filling, tilting and erosion of a shallow sag basin accommodated by the progressive collapse of the diapir. Two other diapirs are investigated using three‐dimensional models built from ‘structure from motion’ drone imagery. Both diapirs are interpreted to have grown predominantly through passive rise (downbuilding). Stratal relationships for all three diapirs indicate that they were uncompacted and fluid‐rich mud beds that became mobilized through soft‐sediment deformation during early burial (i.e. <50 m, likely <10 m depth). Each diapir locally controlled the stratigraphic architecture in the shallow subsurface and potentially influenced local palaeocurrents on the delta. The mud diapirs studied herein are distinct from deeper ‘shale diapirs’ that have been inferred from seismic sections worldwide, now largely disputed.     

Effects of strong network modifiers on Fe<Superscript>3+</Superscript>/Fe<Superscript>2+</Superscript> in silicate melts: an experimental study

The effect of CaO, Na₂O, and K₂O on ferric/ferrous ratio in model multicomponent silicate melts was investigated in the temperature range 1450–1550?°C at 1-atm total pressure in air. It is demonstrated that the addition of these network modifier cations results in an increase of Fe³⁺/Fe²⁺ ratio. The influence of network modifier cations on the ferric/ferrous ratio increases in the order Ca?<?Na?<?K. Some old controversial conceptions concerning the effect of potassium on Fe³⁺/Fe²⁺ ratio in simple model liquids are critically evaluated. An empirical equation is proposed to predict the ferric/ferrous ratio in SiO₂–TiO₂–Al₂O₃–FeO–Fe₂O₃–MgO–CaO–Na₂O–K₂O–P₂O₅ melts at air conditions.     

Microstructural constraints on the mechanisms of the transformation to reidite in naturally shocked zircon

Zircon (ZrSiO₄) is used to study impact structures because it responds to shock loading and unloading in unique, crystallographically controlled manners. One such phenomenon is the transformation of zircon to the high-pressure polymorph, reidite. This study quantifies the geometric and crystallographic orientation relationships between these two phases using naturally shocked zircon grains. Reidite has been characterized in 32 shocked zircon grains (shocked to stages II and III) using a combination of electron backscatter diffraction (EBSD) and focused ion beam cross-sectional imaging techniques. The zircon-bearing clasts were obtained from within suevite breccia from the Nördlingen 1973 borehole, close to the center of the 14.4 Ma Ries impact crater, in Bavaria, Germany. We have determined that multiple sets (up to 4) of reidite lamellae can form in a variety of non-rational habit planes within the parent zircon. However, EBSD mapping demonstrates that all occurrences of lamellar reidite have a consistent interphase misorientation relationship with the host zircon that is characterized by an approximate alignment of a {100}_zircon with a {112}_reidite and alignment of a {112}_zircon with a conjugate {112}_reidite. Given the tetragonal symmetry of zircon and reidite, we predict that there are eight possible variants of this interphase relationship for reidite transformation within a single zircon grain. Furthermore, laser Raman mapping of one reidite-bearing grain shows that moderate metamictization can inhibit reidite formation, thereby highlighting that the transformation is controlled by zircon crystallinity. In addition to lamellar reidite, submicrometer-scale granules of reidite were observed in one zircon. The majority of reidite granules have a topotaxial alignment that is similar to the lamellar reidite, with some additional orientation dispersion. We confirm that lamellar reidite likely forms via a deviatoric transformation mechanism in highly crystalline zircon, whereas granular reidite forms via a reconstructive transformation from low-crystallinity ZrSiO₄ within the reidite stability field. The results of this study further refine the formation mechanisms and conditions of reidite transformation in naturally shocked zircon.     

Indian summer monsoon forcing on the deglacial polar cold reversals

The deglacial transition from the last glacial maximum at \(\sim \)20 kiloyears before present (ka) to the Holocene (11.7 ka to Present) was interrupted by millennial-scale cold reversals, viz., Antarctic Cold Reversal (\(\sim \)14.5–12.8 ka) and Greenland Younger Dryas (\(\sim \)12.8–11.8 ka) which had different timings and extent of cooling in each hemisphere. The cause of this synchronously initiated, but different hemispheric cooling during these cold reversals (Antarctic Cold Reversal \(\sim \)3\(^{\circ }\hbox {C}\) and Younger Dryas \(\sim \)10\(^{\circ }\hbox {C}\)) is elusive because \(\hbox {CO}_{2}\), the fundamental forcing for deglaciation, and Atlantic meridional overturning circulation, the driver of antiphased bipolar climate response, both fail to explain this asymmetry. We use centennial-resolution records of the local surface water \(\delta ^{18}\hbox {O}\) of the Eastern Arabian Sea, which constitutes a proxy for the precipitation associated with the Indian Summer Monsoon, and other tropical precipitation records to deduce the role of tropical forcing in the polar cold reversals. We hypothesize a mechanism for tropical forcing, via the Indian Summer Monsoons, of the polar cold reversals by migration of the Inter-Tropical Convergence Zone and the associated cross-equatorial heat transport.     

Auklet (Charadriiformes: Alcidae, Aethia spp.) chick meals from the Aleutian Islands, Alaska, have a very low incidence of plastic marine debris

The ingestion of plastic marine debris is a chronic problem for some of the world’s seabird species, contributing to reduced chick survival, population declines, and deposition of contaminants via absorption in birds’ gastrointestinal tract. We analysed the frequency of ingested plastic in chick meals delivered by adults in four species of auklet - Crested (Aethia cristatella), Least (A. pusilla), Parakeet (A. psittacula), and Whiskered (A. pygmaea) - from three breeding colonies in the Aleutian Islands, Alaska, USA over a 14-year period from 1993 to 2006. Among 2541 chick meals, we found plastic in only one - from a Whiskered Auklet on Buldir Island in 1993. While adult Parakeet Auklets have a high frequency of plastic ingestion (over 90%), no chick meals contained plastic. Unlike other seabirds, the planktivorous auklets do not appear to offload plastic to their chicks, and we conclude that auklet chicks are probably at a low risk of contamination from plastic debris.     

Variability of baroclinic tidal currents on the Mackenzie Shelf,the Southeastern Beaufort Sea

Semidiurnal tidal currents on the outer shelf of the Mackenzie Shelf in the Beaufort Sea were found to be strongly influenced by the locally generated baroclinic tide. Two primary factors are involved in this process: (1) the sharp shelf break along the northeastern Mackenzie Shelf, promoting the generation of vigorous internal tidal waves; and (2) the proximity to critical latitudes for M₂ and N₂ motions locking these waves and preventing them from leaving the source region. As a result, internal tides are resonantly trapped between the shelf and critical latitudes. The physical properties and temporal variations of tidal motions were examined using current meter measurements obtained from 1987–1988 at four sites (SS1, SS2, SS3, and SS4) offshore of the shelf break at depths of ∼200 m. Each mooring had Aanderaa RCM4s positioned at ∼35 m below the surface and ∼50 m above the bottom. Complex demodulation was used to compute the envelopes (amplitude modulation) of these components. A striking difference in the variability of clockwise (CW) and counterclockwise (CCW) tidal currents was found. The CW tides are highly variable, have greater amplitude, exhibit a burst-like character associated with wind events and contain about 80% of the total energy of the semidiurnal tidal currents. In contrast, the CCW components have a more regular temporal regime with distinct monthly, fortnightly and 10-day modulation at astronomical periodicities associated with frequency differences M₂–N₂ (0.03629 cpd), S₂–M₂ (0.06773 cpd), and S₂–N₂ (0.10402 cpd). Significant horizontal correlation of the CW current envelopes was found only between stations near the northeast Mackenzie Shelf, indicating this to be the main area of baroclinic internal wave generation.     

Uncertainty of moment tensor determination from surface wave analysis for shallow earthquakes

We treat the problem of the non-uniqueness of the moment tensor inversion for shallow earthquakes from long period surface wave data. We give an existence condition for double couples radiating the same long period surface waves as the deviatoric moment tensor obtained by such inversion. We describe the family of such double couples and show that they may provide better estimates of double couple mechanisms than the traditional “best double-couple” solution. We demonstrate the results of application to four large shallow earthquakes with different source signatures.