Fast numerics for the spin orbit equation with realistic tidal dissipation and constant eccentricity

We present an algorithm for the rapid numerical integration of a time-periodic ODE with a small dissipation term that is \(C^1\) in the velocity. Such an ODE arises as a model of spin–orbit coupling in a star/planet system, and the motivation for devising a fast algorithm for its solution comes from the desire to estimate probability of capture in various solutions, via Monte Carlo simulation: the integration times are very long, since we are interested in phenomena occurring on timescales of the order of \(10^6\)–\(10^7\) years. The proposed algorithm is based on the high-order Euler method which was described in Bartuccelli et al. (Celest Mech Dyn Astron 121(3):233–260, 2015), and it requires computer algebra to set up the code for its implementation. The payoff is an overall increase in speed by a factor of about 7.5 compared to standard numerical methods. Means for accelerating the purely numerical computation are also discussed.     

Meteorite impacts in the ocean: the danger of tsunamis on the coast of Buenos Aires Province,Argentina

Natural Hazards - Comets, meteorites, or asteroids impacting against the Earth are not unusual events. Such impacts on the ocean could produce tsunamis which can reach coastal areas. This paper...     

Prediction of inter-storey drifts for regular RC structures with masonry infills based on bare frame modelling

The traditional construction of masonry infills adjacent to RC structural elements is still widely adopted in European countries, including seismically active regions. Given the repeated field observations from damaging earthquakes, pointing to unacceptably high levels of masonry infill damage, the present study is motivated by the need to improve further the European seismic design approach for new RC structures with masonry infills, in order to exclude the poor seismic behaviour probably caused by deficiencies in the verification procedure. Since the in-plane damage to non-structural panels is commonly controlled through the limitation of inter-storey drifts, the possibility to introduce more effective verification criteria, accounting for structural properties, infill layouts and masonry properties is explored. Therefore, starting from the assumption that analyses and verifications in the design of buildings are commonly accomplished neglecting the presence of infills, results of extensive nonlinear numerical analyses for different building configurations are examined. As a result, a simplified procedure for the prediction of expected inter-storey drifts for infilled structures, based on the corresponding demands of bare configurations, in function of a simple parameter accounting for structural properties and the presence of infills, is introduced. Possible implications of the proposed approach aimed at the improvement of the current design provisions are discussed.     

Channelization of a large Alpine river: what is left of its original morphodynamics?

The Adige River drains 12 200 km² of the Eastern Alps and flows for 213 km within this mountain range. Similar to other large rivers in Central Europe, the Adige River was subject to massive channelization works during the 19th century. Thanks to the availability of several historical maps, this river represents a very valuable case study to document the extent to which the morphology of the river changed due to channelization and to understand how much is left of its original morphodynamics. The study was based on the analysis of seven sets of historical maps dating from 1803–1805 to 1915–1927, on geomorphological analysis, on the application of mathematical morphodynamic theories and on the application of bar and channel pattern prediction models. The study concerns 115 km of the main stem and 29 km of its tributaries. In the pre‐channelization conditions, the Adige River presented a prevalence of single‐thread channel planforms. Multi‐thread patterns developed only immediately downstream of the main confluences. During the 19th century, the Adige underwent considerable channel adjustment, consisting of channel narrowing, straightening, and reduction of bars and islands. Multi‐thread and single‐thread reaches evolved through different evolutionary trajectories, considering both the channel width and the bar/vegetation interaction. Bar and channel pattern predictors showed good correspondence with the observed patterns, including the development of multi‐thread morphologies downstream of the confluences. Application of the free‐bar predictor helped to interpret the strong reduction – almost complete loss – of exposed sediment bars after the channelization works, quantifying the riverbed inclination to form alternate bars. This morphological evolution can be observed in other Alpine rivers of similar size and similar massive channelization, therefore, a simplified conceptual model for large rivers subjected to channelization is proposed, showing that a relatively small difference in the engineered channel width may have a strong impact on the river dynamics, specifically on bar formation. Copyright © 2017 John Wiley & Sons, Ltd.     

Pliocene clays from Aegina (Greece): Reference material for chemical provenance studies on bronze age pottery from the Island

Normally, the use of clays as reference materials in chemical provenance studies of ancient ceramics is complicated due to the original clay paste processing. The primary mixing and/or refining of raw materials during pottery production makes a straightforward comparison of archaeological ceramics with extant geological materials difficult if not impossible in many cases. However, in the case of Pliocene clays from Aegina (Greece), which were examined chemically and mineralogically and compared with Bronze Age pottery produced on the island, a successful exception can be presented. The chemical composition of a large group of Aeginetan pottery resembles the chemical composition of clays from a deposit in close vicinity to the main Bronze Age settlement of the island. Clays from specific outcrops exhibit considerable chemical and mineralogical homogeneity, and the suitability of those clays for pottery production apparently made substantial clay paste processing unnecessary. © 2004 Wiley Periodicals, Inc.     

Evidence of multiple thermokarst lake generations from an 11 800‐year‐old permafrost core on the northern Seward Peninsula,Alaska

Permafrost degradation influences the morphology, biogeochemical cycling and hydrology of Arctic landscapes over a range of time scales. To reconstruct temporal patterns of early to late Holocene permafrost and thermokarst dynamics, site‐specific palaeo‐records are needed. Here we present a multi‐proxy study of a 350‐cm‐long permafrost core from a drained lake basin on the northern Seward Peninsula, Alaska, revealing Lateglacial to Holocene thermokarst lake dynamics in a central location of Beringia. Use of radiocarbon dating, micropalaeontology (ostracods and testaceans), sedimentology (grain‐size analyses, magnetic susceptibility, tephra analyses), geochemistry (total nitrogen and carbon, total organic carbon, δ¹³C_org) and stable water isotopes (δ¹⁸O, δD, d excess) of ground ice allowed the reconstruction of several distinct thermokarst lake phases. These include a pre‐lacustrine environment at the base of the core characterized by the Devil Mountain Maar tephra (22 800±280 cal. a BP, Unit A), which has vertically subsided in places due to subsequent development of a deep thermokarst lake that initiated around 11 800 cal. a BP (Unit B). At about 9000 cal. a BP this lake transitioned from a stable depositional environment to a very dynamic lake system (Unit C) characterized by fluctuating lake levels, potentially intermediate wetland development, and expansion and erosion of shore deposits. Complete drainage of this lake occurred at 1060 cal. a BP, including post‐drainage sediment freezing from the top down to 154 cm and gradual accumulation of terrestrial peat (Unit D), as well as uniform upward talik refreezing. This core‐based reconstruction of multiple thermokarst lake generations since 11 800 cal. a BP improves our understanding of the temporal scales of thermokarst lake development from initiation to drainage, demonstrates complex landscape evolution in the ice‐rich permafrost regions of Central Beringia during the Lateglacial and Holocene, and enhances our understanding of biogeochemical cycles in thermokarst‐affected regions of the Arctic.     

Pingüino In-bearing polymetallic vein deposit,Deseado Massif,Patagonia, Argentina: characteristics of mineralization and ore-forming fluids

The Pingüino deposit, located in the low sulfidation epithermal metallogenetical province of the Deseado Massif, Patagonia, Argentina, represents a distinct deposit type in the region. It evolved through two different mineralization events: an early In-bearing polymetallic event that introduced In, Zn, Pb, Ag, Cd, Au, As, Cu, Sn, W and Bi represented by complex sulfide mineralogy, and a late Ag–Au quartz-rich vein type that crosscut and overprints the early polymetallic mineralization. The indium-bearing polymetallic mineralization developed in three stages: an early Cu–Au–In–As–Sn–W–Bi stage (Ps₁), a Zn–Pb–Ag–In–Cd–Sb stage (Ps₂) and a late Zn–In–Cd (Ps₃). Indium concentrations in the polymetallic veins show a wide range (3.4 to 1,184 ppm In). The highest indium values (up to 1,184 ppm) relate to the Ps₂ mineralization stage, and are associated with Fe-rich sphalerites, although significant In enrichment (up to 159 ppm) is also present in the Ps₁ paragenesis associated with Sn-minerals (ferrokesterite and cassiterite). The hydrothermal alteration associated with the polymetallic mineralization is characterized by advanced argillic alteration within the immediate vein zone, and sericitic alteration enveloping the vein zone. Fluid inclusion studies indicate homogenisation temperatures of 308.2–327°C for Ps₁ and 255–312.4°C for Ps₂, and low to moderate salinities (2 to 5 eq.wt.% NaCl and 4 to 9 eq.wt.% NaCl, respectively). δ³⁴S values of sulfide minerals (+0.76‰ to +3.61‰) indicate a possible magmatic source for the sulfur in the polymetallic mineralization while Pb isotope ratios for the sulfides and magmatic rocks (²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb ratios of 17.379 to 18.502; 15.588 to 15.730 and 38.234 to 38.756, respectively) are consistent with the possibility that the Pb reservoirs for both had the same crustal source. Spatial relationships, hydrothermal alteration styles, S and Pb isotopic data suggest a probable genetic relation between the polymetallic mineralization and dioritic intrusions that could have been the source of metals and hydrothermal fluids. Mineralization paragenesis, alteration mineralogy, geochemical signatures, fluid inclusion data and isotopic data, confirm that the In-bearing polymetallic mineralization from Pingüino deposit is a distinct type, in comparison with the well-known epithermal low sulfidation mineralization from the Deseado Massif.     

Simulation of marine controlled source electromagnetic measurements using a parallel fourier <Emphasis Type="Italic">hp</Emphasis>-finite element method

We introduce a new numerical method to simulate geophysical marine controlled source electromagnetic (CSEM) measurements for the case of 2D structures and finite 3D sources of electromagnetic (EM) excitation. The method of solution is based on a spatial discretization that combines a 1D Fourier transform with a 2D self-adaptive, goal-oriented, hp-Finite element method. It enables fast and accurate simulations for a variety of important, challenging and practical cases of marine CSEM acquisition. Numerical results confirm the high accuracy of the method as well as some of the main physical properties of marine CSEM measurements such as high measurement sensitivity to oil-bearing layers in the subsurface. In our model, numerical results indicate that measurements could be affected by the finite oil-bearing layer by as much as 10⁴% (relative difference). While the emphasis of this paper is on EM simulations, the method can be used to simulate different physical phenomena such as seismic measurements.     

Modelling of fractured carbonate reservoirs: outline of a novel technique via a case study from the Molasse Basin, southern Bavaria, Germany

Fluid flow in low-permeable carbonate rocks depends on the density of fractures, their interconnectivity and on the formation of fault damage zones. The present-day stress field influences the aperture hence the transmissivity of fractures whereas paleostress fields are responsible for the formation of faults and fractures. In low-permeable reservoir rocks, fault zones belong to the major targets. Before drilling, an estimate for reservoir productivity of wells drilled into the damage zone of faults is therefore required. Due to limitations in available data, a characterization of such reservoirs usually relies on the use of numerical techniques. The requirements of these mathematical models encompass a full integration of the actual fault geometry, comprising the dimension of the fault damage zone and of the fault core, and the individual population with properties of fault zones in the hanging and foot wall and the host rock. The paper presents both the technical approach to develop such a model and the property definition of heterogeneous fault zones and host rock with respect to the current stress field. The case study describes a deep geothermal reservoir in the western central Molasse Basin in southern Bavaria, Germany. Results from numerical simulations indicate that the well productivity can be enhanced along compressional fault zones if the interconnectivity of fractures is lateral caused by crossing synthetic and antithetic fractures. The model allows a deeper understanding of production tests and reservoir properties of faulted rocks.     

The use of low frequencies in a full‐waveform inversion and impedance inversion land seismic case study

Velocity model building and impedance inversion generally suffer from a lack of intermediate wavenumber content in seismic data. Intermediate wavenumbers may be retrieved directly from seismic data sets if enough low frequencies are recorded. Over the past years, improvements in acquisition have allowed us to obtain seismic data with a broader frequency spectrum. To illustrate the benefits of broadband acquisition, notably the recording of low frequencies, we discuss the inversion of land seismic data acquired in Inner Mongolia, China. This data set contains frequencies from 1.5–80 Hz. We show that the velocity estimate based on an acoustic full‐waveform inversion approach is superior to one obtained from reflection traveltime inversion because after full‐waveform inversion the background velocity conforms to geology. We also illustrate the added value of low frequencies in an impedance estimate.