Application of anisotropy of magnetic susceptibility (AMS) fabrics to determine the kinematics of active tectonics: examples from the Betic Cordillera, Spain, and the Northern Apennines, Italy

The anisotropy of magnetic susceptibility (AMS) technique provides an effective way to measure fabrics and, in the process, interpret the kinematics of actively deforming orogens. We collected rock fabric data of alluvial fan sediments surrounding the Sierra Nevada massif, Spain, and a broader range of Cenozoic sediments and rocks across the Northern Apennine foreland, Italy, to explore the deformation fabrics that contribute to the ongoing discussions of orogenic kinematics. The Sierra Nevada is a regional massif in the hinterland of the Betic Cordillera. We recovered nearly identical kinematics regardless of specimen magnetic mineralogy, structural position, crustal depth, or time. The principal elongation axes are NE–SW in agreement with mineral lineations, regional GPS geodesy, and seismicity results. The axes trends are consistent with the convergence history of the Africa–Eurasia plate boundary. In Italy, we measured AMS fabrics of specimens collected along a NE–SW corridor spanning the transition from crustal shortening to extension in the Northern Apennines. Samples have AMS fabrics compatible only with shortening in the Apennine wedge and have locked in penetrative contractional fabrics, even for those samples that were translated into the actively extending domain. In both regions, we found that specimens have a low degree of anisotropy and oblate susceptibility ellipsoids that are consistent with tectonic deformation superposed on compaction fabrics. Collectively, these studies demonstrate the novel ways that AMS can be combined with structural, seismic, and GPS geodetic data to resolve orogenic kinematics in space and time.


Introduction
A number of circum-Mediterranean orogens are associated with rapid slab rollback, resulting in paired compressional and extensional domains in the orogenic wedge of the retreating upper plate (Elter, 1975;Carminati and Doglioni, 2012). Examples include the Calabria Arc-Tyrrhenian Sea (Beccaluva et al., 1985;Milia et al., 2009), the Hellenic Arc-Aegean Sea (Le Pichon and Angelier, 1979;Papazachos et al., 2000), and the Gibraltar Arc-Alboran Sea (Lonergan and White, 1997;Platt et al., 2006;Fernández-Ibáñez and Soto, 2008). Along these tectonic boundaries, the temporal and spatial relationship between thrust belt contraction, wedge-top basin evolution, hinterland extension, and orogenic uplift are the subjects of continuing controversy.
Finite and incremental strain data provide deformation history and fabric distribution information for kinematic studies of folds, faults, and orogens (e.g., Ramsay and Huber, 1984;Fagereng and Biggs, 2018). However, in orogenic forelands where deformation occurs at shallow depths and low temperatures, ductile penetrative deformation features may be absent and brittle structures may be sparse. Anisotropy of magnetic susceptibility (AMS) results offers an alternative proxy for grain-preferred orientation, and hence rock strain, to determine the tectonic fabric in these orogens where other deformation markers are not available (Borradaile andJackson, 2004, 2010;Borradaile and Henry, 1997;Averbuch et al., 1992;Parés, 2004). In general, comparative studies from siliciclastic rocks show good agreement between both the relative magnitude and orientation of penetrative rock strain determined by traditional geometric methods and AMS principal axes; however, in specimens dominated by diamagnetic mineral abundance, the AMS axes' orientation, and not necessarily their magnitude, correlates to the rock strain (e.g., Latta and Anastasio, 2007;Burmeister et al., 2009). In this paper, we show how AMS can extend the temporal reach of GPS geodesy back in time in orogenic studies of the Betic Cordillera, Spain, and in the Northern Apennines, Italy (e.g., Sagnotti et al., 1998;Mattei et al., 2004;Fig. 1).

Kinematic studies for active tectonic research
Sedimentary rocks acquire a primary depositional fabric, which is bedding parallel. It is measurable with the AMS technique and is further enhanced and modified during burial, compaction, and water loss (e.g., Tarling and Hrouda, 1993;Schwehr et al., 2006). Even unconsolidated rocks record a magnetic fabric that can potentially provide a kinematic record (Mattei et al., 1997;Porreca and Mattei, 2012). The sensitivity of AMS allows its use as a paleogeodetic tool in tectonic studies. Kinematics allow for an assessment of rheology and strain history that are necessary prerequisites for understanding geodynamics, incrementally balancing cross sections, or in paleogeographic reconstructions. We sampled both consolidated sedimentary rocks and unconsolidated sediments in the Betic Cordillera, Spain, and northeastern Apennine ranges, Italy, for AMS analysis. The Betics field sampling was designed to test AMS recovery from unburied and unconsolidated sediments around the Sierra Nevada massif. Here, oriented samples were collected from sites around the Sierra Nevada massif in Plio-Pleistocene terrestrial, siliciclastic deposits (Table A1 in Appendix). The Apennines field sampling was designed to measure the rotation of strain across the foreland as sampling site passes from the actively shortening part of the orogenic wedge near the trench to the actively extending regime further to the southwest. Here, oriented samples were collected from sites along a NE-SW-oriented corridor inclusive of Cenozoic marine and fluvial siliciclastics, marls, and carbonate rocks, and unconsolidated Pleistocene fluvial sediments (Table A1).

The AMS method
The AMS ellipsoid is defined by the principal axes (k 1maximum, k 2 -intermediate, k 3 -minimum) of a specimen. It can be represented by a second-rank tensor that characterizes a material's magnetization response to an applied magnetic field (e.g., Borradaile and Tarling, 1981;Tarling and Hrouda, 1993). The orientation and relative length of the principal anisotropy axes of a specimen are controlled by the preferred alignment of the anisotropy axes of the individual magnetic particles in the specimen and the degree of the individual particle's anisotropy. The anisotropy of individual magnetic grains is controlled by their crystallography and grain shape (Tarling and Hrouda, 1993). For magnetite grains, the anisotropy is controlled by grain shape, whereas for hematite and phyllosilicate particles the anisotropy can be controlled by a grain's crystallography, which, in turn, controls their shape. This does not preclude the possibility of the mutual orientation of different particles (e.g., Housen et al., 1993;Weil and Yonkee, 2009) or the preferential alignment of iron-bearing inclusions within the particles (e.g., Biedermann, 2018;Borradaile and Werner, 1994;Borradaile and Lagroix, 2000;Martín-Hernandez and Hirt, 2003;Parés and van der Pluijm, 2002) controlling the AMS fabric.
Natural processes such as current deposition, lithification, and tectonic deformation all contribute to a specimen's AMS. In deformed rocks, it was shown that the principal susceptibility axis (k 1 ) orientation is typically parallel to the strain long axis and orthogonal to the tectonic shortening direction, whereas the shortest axis (k 3 ) is orthogonal to bedding in orientation (e.g., Kligfield et al., 1981;Hrouda, 1982), regardless of whether the individual particle anisotropy is controlled by crystallography or shape.
The sedimentary rocks and deposits in this study contain enough phyllosilicate minerals to be excellent specimens for AMS studies because of the presence of oblate mineral grains which adjust readily to deposition, lithification, and any subsequent deformation. As grains reorient in response to depositional or tectonic processes, the magnetic fabric will continuously adjust (Parés and van der Pluijm, 2002). Deposition from currents in alluvial fans or rivers like the examples discussed here can cause preferred grain alignment. Because the intermediate and maximum AMS axes of platy grains, such as phyllosilicates, are nearly equal in magnitude, they will be randomly oriented within the bedding plane, with the minimum axes orthogonal to bedding. In mudstones and finegrained sandstones, where both paramagnetism and ferromagnetism contributions were quantified, paramagnetic mineral grains typically dominate the AMS signal (e.g., Martín-Hernández and Hirt, 2001) because of the shape anisotropy of clay minerals, although very fine magnetic particles attached to the clay fabric might also contribute (Kodama and Sun, 1992). Clustering of k 1 axes (magnetic lineation) under progressive deformation has been observed in a number of studies and is occasionally accompanied by a girdle containing k 3 and k 2 axes. This particular axes distribution is thought to be the first piece of evidence for layer parallel shortening in sedimentary rocks (Kissel et al., 1986;Sagnotti and Speranza, 1993;Parés and Dinarès, 1993;Sagnotti D. J. Anastasio et al.: Application of anisotropy of magnetic susceptibility (AMS) fabrics  Parés et al., 1999) in sequences of low to moderately deformed mudstones. Borradaile (1988) had already described that the intersection of two competing fabrics (e.g., bedding and layer parallel shortening) will control the orientation of the k 1 direction, producing an intersection lineation or magnetic lineation. The concept of magnetic lineation as a result of an assemblage of clay minerals (platelets) was originally brought up by Henry (1997), who referred to it as the "zone axis of two planar objects", a concept that was reviewed by Parés and van der Pluijm (2002).
4 Example I: Sierra Nevada massif, Spain

Geologic setting of the Sierra Nevada massif
The Sierra Nevada massif is part of the Betic Cordillera-Rif-Tell orogens that extend along the European-African plate boundary from the southern Iberian Peninsula to northern Africa. These orogens were formed by slab rollback and western migration of the Gibraltar Arc throughout the Neogene (Rosenbaum et al., 2002). Coincident with the translation of the arc, the upper plate experienced shortening, the growth of doubly vergent thrust belts, crustal thickening, and rock uplift (Duggen et al., 2003;Soto et al., 2008;Platt et al., 2013). In the Betics, contraction across the plate boundary was initially directed northward (Sanz De Galdeano, 1990;Lonergan, 1993;Platt et al., 2013). As the contraction continued into the foreland during the late Miocene, it slowed and progressively rotated to the northwest into its present orientation (Mazzoli and Helman, 1994;Rosenbaum et al., 2002). Active tectonics in the Betic Cordillera today is dominated by distributed NW-SE convergence of 4-6 mm yr −1 (Fernández-Ibanez et al., 2007;Koulali et al., 2011;Gutscher et al., 2012;Mancilla et al., 2013) and is accommodated in part on NW-SE-trending normal faults (Martínez-Martínez et al., 2006;Stich et al., 2006;Fernández-Ibáñez and Soto, 2008;Giaconia et al., 2014Giaconia et al., , 2015Fig. 2).
The Sierra Nevada massif is a doubly plunging, actively uplifting (Azañón et al., 2015) elongate dome, characterized by medium to low-grade metamorphic rocks stacked in north-verging thrust sheets (Martínez-Martínez et al., 2002). Previous interpretations are that the Sierra Nevada dome was uplifted following top-to-west extension and isostatic rebound after thrust belt formation (Martínez-Martínez et al., 2006). Alternatively, as many culminations exist in orogenic hinterlands, the massif could have been uplifted during contractional or transpressive strain (e.g., Bernini, 1990;Mitra et al., 1997).
To resolve whether the uplift of the Sierra Nevada dome was the result of extensional exhumation or a compressional orogenic culmination, we collected rock fabric (AMS) data in Plio-Pleistocene deposits around the massif to explore the presence of penetrative tectonic fabrics that can contribute additional constraints to the kinematics of dome emplacement. We focused sampling on unburied alluvial fan deposits in Neogene basins that surround the core of the structure (Fig. 3).

Methods for Example I
We collected samples from six sites distributed around Sierra Nevada, from all structural positions, around the massif in unburied Plio-Pleistocene fan deposits that range from poorly cemented to unconsolidated (Sanz de Galdeano and Vera, 1992; Table A1; Fig. 3). The ages of the deposits sampled were determined from published geologic maps (IGME-1 : 50 000 scale) and bridged the temporal gap between the late Miocene age metamorphic fabrics and the present-day deformation field recorded by GPS geodesy and recent seismicity. At each site, three oriented samples were collected as independent blocks. Before removal from the outcrop, most blocks were hardened with a diluted (∼ 50 %) aqueous solution of sodium silicate (Fig. 4). In the laboratory, two to three oriented cubes (8 cm 3 ) were cut from each block using nonmagnetic Teflon knives and enclosed in standard cubic paleomagnetic boxes. The anisotropy of magnetic susceptibility (AMS) was determined with an Agico KLY-3S Kappabridge at Lehigh University. To determine magnetic mineralogy, a heating stage under the presence of an argon atmosphere and a cold stage accessory to the Kappabridge were used.

Results for Example I
Results from heating and cooling experiments show a complicated magnetic mineralogy composed of nearly 100 % ferromagnetic (magnetite or hematite) to nearly 100 % paramagnetic mineralogy (clays and iron-rich micas; Fig. 5). Since the kinematic interpretation of each of the specimen is the same regardless of magnetic mineralogy, the details of each specimen are not important for subsequent analysis. There is no correlation between the bulk magnetic susceptibility (k m ) and the anisotropy of the magnetic ellipsoid (P j ), so a comparison of the principal axis of susceptibility across the various structural positions around the Sierra Nevada massif specimens can provide useful kinematic information (Fig. 6a). Site 6 has a much higher magnetic susceptibility than the other sites because of possible secondary sulfide minerals at the site indicated by the heating and cooling behavior of the MS vs. T experiments.
Nearly all AMS ellipsoids are characterized by a low anisotropy degree (P j ) and oblate ellipsoid shape (T ) (Jelinek, 1981; Fig. 6b). The AMS axes' determinations record nearly the same axis orientations. At all sites around the Sierra Nevada, the minimum principal axes, k 3 , are nearly orthogonal to bedding. The principal elongation axes means are preferentially oriented NNE-SSW to NE-SW (Fig. 3).  The same is done on a subhorizontal layer of a poorly cemented, fine calcareous sandstone from an upper Middle Pleistocene fluvial terrace exposed in a wine cellar at the Geological Observatory of Coldigioco, Northern Apennines. Both samples were hardened with a dilute sodium silicate solution. Three or four oriented blocks were collected from each sampling site. Samples were oriented with a Brunton compass and located with a handheld GPS receiver, labeled, and photographed. Results from all measurements indicate that the magnetic susceptibility of the Spanish samples varies from being dominated by paramagnetic to ferromagnetic mineral grains. Irrespective of whether the greatest magnetic lineation, k 1 , is controlled by grain shape of crystallography, the kinematic interpretation is the same. (c, d) High temperature (MS vs. T ) measurements showing heating from room temperature (20 • C) to 700 • C and subsequent cooling back to room temperature. All four plots show evidence of the ferromagnetic mineral magnetite (Curie temperature of 580 • C). A lower temperature phase is indicated at Site 3, possibly maghemite. Site 6 shows the formation of additional magnetite during heating because of the much stronger susceptibility upon cooling. Heating curves are in red and cooling curves in blue.
The orientation of the site-mean magnetic susceptibility axes, k 1 , is horizontal or very shallowly plunging to the NE or SW (Fig. 3). In general, k 1 and k 2 are in or near the bedding plane of the specimens, and k 2 and k 3 do not form a girdle pattern in this principle plane.

Discussion of Example I
The AMS principal axes show a consistency between sites (Fig. 3), so we combine the susceptibility axes' orientation data in Fig. 7. These combined data suggest that during deposition the phyllosilicate grains were oriented with their basal planes parallel or slightly imbricated to the depositional surface. Compaction during dewatering and lithification amplified the initial oblate depositional fabric and was coincident with the formation of the tectonic fabric. Regardless of the magnetic mineralogy of the specimens, a well-clustered minimum susceptibility axis (k 3 ) is present, which we interpret as a compaction fabric in these sedimentary deposits. The possibility of a primary depositional current fabric (imbrication) is unlikely because of an independent paleocurrent study on clast imbrication at Site 3 and Site 4, which shows an eastward rather than westward transport direction during deposition (Carrigan et al., 2018).
Irrespective of the structural position around the Sierra Nevada massif, all sites show a preferred orientation of k 1 . The mean principal axis of maximum susceptibility is preferentially oriented at 30-210 • (Figs. 7 and 8). We interpret this as a tectonic fabric due to the tight clustering of k 1 and k 2 , the relationship between k 1 and strike of dipping bedding at sites SN1, SN4, and SN6, and the lack of influence from depositional processes. In specimens dominated by phyllosilicate grains, it is difficult to create a strong lineation by aligning grain crystallographic axes; however, an intersection lineation between slightly rotated clay grains orthogonal to a Jelinek diagram of Sierra Nevada specimens colored by site and consistent with Fig. 3 colors. All AMS measurement have a low anisotropy (less than 12 % P j ) and nearly all specimens are oblate (T >0). T and P j are calculated as follows: if n 1 = ln(t 1 ), n 2 = ln(t 2 ), n 3 = ln(t 3 ), where t 1 , t 2 , and t 3 are the eigenvalues, then T = (2n 2 − n 1 − n 3 )/(n 1 − n 3 ) and P = exp( (2[(n 1 − n mean )2 + (n 2 − n mean )2 + (n 3 − n mean )2])) and n mean = (n 1 + n 2 + n 3 )/3 (Jelinek, 1981). Figure 7. All Sierra Nevada massif AMS data. Lower hemisphere, stereographic projection of the principal axes of susceptibility orientations for all specimens determined from AMS measurements in tilt-corrected coordinates (Fig. 3). Arrows outside the stereonet periphery are parallel to the mean long axis (k 1 ) orientation. k 1 is the maximum axis, k 2 is the intermediate axis, and k 3 is the minimum axis.
shortening direction has been observed (Henry, 1997;Cifelli et al., 2004;Parés et al., 2007;Martín-Hermández and Ferré, 2007;Borradaile and Jackson, 2010). The orientation of k 1 is consistent with the present-day GPS velocity field, being oriented almost perfectly orthogonal to the direction of convergence of the Betic Cordillera to stable Africa (Nubia; Fig. 2; Gutscher et al., 2012), in good agreement with the mineral lineations recorded in the massif's core (Martínez-Martínez et al., 2002), and the Neogene brittle extensional structures and recent seismicity (Mancilla et al., 2013) in the orogen (Fig. 2). Because of the low strains and the orthogonal relationship between contractional and extensional principal directions, it is not possible to distinguish the uplift processes of the Sierra Nevada massif with our results. The AMS ellipsoid orientations, mineralogic stretching lineation from the core of the Sierra Nevada massif, the nearby GPS velocity field, and recent fault slip all have orientations consistent with the same strain field (Fig. 8). The principal elongation direction is interpreted to have persisted across different structural levels from Miocene time to the present (> 10 Myr). The principal elongation direction, k 1 , was collected from only young sediments so this fabric must also be young. The AMS fabric points out that the various phases of deformation affecting the Betic Cordillera were nearly coaxial since the Miocene.

Geologic setting of the Northern Apennines
The Northern Apennines are an accretionary fold and thrust belt (Bally et al., 1986) where crustal deformation, rock uplift, and topographic growth result from the ongoing subduction of Adria beneath Europe (Picotti and Pazzaglia, 2008;Carminati and Doglioni, 2012). The Apennine orogenic wedge initiated ∼ 30 Ma along the southern flank of the Alps (Le Pichon et al., 1971) and has grown at variable rates through the Neogene dependent on the transfer of mass imbricated from the subducting plate (Picotti and Pazzaglia, 2008). Rapid rollback of Adria with respect to Europe results in retreat and stretching of the upper plate, forming a wide zone of back arc crustal extension. The Apennine wedge started to become emergent ∼ 4 Ma (Picotti and Pazzaglia, 2008), uplifting and exposing paired compressional and extensional deformation fronts near the trench and in the forearc, respectively, with the structural transition near the topographic culmination of the range (D' Agostino et al., 2001;Carminati and Doglioni, 2012). Balanced cross sections for the Apennines (Bally et al., 1986;Hill and Hayward, 1988) indicate ∼ 130 to 150 km of subduction over the 30 Myr history of the wedge, which indicates relatively slow long-term rates at ∼ 4 to 5 km Myr −1 (4-5 mm yr −1 ), similar to the GPS geodetic rates (Devoti et al., 2008;Caporali et al., 2011;Bennett et al., 2012).
The northeastern Apennines, including the Umbria-Marche target region of this research, exposes Mesozoicearly Cenozoic carbonates and middle-late Cenozoic mixed carbonate-siliciclastic rocks folded and imbricated into northeast-vergent thrust sheets (Porreca et al., 2018, and references therein; Fig. 9). In Marche, these thrust sheets are located with carbonate ridges and have inferred blind thrusts in their cores (Artoni, 2013). Further west in Umbria, the thrust sheets are dissected by both east-and west-dipping high-angle normal faults Fig. 9). Ongoing thrust earthquakes beneath the Po Plain and Adriatic Sea (Pondrelli et al., 2006;Boccaletti et al., 2011) and normalfault-sense earthquakes beneath the high Apennines (Lavecchia et al., 1994;Doglioni et al., 1999;Ghisetti and Vezzani, 2002;Chiaraluce et al., 2017) speak to concurrent shortening and extension in the wedge.
The paired deformation fronts in the Northern Apennines, Italy, are convolved with an enigmatic, but active, eastdipping (towards Adria), 14-15 km deep detachment called the Alto Tiberina fault that projects to the surface west of the Apennine crest Pialli et al., 1998;Boncio et al., 2004;Chiaraluce et al., 1999;Eva et al., 2014;Lavecchia et al., 2016;Fig. 9). This detachment is one of only a handful of low-angle normal faults globally that are demonstrably seismogenic (Hreinsdóttir and Bennett, 2009;Valoroso et al., 2017), apparently in contradiction to frictional fault reactivation theory that predicts that slip on lowangle normal faults as extremely unlikely (reviewed in Collettini, 2011). Most of the destructive seismicity in the high Apennines tends to nucleate on west-dipping high-angle normal faults that are antithetic to and sole into this east-dipping detachment (Galadini and Galli, 2000;Boncio et al., 2004;Roberts and Michetti, 2004). The most destructive seismicity, including the 2016-2017 earthquake sequence, is tightly focused along the highest crest of the Apennines where it is co-located with young, underfilled, extensional basins, highangle normal faults that rupture the surface (Fig. 9) and geomorphic evidence for an east-marching drainage divide. It is not known if the infrequent but large historic earthquakes east of the divide are indicative of new blind normal faults that have nucleated on the detachment, represent active shortening, or alternatively are responding to a different stress field.
Imbricated foredeep and wedge-top basins contain a timetransgressive range of poorly consolidated deposits that span the compressional and extensional regimes. Conceivably, shortening fabrics could be recorded in lithofacies at the base of one of these basins when it was formed and filled in the shortening part of the wedge, only to be superseded by stretching fabrics in overlying lithofacies as the basin was translated westward and into the extending part of the wedge. Adriatic slope transverse rivers (Alvarez, 1999) traverse both the extending and shortening parts of the wedge and contain Pleistocene alluvial deposits representing an AMS geodetic snapshot of the current crustal strains. Published AMS data from the thrust belt show a strike-parallel (NE-SW and horizontal) extension that is perpendicular to compression and shortening directions (Sagnotti, et al., 1998;Caricchi et al., 2016). To confirm these data towards the southeast and to better locate the kinematic transition region between the contracting and extending regions of the overlying Eurasian plate, we sampled AMS data in Oligocene and younger units, Figure 9. (a) Location map showing the topography, known major faults with black lines, the location of large, historic earthquakes in orange circles (from Boncio et al., 1998), the drainage divide as a red line, and GPS geodetic velocities with uncertainties in red arrows (from Hreinsdóttir and Bennett, 2009)  including Quaternary deposits in a NE-SW-oriented corridor across the thrust belt (Fig. 9).

Methods for Example II
Sampling in the Apennines was designed to identify the location of the modern extensional front. Field collection and specimen preparation occurred as in Example I from Spain, with unconsolidated samples being hardened with sodium silicate before or just after orienting and removal from the outcrop (Fig. 4). We collected samples from 17 sites from sedimentary rocks and poorly consolidated sediments from Late Eocene to late Pleistocene age, with a focus on late Miocene-Pliocene argillaceous marine deposits (Table A1). The Italian specimens were prepared and rock magnetic data were acquired in the Archeomagnetism Laboratory at Centro Nacional de Investigación de la Evolución Humana (CE-NIEH) (Spain). The AMS of the collected specimens was measured on a MFK1-FA Kappabridge (Agico Instruments), a fully automated inductive bridge, at a frequency of 976 Hz Paramagnetic susceptibility clearly dominates all the specimens, as revealed by the slope of the loops. (b) Example of a specimen where the paramagnetic contribution has been removed in order to enhance the ferromagnetic contribution (loop in black). (c) Example of a specimen where diamagnetism dominates the total magnetic susceptibility. and a field of 200 A m −1 . Analysis software (Saphyr6, by Agico) creates a complete susceptibility tensor. Rock magnetic measurements included isothermal remanent magnetization (IRM) acquisition experiments up to 1 Tesla and hysteresis curves to determine the relative contribution of ferromagnetism and paramagnetism to the total susceptibility tensor. These experiments were carried out with a vibrating sample magnetometer (VSM; Micromag 3900).

Results of Example II
Samples from the Apennines have variable magnetic mineralogy and include a wider range of lithologies and ages than the Betics sampling. Samples from sites AP2 and AP7 (Bisciaro Fm.) are dominated by diamagnetic calcite and negative mean susceptibility, which precludes any meaningful analysis of the AMS axes' orientations. At the 1T field, the magnetization was not fully saturated, indicating the presence of hematite in addition to lower coercivity magnetite as the dominant ferromagnetic components (Heller, 1978). Still, the bulk magnetic susceptibility is dominated by paramagnetism, as revealed by the hysteresis curves (Fig. 10). The contribution of paramagnetic grains suggests that the measured magnetic fabric can be used as a proxy for phyllosilicate grains' preferred orientation; therefore, the AMS principal axes are indicators of the orientation of the strain axes' orientation (e.g., Soto et al., 2009).
Representative examples of AMS fabrics are shown in Fig. 11. The mean susceptibility shows no positive correlation with the shape parameter or anisotropy degree (T , P j ; Fig. 12). Similar to the data from Spain, the AMS ellipsoids from the Italian specimens indicate low P j values, revealing a low degree of grain shape preferred orientation and low strains. The AMS axes' distribution is particularly clear in specimens of the argillaceous and semi-consolidated Pliocene Argille Azzurre Fm. At all sites, k 1 axes' orientations are shown as a function of rock formation, as well as the sites at which k 3 is perpendicular to bedding (Fig. 13). All interpretable specimens from the Apennine range samples, including the Pleistocene fluvial deposits, generate a site-mean AMS fabric consistent with contraction and shortening in the wedge.

Discussion for Example II
The k 1 axis orientation is orthogonal to the rock transport and crustal shortening directions, as recorded in GPS geodesy data and seismology (Fig. 9). Irrespective of sample age, we interpret AMS ellipsoids that have the magnetic lineation in a NW-SE orientation as recording contraction, as this is the main trend of the fault traces and strike of bedding and topography (Fig. 13). A few sites do not provide interpretable kinematic results because the axes directions are scattered and suggest inconsistent strain directions. The calcareous marls of the Bisciaro Fm. (AP2, AP7) have a poorly formed AMS fabric. In these specimens, the mean susceptibility is negative and dominated by diamagnetism, most likely calcite. The absence of a compactional fabric in carbonate-dominated specimens (AP2, AP7) likely indicates that these sediments lithified by cementation soon after deposition.
In general, the distribution of the principal axes of the AMS ellipsoid does not significantly vary with stratigraphic age or structural position. For example, the oldest specimens collected from Eocene-middle Miocene marls and Pliocene siliciclastic rocks (AP6, AP14, AP17) uniformly show AMS fabrics consistent with contractional deformation of the orogenic wedge (Fig. 13), similar to the results of Sagnotti et al. (1998). Most importantly, sites collected from thrust structures that are currently in an extending regime (AP11, AP12, AP13) imply that either the AMS fabric was locked after the original deformation due to the high strain required to rotate grain pairs, or that subsequent extension has not affected the previous AMS fabric (e.g., Larrasoña et al., 2004). The same is true for middle and late Miocene siliciclastic deposits next to the Marche ridge (AP3, AP9), where the current orientations of crustal stresses from fault and earthquake data are ambiguous. Pliocene and Pleistocene samples from near the toe of the orogenic wedge show an orientation consistent with ongoing shortening (AP4, AP5, AP8). Wegmann and Pazzaglia (2009) also report ongoing shortening in this region as evidenced by fluvial terrace folding above the Filottrano thrust, which we cross at the location of AP4.
The kinematic transition zone in central Italy aligns with the topography, the seismicity (Pondrelli et al., 2006), and the GPS geodesy (Bennett et al., 2012;Fig. 9). Our AMS data do not improve on the location of the transition zone because of the lack of samples from Plio-Pleistocene deposits directly northeast of the drainage divide (Fig. 13). Unfortunately, the one Pleistocene river terrace deposit northeast of the divide (AP10) has indeterminate axes. As such, our AMS results are not able to support the idea that there is an apparent rotation of the principal compressive stress between the Adriatic coast and the Marche ridge associated with wedge-scale pore-pressure variations (Peacock et al., 2017). Furthermore, the AMS is unable to determine the stress field responsible for the large historic earthquakes in the region between the drainage divide and the Marche ridge. If earthquakes in the region are related to blind normal faults with tips breaking   Rovida et al. (2020). The presence of a tectonic fabric was determined by clustering of k 1 declinations outside of the expected compaction fabric. Axis certainty represents the percentage of specimens of the total used to calculate a mean k 1 vector. Right legend: (1) Holocene fill; (2) first-order Quaternary terrace (Q t 1 ); (3) second-order Quaternary terrace (Q t 2 ); (4) third-order Quaternary terrace (Q t 3 ); (5) Argille Azzurre Fm.; (6) Scaglia Rossa Fm.; (7) Maiolica Fm.; (8) Bisciaro Fm.; (9) thrust fault trace; (10) normal fault trace; (11) Alto Tiberina detachment; (12) drainage divide; (13) large historic but pre-instrument earthquakes (pre-1800) of unknown origin (see Fig. 9). upsection from the Alto Tiberina detachment (Fig. 9), a possible rationale is that, according to extensional critical wedge theory (Davis et al., 1983), a wedge with a taper greater than some critical value is unable to slide over its basal detachment until sufficient wedge thinning on connecting faults reduces the surface slope and wedge taper below the critical value (Xiao et al., 1991). Suitable deposits do outcrop in this critical region, so additional field work and AMS analyses may shed light on this problem.

Conclusions
The AMS technique provides an effective way to identify both modern and paleokinematics from sediments and sedimentary rocks largely independent of the magnetic mineralogy of a specimen. Stratigraphically controlled AMS measurements are a deep-time paleogeodetic technique that can be combined with structural geology, GPS geodesy, and seismic data to collectively describe the kinematics of active orogens and to better understand the nature of seismic hazards. In both the Betic Cordillera (Example I) and Northern Apennines (Example II), weak but well-organized penetrative AMS fabrics were recovered from young unconsolidated and unburied rocks that could not be analyzed with more traditional methods. In the Betic Cordillera, we established a long-term consistency to the strain field from the late Miocene to the present from unburied, young deposits around the Sierra Nevada. For the Northern Apennines, all studied sites, regardless of the sites' stratigraphic age, ubiquitously record NW-SE-oriented k 1 axes' orientations, irrespective of structural position. Contractional strains in the most southwest-located samples are likely locked into the rocks and do not record superposed penetrative extension. In any case, the recovered magnetic fabric orientation successfully determined the kinematics of an area near the synorogenic surface in the still-contracting orogen toe region.  Code availability. The Pmagpy code is available at https:// earthref.org/PmagPy/cookbook/ (Magnetics Information Consortium, 2020). PmagPy code was used to analyze and plot AMS data. The Anisoft code is available at https://www.agico.com/text/ software/anisoft/anisoft.php (AGICO, 2020).
Data availability. The datasets used in this paper are available upon request.
Author contributions. DJA, JMP, and CB conceived the Spanish project and completed sampling, sample preparation, measurement, and analyses. DJA and JMP conceived the Italian project. DJA, FJP, AM, and LKC completed the Italian sampling. DJA and JMP prepared the Italian specimens, measured the samples, and analyzed the results. DJA, FJP, JAF, CB, and KPK analyzed results and drafted figures for the manuscript. DJA and FJP wrote the first draft of the manuscript and edited each subsequent draft. JMP, KPK, CB, and AM edited multiple drafts of the manuscript. DJA completed the final edits.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Tools, data and models for 3-D seismotectonics: Italy as a key natural laboratory". It is a result of the workshop on "Tools, data and models for 3D seismotectonics: the Italian laboratory over time", Perugia, Italy, 9-10 July 2019.