Strike-slip tectonics and development of the Tertiary Queen Charlotte Basin, offshore western Canada: evidence from seismic reflection data

Interpretation of seismic reflection data have led to a new model of the development of the Queen Charlotte Basin. New multi-channel data collected in 1988 and an extensive network of unpublished older single- and multi-channel profiles from industry image a complex network of sub-basins. Structural styles vary along the axis of the basin from broadly spaced mainly N-trending sub-basins in Queen Charlotte Sound, to closely spaced NW-trending sub-basins in Hecate Strait, to an E-W en echelon belt of sub-basins in Dixon Entrance. Transtensional tectonics dominated in the Miocene and transpression dominated in the Pliocene except in Queen Charlotte Sound. The data we present prove that the origin of the basin is extensional and its most recent deformation is compressive. Evidence for the strike-slip origin of tectonism includes along-axis variations in structures, simultaneous extension and compression in adjacent sub-basins, lack of correlations across faults, and mixed normal and reverse faults within structures. We infer that the Pacific-North America plate boundary has been west of the Queen Charlotte Islands since the Miocene when relative plate motions have been dominantly strike-slip. The formation and development of the Queen Charlotte Basin is the result of distributed shear; by which a small percentage of the plate motion has been taken up in a network of faults across the continental margin. As this region of crust deforms it interacts with neighbouring rigid crust resulting in extension dominating in the south of the basin and compression in the north. Continental crust adjacent to some transform plate boundaries can be sheared over a wide region; the network of basins in southwestern California is a good analogue for the Queen Charlotte Basin.     

Late Jurassic subsidence and passive margin evolution in the Vulcan Sub-basin, north-west Australia: constraints from basin modelling

The Vulcan Sub-basin, located in the Timor Sea, north-west Australia, developed during the Late Jurassic extension which ultimately led to Gondwanan plate breakup and the development of the present-day passive continental margin. This paper describes the evolution of upper crustal extension and the development of Late Jurassic depocentres in this subbasin, via the use of forward modelling techniques. The results suggest that a lateral variation in structural style exists. The south of the basin is characterized by relatively large, discrete normal faults which have generated deep sub-basins, whereas more distributed, small-scale faulting further north reflects a collapse of the early basin margin, with the development of a broader, 'sagged' basin geometry. By combining forward and reverse modelling techniques, the degree of associated lithosphere stretching can be quantified. Upper crustal faulting, which represents up to 10% extension, is not balanced by extension in the deeper, ductile lithosphere; the magnitude of this deeper extension is evidenced by the amount of post-Valanginian thermal subsidence. Reverse modelling shows that the lithosphere stretching
factor has a magnitude of up to β=1.55 in the southern Vulcan Sub-basin, decreasing to β=1.2 in the northern Vulcan Sub-basin. It is proposed that during plate breakup, deformation in the Vulcan Sub-basin consisted of depth-dependent lithosphere extension. This additional component of lower crustal and lithosphere stretching is considered to reflect long-wavelength partitioning of strain associated with continental breakup, which may have extended 300–500 km landward of the continent–ocean boundary.     

Regional stratigraphy and subsidence of Orphan Basin near the time of breakup and implications for rifting processes

The stratigraphic, subsidence and structural history of Orphan Basin, offshore the island of Newfoundland, Canada, is described from well data and tied to a regional seismic grid. This large (400 by 400 km) rifted basin is part of the non‐volcanic rifted margin in the northwest Atlantic Ocean, which had a long and complex rift history spanning Middle Jurassic to Aptian time. The basin is underlain by variably thinned continental crust, locally <10‐km thick. Our work highlights the complex structure, with major upper crustal faults terminating in the mid‐crust, while lower crustal reflectivity suggests ductile flow, perhaps accommodating depth‐dependent extension. We describe three major stratigraphic horizons connected to breakup and the early post‐rift. An Aptian–Albian unconformity appears to mark the end of crustal rifting in the basin, and a second, more subdued Santonian unconformity was also noted atop basement highs and along the proximal margins of the basin. Only minor thermal subsidence occurred between development of these two horizons. The main phase of post‐rift subsidence was delayed until post‐Santonian time, with rapid subsidence culminating in the development of a major flooding surface in base Tertiary time. Conventional models of rifting events predict significant basin thermal subsidence immediately following continental lithospheric breakup. In the Orphan Basin, however, this subsidence was delayed for about 25–30 Myr and requires more thinning of the mantle lithosphere than the crust. Models of the subsidence history suggest that extreme thinning of the lithospheric mantle continued well into the post‐rift period. This is consistent with edge‐driven, small‐scale convective flow in the mantle, which may thin the lithosphere from below. A hot spot may also have been present below the region in Aptian–Albian time.     

Tectonic evolution of the North Sea basin: crustal stretching and subsidence

Summary. The lithospheric stretching model for the formation of sedimentary basins was tested in the central North Sea by a combined study of crustal thinning and basement subsidence patterns. A profile of crustal structure was obtained by shooting a long-range seismic experiment across the Central Graben, the main axis of subsidence. A seabed array of 12 seismometers in the graben was used to record shots fired in a line 530 km long across the basin. The data collected during the experiment were interpreted by modelling synthetic seismograms from a laterally varying structure, and the final model showed substantial crustal thinning beneath the graben. Subsidence data from 19 exploration wells were analysed to obtain subsidence patterns in the central North Sea since Jurassic times. Changes in water depth were quantified using foraminiferal assemblages where possible, and observed basement subsidence paths were corrected for sediment loading, compaction and changes in water depth through time. The seismic model is shown to be compatible with the observed gravity field, and the small size of observed gravity anomalies is used to argue that the basin is in local isostatic equilibrium. Both crustal thinning and basement subsidence studies indicate about 70 km of stretching across the Central Graben during the mid-Jurassic to early Cretaceous extensional event. This extension appears to have occurred over crust already slightly thinned beneath the graben, and the seismic data suggest that total extension since the early Permian may have been more than 100km. The data presented here may all be explained using a simple model of uniform extension of the lithosphere.     

The long-term thermal consequences of rifting: implications for basin reactivation

The attenuation of the continental crust during rifting and the subsequent filling of the rift‐related accommodation alter the long‐term thermal and mechanical state of the lithosphere. This is primarily because the Moho is shallowed due to density contrasts between the sediment fill and the crust, but also reflects the attenuation of the pre‐existing crustal heat production and its burial beneath the basin, as well the thermal properties of the basin fill. Moho shallowing and attenuation of pre‐existing heat production contribute to long‐term cooling of the Moho and thus lithospheric strengthening, as has been pointed out in many previous studies. In contrast, basin filling normally contributes to significant Moho heating allowing the possibility of long‐term lithospheric weakening, the magnitude of which is dependent on the thermal properties of the basin‐fill and the distribution of heat sources in the crust. This paper focuses on the thermal property structure of the crust and basin‐fill in effecting long‐term changes in lithospheric thermal regime, with particular emphasis on the distribution of heat producing elements in the crust. The parameter space appropriate to typical continental crust is explored using a formalism for the heat production distributions that makes no priori assumptions about the specific form of the distribution. The plausible parameter space allows a wide range in potential long‐term thermal responses. However, with the proviso that the accommodation created by the isostatic response to rifting is essentially filled, the long‐term thermal response to rift basin formation will generally increase average crustal thermal gradients beneath basins but cool the Moho due to its reduction in depth. The increase in the average crustal thermal gradient induces lateral heat flow that necessarily heats the Moho along basin margins, especially in narrow rift basins. Using coupled thermo‐mechanical models with temperature sensitive creep‐parameters, we show that such heating may be sufficient to localise subsequent deformation in the vicinity of major basin bounding structures, potentially explaining the offset observed in some stacked rift basin successions.     

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.     

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.     

Causes of spatially variable tectonic subsidence in the Miocene Bermejo Foreland Basin, Argentina

ABSTRACT Tectonic subsidence in the 20–9 Ma Bermejo basin resulted from spatially variable crustal loading on a lithosphere of spatially variable strength (e.g. elastic thickness). Reconstruction of the crustal loads added between 20 and 9 Ma, and assessment of the effects of these loads on an elastic, isotropic lithosphere confirm this hypothesis. Elastic models effectively explain tectonic subsidence east of the Iglesia–Calingasta basin, but west of it crustal loads were locally compensated. Elastic models also prove that the 20–9 Ma Frontal Cordillera loading is of no importance in the mechanical system of the Bermejo basin. 2D and 3D elastic models of a uniformly strong lithosphere under 20–9 Ma crustal loads corrected for post‐9 Ma erosion successfully replicate the 9 Ma Bermejo basin's proximal palaeotopography. However, they fail to replicate the 9 Ma basin's medial and distal palaeotopography. A 3D finite element model of a lithosphere with bimodal strength (weak below the Bermejo basin and west of the Precordillera, and strong below the Precordillera and east of the Valle Fértil lineament) successfully replicates the 9 Ma basin's palaeotopography. That variable strength model introduces a southward decrease in the wavelength of flexural deformation, which results in a basin that narrows southward, consistent with the 9 Ma Bermejo basin. The preferred 9 Ma lithospheric strength distribution is similar to the present lithospheric strength field derived from gravity data, suggesting that the bimodal strength signature was retained throughout the entire basin's evolution. Late Miocene flattening of the subducting slab, tectonic change to a broken foreland, or deposition of a thick (～8–10 km) sedimentary cover did not affect the strength of the lithosphere underlying the Bermejo basin. The long‐term bimodal strength field does not correlate with the documented thickness of the seismogenic crust.     

Evidence for underthrusting beneath the Queen Charlotte Margin, British Columbia, from teleseismic receiver function analysis

The Queen Charlotte Fault zone is the transpressive boundary between the North America and Pacific Plates along the northwestern margin of British Columbia. Two models have been suggested for the accommodation of the ∼20 mm yr⁻¹ of convergence along the fault boundary: (1) underthrusting; (2) internal crustal deformation. Strong evidence supporting an underthrusting model is provided by a detailed teleseismic receiver function analysis that defines the underthrusting slab. Forward and inverse modelling techniques were applied to receiver function data calculated at two permanent and four temporary seismic stations within the Queen Charlotte Islands. The modelling reveals a ∼10 km thick low-velocity zone dipping eastward at 28° interpreted to be underthrusting oceanic crust. The oceanic crust is located beneath a thin (28 km) eastward thickening (10°) continental crust.     

Crustal and upper-mantle structure in the Eastern Mediterranean from the analysis of surface wave dispersion curves

The dispersive properties of surface waves are used to infer earth structure in the Eastern Mediterranean region. Using group velocity maps for Rayleigh and Love waves from 7 to 100 s, we invert for the best 1-D crust and upper-mantle structure at a regular series of points. Assembling the results produces a 3-D lithospheric model, along with corresponding maps of sediment and crustal thickness. A comparison of our results to other studies finds the uncertainties of the Moho estimates to be about 5 km. We find thick sediments beneath most of the Eastern Mediterranean basin, in the Hellenic subduction zone and the Cyprus arc. The Ionian Sea is more characteristic of oceanic crust than the rest of the Eastern Mediterranean region as demonstrated, in particular, by the crustal thickness. We also find significant crustal thinning in the Aegean Sea portion of the backarc, particularly towards the south. Notably slower S -wave velocities are found in the upper mantle, especially in the northern Red Sea and Dead Sea Rift, central Turkey, and along the subduction zone. The low velocities in the upper mantle that span from North Africa to Crete, in the Libyan Sea, might be an indication of serpentinized mantle from the subducting African lithosphere. We also find evidence of a strong reverse correlation between sediment and crustal thickness which, while previously demonstrated for extensional regions, also seems applicable for this convergence zone.     

Subsidence in the super-deep Pattani and Malay basins of Southeast Asia: a coupled model incorporating lower-crustal flow in response to post-rift sediment loading

Two Early Cenozoic rifts in Southeast Asia (beneath the Pattani and Malay basins) experienced only limited upper-crustal extension (β≤1.5); yet very thick post-rift sequences are present, with 6–12 km of Late Cenozoic terrestrial and shallow-marine sediment derived from adjacent sources. Conventional post-rift backstripping requires depth-dependent lithospheric thinning by β=2–4 to explain these tremendous thicknesses. We assess an alternative explanation for this post-rift subsidence, involving lower-crustal flow from beneath these basins in response to lateral pressure-gradients induced by the sediment loads and the negative loads arising from the erosion of their sediment sources. We calculate that increased rates of erosion in western Thailand in the Early Miocene placed the crust in a non-steady thermal state, such that the depth (and thus, the pressure) at the base of the brittle upper crust subsequently varied over time. Following such a perturbation, thermal and mass-flux steady-state conditions took millions of years to re-establish. In the meantime, the lateral pressure-gradient caused net outflow of lower crust, thinning the crust beneath the depocentre by several kilometres (mimicking the isostatic effect of greater crustal extension having occurred beforehand) and thickening it beneath the sediment source region. The local combination of hot crust and high rates of surface processes, causing lower-crustal flow to be particularly vigorous and thus making its effects more readily identifiable, means that the Pattani and Malay basins represent a set of conditions different from basins in many other regions. However, lower-crustal flow induced by surface processes will also occur to some extent, but less recognisably, in many other continental crustal provinces, but its effects may be mistaken for those of other processes, such as larger-magnitude stretching and/or depth-dependent stretching.     

Crust‐scale 3D model of the Western Bredasdorp Basin (Southern South Africa): data‐based insights from combined isostatic and 3D gravity modelling

The southern South African continental margin documents a complex margin system that has undergone both continental rifting and transform processes in a manner that its present‐day architecture and geodynamic evolution can only be better understood through the application of a multidisciplinary and multi‐scale geo‐modelling procedure. In this study, we focus on the proximal section of the larger Bredasdorp sub‐basin (the westernmost of the five southern South African offshore Mesozoic sub‐basins), which is hereto referred as the Western Bredasdorp Basin. Integration of 1200 km of 2D seismic‐reflection profiles, well‐logs and cores yields a consistent 3D structural model of the Upper Jurassic‐Cenozoic sedimentary megasequence comprising six stratigraphic layers that represent the syn‐rift to post‐rift successions with geometric information and lithology‐depth‐dependent properties (porosities and densities). We subsequently applied a combined approach based on Airy's isostatic concept and 3D gravity modelling to predict the depth to the crust‐mantle boundary (Moho) as well as the density structure of the deep crust. The best‐fit 3D model with the measured gravity field is only achievable by considering a heterogeneous deep crustal domain, consisting of an uppermost less dense prerift meta‐sedimentary layer [ρ = 2600 kg m^?3] with a series of structural domains. To reproduce the observed density variations for the Upper Cenomanian–Cenozoic sequence, our model predicts a cumulative eroded thickness of ca. 800–1200 m of Tertiary sediments, which may be related to the Late Miocene margin uplift. Analyses of the key features of the first crust‐scale 3D model of the basin, ranging from thickness distribution pattern, Moho shallowing trend, sub‐crustal thinning to shallow and deep crustal extensional regimes, suggest that basin initiation is typical of a mantle involvement deep‐seated pull‐apart setting that is associated with the development of the Agulhas‐Falkland dextral shear zone, and that the system is not in isostatic equilibrium at present day due to a mass excess in the eastern domain of the basin that may be linked to a compensating rise of the asthenospheric mantle during crustal extension. Further corroborating the strike‐slip setting is the variations of sedimentation rates through time. The estimated syn‐rift sedimentation rates are three to four times higher than the post‐rift sedimentation, thereby indicating that a rather fast and short‐lived subsidence during the syn‐rift phase is succeeded by a significantly poor passive margin development in the post‐rift phase. Moreover, the derived lithospheric stretching factors [β = 1.5–1.75] for the main basin axis do not conform to the weak post‐rift subsidence. This therefore suggests that a differential thinning of the crust and the mantle‐lithosphere typical for strike‐slip basins, rather than the classical uniform stretching model, may be applicable to the Western Bredasdorp Basin.     

Timing of depth‐dependent lithosphere stretching on the S. Lofoten rifted margin offshore mid‐Norway: pre‐breakup or post‐breakup?

Depth‐dependent stretching, in which whole‐crustal and whole‐lithosphere extension is significantly greater than upper‐crustal extension, has been observed at both non‐volcanic and volcanic rifted continental margins. A key question is whether depth‐dependent stretching occurs during pre‐breakup rifting or during sea‐floor spreading initiation and early sea‐floor spreading. Analysis of post‐breakup thermal subsidence and upper‐crustal faulting show that depth‐dependent lithosphere stretching occurs on the outer part of the Norwegian volcanic rifted margin. For the southern Lofoten margin, large breakup lithosphere β stretching factors approaching infinity are required within 100 km of the continent–ocean boundary to restore Lower Eocene sediments and flood basalt surfaces (～54 Ma) to interpreted sub‐aerial depositional environments at sea level as indicated by well data. For the same region, the upper crust shows no significant Palaeocene and Late Cretaceous faulting preceding breakup with upper‐crustal β stretching factors <1.05. Further north on the Lofoten margin, reverse modelling of post‐breakup subsidence with a β stretching factor of infinity predicts palaeo‐bathymetries of ～1500 m to the west of the Utrøst Ridge and fails to restore Lower Eocene sediments and flood basalt tops to sea level at ～54 Ma. If these horizons were deposited in a sub‐aerial depositional environment, as indicated by well data to the south, an additional subsidence event younger than 54 Ma is required compatible with lower‐crustal thinning during sea‐floor spreading initiation. For the northern Vøring margin, breakup lithosphere β stretching factors of ～2.5 are required to restore Lower Eocene sediments and basalts to sea level at deposition, while Palaeocene and Late Cretaceous upper‐crustal β stretching factors for the same region are < 1.1. The absence of significant Palaeocene and late Cretaceous extension on the southern Lofoten and northern Vøring margins prior to continental breakup supports the hypothesis that depth‐dependent stretching of rifted margin lithosphere occurs during sea‐floor spreading initiation or early sea‐floor spreading rather than during pre‐breakup rifting.     

The formation of a failed continental breakup basin: The Cenozoic development of the Faroe‐Shetland Basin

Ultra‐large rift basins, which may represent palaeo‐propagating rift tips ahead of continental rupture, provide an opportunity to study the processes that cause continental lithosphere thinning and rupture at an intermediate stage. One such rift basin is the Faroe‐Shetland Basin (FSB) on the north‐east Atlantic margin. To determine the mode and timing of thinning of the FSB, we have quantified apparent upper crustal β‐factors (stretching factors) from fault heaves and apparent whole‐lithosphere β‐factors by flexural backstripping and decompaction. These observations are compared with models of rift basin formation to determine the mode and timing of thinning of the FSB. We find that the Late Jurassic to Late Palaeocene (pre‐Atlantic) history of the FSB can be explained by a Jurassic to Cretaceous depth‐uniform lithosphere thinning event with a β‐factor of ~1.3 followed by a Late Palaeocene transient regional uplift of 450–550 m. However, post‐Palaeocene subsidence in the FSB of more than 1.9 km indicates that a Palaeocene rift with a β‐factor of more than 1.4 occurred, but there is only minor Palaeocene or post‐Palaeocene faulting (upper crustal β‐factors of less than 1.1). The subsidence is too localized within the FSB to be caused by a regional mantle anomaly. To resolve the β‐factor discrepancy, we propose that the lithospheric mantle and lower crust experienced a greater degree of thinning than the upper crust. Syn‐breakup volcanism within the FSB suggests that depth‐dependent thinning was synchronous with continental breakup at the adjacent Faroes and Møre margins. We suggest that depth‐dependent continental lithospheric thinning can result from small‐scale convection that thins the lithosphere along multiple offset axes prior to continental rupture, leaving a failed breakup basin once seafloor spreading begins. This study provides insight into the structure and formation of a generic global class of ultra‐large rift basins formed by failed continental breakup.     

Rifting and pre‐rift lithosphere variability in the Orphan Basin,Newfoundland margin,Eastern Canada

The Orphan Basin, lying along the Newfoundland rifted continental margin, formed in Mesozoic time during the opening of the North Atlantic Ocean and the breakup of Iberia/Eurasia from North America. To investigate the evolution of the Orphan Basin and the factors that governed its formation, we (i) analysed the stratigraphic and crustal architecture documented by seismic data (courtesy of TGS), (ii) quantified the tectonic and thermal subsidence along a constructed geological transect, and (iii) used forward numerical modelling to understand the state of the pre‐rift lithosphere and the distribution of deformation during rifting. Our study shows that the pre‐rift lithosphere was 200‐km thick and rheologically strong (150‐km‐thick elastic plate) prior to rifting. It also indicates that extension in the Orphan Basin occurred in three distinct phases during the Jurassic, the Early Cretaceous and the Late Cretaceous. Each rifting phase is characterized by a specific crustal and subcrustal thinning configuration. Crustal deformation initiated in the eastern part of the basin during the Jurassic and migrated to the west during the Cretaceous. It was coupled with a subcrustal thinning which was reduced underneath the eastern domain and very intense in the western domains of the basin. The spatial and temporal distribution of thinning and the evolution of the lithosphere rheology through time controlled the tectonic, stratigraphic and crustal architecture that we observe today in the Orphan Basin.     

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.     

Interplay between lithospheric flexure and river transport in foreland basins

ABSTRACT Foreland basins form by lithospheric flexure under orogenic loading and are filled by surface transport of sediment. This work readdresses the interplay between these processes by integrating in a 3D numerical model: the mechanisms of thrust stacking, elastic flexural subsidence and sediment transport along the drainage network. The experiments show that both crustal tectonic deformation and vertical movements related to lithospheric flexure control and organise the basin-scale drainage pattern, competing with the nonlinear, unpredictable intrinsic nature of river network evolution. Drainage pattern characteristics are predicted that match those observed in many foreland basins, such as the axial drainage, the distal location of the main river within the basin, and the formation of large, long-lasting lacustrine systems. In areas where the river network is not well developed before the formation of the basin, these lithospheric flexural effects on drainage patterns may be enhanced by the role of the forebulge uplift as drainage divide. Inversely, fluvial transport modifies the flexural vertical movements differently than simpler transport models (e.g. diffusion): Rivers can drive erosion products far from a filled basin, amplifying the erosional rebound of both orogen and basin. The evolution of the sediment budget between orogen and basin is strongly dependent on this coupling between flexure and fluvial transport: Maximum sediment accumulations on the foreland are predicted for a narrow range of lithospheric elastic thickness between 15 and 40 km, coinciding with the T _e values most commonly reported for foreland basins.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

Sediment blanketing and the flexural strength of extended continental lithosphere

The flexural rigidity of the oceanic lithosphere is strongly dependent on its temperature structure at the time of loading. It is commonly assumed that the depth to the 450°C isotherm defines the effective elastic thickness T_e of the lithosphere. However, recent gravity studies across the Baltimore Canyon and Nova Scotian margins suggest that temperature may play a more complicated role in controlling the mechanical strength of extended continental lithosphere. For example, the flexural strength of the Baltimore Canyon margin (with sediment thicknesses of ? 15 km) appears to be controlled by the depth to the 150°C isotherm whereas the strength of the Nova Scotian margin (with sediment thicknesses cf ? 10 km) is controlled by the depth to the 250°C isotherm. The apparent correlation between sediment thickness and controlling isotherm suggests that sediment blanketing may play a role in modifying the flexural strength of extended continental lithosphere. This hypothesis was investigated by simulating the sedimentation history of a margin as a Gaussian function in which sedimentation peak and rate are determined by the mean and standard deviation of the function. The temperature structure of the lithosphere is continually modified as sediments are deposited on, and incorporated into the temperature structure of, the underlying lithosphere. Given a ‘starting’ value of T_e defined by the degree of extension of the lithosphere, the modification of T_e appears to be directly proportional to the sedimentation rate and cumulative sediment thickness, and inversely proportional to the time at which the sedimentation rate is a maximum. The first-order consequence of sediment blanketing is to reduce the cooling rate of the lithosphere relative to cooling in the absence of sediments. At thermal equilibrium, the initial value of T_e is reduced by the cumulative sediment thickness. Local isostatic conditions (i. e. T_e? 0) can only be approached when the sedimentation rate is unrealistically high (> 1000 m/Myr) during the rift or early post-rift phase of basin development. However, while these early loads may be locally compensated, any subsequent loads will be regionally compensated. Thus, it is unlikely that the low present-day flexural strengths interpreted from the Baltimore Canyon and Nova Scotian passive continental margins are a consequence of sediment blanketing.