Subsidence, sedimentation and sea-level changes in the Eromanga Basin, Australia

Abstract The Jurassic-Cretaceous subsidence history of the Eromanga Basin, a large intracratonic sedimentary basin in central eastern Australia, has been examined using standard backstripping techniques, allowing for porosity reduction by compaction and cementation. Interpretation of the results suggests that during the Jurassic the basin was subsiding in a manner consistent with the exponentially decreasing form predicted by simple thermally based tectonic models. By the Early Cretaceous, the rate of subsidence was considerably higher than that expected from such models and nearly half of the total sediment thickness was deposited over the final 20 Myr of the basin's 95 Myr Mesozoic depositional history. The Early Cretaceous also marks the first marine incursion into the basin, consistent with global sea-level curves. Subsequently, however, the sediments alternate between marine and non-marine, with up to 1200 m of fluvial sediments being deposited, and this was followed by a depositional hiatus of about 50 Myr in the Late Cretaceous. This occurred at a time when global sea-level was rising to its peak. A model is presented which is consistent with the rapid increase in tectonic subsidence rate and the transgressive-regressive nature of the sediments. The model incorporates a sediment influx which is greater than that predicted by the thermally based tectonic models implied by the Jurassic subsidence history. The excess sedimentation results in the basin region attaining an elevation which exceeds that of the contemporary sea-level, and thereby giving the appearance of a regression. The present day elevation of the region predicted by the model is about 100–200 m above that observed. This discrepancy may arise because the primary tectonic subsidence is better represented by a linear function of time rather than an exponentially decreasing form.     

Feedbacks of sedimentation on crustal heat flow: New insights from the Vøring Basin,Norwegian Sea

Basement heat flow is one of the key unknowns in sedimentary basin analysis. Its quantification is challenging not in the least due to the various feedback mechanisms between the basin and lithosphere processes. This study explores two main feedbacks, sediment blanketing and thinning of sediments during lithospheric stretching, in a series of synthetic models and a reconstruction case study from the Norwegian Sea. Three types of basin models are used: (1) a newly developed one‐dimensional (1D) forward model, (2) a decompaction/backstripping approach and (3) the commercial basin modelling software TECMOD2D for automated forward basin reconstructions. The blanketing effect of sedimentation is reviewed and systematically studied in a suite of 1D model runs. We find that even for moderate sedimentation rates (0.5 mm year^?1), basement heat flow is depressed by ～25% with respect to the case without sedimentation; for high sedimentation rates (1.5 mm year^?1), basement heat flow is depressed by ～50%. We have further compared different methods for computing sedimentation rates from the presently observed stratigraphy. Here, we find that decompaction/backstripping‐based methods may systematically underestimate sedimentation rates and total subsidence. The reason for this is that sediments are thinned during lithosphere extension in forward basin models while there are not in backstripping/decompaction approaches. The importance of sediment blanketing and differences in modelling approaches is illustrated in a reconstruction case study from the Norwegian Sea. The thermal and structural evolution of a transect across the Vøring Basin has been reconstructed using the backstripping/decompaction approach and TECMOD2D. Computed total subsidence curves differ by up to ～3 km and differences in computed basement heat flows reach up to 50%. These findings show that strong feedbacks exist between basin and lithosphere processes and that resolving them require integrated lithosphere‐scale basin models.     

Subsidence history from a backstripping analysis of the Permo‐Mesozoic succession of the Central Southern Alps (Northern Italy)

Seven tectonic subsidence curves, based on outcrop data, have been calculated in order to constrain the geodynamic evolution of the Permian–Mesozoic sedimentary succession (up to 10 km thick) of the Central Southern Alps basin (Italy). The analysis of the tectonic subsidence curves, covering a time span of about 200 Ma, allowed us to quantify the subsidence rates, to document the activity of syndepositional fault systems and calculate their slip rates. Different stages, in terms of duration and magnitude of subsidence‐uplift trends, have been identified in the evolution of the basin. The fault activity, reconstructed by comparing subsidence curves from adjacent sectors, resulted as highly variable both temporally and spatially. Strike‐slip tectonics was coeval to Permian sedimentation, as suggested by the strong differences in the subsidence rates in the sections. The evolution and subsidence rates suggest a continental shelf deposition from Early Triassic to Carnian, when subsidence came to a stop. A rapid resumption of subsidence is observed from the Norian, with a subsidence pulse in the Late Norian, followed by the regional uplift, in the Late Rhaetian. The following Early Jurassic subsidence is characterized by tectonic subsidence similar to that of the Norian. The Norian and Early Jurassic pulses were characterized by the highest slip rates along growth faults and are identified as two distinct tectonic events. The Norian–Rhaetian event is tentatively related to transtensional tectonics whereas the Early Jurassic event is related to crustal extension. The Early Jurassic subsidence records a shift in space an time of the beginning of the extensional stage, from Late Hettangian to the east to Late Pliensbachian–Toarcian to the west. From the Toarcian to the Aptian, the curves are compatible with regional thermal subsidence, later followed (Albian–Cenomanian) by uplift pulses in a retrobelt foreland basin (from Cenomanian onward).     

Location,extent, and magnitude of dynamic topography in the Late Cretaceous Cordilleran Foreland Basin,USA: New insights from 3D flexural backstripping

Mantle-induced dynamic topography (i.e., subsidence and uplift) has been increasingly recognized as an important process in foreland basin development. However, characterizing and distinguishing the effects (i.e., location, extent and magnitude) of dynamic topography in ancient foreland basins remains challenging because the spatio-temporal footprint of dynamic topography and flexural topography (i.e., generated by topographic loading) can overlap. This study employs 3D flexural backstripping of Upper Cretaceous strata in the central part of the North American Cordilleran foreland basin (CFB) to better quantify the effects of dynamic topography. The extensive stratigraphic database and good age control of the CFB permit the regional application of 3D flexural backstripping in this basin for the first time. Dynamic topography started to influence the development of the CFB during the late Turonian to middle Campanian (90.2–80.2 Ma) and became the dominant subsidence mechanism during the middle to late Campanian (80.2–74.6 Ma). The area influenced by >100 m dynamic subsidence is approximately 400 by 500 km, within which significant (>200 m) dynamic subsidence occurs in an irregular-shaped (i.e., lunate) subregion. The maximum magnitude of dynamic subsidence is 300 ± 100 m based on the 80.2–74.6 Ma tectonic subsidence maps. With the maximum magnitude of dynamic uplift being constrained to be 200–300 m, the gross amount of dynamic topography in the Late Cretaceous CFB is 500–600 m. Although the location of dynamic subsidence revealed by tectonic subsidence maps is generally consistent with isopach map trends, tectonic subsidence maps developed through 3D flexural backstripping provide more accurate constraints of the areal extent, magnitude and rate of dynamic topography (as well as flexural topography) in the CFB through the Late Cretaceous. This improved understanding of dynamic topography in the CFB is critical for refining current geodynamic models of foreland basins and understanding the surface expression of mantle processes.     

Late Cretaceous to Miocene sea-level estimates from the New Jersey and Delaware coastal plain coreholes: an error analysis

Sea level has been estimated for the last 108 million years through backstripping of corehole data from the New Jersey and Delaware Coastal Plains. Inherent errors due to this method of calculating sea level are discussed, including uncertainties in ages, depth of deposition and the model used for tectonic subsidence. Problems arising from the two-dimensional aspects of subsidence and response to sediment loads are also addressed. The rates and magnitudes of sea-level change are consistent with at least ephemeral ice sheets throughout the studied interval. Million-year sea-level cycles are, for the most part, consistent within the study area suggesting that they may be eustatic in origin. This conclusion is corroborated by correlation between sequence boundaries and unconformities in New Zealand. The resulting long-term curve suggests that sea level ranged from about 75–110 m in the Late Cretaceous, reached a maximum of about 150 m in the Early Eocene and fell to zero in the Miocene. The Late Cretaceous long-term (10⁷ years) magnitude is about 100–150 m less than sea level predicted from ocean volume. This discrepancy can be reconciled by assuming that dynamic topography in New Jersey was driven by North America overriding the subducted Farallon plate. However, geodynamic models of this effect do not resolve the problem in that they require Eocene sea level to be significantly higher in the New Jersey region than the global average.     

Calculating rates of syndepositional normal faulting in the western margin of the Mesozoic Subalpine Basin (south-east France)

A method is developed to quantify the rate of fault movement, with a very fine time-resolution, so that relevant histories of fault movements can be obtained. The study subject is a Triassic–Jurassic syndepositional normal fault located at the margin of an intracratonic deep basin, the Subalpine basin of south-eastern France. The fault has recently been identified and specifically investigated by a seismic survey along with drilling (Géologie profonde de la France Program). The investigation is based on correlation of time-lines on both sides of the structure through a period of about 70 Myr. Correlations have been made using variable approaches depending on the stratigraphic interval: recognition of laterally continuous marker-beds, biostratigraphic information and application of genetic stratigraphy concepts. In the case of biostratigraphic data, absolute ages are assigned to selected lines of correlation to determine time lengths and calculate velocities of fault movements. A specific backstripping procedure is established. The differential subsidence history between the two sites is restored not as a simple subtraction made after conventional backstripping on each site but as the sum of discrete differential subsidence increments calculated for each chronostratigraphic interval. The originality of the work lies in the completion of the supporting data base, implementation of high-resolution correlations within the large time-span of the investigation and development of a method to calculate the differential subsidence. Even though unassessable errors and uncertainties are still associated with the stratigraphic correlations, the backstripping procedure and the chronological control, the overall method offers a certain validity because the calculated and the observed differential subsidences are close. Beyond the obvious control on depositional geometries and localization of some reservoirs at the toe of the fault, the kinetic regime of the normal fault played an indirect but crucial part in the differential burial-related alteration of the reservoirs recorded on both sides of the fault. The high accuracy of the calculation has revealed that: (1) the growth pattern of the fault does not result from a continuous thermomechanical process but is composed of a series of rifting and sliding events related to gravity-driven extension; (2) the spectacular differential stratigraphic record on both sides of the fault is associated with fairly low values of the fault growth rate (maximum of 165 m Myr^?1). The method for measuring the growth of structures can be applied to any tectonic and sedimentary environment and offers a wide range of applications.     

Late Cretaceous basement foundering of the Rosario embayment, Peninsular Ranges forearc basin: backstripping of the El Gallo Formation, Baja California Norte, Mexico

Excellent exposure, well-controlled palaeobathymetry, and tightly-spaced, high-precision radiometric age control in the El Gallo Fm. permit rigorous quantitative analysis. Backstripping of these proximal nonmarine, forearc basin deposits reveals that, during the Late Cretaceous, the Rosario embayment of the Peninsular Ranges forearc was undergoing an episode of rapid tectonic subsidence. This subsidence had several marked effects on the sedimentology of the Rosario embayment: formation of a broad alluvial plain consisting of coarse-grained clastics; rapid (∼ 600 m Myr^-1) aggradation of sediments; and a retrogradational succession of facies, capped by a marine transgression, as deposition failed to keep pace with eustatic rise and subsidence.
Long-term sedimentation is driven by some combination of two allocyclic mechanisms: tectonic subsidence and eustatic sea-level rise. In order to evaluate which force predominated during deposition of the El Gallo Fm., the processes of sedimentation, compaction, and isostasy are evaluated through the interval in question. A sensitivity analysis is performed, in which the maximum tectonic and maximum eustatic contributions are estimated, along with the best-fit model. These results are qualitatively the same: tectonic subsidence was the major driving force of sedimentation in the Rosario embayment in late Campanian time. Regional sedimentological similarities suggest that this tectonic subsidence may have characterized the Peninsular Ranges forearc margin at this time, reflecting an episode of active down-faulting during the Late Cretaceous.     

Modelling rates and distribution of subsidence due to dynamic topography over subducting slabs: is it possible to identify dynamic topography from ancient strata?

Dynamic topography formed over subducting oceanic lithosphere has been proposed as a mechanism to explain certain otherwise anomalous long-wavelength patterns of subsidence inferred from ancient strata. Forward modelling of mantle flow in response to a subducting slab predicts amplitudes and distributions of dynamic topography that may occur due to various subducting slab geometries and histories. Plotting calculated dynamic topographies at a point against time produces tectonic subsidence curves. These subsidence curves show features such as evolution from convex to concave shape, amplitudes up to ~2000 m, subsidence rates up to ~220 m Myr⁻¹, and a general decrease in subsidence amplitude away from the subduction zone, over a distance of ~2000 km. On many convergent continental margins, dynamic topography is likely to be superimposed on other subsidence mechanisms. In back-arc basins, subsidence due to dynamic topography should be distinguishable from that due to extensional tectonics based simply on the temporal subsidence evolution expressed in the subsidence curve shapes. In a foreland basin setting, comparing dynamic topography models with forward models of flexural loading suggest the two processes can generate similar temporal subsidence patterns, but that dynamic topography causes subsidence over significantly greater wavelengths. Matches between calculated subsidence due to dynamic topography and backstripped subsidence patterns from Upper Cretaceous strata of the Western Interior Basin, USA, support the hypothesis that a long-wavelength 'background subsidence' was caused by dynamic topography.     

Gravity anomalies, subsidence history and the tectonic evolution of the Malay and Penyu Basins (offshore Peninsular Malaysia)

The tectonic subsidence and gravity anomalies in the Malay and Penyu Basins, offshore Peninsular Malaysia, were analysed to determine the isostatic compensation mechanism in order to investigate their origin. These continental extensional basins contain up to 14  km of sediment fill which implies that the crust had been thinned significantly during basin development. Our results suggest, however, that the tectonic subsidence in the basins cannot be explained simply by crustal thinning and Airy isostatic compensation.
The Malay and Penyu Basins are characterized by broad negative free-air gravity anomalies of between −20 and −30  mGal. To determine the cause of the anomaly, we modelled four gravity profiles across the basins using a method that combines two-dimensional flexural backstripping and gravity modelling techniques. We assumed a model of uniform lithospheric stretching and Airy isostasy in the analysis of tectonic subsidence. Our study shows that the basins are probably underlain by relatively thinned crust, indicating that some form of crustal stretching was involved. To explain the observed gravity anomalies, however, the Moho depth that we calculated based on the free-air gravity data is about 25% deeper than the Moho predicted by assuming Airy isostasy (Backstrip Moho). This suggests that the Airy model overestimates the compensation and that the basins are probably undercompensated isostatically. In other words, there is an extra amount of tectonic subsidence that is not compensated by crustal thinning, which has resulted in the discrepancy between the gravity-derived Moho and the Backstrip Moho. We attribute this uncompensated or anomalous tectonic subsidence to thin-skinned crustal extension that did not involve the mantle lithosphere. The Malay and Penyu Basins are interpreted therefore as basins that formed by a combination of whole-lithosphere stretching and thin-skinned crustal extension.     

洞庭湖演变趋势探讨

洞庭湖的演变主要受构造沉降、泥沙淤积和人类活动影响三大因素的控制。目前 ,洞庭湖盆的构造沉降速率虽然较低 (3～ 10mm/a) ,但构造沉降量仍抵消了部分泥沙淤积量 ,在一定程度上抑制了洞庭湖日益萎缩的趋势。由于湖盆中泥沙淤积速率大于构造沉降速率 ,洲滩会继续发育和扩展 ,洞庭湖仍保持淤高的趋势。在“4 35 0工程”完成和三峡建坝后 80年内 ,湖盆泥沙淤积速率将降低到 1 79mm/a ,洞庭湖不断淤高的趋势将得到减缓。     

Natural subsidence of the Venice area during the last 60 Myr

Integrated geohistory analysis performed on high‐resolution stratigraphy of Venezia 1 and Lido 1 wells (Quaternary–Pliocene interval) and low‐resolution stratigraphy of a simulated well extending Lido 1 down to the base of Cenozoic (Palaeocene–Miocene interval) is used to reconstruct the interplay between subsidence and sedimentation that occurred in the Venice area (eastern Po Plain) during the last 60 Myr, and to discuss the relationships between calculated subsidence rates and time resolution of stratigraphic data. Both subsidence and sedimentation are mostly related to the tectonic evolution of the belts that surround the Venice basin, influencing the lithosphere vertical motions and the input of clastic sediments through time. In particular, two subsidence phases are recorded between 40–33.5 and 32.5–24 Myr (0.13 and 0.14 mm year^?1, respectively), coeval with tectonic phases in the Dinaric belt. Vice versa, during the main South‐Alpine orogenic phase (middle–late Miocene), quiescence or little uplift (?0.03 mm year^?1) reflects the location of the Venice area close to the peripheral bulge of the South‐Alpine foreland system. Early Pliocene evolution is characterised by a number of subsidence/uplift events, among which two uplifts occurred between 5–4.5 and 3–2.2 Myr (at ?0.4 and ?0.2 mm year^?1, respectively) and can be correlated with tectonic motions in the Apennines. During the last million years, the Venice area was initially characterised by uplift (?0.6 mm year^?1 rising to ?1.5 mm year^?1 between 0.4 and 0.38 Myr), eventually replaced by subsidence at a rate ranging between 1.6 and 1.0 mm year^?1 up to 0.12 Myr and then decreased to 0.4 mm year^?1, as an average, up to present. Our results highlight that time resolution of the stratigraphic dataset deeply influences the order of magnitude obtained for the calculated subsidence rate. This is because subsidence seems to have worked through short‐lived peaks (in the order of 10⁵ years), alternating with long relatively quiescent intervals. This suggests caution when components of subsidence are deduced by subtracting long‐term to short‐term subsidence rate.     

Subsidence analyses from the Betic Cordillera, southeast Spain

Fifty‐four Mesozoic–Cenozoic stratigraphic sections from the Betic Cordillera of southeast Spain have been analysed in order to estimate the timing and amount of lithospheric stretching that occurred at the western end of the Tethyan Ocean since the Hercynian Orogeny. The standard backstripping technique has been used in order to calculate the water‐loaded subsidence of basement for each section. Water‐loaded subsidence curves were then inverted in order to determine the variation of lithospheric strain rate as a function of time, which yields estimates of timing, magnitude and intensity of stretching. Rifting commenced during the Late Permian/Early Triassic times and continued intermittently throughout the Mesozoic in response to the opening of the Tethyan Ocean to the east and the opening of the Atlantic Ocean to the west. Two major events in the Permo‐Triassic/Early Jurassic and the Late Jurassic/Early Cretaceous can be clearly identified. Stretching factors are generally small (1.1–1.25) probably because the Betic Cordillera was located at the westernmost end of the Tethys. Peak strain rates of ～10^?15 s^?1 were obtained for Mesozoic rift events and these values are in broad agreement with those obtained throughout the Tethyan Realm. We have also analysed the Neogene extensional event, which played an important role in forming the existing Mediterranean Sea. A combination of well‐log information and calibrated seismic reflection data was modelled. Peak strain rates in these younger basins are almost one order of magnitude faster than those estimated for the Mesozoic basins. These higher values appear to be typical of back‐arc extensional basins elsewhere. To the west and north of the Betic Cordillera, the Guadalquivir foreland basin developed as extension took place further east. Backstripped sections from this basin clearly record the northward migration of foreland basin subsidence through time.     

The large-scale dynamics of grain-size variation in alluvial basins, 2: Application to syntectonic conglomerate

The concept of‘syntectonic’ conglomerate is based on the idea that gravel progradation is mainly generated by an increase in tectonic uplift and erosion of a source area with attendant increase in sediment flux supplied to a basin. However, other mechanisms, such as changes in basin subsidence rates, sorting of supplied sediment, and capability of transporting streams, can also lead to progradation and be difficult to distinguish from a syntectonic origin. Here we use our previously developed model to help understand the origin of gravel progradation in three Neogene alluvial basins - the Bermejo Basin of Argentina, the Himalayan Foreland Basin, and the San Pedro Basin of southern Arizona - all of which have available high-resolution magnetostratigraphy. Interpretation of the origin of gravel progradation in these basins begins with calculation of basin equilibrium time, which is the time-scale required for the streams to reach a steady-state profile, assuming constant conditions. We then compare the time-scale of the observed changes in the basin with the equilibrium time to determine if and how the model can be applied to the stratigraphic record. Most of the changes we have studied occur on time scales longer than the equilibrium time (‘slow variations’), in which case the key to interpretation is the relationship between overall grain-size change and sedimentation rate in vertical sections. Of the three examples studied only one, the Bermejo Basin, is consistent with the traditional model of syntectonic progradation. Overall progradation in the two other basins is most consistent with a long-term reduction in basin subsidence rates. In addition, short-term variation in diffusivity or sediment flux, probably climatically driven, is the most likely control of small-scale progradation of gravel tongues in the San Pedro Basin. These results, along with observations from other basins, suggest that subsidence is clearly an important control on clastic progradation on ‘slow’ time scales (i.e. generally a million years or more). If subsidence rates are directly linked to tectonic events, then subsidence-driven progradation marks times of tectonic quiescence and is clearly not syntectonic in the traditional sense. These examples show that the model can be useful in interpreting the rock record, particularly when combined with other traditional basin-analysis techniques. In particular, our results can be used to help discriminate between clastic progradation due to tectonic origin and progradation resulting from other mechanisms in alluvial basins.     

Tectonic evolution of the Alboran Sea basin

The Alboran Sea is an extensional basin of Neogene age that is surrounded by highly arcuate thrust belts. Multichannel seismic (MCS) reflection profile data suggest the basin has a complex tectonic fabric that includes extensional, compressional and strike-slip structures. The early Miocene history appears to be dominated by graben formation with border faults that are in large part contemporaneous with thrust movements in the external zones of the Betic and Rif mountains. Extension appears to have continued into the late Miocene although the main movements were probably completed by the time of the Messinian ‘salinity crisis’. The Pliocene and younger history of the basin is dominated by infilling of the Messinian topography, gentle subsidence, and extensional, compressional and strike-slip movements. There is evidence from the sea-floor morphology and seismicity patterns that the basin is actively deforming in response to present-day plate motions. Backstripping of well data in the basin margin suggests that the initial extensional event was accompanied by crustal and lithospheric thinning. The depth to Moho inferred from backstripping is greater than the depth expected based on seismic and gravity modelling, suggesting that backstripping underestimates the true amount of thinning. One explanation is that some of the thinning occurred while the crust was above sea level, perhaps as a result of either crustal thickening, or a period of lithospheric heating and thinning, prior to rifting. We found that a model with a ‘normal’ crustal thickness of 31.2 km, a lithospheric thickness of 50 km, and β= 1.4 predicts 0.8 km of initial uplift. These parameters fit the well subsidence data and bring the backstripped Moho into better agreement with the seismic and gravity Moho. The origin of such a thin lithosphere is not constrained by the data, but we believe that it may be a result of the detachment of a cold lithospheric ‘root’ that formed during pre-Neogene collisional orogeny in the region.     

Simple equations of sedimentation: applications to sequence stratigraphy

Abstract An equation to relate the thickness of sediment deposited (ΔSed), eustatic sea-level change (ΔE), and subsidence (ΔSub), to changes in depth of water (ΔD) is: ΔSub +ΔE-ΔSed =ΔD.
Using existing sea-level curves, the equation shows that some transgressive-regressive sequences in a foreland basin and a composite seismic facies sequence on a passive margin cannot result solely from eustatic variation. In each case, the space created by subsidence is greater than that provided by eustatic rise. However, eustatic variation could have triggered sequence development if superimposed on a basin with relatively constant values of the other parameters. Short-period sea-level fluctuations with high rates of change, exceeding 70–100 m Myr^-1 for periods less than 2–3 Myr, affect the stratigraphy and sedimentology more than longer period, higher amplitude variations.
Clinoforms are generated because of lateral variations in sedimentation rate compared to the rate of creation of accommodation space. These variations may result from differing sedimentation rates, subsidence rates, or rates of eustatic change, superimposed on a basin with lateral sediment supply. Clinoform slopes and curvatures are interpre table in terms of these variables as well as the type of sediment supplied and the energy distribution in the basin.
These equations put some well-known geological principles on a simple quantitative basis. They force precision in definition of variables, and may lead to further development of quantitative techniques in stratigraphy and sedimentology.     

Stratigraphic and structural expression of the lateral growth of thrust fault-propagation folds: results and implications from kinematic modelling

In order to better understand the development of thrust fault‐related folds, a 3D forward numerical model has been developed to investigate the effects that lateral slip distribution and propagation rate have on the fold geometry of pre‐ and syn‐tectonic strata. We consider a fault‐propagation fold in which the fault propagates upwards from a basal decollement and along‐strike normal to transport direction. Over a 1 Ma runtime, the fault reaches a maximum length of 10 km and accumulates a maximum displacement of 1 km. Deformation ahead of the propagating fault tip is modelled using trishear kinematics while backlimb deformation is modelled using kink‐band migration. The applicability of two different lateral slip distributions, namely linear‐taper and block‐taper, are firstly tested using a constant lateral propagation rate. A block‐taper slip distribution replicates the geometry of natural fold‐thrusts better and is then used to test the sensitivity of thrust‐fold morphology to varied propagation rates in a set of fault‐propagation folds that have identical final displacement to length (D_max/L_max) ratios. Two stratigraphic settings are considered: a model in which background sedimentation rates are high and no topography develops, and a model in which a topographic high develops above the growing fold and local erosion, transport and deposition occur. If the lateral propagation rate is rapid (or geologically instantaneous), the fault tips quickly become pinned as the fault reaches its maximum lateral extent (10 km), after which displacement accumulates. In both stratigraphic settings, this leads to strike‐parallel rotation of the syn‐tectonic strata near the fault tips; high sedimentation rates relative to rates of uplift result in along‐strike thinning over the structural high, while low sedimentation rates result in pinchout against it. In contrast, slower lateral propagation rates (i.e. up to one order of magnitude greater than slip rate) lead to the development of along‐strike growth triangles when sedimentation rates are high, whereas when sedimentation rates are low, offflap geometries result. Overall we find that the most rapid lateral propagation rates produce the most realistic geometries. In both settings, time‐equivalent units display both nongrowth and growth stratal geometries along‐strike and the transition from growth to nongrowth has the potential to delineate the time of fault/fold growth at a given location. This work highlights the importance of lateral fault‐propagation and fault tip pinning on fault and fold growth in three dimensions and the complex syn‐tectonic geometries that can result.     

Late orogenic intramontane basin development: the Granada basin, Betics (southern Spain)

The quantitative study of subsidence in the Granada basin, using decompaction and backstripping techniques, and contemporaneous relief development in the surrounding areas, especially in the Sierra Nevada, provides a good case example of the development of an intramontane basin. In the Granada basin, according to the interpretation of the seismic profiles and results of the backstripping analysis, subsidence and sedimentation rates were at a maximum in the late Tortonian and decreased progressively; meanwhile, the neighbouring areas were uplifted forming important relief. Chronostratigraphical revisions of the marine sediments show that the marine incursion that deposited sediments in the Granada basin lasted only 1.3 Ma, between 8.5 and 7.2 Ma. The gradual retreat of the sea in the Granada basin is not attributable to global eustatic fluctuations, but rather to uplift in the Sierra Nevada and its adjacent areas. From latest Tortonian to early Messinian times, the region became continental and the Granada basin acquired its present physiography and was differentiated as such. From the late Tortonian onwards, NNW–SSE compression combined with ENE–WSW extension affected the cordillera. In the Granada basin, extension controlled fault movements. There are two well-defined fault sets: the first trends 70°N–90°E, with low angle faults (less than 30°) dipping towards the north and south, defining the subsiding areas which have approximately E–W direction; whereas the second set has a NW–SE direction, and cuts and displaces the previous ones, defining the main subsiding areas in the eastern part of the basin. The reinterpretation of seismic profiles reveals that the subsiding axes within the Granada basin persisted from the Tortonian to the present because of continued displacements of the main faults.     

Joint inversion of temperature,vitrinite reflectance and fission tracks in apatite with examples from the eastern North Sea area

As sediment accumulation indicates basin subsidence, erosion often is understood as tectonic uplift, but the amplitude and timing may be difficult to determine because the sedimentary record is missing. Quantification of erosion therefore requires indirect evidence, for example thermal indicators such as temperature, vitrinite reflectance and fission tracks in apatite. However, as always, the types and quality of data and the choice of models are important to the results. For example, considering only the thermal evolution of the sedimentary section discards the thermal time constant of the lithosphere and essentially ignores the temporal continuity of the thermal structure. Furthermore, the types and density of thermal indicators determine the solution space of deposition and erosion, the quantification of which calls for the use of inverse methods, which can only be successful when all models are mutually consistent. Here, we use integrated basin modelling and Markov Chain Monte Carlo inversion of four deep boreholes to show that the erosional pattern along the Sorgenfrei–Tornquist Zone (STZ) in the eastern North Sea is consistent with a tectonic model of tectonic inversion based on compression and relaxation of an elastic plate. Three wells in close proximity SW of the STZ have different data and exhibit characteristic differences in erosion estimates but are consistent with the formation of a thick chalk sequence, followed by minor Cenozoic erosion during relaxation inversion. The well on the inversion ridge requires ca. 1.7 km Jurassic-Early Cretaceous sedimentation followed by Late Cretaceous–Palaeocene erosion during inversion. No well demands thick Cenozoic sedimentation followed by equivalent significant Neogene exhumation. When data are of high quality and models are consistent, the thermal indicator method yields significant results with important tectonic and geodynamic implications.     

Tectonic and climatic controls on the sequential arrangement of an alluvial fan/fan‐delta complex (Montserrat,Eocene, Ebro Basin,NE Spain)

A magnetostratigraphy‐based chronological framework has been constructed in the Eocene sediments of the Montserrat alluvial fan/fan‐delta complex (southeast Ebro Basin), in order to unravel forcing controls on their sequential arrangement and to revise the tectonosedimentary history of the region. The palaeomagnetic study is based on 403 sites distributed along an 1880‐m‐thick composite section, and provides improved temporal constraints based on an independent correlation to the geomagnetic polarity time scale. The new chronological framework together with sequence stratigraphy and geohistory analysis allow us to investigate the interplay between factors controlling the sequential arrangement of the Montserrat complex at the different temporal scales and to test for orbitally driven climate forcing. The results suggest that the internal stacking pattern in transgressive and regressive sequences sets within the more than 1000‐m‐thick Milany Composite Megasequence can be explained as the result of subsidence‐driven accommodation changes under a general increase of sediment supply. Composite sequences (tens to hundreds of metres thick) likely reflect orbitally forced cyclicity related to the 400‐kyr eccentricity cycle, possibly controlled by climatically induced sea‐level fluctuations. This study also provides new insights on the deformational history of the area, and shows a correlation between (tectonic) subsidence and forelimb rotation measured on basin‐margin deformed strata. Integration of subsidence curves from different sectors of the eastern Ebro Basin allows us to estimate the variable contribution of tectonic loads from the two active basin margins: the Catalan Coastal Ranges and the Pyrenees. The results support the presence of a double flexure from Late Lutetian to Late Bartonian, associated with the two tectonically active margins. From Late Bartonian to Early Priabonian the homogenization of subsidence values is interpreted as the result of the coupling of the two sources of tectonic load.     

Evolution of the intracratonic Officer Basin, central Australia: implications from subsidence analysis and gravity modelling

ABSTRACT The intracratonic basins of central Australia are distinguished by their large negative Bouguer gravity anomalies, despite the absence of any significant topography. Over the Neoproterozoic to Palaeozoic Officer Basin, the anomalies attain a peak negative amplitude in excess of 150 mGal, amongst the largest of continental anomalies observed on Earth. Using well data from the Officer and Amadeus basins and a data grid of sedimentary thicknesses from the eastern Officer Basin, we have assessed the evolution of these intracratonic basins. One-dimensional backstripping analysis reveals that Officer and Amadeus basin tectonic subsidence was not entirely synchronous. This implies that the basins evolved as discrete geological features once the Centralian Superbasin was dismembered into its constituent basins. Two- and three-dimensional backstripping and gravity modelling suggest that the eastern Officer Basin evolved from a broad continental sag into a region of intracratonic flexural subsidence from the latest Neoproterozoic, when flexure of the lithosphere deepened the northern basin. The results from gravity modelling improve when the crust is thickened beneath the northern margin of the basin and thinned at the southern margin, as has been suggested by recent deep seismic data. The crustal thickening beneath the basin's northern margin abuts the region of greatest topographic relief and is consistent with the observed structure at the edges of many orogenic belts. If the Officer Basin evolved as a foreland-type basin from the late Proterozoic and has retained those features to the present, then one implication is that in the absence of any significant topography, cratonic lithosphere must be able to support stresses over very long periods of geological time.