Spatial variations in the similarity of earthquake populations: The case of the 1997 Colfiorito–Sellano (northern Apennines, Italy) seismic sequence

The spatial properties of events in the 1997 Colfiorito–Sellano seismic sequence (Northern Apennines, Italy) were investigated using coherence, a parameter derived from seismic moment tensors that quantifies the kinematic similarity between focal mechanisms. The 1997 Colfiorito–Sellano seismic sequence predominantly consists of normal faulting earthquakes, with a few strike-slip and reverse faulting episodes. This kinematic heterogeneity is possibly related to the contemporaneous activity of two different sets of faults: NW–SE normal faults and NNE–SSW sub-vertical faults, the latter inherited from the previous Miocene compressional phase. The study used two independently-derived data sets of the same seismic sequence characterized by a different number of events and by different precision of spatial localisation. Their statistical significances, assessed through a reshuffling procedure, reveal that data sets with at least some hundreds of events and good positional precision are required to obtain significant results through coherence analysis. Results from the better quality data set indicate that this seismic sequence is characterized by a rapid decrease in the kinematic similarity between earthquake pairs within 2 km of separation, particularly along directions sub-perpendicular to the normal fault strike. The decrease rate seems to be controlled by the geometric characteristics of the normal faults, given that the mean along-dip distance between fault segments is 2 km. In proximity to pre-existing tectonic lineaments the relative abundance of strike-slip and reverse faults tends to decrease the kinematic similarity between events but does not influence the coherence decrease rate. The presence of mixed focal mechanisms (normal, reverse and strike-slip) in a single seismic phase implies that mixed fault types are not restricted to polyphase tectonic histories: such heterogeneous kinematics during a single phase may be induced by the presence of inherited discontinuities.     

Effect of Aftershocks on Earthquake Hazard Estimation: An Example from the North Anatolian Fault Zone

In order to investigate the effect of aftershocks on earthquake hazard estimation, earthquake hazard parameters (_m, b and M_max) have been estimated by the maximum likelihood method from the main shocks catalogue and the raw earthquakes catalogue for the North Anatolian Fault Zone (NAFZ). The main shocks catalogue has been compiled from the raw earthquake catalogue by eliminating the aftershocks using the window method. The raw earthquake catalogue consisted of instrumentally detected earthquakes between 1900 and 1992, and historical earthquakes that occurred between 1000–1900. For the events of the mainshock catalogue the Poisson process is valid and for the raw earthquake catalogue it does not fit. The paper demonstrates differences in the hazard outputs if on one hand the main catalogues and on the other hand the raw catalogue is used. The maximum likelihood method which allows the use of the mixed earthquake catalogue containing incomplete (historical) and complete (instrumental) earthquake data is used to determine the earthquake hazard parameters. The maximum regional magnitude (M_max, the seismic activity rate (_m), the mean return period (R) and the b value of the magnitude-frequency relation have been estimated for the 24°–31° E, 31°–41° E, 41°–45° E sections of the North Anatolian Fault Zone from the raw earthquake catalogue and the main shocks catalogue. Our results indicate that inclusion of aftershocks changes the b value and the seismic activity rate _m depending on the proportion of aftershocks in a region while it does not significantly effect the value of the maximum regional magnitude since it is related to the maximum observed magnitude. These changes in the earthquake hazard parameters caused the return periods to be over- and underestimated for smaller and larger events, respectively.     

Probabilistic assessment of seismic hazard of dam sites in Jordan

A Probabilistic method is used to evaluate the seismic hazard of nineteen embankment dam sites in Jordan. A line source model developed by McGuire (1978) is used in this study. An updated earthquake catalogue covering the period from 1 A.D. to 1991 A.D. is used for this purpose. This catalogue includes all earthquakes that occurred in Jordan and adjacent areas, more specifically between latitudes 27.0°–35.5° N and longitudes 32.0°–39.0° E.Nine distinct seismic sources of potential seismic activities are identified. The seismic hazard parameters are determined using the method suggested by Kijko and Sellevoll (1989).The Peak Ground Acceleration (PGA) is selected as a measure of ground motion severity. Esteva (1974) attenuation relationship is used in evaluating PGA values at each dam site. Analysis is carried out for 50%, 90%, and 95% probability that is not being exceeded in a life time of 50, 100, and 200 years.Results of analysis indicate that PGA values are higher for dam sites closer to the Dead Sea Fault. This fault is believed to be responsible for most earthquake activities in Jordan and vicinity. The highest PGA value is found to be for Al-Karama dam site.     

Analysis of input data and seismic hazard assessment for the low seismicity area ‘Belgium,The Netherlands and NW Germany’

From the earthquake catalogue of the region of interest 2–8°E, 49–52°N, only events with intensity > 3 are considered. The fore/aftershocks and the induced seismicity are removed. The completeness and other properties of the remaining data set are investigated. The seismic data, as well as the geological, tectonic, and other geophysical data of the area are used for the determination of the seismotectonic zones. For each zone, the cumulative intensity-frequency relation (taking into account the completeness of the catalogue), the attenuation depending on the direction, the upper bound of intensity and the average depth are calculated. When the seismic activity is changing within a zone, this zone is divided into subzones with an almost constant activity and a b value equal to the value of the whole zone. When necessary for each subzone, a different attenuation coefficient can be used, to take into account more regional effects. Then the seismic hazard is assessed with the modified McGuire program. Annual probability versus intensity for several sites and hazard intensity and probability maps are calculated. Finally, an error discussion of the whole procedure for one test site is given.     

Density distribution functions and correlations of earthquake parameters

The statistical analysis of the source parameters of 9 earthquake sets of different types (aftershocks, scattered events, swarms) and of different seismic regions shows that the density distribution function (ddf) of the linear dimension l of a fault is represented by a negative power law, as well as the ddf of the static stress drop σ and of the scalar seismic moment M_o. It is then suggested, and tentatively verified, that also the ddf of the root mean square ground acceleration, defined as a function of l and σ, may be represented by a negative power law and that, at least in the cases examined, it scales like the ddf of σ. It is seen that the variability of the static stress drop is significant from one region to another, as is well known, but it seems remarkable also in the same seismic region (in particular in California, σ varies by several orders of magnitude) and in the different sets of events of a given region (as observed again for California). It is hypothesized that a correlation, although weak, between the stress drop and the linear dimension of a fault exists and the analyses seem not to contradict that σ may be a decreasing function of l. Finally, it is suggested that the seismicity of a region may be represented two-dimensionally as a function of the ddf of the stress drop and of the linear dimension of a fault instead of the classic b and b_o values.     

Spatial relationships between natural seismicity and faults, southeastern Ontario and north-central New York state

A statistical analysis was carried out to investigate spatial associations between natural seismicity and faults in southeastern Ontario and north-central New York State (between 73°18′ and 77°00′W and 43°30′ and 45°18′N). The study area is situated to the west of the seismically active St. Lawrence fault zone, and to the east of the Lake Ontario basin where recently documented geological and geophysical evidence points to possible neotectonic faulting. The weights of evidence method was used to judge the spatial associations between seismic events and populations of faults in eight arbitrarily defined orientation groups. Spatial analysis of data sets for seismic events in the periods 1930–1970 and post-1970 suggest stronger spatial associations between earthquake epicentres and faults with strikes that lie in the NW–SE quadrants, and weaker spatial associations of epicentres with faults that have strikes in the NE–SW quadrants. The strongest spatial associations were determined for groups of faults with strikes between 101° and 146°. The results suggest that faults striking broadly NW–SE, at high angles to the regional maximum horizontal compressive stress, are statistically more likely to be spatially associated with seismic events than faults striking broadly NE–SW. If the positive spatial associations can be interpreted as indicating genetic relationships between earthquakes and mapped faults, then the results may suggest that, as a population, NW–SE trending faults are more likely to be seismically active than NE–SW striking faults. Detailed geological studies of faults in the study area would be required to determine possible neotectonic displacements and the kinematics of the displacements.     

Estimation of earthquake hazard parameters for incomplete and uncertain data files

The maximum likelihood estimation of earthquake hazard parameters (maximum regional magnitudem _max, activity rate , and theb parameter in the Gutenberg-Richter distribution) is extended to the cases of incomplete and uncertain data. The method accepts mixed data containing only large (extreme) events and a variable quality of complete data with different threshold magnitude values. Uncertainty of earthquake magnitude is specified by two values, the lower and upper magnitude limits. It is assumed that such an interval contains the real unknown magnitude. The proposed approach allows the combination of different quality catalog parts, e.g. those where the assignment of magnitude is questionable and those with magnitudes precisely determined.As an illustration of the method, the seismic hazard analysis for western Norway and adjacent sea area (4–8°E, 58–64°N) is presented on the basis of the strongest earthquakes felt during the period 1831–1889 and three complete catalog parts, covering the period 1890–1987.Paper presented at the 21st General Assembly of the European Seismological Commission held in Sofia, 1988.On leave from Institute of Geophysics, Polish Academy of Sciences, Warsaw, Poland.     

Local seismic array observations at north Evoikos, central Greece, delineate crustal deformation between the North Aegean Trough and Corinthiakos Rift

In order to better constrain and define the microseismic activity at the north Evoikos Gulf and its surrounding area we deployed an onshore/offshore seismic array consisting of 31 three-component seismic digital stations. The array was active from 30 June to 24 October 2003, and covered an area of 2500 km². We located more than 2000 seismic events ranging from 0.7 to 4.5 M_L by using six stations as a minimum in order to define the foci parameters. Recorded seismicity delineated three major zones of deformation: from south to north, the Eretria–Parnis–eastern Corinthiakos zone, the Psachna–Viotia zone, and the Northern Sporades–North Evia–Bralos zone. Alignments of the recorded seismicity follow the tectonic trends and their orientation in the above zones. The whole area accommodates the stress field between the North Aegean Trough and the Corinthiakos Gulf. Rate of deformation intensifies from north to south, as revealed also by historical and instrumental seismicity. The successive change of orientation between the two stress fields fragments the crust in relatively small units and the fault systems developed do not permit the generation of major earthquakes in the north Evoikos area and its immediate vicinity. This is also supported by the instrumental seismicity of the last century. Larger events reported in historical times are probably overestimated.Most seismic activity is crustal. Subcrustal events were recorded mainly below the Lichades area and are interpreted as the consequence of the subduction of the Ionian oceanic lithosphere below the Hellenides. The Lichades volcano is the most northern end of the Hellenic volcanic arc.At present the highest seismic activity is associated with the Psachna region of north Evia that has been continuously active since 2001. Considering, however, the development of the seismic activity during the last decade, there has been a sequence of large events, i.e., Parnis in 1999, Skyros in 2001 and Psachna in 2001–2003. This demonstrates the fact that the tectonic deformation in all this area is intense and important for the accommodation of the stress field of the North Aegean Trough to that of the Corinthiakos Rift.     

Integration of the THRUST project into the Chile Tsunami Warning System

Several seismic gaps are identified along the coast of Chile. The main one is located in the northern part of the country (18° S–20° S), where a THRUST seismic trigger has been installed; two more triggers will be placed in other seismic gaps. The Standard Operations Plan (SOP) was tested during a disaster simulation exercise based on a major earthquake and tsunami situation. Response actions by governmental and civil agencies were monitored and a performance evaluation was done. Modifications of the SOP were found necessary to adjust the interactions between agencies.     

Double seismic zone beneath the middle of Hokkaido, Japan, in the southwestern side of the Kurile Arc

The vertical section of microearthquakes, determined accurately by using the Hokkaido University network, shows two dipping zones (the double seismic zone) 25–30 km apart in the depth range of 80–150 km beneath the middle of Hokkaido in the southwestern side of the Kurile arc. Hypocentral distribution of large earthquakes (m_b > 4) based on the ISC (International Seismological Centre) bulletin also shows the double seismic zone beneath the same region. The hypocentral distribution indicates that the frequency of events occurring in the lower zone is four times greater than that in the upper zone. The difference in seismic activity between the two zones beneath Hokkaido is in contrast with the region beneath northeastern Honshu in the northeastern Japan arc.Composite focal mechanisms of microearthquakes and individual mechanisms of large events mainly characterize the down-dip extension for the lower zone as is observed beneath northeastern Honshu. For the upper zone, however, the stress field is rather complex and not necessarily similar to that beneath northeastern Honshu. This may be considered to indicate the influence of slab contortion or transformation in the Hokkaido corner between the Kurile and the northeastern Japan arcs.     

An interpretation and catalogue of paleoseismicity in Sweden

Paleoseismic data provide a long-term record of seismic activity to predict hazards for periods longer than one to a few centuries. In Sweden, the analysis reveals there was a drastic change in dominant seismic mode from a high to super-high deglacial mode to a low to moderately low mode in present and Late Holocene time. Paleoseismic criteria and characteristics include numerous different sources of information; viz. primary faults, bedrock deformation, sedimentary deformation, rock and sediment slides, liquefaction, sorting by shaking, tsunamis, differing geomorphic expressions, disordering and ordering of magnetic particles. By applying multiple criteria, it was possible to identify 44 paleoseismic events, including 23 events of estimated M 6–7, 12 events of M 7–8 and 6 events of M > 8. Varve-dating often allows a precision as to a single year, in one case even to the season of a year. The key for paleoseismic reconstruction and testing is the application of multiple criteria.     

Seismic potential of Southern Italy

To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.     

A New Model for Quantifying Anisotropic Scale Invariance and for Decomposition of Mixing Patterns

A new power–law function has been derived to represent the relationship between area of the set consisting of wave numbers with spectral energy density above S (A(>S)) on the two-dimensional frequency plane and S. The power–law relation holds if the field concerned possessing isotropic scale invariance or generalized scaling invariance involves rotational and ratio-scale changing transforms. The equation is valid for dealing with common exploration geophysical and geochemical fields encountered in mineral exploration and environmental assessment. This power–law function not only provides a new model for characterizing anisotropic scaling invariance for generalized scaling field, for example, estimating the power exponent of power spectrum of generalized scale invariance measure in frequency domain, but also forms a theoretical base for the S–A filtering technique developed for decomposing a mixing field into components on the basis of distinct scaling properties in the frequency domain. It is demonstrated that the method has potential to become a general technique for image processing and pattern recognition.     

Dependence of earthquake recurrence times and independence of magnitudes on seismicity history

The fulfillment of a scaling law for earthquake recurrence–time distributions is a clear indication of the importance of correlations in the structure of seismicity. In order to characterize these correlations we measure conditional recurrence–time and magnitude distributions for worldwide seismicity as well as for Southern California during stationary periods. Disregarding the spatial structure, we conclude that the relevant correlations in seismicity are those of the recurrence time with previous recurrence times and magnitudes; in the latter case, the conditional distribution verifies a scaling relation depending on the difference between the magnitudes of the two events defining the recurrence time. In contrast, with our present resolution, magnitude seems to be independent on the history contained in the seismic catalogs (except perhaps for Southern California for very short time scales, less than about 30 min for the magnitude ranges analyzed).     

Fluid inclusion study of Xihuashan tungsten deposit in the southern Jiangxi province,China

The Xihuashan tungsten deposit is closely related to a small highly evolved granitic intrusion. The fluid phases associated with the wolframite-bearing quartz veins have been investigated using microthermometry and the Raman microprobe; they are highly variable in density and composition. The earlier fluids are low-density and low-salinity CO₂-bearing aqueous solutions circulating at temperatures up to 420 °C, and low-salinity (2–3 equiv. wt% NaCl) aqueous solutions without traces of CO₂ circulating at high temperatures 280°–400 °C) involved in a specific hydrothermal fracturing event; limited unmixing occurs at 380 °C and 200–100 bar in response to a sudden pressure drop. The second types of fluids related to deposition of idiomorphic drusy quartz are typical CO₂-bearing aqueous solutions with low salinity (2.5 equiv. wt% NaCl) homogenizing at low to moderate temperatures (180°–340 °C). The late fluids characterize the sulfide deposition stage; they are aqueous fluids with variable salinities homogenizing in the liquid phase between 100° and 275 °C. The Xihuashan hydrothermal evolution resulted from a discontinuous sequence of specific events occurring between 420° and 150 °C and during a continuous hydrothermal evolution of the system during cooling. The role played by the CO₂-rich fluids in the transport and deposition of tungsten in the hydrothermal environment is discussed.     

Interpretation of tectonic versus magmatic processes for resetting apatite fission track ages in the Niğde Massif, Turkey

A network of 10 OBS operating in continuous recording mode was deployed in June 2003 in the Galicia passive Margin, surrounding the sinking zone of the ‘Prestige’ oil tanker. The OBSs recorded very different signals, including air-gun shots, seismic events, whale vocalizations or noise associated to ships, as well as two sets of atypical seismic signals that are the object of this contribution and can be described as harmonic tremors and short-duration events. Harmonic tremors appear in most sites as intense monochromatic signals lasting for up to 4 h, with coherent polarization and a wide range of amplitudes. Even if they may be recorded at different sites within a time interval, their envelopes are not correlated, suggesting a source origin at the vicinity of each receiver. A time periodicity close to 6 h is observed, suggesting a correlation with the tidal variation. Most OBSs show also a large number of impulsive events of short-duration, between 1 and 2 s, a very regular decrease of amplitude in the signal coda and characteristic frequency spectrum with one or two narrow peaks. As for the tremors, those events are detected only in individual sites, suggesting very local source zones. We favor the interpretation of both seismic features as coming from the resonance of fluid-filled cracks induced by impulsive pressure transients, in agreement with the theoretical model of Chouet [Chouet, B. 1988. Resonance of a fluid-driven crack: radiation properties and implications for the source of long-period events and harmonic tremor. J. Geophys. Res., 93, 4375–4400.]. The observed harmonic tremors would be the result of sustained pressure fluctuations, probably related to the stress variations induced by the tidal change of the oceanic load, while the short duration events are interpreted as the impulse response to the tremor-generating system.     

Seismic activity at the western Pyrenean edge

The present-day seismicity at the westernmost part of the Pyrenean domain reported from permanent networks is of low to moderate magnitude. However, it is poorly constrained due to the scarce station coverage of the area. We present new seismic data collected from a temporary network deployed there for 17 months that provides an enhanced image of the seismic activity and its tectonic implications. Our results delineate the westward continuity of the E–W Pyrenean band of seismicity, through the Variscan Basque Massifs along the Leiza Fault, ending up at the Hendaya Fault. This seismicity belt is distributed on a crustal scale, dipping northward to almost 30 km depth. Other relevant seismic events located in the area can be related to the central segment of the Pamplona fault, and to different E–W thrust structures.     

Evaluation of Strong Motion Acceleration for Embankment Dam Design Considering Local Seismotectonics

A probabilistic method is used to evaluate the seismichazard of Adassiya dam site on the Yarmouk river in Jordan. A line source model developedby McGuire (1978) is used in this study. An updated earthquake catalogue coveringthe period from 1 A.D. to 1996 A.D. is used for this purpose. This catalogue includesall earthquakes that occurred in Jordan and adjacent areas, more specifically between latitudes27.0°–35.5°N and longitudes 32.0°–39.0°E.Nine distinct seismic sources of potential seismic activitiesare identified. The seismic hazard parameters are determined using the method suggested by Kijko and Sellevoll (1989).The Peak Ground Acceleration (PGA) is selected as a measure of ground motion severity. Esteva (1974) attenuation relationship is used in evaluating PGA values at each dam site. Analysis is carried out for 50%, 90%, and 95% probability that is not being exceeded in a life time of 50, 100, and 200 years.Results of analysis indicate that PGA values at the dam site are as follows:[] Operating Basis Earthquake (OBE) (50% probabilityof non-exceedance for a design life of 100 years – corresponding to a return period of 145 years) is 133.6 cm/sec².[] An earthquake with 90% probability of non-exceedancefor a design life of 50 years – corresponding to a return period of 475 years is 214.9 cm/sec².[] Maximum Credible Earthquake (MCE) (Return period of900 years) is 283.0 cm/sec².Strong motion acceleration time history of these earthquakes are givenbased on strong motion records of the November 1995 Gulf of Aqaba earthquake.Local site effect analysis for Adassiya Dam site using SHAKE program showed no amplification. Normalized site-specific acceleration response spectra for OBE and MCE design earthquakes is also given.     

Structure of the crust and uppermost mantle beneath southern Finland revealed by analysis of local events registered by the SVEKALAPKO seismic array

In 1998–1999, a large-scale seismic array was deployed in Finland as a part of the EUROPROBE/SVEKALAPKO subproject, involving 14 European universities and research institutes. The objective of the project was to map the deep lithosphere structure and thickness beneath the Fennoscandian Shield by means of teleseismic events. In addition, about 580 local seismic events were registered during the data acquisition period. Among them, only eight local earthquakes were recorded, the rest being quarry blasts from mining sites in Russia, Finland, Estonia and Sweden. In this study, we present the analysis of the seismic wave field from the strongest local events registered by the majority of the stations of the SVEcofennian–KArelian–LAPland–KOla Transect (SVEKALAPKO) array with the aim of mapping the structure of the upper mantle beneath the array. For this purpose, we selected the events corresponding to a single source type and compared these recordings with those from wide-angle reflection and refraction experiments in the area to identify the regional phases. The record sections of selected events demonstrate strong reflections (PmP) from the Moho boundary. The refracted Pn phases can be seen as first arrivals at distances of about 200–400 km from the source. At offsets of about 400–800 km, phases reflected from inhomogeneities in the uppermost mantle (P1) and double reflections from the Moho boundary (PmPPmP) were recorded.Results from 2D forward ray trace modeling of reflected and refracted P-waves along four profile swathes from SVEKALAPKO stations demonstrate that the mantle reflections originate from two different groups of boundaries beneath the array: one group of phases arrive from subhorizontal and gently dipping reflectors below the Moho boundary at a depth of 70–90 km, while the other group are phases originating from a depth of 100 to 130 km. Based on the irregular character of the first group of reflections, their different spatial orientation and correlation with the Moho offsets, we interpret the boundaries of this group as relicts of ancient subduction and collision processes. The second group of reflections can be explained by a transition from mechanically strong to mechanically weak lithosphere.     

The Sanfandila earthquake sequence of 1998, Queretaro, Mexico: activation of an undocumented fault in the northern edge of central Trans-Mexican Volcanic Belt

A sequence of small earthquakes occurred in Central Mexico, at the northern edge of the Trans-Mexican Volcanic Belt (TMVB) in the State of Queretaro, during the first 3 months of 1998. Medium to large events in the continental regime of central Mexico are not common, but the seismic history of the region demonstrates that faults there are capable of generating destructive events. The sequence was analyzed using data from a temporary network with the goals of identifying the causative fault and its relation to regional tectonics. Employing a waveform inversion scheme adapted from a method used for regional studies, we found that the source mechanisms conform to the style of faulting (i.e. extension in the E–W direction) representative of the Taxco–San Miguel Allende Fault system. This system has been proposed as the southernmost extension of the Basin and Range (BR) Province. The spatial distribution of hypocenters and source mechanisms indicate that the seismogenic segment was a fault with an azimuth of approximately 334° with almost pure dip slip. Since events which occurred just south from this region show features which are consistent with TMVB tectonics (i.e. extension in an N–S direction), the sequence may mark the boundary between the TMVB and BR stress domains.