Cretaceous deformation of the southern Central Andes: synorogenic growth strata in the Neuquén Group (35° 30′–37° S)

The Neuquén Group is an Upper Cretaceous continental sedimentary unit exhumed during the latest Miocene contractional phase occurred in the southern Central Andes, allowing a direct field observation and study of the depositional geometries. The identification of growth strata on these units surrounding the structures of the frontal parts of the Andes, sedimentological analyses and U–Pb dating of detrital components, allowed the definition of a synorogenic unit that coexisted with the uplift of the early Andean orogen since ca. 100 Ma, maximum age obtained in this work, compatible with previous assignments and constrained in the top by the deposition of the Malargüe Group, in the Maastrichtian (ca. 72 Ma). The definition of a wedge top area in this foreland basin system, where growth strata were described, permitted to identify a Late Cretaceous orogenic front and foredeep area, whose location and amplitude contrast with previous hypotheses. This wedge top area was mostly fed from the paleo‐Andes with small populations coming from sources in the cratonic area that are interpreted as a recycling in Jurassic and Lower Cretaceous sections, which contrasts with other analyses performed at the foredeep zone that have mixed sources. In particular, Permian sources are interpreted as coming directly from the cores of the basement structures, where Neopaleozoic sections are exposed, next to the synorogenic sedimentation, implying a strong incision in Late Cretaceous times with an exhumation structural level similar to the present. The maximum recognised advance for this Late Cretaceous deformation in the study area is approximately 500 km east of the Pacific trench, which constitutes an anomaly compared with neighbour segments where Late Cretaceous deformations were found considerably retracted. The geodynamic context of the sedimentation of this unit is interpreted as produced under the westward fast moving of South America, colliding with two consecutive mid‐ocean ridges during a period of important plate reorganisation. The subduction of young, anhydrous, buoyant lithosphere would have produced changes in the subduction geometry, reflected first by an arc waning/gap and subsequently by an arc migration that coexisted with synorogenic sedimentation. These magmatic and deformational processes would be the product of a shallow subduction regime, following previous proposals, which occurred in Late Cretaceous times, synchronous to the sedimentation of the Neuquén Group.     

Nonlinear FE model updating and reconstruction of the response of an instrumented seismic isolated bridge to the 2010 Maule Chile earthquake

Nonlinear finite element (FE) modeling has been widely used to investigate the effects of seismic isolation on the response of bridges to earthquakes. However, most FE models of seismic isolated bridges (SIB) have used seismic isolator models calibrated from component test data, while the prediction accuracy of nonlinear FE models of SIB is rarely addressed by using data recorded from instrumented bridges. In this paper, the accuracy of a state‐of‐the‐art FE model is studied through nonlinear FE model updating (FEMU) of an existing instrumented SIB, the Marga‐Marga Bridge located in Viña del Mar, Chile. The seismic isolator models are updated in 2 phases: component‐wise and system‐wise FEMU. The isolator model parameters obtained from 23 isolator component tests show large scatter, and poor goodness of fit of the FE‐predicted bridge response to the 2010 Mw 8.8 Maule, Chile Earthquake is obtained when most of those parameter sets are used for the isolator elements of the bridge model. In contrast, good agreement is obtained between the FE‐predicted and measured bridge response when the isolator model parameters are calibrated using the bridge response data recorded during the mega‐earthquake. Nonlinear FEMU is conducted by solving single‐ and multiobjective optimization problems using high‐throughput cloud computing. The updated FE model is then used to reconstruct response quantities not recorded during the earthquake, gaining more insight into the effects of seismic isolation on the response of the bridge during the strong earthquake.     

Climatic signals in lacustrine deposits of the Upper Yacoraite Formation,Western Argentina: Evidence from clay minerals,analcime, dolomite and fibrous calcite

Authigenic clays are an important control on reservoir quality in lacustrine carbonates but remain challenging to predict. Lacustrine depositional systems respond to climatic variations in rainfall, surface runoff and groundwater input, and evaporation, and result in rapid and frequent changes in lake volume; this is expressed through changing water depth and shoreline position. In the upper portion of the Early Palaeocene Yacoraite Formation of the Salta Basin in Argentina, extensive lacustrine deposits were deposited during the sag phase of rifting. Prior high-resolution stratigraphic studies have suggested that climatic factors control microbial carbonate sequences within a ‘balanced fill’ lake, with variation in the lake level having a major influence on facies association changes. This study characterizes the evolution of facies and mineralogy within the Yacoraite Formation, focusing on the distribution of clay minerals, making a link between the high, medium and low-frequency sequence stratigraphic cycles. The low-frequency transgressive hemicycle of the upper portion of the Yacoraite Formation is comprised of abundant siliciclastic facies, suggesting a wetter period. Microbialites occurring in this interval are coarse-grained and agglutinated. Detrital clay minerals such as illite and chlorite and associated siliciclastic sediments were input to the lake during high-frequency transgressive periods. During high-frequency regressive hemicycles, sedimentation was dominated by carbonate facies with Ca-rich dolomite and the authigenic clays are comprised of chlorite/smectite mixed-layers. By contrast, the low frequency regressive hemicycle records fine-grained agglutinated microbialite with horizons of fibrous calcite, more stoichiometric dolomite, barite and authigenic magnesian smectite. This indicates elevated ion concentrations in the lake under intense evaporation during an arid period. Understanding the conditions that are favourable for formation and preservation of authigenic clays within the lacustrine environment can improve understanding of reservoir quality in comparable economically important deposits.     

Digital mapping of coastal boulders – high-resolution data acquisition to infer past and recent transport dynamics

Coastal boulder fields provide clues to long-term frequency-magnitude patterns of coastal flooding events and have the potential to play an important role in coastal hazard assessment. Mapping boulders in the field is time and labour-intensive, and work on intertidal reef platforms, as in the present study, is physically challenging. By addressing coastal scientists who are not specialists in remote sensing, this contribution reports on the possibilities and limitations of digital applications in boulder mapping in Eastern Samar, Philippines, where recent supertyphoons Haiyan and Hagupit induced high waves, coastal flooding and boulder transport. It is demonstrated how satellite imagery of sub-metre resolution (from Pléiades and WorldView-3 imagery) enables efficient analysis of transport vectors and distances of larger boulders, reflecting variation in latitudes of both typhoon tracks and approaching angles of typhoon-generated waves. During the investigated events, boulders with a-axes of up to 8 m were clearly identified to have been shifted for up to 32 m, mostly along the seaward margin of the boulder field. It is, however, hard to keep track of smaller boulders, and the length of a-axes and b-axes including their orientation is often impossible to map with sufficient accuracy. Orthophotographs and digital surface models created through the application of an unmanned aerial vehicle and the ‘Structure from Motion’ technique provide ultra-high-resolution data, and have the potential to not only improve the results of satellite image analysis, but also those from field mapping and may significantly reduce overall time in the field. Orthophotographs permit unequivocal mapping of a-axes and b-axes including their orientation, while precise values for c-axes can be derived from the respective digital surface models. Volume of boulders is best inferred from boulder-specific Structure from Motion-based three-dimensional models. Battery power, flight speed and altitude determine the limits of the area covered, while patches shielded by the boulders are difficult to resolve. For some tasks, field mapping remains mandatory and cannot be replaced by currently available remote sensing tools: for example, sampling for rock type, density and age dating, recording of lithological separation of boulders from the underlying geological unit and of geomorphic features on a millimetre to decimetre-scale, or documentation of fine-grained sediment transport in between the boulders in supratidal settings. In terms of future events, the digital products presented here will provide a valuable reference to track boulder transport on a centimetre to decimetre-scale and to better understand the hydrodynamics of extreme-wave events on a fringing reef coastline.     

Adjustability of exposed elements by updating their capacity for resistance after a damaging event: application to an earthquake–tsunami cascade scenario

Natural Hazards - The 2011 Great East Japan earthquake and tsunami remind us once again that these types of cascade event can occur and cause considerable damage. The scientific community realizes...     

Montecarlo DLA-type simulations of wetting effects in fluid displacement in porous media

In this work, we report diffusion-limited aggregation (DLA)-type Montecarlo computations of a stochastic model of displacement of a viscous fluid by another that preferentially wets a porous medium, for the case when both fluids are immiscible in the absence of buoyancy forces. The model has the aim to simulate cooperative invasion processes found in experiments of immiscible wetting displacement. The model considers the nonlocal effects of the Laplacian pressure field and the capillary forces via hydrodynamic equations in the Darcy regime with a boundary condition for the pressure at the interface. The boundary condition contains two different types of disorder: the capillary term, which constitutes an additive random disorder, and a term containing an effective random surface tension, which couples to a curvature (it constitutes a multiplicative random term that carries nonlocal information of the whole pressure). We generate different displacement patterns for different setting of the parameters of the model. We analyze these patterns by studying the scaling properties of the interface that separate the two fluids and calculating the fractal dimension of the interface. The results show the existence of three distinct regimes of scaling. One regime at the smallest-length scales is due to the multiplicative random disorder together with the nonlocal coupling; it reveals itself in a roughness exponent α ≈ 0.80. Additionally, we find a DLA-type scaling regime with a roughness exponent α ≈ 0.60 at the largest scales and intermediate scaling regime with α ≈ 0.70 corresponding to invasion percolation with trapping. Each regime has definite scaling ranges that depend on the capillary number and the relative wetting tendency of the fluids. The behavior of the fractal dimensions of the interfaces of the aggregates constitutes a further confirmation of the existence of three scaling regimes and the multi-self-affinity of the perimeter of the interface boundaries.     

Effective sea ice area based on a thickness threshold

As Arctic sea ice declines in response to climate change, a shift from thick multiyear ice to a thinner ice cover is occurring. With this transition, ice thicknesses approach a threshold below which ice no longer insulates the atmosphere from oceanic surface fluxes. While this is well known, there are no estimates of the magnitude of this threshold, nor of the proportion of sea ice area that is below this threshold as ice thins. We determine this threshold by simulating the atmospheric response to varying thicknesses, ranging from 0.0 to 2.0 m and determine that threshold to be 0.40–0.50 m. The resulting “effective” ice area is 4–14% lower than reported total ice area, as 0.39–0.97 × 10⁶ km² of the total ice area falls below the threshold throughout the twentieth century, including during notable ice minima. The atmosphere above large non-insulating ice-covered regions is susceptible to more than 2 °C of warming despite ice presence. Observed mean Arctic Ocean ice thickness is projected to fall below this threshold as early as the mid-2020s. Studies on ocean–atmosphere interactions in relation to sea ice area should focus on this insulating sea ice area, where ice is at least 0.40–0.50 m thick, and treat ice regions below 0.40–0.50 m thickness with caution.