Next Article in Journal
Application of Experimental Configurations of Seismic and Electric Tomographic Techniques to the Investigation of Complex Geological Structures
Next Article in Special Issue
Investigation of the Effect of Integrated Offset, GPS, and InSAR Data in the Stochastic Source Modeling of the 2002 Denali Earthquake
Previous Article in Journal
Application of the Coastal and Marine Ecological Classification Standard (CMECS) to Create Benthic Geologic Habitat Maps for Portions of Acadia National Park, Maine, USA
Previous Article in Special Issue
Conservative Evaluation of Fault Displacement Hazard for a Nuclear Site in Case of Insufficient Data on the Fault Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seismotectonic Setting of the Andes along the Nazca Ridge Subduction Transect: New Insights from Thermal and Finite Element Modelling

1
Dipartimento di Scienze Pure e Applicate (DiSPeA), Università di Urbino “Carlo Bo”, Via Aurelio Saffi, 2, 61029 Urbino, Italy
2
Scuola di Scienze e Tecnologie, Sezione di Geologia, Università degli Studi di Camerino, Via Gentile III da Varano, 7, 62032 Camerino, Italy
3
Istituto Nazionale di Geofisica e Vulcanologia, Via di Vigna Murata, 605, 00143 Rome, Italy
*
Author to whom correspondence should be addressed.
Geosciences 2024, 14(10), 257; https://doi.org/10.3390/geosciences14100257
Submission received: 20 July 2024 / Revised: 22 September 2024 / Accepted: 25 September 2024 / Published: 28 September 2024
(This article belongs to the Special Issue New Trends in Earthquake Engineering and Seismotectonics)

Abstract

:
The structural evolution of Andean-type orogens is strongly influenced by the geometry of the subducting slab. This study focuses on the flat-slab subduction of the Nazca Ridge and its effects on the South American Plate. The process of flat slab subduction impacts the stress distribution within the overriding plate and increases plate coupling and seismic energy release. Using the finite element method (FEM), we analyse interseismic and coseismic deformation along a 1000 km transect parallel to the ridge. We examine stress distribution, uplift patterns, and the impact of megathrust activity on deformation. To better define the crust’s properties for the model, we developed a new thermal model of the Nazca Ridge subduction zone, reconstructing the thermal structure of the overriding plate. The results show concentrated stress at the upper part of the locked plate interface, extending into the Coastal and Western Cordilleras, with deeper stress zones correlating with seismicity. Uplift patterns align with long-term rates of 0.7–1 mm/yr. Cooling from flat-slab subduction strengthens the overriding plate, allowing far-field stress transmission and deformation. These findings provide insights into the tectonic processes driving stress accumulation, seismicity, and uplift along the Peruvian margin.

1. Introduction

The structural evolution of Andean-type orogens is strongly influenced by the geometry of the subducting slab. Several researchers [1,2,3,4,5,6] have widely investigated how anomalies in the thickness of the oceanic crust, such as ridges and plateaus, have significant implications for the geodynamics of subduction margins and the overriding plate [3,6]. In this study, we focus on the Nazca Ridge and the effects of its subduction on the overriding South American Plate.
The Nazca Ridge, shown in Figure 1, is a portion of oceanic crust with an estimated thickness of 18 ± 3 km [7] that is currently subducting beneath the South American Plate between 14° and 16° S [7]. The origin of this aseismic ridge is related to the spreading centre of the Pacific–Farallon/Nazca plates, and it was formed during the early Cretaceous [4,8,9].
The collision and the consequent subduction of the ridge with the South American Plate occurred around 11.2 Ma at a latitude of 11° S. Having an orientation of N42° E, the ridge tends to migrate southward, and it reached the current subduction zone around 4 to 5 Ma ago [10,11]. The ridge extends approximately 200 km in width and over 1100 km in length, with an additional estimated 900 km of the ridge already subducted beneath the South American Plate [9]. The presence of this thickened crust has probably induced the geological evolution of the Peruvian margin and a flat-slab subduction phenomenon beneath the South American continental margin [12]. According to tomographic data, the flat slab extends for nearly 300 km at depths between 80 and 100 km [13,14,15]. The presence of a flat slab in this region leads to several significant consequences. The stress distribution within the overriding plate is altered, resulting in a reduction of extensional stresses and localised shortening in the forearc region. In these conditions, the coupling between the overriding and subducting plates is strengthened and induces a greater release of seismic energy [16]. Despite this, there are notably few intermediate-depth earthquakes within the lower plate and along its upper boundary [17]. Additionally, the flat slab causes the asthenospheric wedge to withdraw from the trench, producing a westward pinching out of the asthenosphere beneath the overriding plate, and significantly impacting the thermal structure of the upper plate. Indeed, the withdrawal produces the eastward migration of arc magmatism and, sometimes, the complete deletion of volcanic activity [16]. Nevertheless, the impact of the Nazca Ridge’s flat subduction on the geological evolution of the Peruvian Cordillera, as well as the geodynamic behaviour of the Peruvian margin and the South American Plate, remains a subject of ongoing debate. In particular, the factors driving crustal shortening are still uncertain, with interplate coupling [16,18,19] and structural reactivation [20,21] both proposed as potential influences.
In this study, we investigated the stress distribution on the overriding plate and, in particular, during the interseismic and coseismic stages of the subduction zone and the megathrust. With the aim of simulating strain accumulation and interseismic stress, we discretized a 2D FEM model along a transect of 1000 km that runs almost parallel to the Nazca Ridge and the flat slab (Figure 1). The transect intersects the main features of the South Peruvian Cordillera until the undeformed area of the Amazonian Basin. The coseismic slip was constrained, taking into consideration the earthquakes with Mw ≥ 6 that occurred during the last 50 years in a buffer of 100 km around the study area [22]. These events are related to the megathrust and the distribution describes the subduction interface. To perform the FEM modelling of the megathrust, we used a modern and asserted methodology already employed in several studies [23] and in different geological contexts, (i.e., Chile [24], China [25]; Val D’Agri (Italy) [26]; Zagros (Iran) [27]. The procedure consists of an analysis of interseismic deformation in zones of plate convergence, which results from the elastic strain in a brittle layer above a ductile half space [28]. Understanding this deformation is crucial for comprehending the tectonic processes driving seismic activity and surface uplift. This model assumes that the brittle layer encompasses the entire model thickness, exhibiting stick-slip behaviour. Furthermore, we modelled coseismic deformation linked with certain seismic occurrences using the same numerical technique. To constrain our model, we utilized available geological information derived from both published works [11,12,13,14,15] and our fieldwork in the study area. Moreover, investigating the stress distribution and geodynamics of the area requires understanding the thermal structure of the overriding plate [29,30,31,32,33]. Using this approach, data about rheology are crucial for performing an accurate analysis of the co-seismic field, identifying areas with high strain rates, and determining where stress tends to accumulate. Therefore, in this study, we propose an upgrade to the thermal model proposed by Ciattoni et al. [34] for the same area.
The new thermal model considers new parameters of density and a more accurate geometry of the overriding plate. This combined approach enabled us to study recent uplift through topographic analysis and lithologic boundary conditions. By using FEM methodology, we aim to explore the interseismic and coseismic deformation processes and compare them with the observed uplift patterns. This comprehensive analysis not only enhances our understanding of crustal deformation but also provides insights into the underlying tectonic mechanisms.
Figure 1. Map of the study area (Esri World Imagery basemap). The red line (A–A’) is the trace of the studied transect; the white lines are the 100 km contour line indicating the depth of the slab from Hayes et al. [15]. The grey area shows the Nazca Ridge, while the light-grey area represents the subducted portion of the Nazca Ridge.
Figure 1. Map of the study area (Esri World Imagery basemap). The red line (A–A’) is the trace of the studied transect; the white lines are the 100 km contour line indicating the depth of the slab from Hayes et al. [15]. The grey area shows the Nazca Ridge, while the light-grey area represents the subducted portion of the Nazca Ridge.
Geosciences 14 00257 g001

2. Geological Setting

The Andean Cordillera is a continental mountain belt that ranges from 8° N to 52° S and is the result of the effect of the Nazca Plate subduction beneath the South American Plate [35]. The Peruvian Andes are a part of the Central Cordillera and, more specifically, belong to the southern portion of the Northern Central Andes and the northern branch of The Central Andes Orocline. The current configuration of this portion of the Cordillera is the expression of a compelling geological history [36,37]. During the Paleozoic, the first deformation step occurred during the Hercynian orogenic phase (Late Devonian and Early Mississippian), which strongly deformed the high-grade metamorphic basement (Precambrian) and the overlying marine sediments (Middle Ordovician to the Upper Devonian) [36]. The second crucial deformation stage, called the Andean orogenic cycle, began during the Cenozoic, and it is still active. In this tectonic phase, we can identify three compressional periods: (i) Peruvian (Late Cretaceous); (ii) Incaic (Eocene); and (iii) Quechua (Miocene). The Quechua compressional phase is subdivided into three subphases: (i) Quechua 1 (Lower Miocene), (ii) Quechua 2 (Tortonian), and (iii) Quechua 3 (Messinian). All of these compressional periods are combined with plutonic activity, volcanism, and periods of sedimentation, and part of the deformation occurred along the inherited structures of the Hercynian phase [36]. The most recent history of this part of the Cordillera is strongly linked to the Nazca Ridge subduction, and to its southwest migration. Since the beginning of its subduction 11 Ma at 11° S, the presence of the Nazca Ridge considerably affected the development of the Peruvian forearc system and the deformation of the innermost part of the overriding plate along the whole Peruvian margin until it came to its present position (14°–16° S) [9,38]. The most evident effect of the Nazca Ridge subduction is the evolution and uplift of the East Pisco forearc Basin [39,40] and the deformation of the Amazonian Basin with the formation of the Fitzcarrald Arch [9,19,41]. In this study, we focused our attention on the current subduction zone of the Nazca Ridge. The study area (Figure 2) goes from 78° W 16° S to 71° W 10° S and includes (i) the East Pisco Basin, which represents the subaerial portion of the Peruvian Forearc System; (ii) the Coastal Cordillera and the Marañón fold-and-thrust belt, which compose the Western Cordillera; (iii) the Eastern Cordillera; (iv) the Sub Andean Zone; and (v) the Fitzcarrald arch that is part of the Amazonian Basin.
Figure 2. Geological map of the study area based on the 1:100,000 scale geological maps from INGEMMET database. Modified after Ciattoni et al. [34]. The black line (A–A’) is the trace of the studied transect, and the white dots indicate the pseudowell location.
Figure 2. Geological map of the study area based on the 1:100,000 scale geological maps from INGEMMET database. Modified after Ciattoni et al. [34]. The black line (A–A’) is the trace of the studied transect, and the white dots indicate the pseudowell location.
Geosciences 14 00257 g002

3. Methods

In this study, FEM analysis is used to examine both interseismic and coseismic deformation along the section. This procedure consists of the geometric modelling of the fault and the crustal geological section, discretizing the model by 2D finite-element method (FEM) by MSC-Marc software (2012) [42] and recreating the main geological features of the tectono-stratigraphy units of the study area. Since the strength and the rheology of the lithosphere are functions of the temperature and the heat flow, we also analysed the thermal structure of the overriding plate. In particular, the model proposed by Ciattoni et al. [34] has been improved and integrated with new available data. Density data for the crustal layers were obtained from the models provided by Ciattoni et al. [43].

3.1. Thermal Model

The method used to develop the model is the same as that applied by Ciattoni et al. [34], where they used the one-dimensional calculation approach to 11 pseudo-wells placed along the section. In this research, we added a further pseudo-well to improve the model’s constraints in the Western Cordillera. The pseudowells reach the Moho interface, so their depths range from a minimum of 21 km (PW_2) to a maximum of 51 km (PW_5). Next, the heat flow obtained was interpolated through a sixth-order polynomial equation. Although the procedure and the zone remained unchanged from Ciattoni et al. [34], adjustments were made to the density and geometry of the layers based on new insights suggested by Ciattoni et al. [43]. These authors [43] produced a new integrated model of the crustal structure up to depth of 130 km along the transect across the Peruvian Forearc and Andean Cordillera. They integrated geological and geophysical data and performed forward modelling of the Bouguer anomaly to develop a 2D density model. They propose a more detailed geometry of the geological structures present in the area, including the slab and the Moho. Moreover, the density data are validated by a gravimetric model that takes into account the petrological properties of the materials and the P-T conditions. The adjustments provided in this study to the thermal model consist of changes to the basement thicknesses by 10% and values of the densities that vary by less than 5%. A thickness change can affect the heat flow amount, while a different density can affect the frictional heat calculation. To calculate: (i) heat flow at the surface, (ii) temperature gradients, and (iii) lines of equal temperature within the continental crust of South America, we utilized the method previously applied in various regions by Basilici et al. [44], Santini et al. [45], and Ciattoni et al. [34]. The thermal model described here accounts for four layers that represent the primary elements of the tectono-stratigraphy in the studied region: (i) sedimentary units from Mesozoic to Cenozoic, (ii) sedimentary units from Paleozoic, (iii) upper basement, and (iv) lower basement. The Paleozoic sequence is primarily made up of a mix of schists, sandstones, and limestones, while the Cenozoic–Mesozoic succession is mostly composed of micritic limestones and calcareous sandstones, with occasional quartz-rich sandstones, conglomerates, and volcanic materials like andesitic lavas and tuffs [46,47,48,49]. The upper basement consists of high-grade metamorphic rocks like granulites and amphibolites, while the lower basement is made up of mafic rocks. Figure 3 represents the entire plan.
Calculating surface heat flow Q s for the model in Figure 3, we considered the mantle heat flow Q m , radiogenic heat in the sedimentary cover ( Q S C ), and in the basement ( Q B ), as well as frictional heat Q F in the megathrust zone, resulting in the following equation:
Q S = Q + i = 1 2 Q S C i + Q B    
In the formula, Q represents the heat transfer at the bottom of the model caused by the heat transfer from the mantle and the frictional heat from the thrust if present. In the megathrust area, we took into account Q = ( Q 0 + Q F )   S 1 , where Q 0 represents the heat flow of the oceanic mantle and S = 1 + b   z f v s i n δ / κ 1 / 2 , with b = 1 , is the decreasing factor close to the megathrust fault [4]. For the continental zone, we took Q to be equal to Q m . Instead, Q S C i = H S C i h S C i represents the contributions from individual sedimentary layers with a constant radiogenic source, as mentioned in references [50,51,52,53], while Q B = H B D 1 e h U B + h L B D denotes the contribution of the basement with an exponential trend due to a radiogenic source, as discussed in references [54,55].
With a constant radiogenic source, the temperature in sedimentary cover layers is determined as follows:
T S C i z = T S C i + Q S i = 1 2 Q S C i 1 k S C i z z t i = 1 2 h S C i 1 H S C i 2 k S C i z z t i = 1 2 h S C i 1 2
The temperature at the top of layer i is denoted as T S C i .
Taking into account a radiogenic source exhibiting an exponential pattern, we calculate the temperature found within the basement as follows:
T B z = T B + Q   k B z z t i = 1 2 h S C i + H B D 2 k B 1 e z z t i = 1 2 h S C i D
T B represents the temperature at the top of the basement, while D represents the depth scale. The definitions of the terms in the equations can be found in Table 1.
Table 1. Definitions of terms considered in the thermal model [56,57,58,59,60,61].
Table 1. Definitions of terms considered in the thermal model [56,57,58,59,60,61].
Equation ParametersModel Data
z t (m) Ground Elevation (Negative) or Sea Depth (Positive) related to m.s.l.
h S C _ 1 (m) Sedimentary cover (Cenozoic–Mesozoic) thickness
h S C _ 2 (m) Sedimentary cover (Paleozoic) thickness
h U B (m) Upper basement thickness
h L B (m) Lower basement thickness
H S C _ 1 = 0.83   µ W m 3 RHP * of sedimentary cover (Cenozoic–Mesozoic) [56]
H S C _ 2 = 1.2   µ W m 3 RHP * of sedimentary cover (Paleozoic) [56]
H B = 1.6   µ W m 3 RHP * of basement thickness [56]
k S C _ 1 = 2.1   W m 1 K 1 Thermal conductivity of sedimentary cover (Cenozoic–Mesozoic) [57]
k S C _ 2 = 2.4   W m 1 K 1 Thermal conductivity of sedimentary cover (Paleozoic) [57]
k B = 2.7   W m 1 K 1 Thermal conductivity of the basement [57]
κ = 10 6   m 2 s 1 Thermal diffusivity [58]
µ = 0.7 Coefficient of static friction [58]
ρ S C _ 1 = 2.4 10 3   k g m 3 Density of sedimentary cover (Cenozoic–Mesozoic) [58]
ρ S C _ 2 = 2.4 10 3   k g m 3 Density of sedimentary cover (Paleozoic) [59]
ρ U B = 2.83 10 3   k g m 3 Density of the upper basement [59]
ρ L B = 2.98 10 3   k g m 3 Density of the lower basement [59]
x f   ( k m ) Horizontal coordinate for points on the megathrust fault
z f   ( k m ) Vertical coordinate for points on the megathrust fault
δ = a r c t g z f x f Dip angle of the megathrust fault
λ = 0.96 Pore fluid factor [60]
v = 7.7   c m y r 1 Relative plate velocity [60]
Q F = µ ρ g h 1 λ v Frictional heat flow density
Q o = 60   m W m 2 Oceanic heat flow density from the mantle [61]
Q m = 20   m W m 2 Continental heat flow density from the mantle [61]
Q S   ( m W m 2 ) Surface heat flow density
T S = 15   ° C Surface temperature
* Radiogenic heat productivity.

3.2. Finite-Element Modelling

This study examines the coseismic and interseismic deformation along the segment using finite element analysis. Using MSC-Marc software [42], the crustal geological section has been discretized in order to create a 2D FEM model that replicates the major tectono-stratigraphic features of the study region. As represented in Figure 4, in the model, the crust was considered to be composed of four main layers: (i) sedimentary cover, which includes the entire sedimentary succession without any distinction in age or composition; (ii) basement, composed of rocks of granitoid composition (including gneiss and granulites); (iii) oceanic crust; and (iv) lithospheric mantle. To define the rheological behaviour of the basement, we computed its strain rate using the empirical flow law in the well-known power law from Kirby [62]:
ϵ ˙ = A σ n e Q / R T
where σ is differential stress, A is a material property, Q is activation energy, R is the molar gas constant, and T is absolute temperature. The most appropriate values for the parameters included in the equation above were chosen by comparing both the observed Vp (recorded by seismic stations) and lithology type with published Vp values and available crustal rocks, as discussed by Ranalli and Murphy [63] and Ranalli [64]. When information about rock density is unavailable, it is often estimated from P-wave velocity (Vp) using an empirical relationship, such as Gardner’s equation [65]. For the basement layer, we selected a value of 1.8   ·   10 9   M P a n   s 1 for Parameter A, 123 kJ m o l 1 for the activation energy (Q), and 3.2 for the power-law exponent. By defining both temperature and differential stress as a function of depth in the calculation, we can describe the variations of strain rate with depth.
The proportion between the viscous and elastic behaviour of a particular element is governed by the local crustal geotherm derived from heat flow measurements. Equation (4) can be expressed in terms of a simple viscous relation:
ϵ ˙ = 1 2   η   σ
This is defined by an effective viscosity:
η = 1 2 A   σ 1 n   e Q / R T
which varies with depth as a function of temperature T [66,67].
The friction coefficient was set as μ = 0.7 according to Cermak et al. [58]. The model was divided into 65,670 quadrangular elementary cells and 66,740 nodes. Numerical modelling included two independent FEM simulations conducted using the same mesh: a first procedure was used to analyse interseismic stress and strain accumulations, while a second procedure was applied to investigate surface vertical motions associated with the megathrust coseismic stage (by modelling the characteristic earthquake). Modelling of the interseismic stage was divided into two different steps: the first step consisted of setting the boundary conditions of the model in order to observe the effects of gravity alone. In the model, the following boundary conditions were applied: (i) the surface was considered free to move in all directions; (ii) the SW and NE boundaries were considered free to move in the vertical direction and locked in the horizontal direction; and (iii) the base was treated as a Winkler’s foundation [68]. In the simulation of the hydrostatic pressure of the Earth’s mantle, the model was set in order to allow free horizontal movement, while the vertical motion was subject to the influence of an elastic spring with stiffness coefficient K equal to K = A/L · E, where A is the base length, L is the thickness of the model, and E is the average of Young’s modulus of the rocks included in the model.
Using this method, we were able to estimate the amount of equivalent Von Mises stress (σVM):
σ V M = 3 2 i j σ i j   σ i j ,   with     σ i j =   σ i j 1 3 k δ i j   σ k k
where σij is the stress component and δij is the Kronecker delta, and of equivalent strain (εeq):
ε e q = 2 3 i j ε i j   ε i j ,   with     ε i j =   ε i j 1 3 k δ i j   ε k k
where εij is the strain component [69]. This procedure also yields values of surface motion in the vertical direction, accumulated over a particular time interval during the interseismic stage. Modelling of the coseismic stage shares the first step with the interseismic simulation procedure. However, in this case, the second step consisted of the sudden movement of a part of the fault plane to simulate stick-slip behaviour. This behaviour was achieved by imposing a slip value consistent with the magnitude of the seismic event to be modelled (characteristic earthquake in our instance).
Studies on the historical and recent seismicity of the area show that the most significant events do not exceed Mw 8, in particular, an event of Mw 7.7 occurred on 12 November 1996, along our section at a depth of 33 km. In addition, other studies indicate a period of 100 years as seismic periodicity [70]. Based on these two indications, we created our FEM model (Figure 5) and calculated surface displacement, stress, and strain. For the coseismic phase, we considered a fault with a width of about 50 km, imposing a constant slip of 70 cm on it to simulate, as the FEM model was made, an event with a Mw 7.5. To define the fault width, we considered the portion of the slab where most events with magnitudes greater than 6 occur. The depth range for our target fault is between 30 km and 50 km. During the interseismic phase, we applied the plate convergence rate of 7.7 cm/yr [61], which, when projected along our section, provides a velocity of 6.5 cm/yr. This was applied to the southwestern edge of the FEM model for a duration of 100 years, locking the upper part of the slab down to the simulated fault’s base while allowing the slab beneath to move freely, representing aseismic movement.

4. Results

4.1. Thermal Model

The computed geotherms (Figure 6) are characterised by lower temperature values for Wells PW1 to PW4, while higher temperatures are obtained for Well PW5 at the same depth. Notably, there is a similar trend in the geothermal profiles for Wells PW7 to PW11. The variation of temperature as a function of depth along the entire geological profile leads to a complex trend of the isotherms. In Figure 7, moving from SW to NE, we observe: (i) a deepening of the isotherms near the cold plate beneath the megathrust; (ii) an arching of the curves in the central part of the profile, where there is a thickening of the sedimentary layers and the basement; and finally (iii) a trend subparallel to the surface for the remaining part.
The trend of the heat flux at the surface was defined by calculating the Qs of the 12 pseudo-wells through the analytical treatment described in Section 3.1, and the results were compared with the observed data available in the literature by projecting them orthogonally on the profile (Figure 7). In particular, the results show low values of Qs along the southwestern shallow part of the megathrust ranging from about 24 mW/m2 to 36 mW/m2, and then it continues to increase to about 47 mW/m2 at Pseudo-well 5. In the central part of the section, the heat flux decreases slightly and then remains constant between about 40 mW/m2 and 44 mW/m2. The low values of Qs (30 ± 5 mW/m2) are due to the presence of oceanic plate subduction, between the trench and the volcanic arc [34]. The higher values of the heat flux (45 − 50 mW/m2) depend mostly on the thickening of the crustal structure since the contribution of mantle flux is constant.

4.2. Finite-Element Model

Figure 8 shows the surface displacements for the coseismic and interseismic phase, as well as the cumulative coseismic–interseismic phase. The cumulative displacement has a maximum of about 0.3 m, over a period of 100 years (0.3 mm/yr), at the top of the fault similar to the coseismic displacement but decreases less rapidly than the coseismic vertical displacement due to contribution from the interseismic displacement, which ranges from 0.05 to 0.2 m. The interseismic results of the equivalent Von Mises stress (Figure 9) and the equivalent strain (Figure 10a) show an accumulation of stress and strain in the locked part of the slab, as well as in the deepest part. In order to compare stress and strain accumulation with the pattern of seismicity, the hypocenters of Mw ≥ 5 seismic events that occurred within 100 km over the last 40 years are shown in the geological section (Figure 10b). Further higher values of stress and strain are observed in the accretionary prism and in some areas of the interface between the basement and the sedimentary cover.

5. Discussion and Conclusions

The integration of recent data from Ciattoni et al. [43] has resulted in variations in both the heat flow and the geothermal trend, predominantly observed in the initial segment of the section. Notably, in contrast to the model proposed by Ciattoni et al. [34], heat flow decreases from PW2 to PW5, with the most significant variation reaching about 4 mW/m2 at PW3. Additionally, the isotherms are positioned at greater depths due to the reduction in temperature at those depths. In general, the thermal structure of the South American continental crust is characterised by unusually low temperatures at the Moho depth. This is consistent with the abrupt cooling of the upper plate, which is believed to accompany flat-slab subduction [16,72]. Therefore, the recent (i.e., Pliocene) establishment of flat-slab subduction along the study transect bears major implications for upper-plate rheology, as cooling produces a strengthening of the overriding plate. Upper-plate strengthening may allow far-field stress transmission, which in turn leads to foreland deformation. Therefore, it may be envisaged that “young” flat-slab subduction in the study area controls deformation, uplift, and recent (<4 Ma) [9,19,41] exhumation as far as the Fitzcarrald Arc foreland domain to the NE. Our study shows for the first time the marked effect of the intensely debated, recent (<5 Ma) Nazca ridge subduction in perturbing the thermal structure of the overriding South American crust. The thermal model conducted in this study was also instrumental in subsequent finite element modelling, as temperature controls effective viscosity in the FEM and, hence, have a significant impact on modelled rheology.
Finite-element modelling conducted in this study was focussed on permanent deformation and uplift patterns, neglecting any coseismic elastic rebound effect, which is not modelled in our FEM. The accumulation of elastic strain is well known to occur where plate motion is hindered by frictional coupling associated with asperities along the plate interface, representing a fundamental aspect of the well-established earthquake cycle model for megathrusts located at subduction plate boundaries [73]. Rupturing of the asperities and recovery of the accumulated strain occurs by elastic rebound of both subducting and overriding plates during megathrust earthquakes, as testified by the fact that most directly observable coseismic strain recovery occurs elastically during the seismic event. As a matter of fact, measured crustal displacements associated with megathrust earthquakes that originated at various subduction zones are consistent with such elastic rebound models, the region of maximum coseismic surface displacement being commonly located close to the trench (and offshore in most of the cases) [74].
The uplift pattern depicted by our FEM results (Figure 8), not taking into account coseismic elastic strain recovery, is consistent with the long-term uplift observed along the coastal area of the active margin of Peru. Macharé and Ortlieb [10] determined a maximum uplift rate of around 0.7 mm/yr, while Zeumann and Hampel [11] suggested that long-term surface maximum uplift in the studied sector of the Andes during the Quaternary is of ca. 1 mm/yr. These values agree with the values estimated by Sillard et al. [75], who, using the technique of cosmogenic 10Be exposure dating, estimated that cumulative uplift rates have progressively increased since approximately the middle Pleistocene from about 0.4 to 0.9 mm/yr. Moreover, these authors suggested that the increasing uplift rate for 800.000 years is related to the subduction of the Nazca Ridge.
However, such surface uplift is affected by erosion, which is not considered in our rock uplift model. As erosion rates obtained by Wipf et al. [39] for the study area range from 0.06 mm/yr to 4 mm/yr, our model satisfactorily compares with measured uplift rates. Due to the locked plate interface, modelled interseismic deformation produces an uplift pattern that extends inboard to the Coastal Cordillera and, to a lesser extent, even to the Western Cordillera. On the other hand, coseismic deformation is characterised by a distinct maximum of focussed ground motion in the offshore forearc basin domain. This latter uplift component is consistent with the Pliocene–Quaternary vertical motion and emergence of the eastern sector of the forearc basin system (i.e., the East Pisco basin), which is currently located at a maximum elevation of about 650 m above sea level. Within this framework, the uncertainties involved in earthquake recurrence intervals and erosion rates in the extrapolation of the coupled coseismic–interseismic deformation do not represent a major issue, since our aim is to obtain information on the location of zones of larger uplift and not on the absolute amounts of surface motion.
FEM results, besides showing predictable stress and strain accumulation at the upper tip of the unlocked (i.e., free to slip) plate interface, also show zones of stress concentration (Figure 9) and strain accumulation at deeper levels along the plate boundary (Figure 10a). This is consistent with transient seismic behaviour at much greater depths, which is testified by the seismicity (Figure 10b).
In conclusion, integrated thermal and mechanical modelling conducted in this study provide useful insights into the active tectonic setting and plate-scale geodynamics of the Peruvian margin along the axis of the subducting aseismic Nazca Ridge.

Author Contributions

Conceptualization, S.C., S.M., A.M. and S.S.; methodology, S.C., S.M., A.M. and S.S.; software, A.M. and S.C.; validation, S.C., A.M., S.M. and S.S.; formal analysis, S.C., A.M. and S.S.; investigation, S.C., S.M., A.M. and S.S.; resources, S.C., A.M. and S.S.; data curation, A.M., S.C., S.M. and S.S.; writing—original draft preparation, S.C.; writing—review and editing, S.M.; visualization, S.C.; supervision, S.M. and S.S.; project administration, S.S.; funding acquisition, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are provided in the paper.

Acknowledgments

The authors are grateful to the editor for his helpful support in organising the manuscript, and to three anonymous reviewers for useful and constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, C.; Clark, A.L. Tectonic Effects of the Subducting Nazca Ridge on the Southern Peru Continental Margin. Mar. Pet. Geol. 1994, 11, 575–586. [Google Scholar]
  2. Hagen, R.A.; Moberly, R. Tectonic Effects of a Subducting Aseismic Ridge: The Subduction of the Nazca Ridge at the Peru Trench. Mar. Geophys. Res. 1994, 16, 145–161. [Google Scholar] [CrossRef]
  3. Gutscher, M.A. Andean Subduction Styles and Their Effect on Thermal Structure and Interplate Coupling. J. South Am. Earth Sci. 2002, 15, 3–10. [Google Scholar] [CrossRef]
  4. Hampel, A.; Kukowski, N.; Bialas, J.; Huebscher, C.; Heinbockel, R. Ridge subduction at an erosive margin: The collision zone of the Nazca Ridge in southern Peru. J. Geophys. Res. Solid Earth 2004, 109, B02101. [Google Scholar] [CrossRef]
  5. Antonijevic, S.K.; Wagner, L.S.; Kumar, A.; Beck, S.L.; Long, M.D.; Zandt, G.; Tavera, H.; Condori, C. The role of ridges in the formation and longevity of flat slabs. Nature 2015, 524, 212–215. [Google Scholar] [CrossRef] [PubMed]
  6. Flórez-Rodríguez, A.G.; Schellart, W.P.; Strak, V. Impact of aseismic ridges on subduction systems: Insights from analog modeling. J. Geophys. Res. Solid Earth 2019, 124, 5951–5969. [Google Scholar] [CrossRef]
  7. Woods, M.T.; Okal, E.A. The structure of the Nazca ridge and Sala y Gomez seamount chain from the dispersion of Rayleigh waves. Geophys. J. Int. 1994, 117, 205–222. [Google Scholar] [CrossRef]
  8. Pilger, R.H., Jr. Plate Reconstructions, Aseismic Ridges, and Low-Angle Subduction beneath the Andes. Geol. Soc. Am. Bull. 1981, 92, 448–456. [Google Scholar] [CrossRef]
  9. Hampel, A. The Migration History of the Nazca Ridge along the Peruvian Active Margin: A Re-Evaluation. Earth Planet. Sci. Lett. 2002, 203, 665–679. [Google Scholar] [CrossRef]
  10. Macharé, J.; Ortlieb, L. Plio-Quaternary vertical motions and the subduction of the Nazca Ridge, central coast of Peru. Tectonophysics 1992, 205, 97–108. [Google Scholar] [CrossRef]
  11. Zeumann, S.; Hampel, A. Deformation of erosive and accretive forearcs during subduction of migrating and non-migrating aseismic ridges: Results from 3-D finite element models and application to the Central American, Peruvian, and Ryukyu margins. Tectonics 2015, 34, 1769–1791. [Google Scholar] [CrossRef]
  12. Bishop, B.T.; Beck, S.L.; Zandt, G.; Wagner, L.; Long, M.; Antonijevic, S.K.; Kumar, A.; Tavera, H. Causes and Consequences of Flat-Slab Subduction in Southern Peru. Geosphere 2017, 13, 1392–1407. [Google Scholar] [CrossRef]
  13. Scire, A.; Zandt, G.; Beck, S.; Long, M.; Wagner, L.; Minaya, E.; Tavera, H. Imaging the transition from flat to normal subduction: Variations in the structure of the Nazca slab and upper mantle under southern Peru and northwestern Bolivia. Geophys. J. Int. 2016, 204, 457–479. [Google Scholar] [CrossRef]
  14. Scire, A.; Zandt, G.; Beck, S.; Long, M.; Wagner, L. The deforming Nazca slab in the mantle transition zone and lower mantle: Constraints from teleseismic tomography on the deeply subducted slab between 6 S and 32 S. Geosphere 2017, 13, 665–680. [Google Scholar] [CrossRef]
  15. Hayes, G.P.; Moore, G.L.; Portner, D.E.; Hearne, M.; Flamme, H.; Furtney, M.; Smoczyk, G.M. Slab2, a comprehensive subduction zone geometry model. Science 2018, 362, 58–61. [Google Scholar] [CrossRef]
  16. Gutscher, M.A.; Spakman, W.; Bijwaard, H.; Engdahl, E.R. Geodynamics of Flat Subduction: Seismicity and Tomographic Constraints from the Andean Margin. Tectonics 2000, 19, 814–833. [Google Scholar] [CrossRef]
  17. Kumar, A.; Wagner, L.S.; Beck, S.L.; Long, M.D.; Zandt, G.; Young, B.; Tavera, H.; Minaya, E. Seismicity and State of Stress in the Central and Southern Peruvian Flat Slab. Earth Planet. Sci. Lett. 2016, 441, 71–80. [Google Scholar] [CrossRef]
  18. Ramos, V.A. Seismic Ridge Subduction and Topography: Foreland Deformation in the Patagonian Andes. Tectonophysics 2005, 399, 73–86. [Google Scholar] [CrossRef]
  19. Bishop, B.T.; Beck, S.L.; Zandt, G.; Wagner, L.S.; Long, M.D.; Tavera, H. Foreland Uplift during Flat Subduction: Insights from the Peruvian Andes and Fitzcarrald Arch. Tectonophysics 2018, 731, 73–84. [Google Scholar] [CrossRef]
  20. McQuarrie, N. The Kinematic History of the Central Andean Fold-Thrust Belt, Bolivia: Implications for Building a High Plateau. Geol. Soc. Am. Bull. 2002, 114, 950–963. [Google Scholar] [CrossRef]
  21. Gautheron, C.; Espurt, N.; Barbarand, J.; Roddaz, M.; Baby, P.; Brusset, S.; Douville, E. Direct Dating of Thick-and Thin-skin Thrusts in the Peruvian Subandean Zone through Apatite (U–Th)/He and Fission Track Thermochronometry. Basin Res. 2013, 25, 419–435. [Google Scholar] [CrossRef]
  22. United States Geological Survey (USGS) Catalog. Available online: https://earthquake.usgs.gov/earthquakes/search/ (accessed on 23 May 2024).
  23. Megna, A.; Barba, S.; Santini, S.; Dragoni, M. Effects of Geological Complexities on Coseismic Displacement: Hints from 2D Numerical Modelling. Terra Nova 2008, 20, 173–179. [Google Scholar] [CrossRef]
  24. Vigny, C.; Rudloff, A.; Ruegg, J.C.; Madariaga, R.; Campos, J.; Alvarez, M. Upper Plate Deformation Measured by GPS in the Coquimbo Gap, Chile. Phys. Earth Planet. Inter. 2009, 175, 86–95. [Google Scholar] [CrossRef]
  25. Zhu, S.; Zhang, P. FEM Simulation of Interseismic and Coseismic Deformation Associated with the 2008 Wenchuan Earthquake. Tectonophysics 2013, 584, 64–80. [Google Scholar] [CrossRef]
  26. Candela, S.; Mazzoli, S.; Megna, A.; Santini, S. Finite Element Modelling of Stress Field Perturbations and Interseismic Crustal Deformation in the Val d’Agri Region, Southern Apennines, Italy. Tectonophysics 2015, 657, 245–259. [Google Scholar] [CrossRef]
  27. Basilici, M.; Ascione, A.; Megna, A.; Santini, S.; Tavani, S.; Valente, E.; Mazzoli, S. Active Deformation and Relief Evolution in the Western Lurestan Region of the Zagros Mountain Belt: New Insights from Tectonic Geomorphology Analysis and Finite Element Modeling. Tectonics 2020, 39, e2020TC006402. [Google Scholar] [CrossRef]
  28. Savage, J.C. A Dislocation Model of Strain Accumulation and Release at a Subduction Zone. J. Geophys. Res. Solid Earth 1983, 88, 4984–4996. [Google Scholar] [CrossRef]
  29. Jaupart, C.; Mareschal, J.C. The thermal structure and thickness of continental roots. Lithos 1999, 48, 93–114. [Google Scholar] [CrossRef]
  30. Cheng, L.Z.; Mareschal, J.C.; Jaupart, C.; Rolandone, F.; Gariépy, C.; Radigon, M. Simultaneous inversion of gravity and heat flow data: Constraints on thermal regime, rheology and evolution of the Canadian Shield crust. J. Geodyn. 2002, 34, 11–30. [Google Scholar] [CrossRef]
  31. McKenzie, D.; Jackson, J.; Priestley, K. Thermal structure of oceanic and continental lithosphere. Earth Planet. Sci. Lett. 2005, 233, 337–349. [Google Scholar] [CrossRef]
  32. Chen, W. Rheology of the continental lithosphere: Progress and new perspectives. Gondwana Res. 2012, 21, 4–18. [Google Scholar] [CrossRef]
  33. Moeck, I.S. Catalog of geothermal play types based on geologic controls. Renew. Sustain. Energy Rev. 2014, 37, 867–882. [Google Scholar] [CrossRef]
  34. Ciattoni, S.; Mazzoli, S.; Megna, A.; Basilici, M.; Santini, S. Two-Dimensional Geothermal Model of the Peruvian Andes above the Nazca Ridge Subduction. Energies 2023, 16, 7697. [Google Scholar] [CrossRef]
  35. Ramos, V.A.; Aleman, A. Tectonic evolution of the Andes. In Tectonic Evolution of South America, Proceedings of the 31st International Geological Congress, Rio de Janeiro, Brazil, 6–17 August 2000; Cordani, U.G., Milani, E.J., Thomaz Filho, A., Campos, D.A., Eds.; Geological Society: London, UK, 2000; pp. 635–685. [Google Scholar]
  36. Mégard, F.; Schaer, J.; Rodgers, J. Structure and evolution of the Peruvian Andes. In The Anatomy of Mountain Ranges; Schaer, J.-P., Rodgers, J., Eds.; Princeton University Press: Princeton, NJ, USA, 1987; pp. 179–210. [Google Scholar]
  37. Sempere, T.; Folguera, A.; Gerbault, M. New Insights into Andean Evolution: An Introduction to Contributions from the 6th ISAG Symposium (Barcelona, 2005). Tectonophysics 2008, 459, 1–13. [Google Scholar] [CrossRef]
  38. Krabbenhöft, A.; Bialas, J.; Kopp, H.; Kukowski, N.; Hübscher, C. Crustal Structure of the Peruvian Continental Margin from Wide-Angle Seismic Studies. Geophys. J. Int. 2004, 159, 749–764. [Google Scholar] [CrossRef]
  39. Wipf, M.; Zeilinger, G.; Seward, D.; Schlunegger, F. Focused Subaerial Erosion during Ridge Subduction: Impact on the Geomorphology in South-central Peru. Terra Nova 2008, 20, 1–10. [Google Scholar] [CrossRef]
  40. Di Celma, C.; Pierantoni, P.P.; Volatili, T.; Molli, G.; Mazzoli, S.; Sarti, G.; Ciattoni, S.; Bosio, G.; Malinverno, E.; Collareta, A.; et al. Towards Deciphering the Cenozoic Evolution of the East Pisco Basin (Southern Peru). J. Maps 2022, 18, 397–412. [Google Scholar] [CrossRef]
  41. Espurt, N.; Baby, P.; Brusset, S.; Roddaz, M.; Hermoza, W.; Barbarand, J. The Nazca Ridge and Uplift of the Fitzcarrald Arch: Implications for Regional Geology in Northern South America. In Amazonia: Landscape and Species Evolution: A Look into the Past; Wiley-Blackwell: Hoboken, NJ, USA, 2010; pp. 89–100. [Google Scholar]
  42. MSC Software Corporation, Home Page (2012). MSC-Marc Software. Available online: http://www.mscsoftware.com/ (accessed on 25 June 2024).
  43. Ciattoni, S.; Cella, F.; Mazzoli, S.; Zambrano, M.; Megna, A.; Santini, S.; Butler, R.; Pierantoni, P.P.; Di Celma, C. Ciattoni_2024.Zip. Figshare. 2024. Available online: https://figshare.com/articles/dataset/Ciattoni_2024_zip/26029138/1?file=47029078 (accessed on 18 June 2024).
  44. Basilici, M.; Mazzoli, S.; Megna, A.; Santini, S.; Tavani, S. 3-D Geothermal Model of the Lurestan Sector of the Zagros Thrust Belt, Iran. Energies 2020, 13, 2140. [Google Scholar] [CrossRef]
  45. Santini, S.; Basilici, M.; Invernizzi, C.; Jablonska, D.; Mazzoli, S.; Megna, A.; Pierantoni, P.P. Controls of Radiogenic Heat and Moho Geometry on the Thermal Setting of the Marche Region (Central Italy): An Analytical 3D Geothermal Model. Energies 2021, 14, 6511. [Google Scholar] [CrossRef]
  46. DeVries, T.J.; Urbina, M.; Jud, N.A. The Eocene-Oligocene Otuma Depositional Sequence (East Pisco Basin, Peru): Paleogeographic and Paleoceanographic Implications of New Data. Boletín Soc. Geológica Perú 2017, 112, 14–38. [Google Scholar]
  47. Di Celma, C.; Malinverno, E.; Bosio, G.; Collareta, A.; Gariboldi, K.; Gioncada, A.; Molli, G.; Basso, D.; Malca, V.; Varas-Malca, R.M.; et al. Sequence Stratigraphy and Paleontology of the Upper Miocene Pisco Formation along the Western Side of the Lower Ica Valley (Ica Desert, Peru). Riv. Ital. Paleontol. Stratigr. 2017, 123, 255–273. [Google Scholar]
  48. Di Celma, C.; Malinverno, E.; Bosio, G.; Gariboldi, K.; Collareta, A.; Gioncada, A.; Landini, W.; Pierantoni, P.P.; Bianucci, G. Intraformational Unconformities as a Record of Late Miocene Eustatic Falls of Sea Level in the Pisco Formation (Southern Peru). J. Maps 2018, 14, 607–619. [Google Scholar] [CrossRef]
  49. Pfiffner, O.A.; Gonzalez, L. Mesozoic–Cenozoic Evolution of the Western Margin of SouthAmerica: Case Study of the Peruvian Andes. Geosciences 2013, 3, 262–310. [Google Scholar] [CrossRef]
  50. Molnar, P.; England, P. Temperatures, heat flux, and frictional stress near major thrust faults. J. Geophys. Res. Solid Earth 1990, 95, 4833–4856. [Google Scholar] [CrossRef]
  51. Santini, S.; Basilici, M.; Invernizzi, C.; Mazzoli, S.; Megna, A.; Pierantoni, P.P.; Spina, V.; Teloni, S. Thermal Structure of the Northern Outer Albanides and Adjacent Adriatic Crustal Sector, and Implications for Geothermal Energy Systems. Energies 2020, 13, 6028. [Google Scholar] [CrossRef]
  52. Valdenegro, P.; Muñoz, M.; Yáñez, G.; Parada, M.A.; Morata, D. A model for thermal gradient and heat flow in central Chile: The role of thermal properties. J. South Am. Earth Sci. 2019, 91, 88–101. [Google Scholar] [CrossRef]
  53. Luo, T.; Leng, W. Thermal structure of continental subduction zone: High temperature caused by the removal of the preceding oceanic slab. Earth Planet. Phys. 2021, 5, 290–295. [Google Scholar] [CrossRef]
  54. Dragoni, M.; Doglioni, C.; Mongelli, F.; Zito, G. Evaluation of Stresses in Two Geodynamically Different Areas: Stable Foreland and Extensional Backarc. Pure Appl. Geophys. 1996, 146, 319–341. [Google Scholar] [CrossRef]
  55. Molnar, P.; Chen, W.P.; Padovani, E. Calculated temperatures in overthrust terrains and possible combinations of heat sources responsible for the tertiary granites in the greater Himalaya. J. Geophys. Res. 1983, 88, 6415–6429. [Google Scholar]
  56. Megna, A.; Candela, S.; Mazzoli, S.; Santini, S. An analytical model for the geotherm in the Basilicata oil fields area (southern Italy). Ital. J. Geosci. 2014, 133, 204–213. [Google Scholar] [CrossRef]
  57. Vilà, M.; Fernández, M.; Jiménez-Munt, I. Radiogenic heat production variability of some common lithological groups and its significance to lithospheric thermal modeling. Tectonophysics 2010, 490, 152–164. [Google Scholar] [CrossRef]
  58. Cermak, V.; Rybach, L. Thermal Conductivity and Specific Heat of Mineral and Rocks. In Landolt–Bornstein: Numerical Data and Functional Relationships in Science and Technology, Physical Properties of Rocks; Angenheister, G., Ed.; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1982; Volume 1, pp. 305–343. [Google Scholar]
  59. Dragoni, M.; Santini, S. Contribution of the 2010 Maule Megathrust Earthquake to the Heat Flow at the Peru-Chile Trench. Energies 2022, 15, 2253. [Google Scholar] [CrossRef]
  60. Rodriguez Piceda, C.; Scheck Wenderoth, M.; Gomez Dacal, M.L.; Bott, J.; Prezzi, C.B.; Strecker, M.R. Lithospheric density structure of the southern Central Andes constrained by 3D data-integrative gravity modeling. Int. J. Earth Sci. 2021, 110, 2333–2359. [Google Scholar] [CrossRef]
  61. Seno, T. Determination of the pore fluid pressure ratio at seismogenic megathrusts in subduction zones: Implications for strength of asperities and Andean-type Mountain building. J. Geophys. Res. 2009, 114, B05405. [Google Scholar] [CrossRef]
  62. Kirby, S.H. Rheology of the Lithosphere. Rev. Geophys. 1983, 21, 1458–1487. [Google Scholar] [CrossRef]
  63. Ranalli, G.; Murphy, D.C. Rheological Stratification of the Lithosphere. Tectonophysics 1987, 132, 281–295. [Google Scholar] [CrossRef]
  64. Ranalli, G. Rheology of the Earth; Chapman & Hall: London, UK, 1995. [Google Scholar]
  65. Miller, S.L.M.; Stewart, R.R. The Relationship between Elastic-Wave Velocities and Density in Sedimentary Rocks: A Proposal. CREWES Res. Rep. 1991, 3, 14. [Google Scholar]
  66. Kirby, S.H.; Kronenberg, A.K. Rheology of the Lithosphere: Selected Topics. Rev. Geophys. 1987, 25, 1219–1244. [Google Scholar] [CrossRef]
  67. Parsons, T. Post-1906 Stress Recovery of the San Andreas Fault System Calculated from Three-dimensional Finite Element Analysis. J. Geophys. Res. Solid Earth 2002, 107, ESE 3-1–ESE 3-13. [Google Scholar] [CrossRef]
  68. Williams, C.A.; Richardson, R.M. A Rheologically Layered Three-dimensional Model of the San Andreas Fault in Central and Southern California. J. Geophys. Res. Solid Earth 1991, 96, 16597–16623. [Google Scholar] [CrossRef]
  69. Zhuang, Z.; Liu, Z.; Cui, Y. Equivalent Stress and Equivalent Strain. In Dislocation Mechanism-Based Crystal Plasticity; Academic Press: Cambridge, MA, USA, 2019; pp. 67–69. [Google Scholar]
  70. Dorbath, L.; Cisternas, A.; Dorbath, C. Assessment of the Size of Large and Great Historical Earthquakes in Peru. Bull. Seismol. Soc. Am. 1990, 80, 551–576. [Google Scholar]
  71. International Seismological Centre (20XX), On-Line Bulletin. Available online: http://www.isc.ac.uk/iscbulletin/search/ (accessed on 18 July 2024).
  72. Liu, X.; Wagner, L.S.; Currie, C.A.; Caddick, M.J. Implications of flat-slab subduction on hydration, slab seismicity, and arc volcanism in the Pampean region of Chile and Argentina. Geochem. Geophys. Geosystems 2024, 25, e2023GC011317. [Google Scholar] [CrossRef]
  73. Hayes, G.P. The finite, kinematic rupture properties of great-sized earthquakes since 1990. Earth Planet. Sci. Lett. 2017, 468, 94–100. [Google Scholar] [CrossRef]
  74. Govers, R.; Furlong, K.P.; Van de Wiel, L.; Herman, M.W.; Broerse, T. The Geodetic Signature of the Earthquake Cycle at Subduction Zones: Model Constraints on the Deep Processes. Rev. Geophys. 2018, 56, 6–49. [Google Scholar] [CrossRef]
  75. Saillard, M.; Hall, S.R.; Audin, L.; Farber, D.L.; Regard, V.; Hérail, G. Andean Coastal Uplift and Active Tectonics in Southern Peru: 10Be Surface Exposure Dating of Differentially Uplifted Marine Terrace Sequences (San Juan de Marcona, ~15.4 S). Geomorphology 2011, 128, 178–190. [Google Scholar] [CrossRef]
Figure 3. Sketch of the layers used in the analytical procedure (the parameters are listed in Table 1).
Figure 3. Sketch of the layers used in the analytical procedure (the parameters are listed in Table 1).
Geosciences 14 00257 g003
Figure 4. Sketch of the FEM model, showing the layers used in modelling.
Figure 4. Sketch of the FEM model, showing the layers used in modelling.
Geosciences 14 00257 g004
Figure 5. Detail of the FEM model, showing the constructed mesh pattern. The cells are of smaller size in areas that are more critical for modelling deformation and stress accumulation.
Figure 5. Detail of the FEM model, showing the constructed mesh pattern. The cells are of smaller size in areas that are more critical for modelling deformation and stress accumulation.
Geosciences 14 00257 g005
Figure 6. Computed geotherms for the 12 pseudo-wells.
Figure 6. Computed geotherms for the 12 pseudo-wells.
Geosciences 14 00257 g006
Figure 7. The values of heat flow Qs (blue dots) are interpolated using a sixth-degree polynomial (solid blue line) to discern the trend across the crustal cross-section. A graphical comparison is made with the observed data (black dots) and their associated relative errors (black bars), taking into account the maximum and minimum variations in radiogenic heat production, heat flow, and slab dip during the heat-flow calculations (dashed blue lines). The isotherms are displayed in the section at intervals of 50 °C (dashed and solid black lines). For comparative purposes, the lithospheric section, provided by Ciattoni et al. [43], is included in the background.
Figure 7. The values of heat flow Qs (blue dots) are interpolated using a sixth-degree polynomial (solid blue line) to discern the trend across the crustal cross-section. A graphical comparison is made with the observed data (black dots) and their associated relative errors (black bars), taking into account the maximum and minimum variations in radiogenic heat production, heat flow, and slab dip during the heat-flow calculations (dashed blue lines). The isotherms are displayed in the section at intervals of 50 °C (dashed and solid black lines). For comparative purposes, the lithospheric section, provided by Ciattoni et al. [43], is included in the background.
Geosciences 14 00257 g007
Figure 8. Modelled results for coseismic, interseismic, and cumulative vertical displacement.
Figure 8. Modelled results for coseismic, interseismic, and cumulative vertical displacement.
Geosciences 14 00257 g008
Figure 9. Von Mises stress results. White fault segment is unlocked (i.e., free to slip).
Figure 9. Von Mises stress results. White fault segment is unlocked (i.e., free to slip).
Geosciences 14 00257 g009
Figure 10. (a) Equivalent total strain results. White fault segment is unlocked (i.e., free to slip). (b) Earthquake hypocenters of events with Mw ≥ 5 that occurred in the last 40 years, projected from a 100 km wide region along the study section. Data from the USGS catalogue [22] and ISC catalogue [71].
Figure 10. (a) Equivalent total strain results. White fault segment is unlocked (i.e., free to slip). (b) Earthquake hypocenters of events with Mw ≥ 5 that occurred in the last 40 years, projected from a 100 km wide region along the study section. Data from the USGS catalogue [22] and ISC catalogue [71].
Geosciences 14 00257 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ciattoni, S.; Mazzoli, S.; Megna, A.; Santini, S. Seismotectonic Setting of the Andes along the Nazca Ridge Subduction Transect: New Insights from Thermal and Finite Element Modelling. Geosciences 2024, 14, 257. https://doi.org/10.3390/geosciences14100257

AMA Style

Ciattoni S, Mazzoli S, Megna A, Santini S. Seismotectonic Setting of the Andes along the Nazca Ridge Subduction Transect: New Insights from Thermal and Finite Element Modelling. Geosciences. 2024; 14(10):257. https://doi.org/10.3390/geosciences14100257

Chicago/Turabian Style

Ciattoni, Sara, Stefano Mazzoli, Antonella Megna, and Stefano Santini. 2024. "Seismotectonic Setting of the Andes along the Nazca Ridge Subduction Transect: New Insights from Thermal and Finite Element Modelling" Geosciences 14, no. 10: 257. https://doi.org/10.3390/geosciences14100257

APA Style

Ciattoni, S., Mazzoli, S., Megna, A., & Santini, S. (2024). Seismotectonic Setting of the Andes along the Nazca Ridge Subduction Transect: New Insights from Thermal and Finite Element Modelling. Geosciences, 14(10), 257. https://doi.org/10.3390/geosciences14100257

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop