@article{BilbaoLasaJaraMunozPedojaetal.2020, author = {Bilbao-Lasa, Peru and Jara-Mu{\~n}oz, Julius and Pedoja, Kevin and {\´A}lvarez, Irantzu and Aranburu, Arantza and Iriarte, Eneko and Galparsoro, Ibon}, title = {Submerged marine terraces identification and an approach for numerical modeling the sequence formation in the Bay of Biscay (Northeastern Iberian Peninsula)}, series = {Frontiers in Earth Science}, volume = {8}, journal = {Frontiers in Earth Science}, number = {47}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2020.00047}, pages = {1 -- 20}, year = {2020}, abstract = {Submerged sequences of marine terraces potentially provide crucial information of past sea-level positions. However, the distribution and characteristics of drowned marine terrace sequences are poorly known at a global scale. Using bathymetric data and novel mapping and modeling techniques, we studied a submerged sequence of marine terraces in the Bay of Biscay with the objective to identify the distribution and morphologies of submerged marine terraces and the timing and conditions that allowed their formation and preservation. To accomplish the objectives a high-resolution bathymetry (5 m) was analyzed using Geographic Information Systems and TerraceM(R). The successive submerged terraces were identified using a Surface Classification Model, which linearly combines the slope and the roughness of the surface to extract fossil sea-cliffs and fossil rocky shore platforms. For that purpose, contour and hillshaded maps were also analyzed. Then, shoreline angles, a geomorphic marker located at the intersection between the fossil sea-cliff and platform, were mapped analyzing swath profiles perpendicular to the isobaths. Most of the submerged strandlines are irregularly preserved throughout the continental shelf. In summary, 12 submerged terraces with their shoreline angles between approximately: -13 m (T1), -30 and -32 m (T2), -34 and 41 m (T3), -44 and -47 m (T4), -49 and 53 m (T5), -55 and 58 m (T6), -59 and 62 m (T7), -65 and 67 m (T8), -68 and 70 m (T9), -74 and -77 m (T10), -83 and -86 m (T11) and -89 and 92 m (T12). Nevertheless, the ones showing the best lateral continuity and preservation in the central part of the shelf are T3, T4, T5, T7, T8, and T10. The age of the terraces has been estimated using a landscape evolution model. To simulate the formation and preservation of submerged terraces three different scenarios: (i) 20-0 ka; (ii) 128-0 ka; and (iii) 128-20 ka, were compared. The best scenario for terrace generation was between 128 and 20 Ka, where T3, T5, and T7 could have been formed.}, language = {en} } @article{JaraMunozMelnickPedojaetal.2019, author = {Jara-Munoz, Julius and Melnick, Daniel and Pedoja, Kevin and Strecker, Manfred}, title = {TerraceM-2: A MatlabR (R) Interface for Mapping and Modeling Marine and Lacustrine Terraces}, series = {Frontiers in Earth Science}, volume = {7}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Research Foundation}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2019.00255}, pages = {18}, year = {2019}, abstract = {The morphology of marine and lacustrine terraces has been largely used to measure past sea- and lake-level positions and estimate vertical deformation in a wealth of studies focused on climate and tectonic processes. To obtain accurate morphometric assessments of terrace morphology we present TerraceM-2, an improved version of our MatlabR (R) graphic-user interface that provides new methodologies for morphometric analyses as well as landscape evolution and fault-dislocation modeling. The new version includes novel routines to map the elevation and spatial distribution of terraces, to model their formation and evolution, and to estimate fault-slip rates from terrace deformation patterns. TerraceM-2 has significantly improves its processing speed and mapping capabilities, and includes separate functions for developing customized workflows beyond the graphic-user interface. We illustrate these new mapping and modeling capabilities with three examples: mapping lacustrine shorelines in the Dead Sea to estimate deformation across the Dead Sea Fault, landscape evolution modeling to estimate a history of uplift rates in southern Peru, and dislocation modeling of deformed marine terraces in California. These examples also illustrate the need to use topographic data of different resolutions. The new modeling and mapping routines of TerraceM-2 highlight the advantages of an integrated joint mapping and modeling approach to improve the efficiency and precision of coastal terrace metrics in both marine and lacustrine environments.}, language = {en} } @article{JaraMunozMelnickZambranoetal.2017, author = {Jara-Munoz, Julius and Melnick, Daniel and Zambrano, Patricio and Rietbrock, Andreas and Gonzalez, Javiera and Argandona, Boris and Strecker, Manfred}, title = {Quantifying offshore fore-arc deformation and splay-fault slip using drowned Pleistocene shorelines, Arauco Bay, Chile}, series = {Journal of geophysical research : Solid earth}, volume = {122}, journal = {Journal of geophysical research : Solid earth}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1002/2016JB013339}, pages = {4529 -- 4558}, year = {2017}, abstract = {Most of the deformation associated with the seismic cycle in subduction zones occurs offshore and has been therefore difficult to quantify with direct observations at millennial timescales. Here we study millennial deformation associated with an active splay-fault system in the Arauco Bay area off south central Chile. We describe hitherto unrecognized drowned shorelines using high-resolution multibeam bathymetry, geomorphic, sedimentologic, and paleontologic observations and quantify uplift rates using a Landscape Evolution Model. Along a margin-normal profile, uplift rates are 1.3m/ka near the edge of the continental shelf, 1.5m/ka at the emerged Santa Maria Island, -0.1m/ka at the center of the Arauco Bay, and 0.3m/ka in the mainland. The bathymetry images a complex pattern of folds and faults representing the surface expression of the crustal-scale Santa Maria splay-fault system. We modeled surface deformation using two different structural scenarios: deep-reaching normal faults and deep-reaching reverse faults with shallow extensional structures. Our preferred model comprises a blind reverse fault extending from 3km depth down to the plate interface at 16km that slips at a rate between 3.0 and 3.7m/ka. If all the splay-fault slip occurs during every great megathrust earthquake, with a recurrence of similar to 150-200years, the fault would slip similar to 0.5m per event, equivalent to a magnitude similar to 6.4 earthquake. However, if the splay-fault slips only with a megathrust earthquake every similar to 1000years, the fault would slip similar to 3.7m per event, equivalent to a magnitude similar to 7.5 earthquake.}, language = {en} } @article{JaraMunozMelnickBrilletal.2015, author = {Jara-Munoz, Julius and Melnick, Daniel and Brill, Dominik and Strecker, Manfred}, title = {Segmentation of the 2010 Maule Chile earthquake rupture from a joint analysis of uplifted marine terraces and seismic-cycle deformation patterns}, series = {Quaternary science reviews : the international multidisciplinary research and review journal}, volume = {113}, journal = {Quaternary science reviews : the international multidisciplinary research and review journal}, publisher = {Elsevier}, address = {Oxford}, issn = {0277-3791}, doi = {10.1016/j.quascirev.2015.01.005}, pages = {171 -- 192}, year = {2015}, abstract = {The segmentation of major fault systems in subduction zones controls earthquake magnitude and location, but the causes for the existence of segment boundaries and the relationships between long-term deformation and the extent of earthquake rupture, are poorly understood. We compare permanent and seismic-cycle deformation patterns along the rupture zone of the 2010 Maule earthquake (M8.8), which ruptured 500 km of the Chile subduction margin. We analyzed the morphology of MIS-5 marine terraces using LiDAR topography and established their chronology and coeval origin with twelve luminescence ages, stratigraphy and geomorphic correlation, obtaining a virtually continuous distribution of uplift rates along the entire rupture zone. The mean uplift rate for these terraces is 0.5 m/ka. This value is exceeded in three areas, which have experienced rapid emergence of up to 1.6 m/ka; they are located at the northern, central, and southern sectors of the rupture zone, referred to as Topocalma, Carranza and Arauco, respectively. The three sectors correlate with boundaries of eight great earthquakes dating back to 1730. The Topocalma and Arauco sectors, located at the boundaries of the 2010 rupture, consist of broad zones of crustal warping with wavelengths of 60 and 90 km, respectively. These two regions coincide with the axes of oroclinal bending of the entire Andean margin and correlate with changes in curvature of the plate interface. Rapid uplift at Carranza, in turn, is of shorter wavelength and associated with footwall flexure of three crustal-scale normal faults. The uplift rate at Carranza is inversely correlated with plate coupling as well as with coseismic slip, suggesting permanent deformation may accumulate interseismically. We propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust resulting in barriers to the propagation of great earthquakes. Slip during the 1960 (M9.5) and 2010 events overlapped with the similar to 90-km-long zone of rapid uplift at Arauco; similarly, slip in 2010 and 1906 extended across the similar to 60-km-long section of the megathrust at Topocalma, but this area was completely breached by the 1730 (M similar to 9) event, which propagated southward until Carranza. Both Arauco and Topocalma show evidence of sustained rapid uplift since at least the middle Pleistocene. These two sectors might thus constitute discrete seismotectonic boundaries restraining most, but not all great earthquake ruptures. Based on our observations, such barriers might be breached during multi-segment super-cycle events. (C) 2015 Elsevier Ltd. All rights reserved.}, language = {en} }