Rift-phase extensional fabrics of the back-arc Drummond Basin, eastern Australia

The base of the Late Devonian–Early Carboniferous Drummond Basin, a major backarc extensional feature in eastern Australia which formed in response to detachment faulting, is extensively exposed in central Queensland. Here a crystalline basin floor is overlain by the Silver Hills Volcanics, a synrift sequence of predominantly silicic ash flow tuffs and lavas ranging to over 2 km in thickness. Detailed mapping of faults and stratigraphic logging of thickness changes within the Silver Hills Volcanics have allowed the rift-phase structural architecture that accompanied initial subsidence near the basin margin to be resolved. A complex mosaic of block faults with throws of up to 1 km is indicated. Locally developed mosaics may conform to, or depart from, the configuration predicted by the detachment faulting model. Structural fabric of the basement was a critical determinant of the extensional geometry. Distributed shear along pre-existing penetrative planar fabrics is considered to have accommodated hangingwall extension at lower strain rates whereas the propagation of tension fractures and the development of block faults by failure on pre-existing, brittle, basement dislocations facilitated extension at higher strain rates. The detachment fault inferred to lie beneath the extended hangingwall carapace has not been mapped at the surface and is thought to dissipate into a broad zone of distributed shear within basement to the east of the basin. Volcanism coincided with the initiation of extensional movements at which time deep crustal repositories for evolved magma were tapped by extensional fractures. The main extensional faults cutting the basinal succession were not used as conduits for magmatic products which were sourced from the basin margin and from extended hinterland to the east.     

Miocene-Pliocene half-graben evolution, detachment faulting and late-stage core complex uplift from reflection seismic data in south-east Arizona

Reflection seismic data show that the late Cenozoic Safford Basin in the Basin and Range of south-eastern Arizona, is a 4.5-km-deep, NW-trending, SW-dipping half graben composed of middle Miocene to upper Pliocene sediments, separated by a late Miocene sequence boundary into lower and upper basin-fill sequences. Extension during lower basin-fill deposition was accommodated along an E-dipping range-bounding fault comprising a secondary breakaway zone along the north-east flank of the Pinaleño Mountains core complex. This fault was a listric detachment fault, active throughout the mid-Tertiary and late Cenozoic, or a younger fault splay that cut or merged with the detachment fault. Most extension in the basin was accommodated by slip on the range-bounding fault, although episodic movement along antithetic faults temporarily created a symmetric graben. Upper-plate movement over bends in the range-bounding fault created rollover structures in the basin fill and affected deposition within the half graben. Rapid periods of subsidence relative to sedimentation during lower basin-fill deposition created thick, laterally extensive lacustrine or alluvial plain deposits, and restricted proximal alluvian-fan deposits to the basin margins. A period of rapid extension and subsidence relative to sediment influx, or steepening of the upper segment of the range-bounding fault at the start of upper basin-fill deposition resulted in a large downwarp over a major fault bend. Sedimentation was restricted to this downwarp until filled. Episodic subsidence during upper basin-fill deposition caused widespread interbedding of lacustrine and fluvial deposits. Northeastward tilting along the south-western flank of the basin and north-eastward migration of the depocentre during later periods of upper basin-fill deposition suggest decreased extension rates relative to late-stage core complex uplift.     

A New Model for Heat Flow in Extensional Basins: Estimating Radiogenic Heat Production

Radiogenic heat production (RHP) represents a significant fraction of surface heat flow, both on cratons and in sedimentary basins. RHP within continental crust—especially the upper crust—is high. RHP at any depth within the crust can be estimated as a function of crustal age. Mantle RHP, in contrast, is always low, contributing at most 1 to 2 mW/m² to total heat flow. Radiogenic heat from any noncrystalline basement that may be present also contributes to total heat flow. RHP from metamorphic rocks is similar to or slightly lower than that from their precursor sedimentary rocks. When extension of the lithosphere occurs—as for example during rifting—the radiogenic contribution of each layer of the lithosphere and noncrystalline basement diminishes in direct proportion to the degree of extension of that layer. Lithospheric RHP today is somewhat less than in the distant past, as a result of radioactive decay. In modeling, RHP can be varied through time by considering the half lives of uranium, thorium, and potassium, and the proportional contribution of each of those elements to total RHP from basement. RHP from sedimentary rocks ranges from low for most evaporites to high for some shales, especially those rich in organic matter. The contribution to total heat flow of radiogenic heat from sediments depends strongly on total sediment thickness, and thus differs through time as subsidence and basin filling occur. RHP can be high for thick clastic sections. RHP in sediments can be calculated using ordinary or spectral gamma-ray logs, or it can be estimated from the lithology.     

Stratigraphy and 40Ar/39Ar geochronology of the Santa Rosa basin,Baja California: Dynamic evolution of a constrictional rift basin during oblique extension in the Gulf of California

The Santa Rosa basin of northeastern Baja California is one of several transtensional basins that formed during Neogene oblique opening of the Gulf of California. The basin comprises Late Miocene to Pleistocene sedimentary and volcanic strata that define an asymmetric half‐graben above the Santa Rosa detachment, a low‐angle normal fault with ca. 4–5 km of SE‐directed displacement. Stratigraphic analysis reveals systematic basin‐scale facies variations both parallel and across the basin. The basin‐fill exhibits an overall fining‐upward cycle, from conglomerate and breccia at the base to alternating sandstone‐mudstone in the depocentre, which interfingers with the fault‐scarp facies of the detachment. Sediment dispersal was transverse‐dominated and occurred through coalescing alluvial fans from the immediate hanging wall and/or footwall of the detachment. Different stratigraphic sections reveal important lateral facies variations that correlate with major corrugations of the detachment fault. The latter represent extension‐parallel folds that formed largely in response to the ca. N‐S constrictional strain regime of the transtensional plate boundary. The upward vertical deflection associated with antiformal folding dampened subsidence in the northeastern Santa Rosa basin, and resulted in steep topographic gradients with a high influx of coarse conglomerate here. By contrast, the downward motion in the synform hinge resulted in increased subsidence, and led to a southwestward migration of the depocentre with time. Thus, the Santa Rosa basin represents a new type of transtensional rift basin in which oblique extension is partitioned between diffuse constriction and discrete normal faulting. ⁴⁰Ar/³⁹Ar geochronology of intercalated volcanic rocks suggests that transtensional deformation began during the Late Miocene, between 9.36 ± 0.14 Ma and 6.78 ± 0.12 Ma, and confirms previous results from low‐temperature thermochronology (Seiler et al., 2011). Two other volcanic units that appear to be part of a conformable syn‐rift sequence are, in fact, duplicates of pre‐rift volcanics and represent allochthonous, gravity‐driven slide blocks that originated from the hanging wall.     

Tectono-thermal modelling of sedimentary basins with episodic extension and inversion, a case history of the Jiyang Basin, North China

A two‐dimensional kinematic model is presented for superimposed basins. It is based on a finite‐element algorithm in the Lagrangian system, which incorporates different stages of lithosphere stretching and shortening to simulate alternating extension and inversion. The Jiyang Basin, developed in the North China, is a superimposed basin comprising four proto‐type basins separated by several unconformities. Four‐phase extension and two‐phase inversion have developed in this basin since the Late Mesozoic era. The thermal history of the basin is modelled based on a seismic cross‐section across the basin. Tectonic subsidence (or uplift) histories from backstripping serve as the objective functions, and crustal thickness, as well as heat flow, provides additional constraints. Effects of different Mesozoic erosion on the thermal history are discussed. Modelling results show that the thermal history of the Jiyang Basin since the Late Mesozoic can be divided into six stages, including four phases of heating accompanied by following thermal attenuation, and two phases of cooling with following thermal recovery. The model also implies a variant pattern of thermal regime in the basin. In the deepest centres of the depressions, the maximum heat flow occurred during the Late Mesozoic, but in the slopes of the depressions, the maximum heat flow appeared in the Cenozoic era.     

Subsidence and thermal history of an inverted Late Jurassic‐Early Cretaceous extensional basin (Cameros,North‐central Spain) affected by very low‐ to low‐grade metamorphism

The Cameros Basin (North Spain) is a Late Jurassic‐Early Cretaceous extensional basin, which was inverted during the Cenozoic. It underwent a remarkable thermal evolution, as indicated by the record of anomalous high temperatures in its deposits. In this study, the subsidence and thermal history of the basin is reconstructed, using subsidence analysis and 2D thermal modelling. Tectonic subsidence curves provide evidence of the occurrence of two rapid subsidence phases during the syn‐extensional stage. In the first phase (Tithonian‐Early Berriasian), the largest accommodation space was formed in the central sector of the basin, whereas in the second (Early Barremian‐Early Albian), it was formed in the northern sector. These rapid subsidence phases could correspond to relevant tectonic events affecting the Iberian Plate at that time. By distinguishing between the initial and thermal subsidence and defining their relative magnitudes, Royden's (1986) method was used to estimate the heat flow at the end of the extensional stage. A maximum heat flow of 60–65 mW/m² is estimated, implying only a minor thermal disturbance associated with extension. In contrast with these data, very high vitrinite reflectance, anomalously distributed in some case with respect to the typical depth‐vitrinite reflectance relation, was measured in the central‐northern sector of the basin. Burial and thermal data are used to construct a 2D thermal basin model, to elucidate the role of the processes involved in sediment heating. Calibration of the thermal model with the vitrinite reflectance (%Ro) and fluid inclusion (FI) data indicates that in the central and northern sectors of the basin, an extra heat source, other than a typical rift, is required to explain the observed thermal anomalies. The distribution of the %Ro and FI values in these sectors suggests that the high temperatures and their distribution are related to the circulation of hot fluids. Hot fluids were attributed to the hydrothermal metamorphic events affecting the area during the early post‐extensional and inversion stages of the 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.     

Testing models of rift basins: structure and stratigraphy of an Eocene–Oligocene supradetachment basin, Muddy Creek half graben, south-west Montana

We examine the basin geometry and sedimentary patterns in the Muddy Creek half graben of south-west Montana by integrating geological mapping, structural and basin analysis, ⁴⁰Ar/³⁹Ar geochronology, biostratigraphy and reflection seismic data. The half graben formed in late Middle Eocene to early Oligocene (?) time at the breakaway of a regional, WSW-dipping detachment system. Although the structure of the half graben is that of a supradetachment basin, facies patterns and basin architecture do not conform to a recent model for extensional basins above detachment faults. The border fault, the Muddy Creek fault system, consists of three en echelon, left-stepping normal faults separated by two relay ramps. The fault steepens southward toward each en echelon step, ranges in dip from 8 to 60° near the surface, but flattens at depths between 0 and 3 km. A broad ENE-plunging displacement-gradient syncline defines the central part of the half graben and is flanked by narrow SE-and NE-plunging anticlines to the north and south. Fine-grained deposits of the syntectonic basin-fill are thickest in the central syncline and interfinger with footwall-derived conglomerate near the adjacent anticlines. These facies patterns suggest that folding was coeval with extension and sedimentation in the half graben. Pre-extensional volcanic rocks and interbedded conglomerate filled a major ESE-trending palaeovalley along the future axis of the Muddy Creek half graben. Synextensional sedimentary deposits include lacustrine and paludal shale, mudstone and sandstone ponded in the centre of the half graben, and a narrow (typically <1.5 km wide) fringe of coarse alluvial-fan and fan-delta conglomerate and sandstone derived from the footwall. Angular unconformities and rock-slide deposits occur only locally within the syntectonic sequence. These facies patterns agree well with the half-graben depositional model of Leeder & Gawthorpe but not with a more recent supradetachment basin model of Friedmann & Burbank despite the demonstrably low dip-angle of the basin-bounding normal fault. These data show that it may not be possible to differentiate between supradetachment basins and half graben with steeper border faults using the architecture of the associated basin-fill deposits.     

Numerical models of a subsidence mechanism in intracratonic basins: application to North American basins

McKenzie's model of sedimentary basin evolution and its modification, widely used in geophysics, sometimes fails to explain discrepancies between predicted and observed values of extension, thinning and subsidence of the Earth's crust, as for the North Sea. We develop a numerical model of sedimentary basin evolution based on the mechanism suggested by Lobkovsky. In the course of rifting, accompanied by thinning of lower parts of the lithosphere, the roof of the underlying asthenosphere moves upward. the material of the mantle lifts and partially melts owing to the reduction of pressure. the density difference between the melt and the crystalline skeleton results in the filtration of the lighter melt and its accumulation in the form of a magmatic lens. Due to changed P-T conditions, the material of the lens undergoes the gabbro-eclogite phase transformation. the resultant anomalously heavy eclogite lens sinks in the surrounding material. This induces a viscous flow, changing the surface topography and forming a sedimentary basin. We construct a 2-D numerical model describing a viscous flow induced by subsidence of a heavy body and compute changes of surface topography. to compute the flow we employ the Galerkin-spline approach, with modifications allowing for density discontinuities and time dependence of the phase transformation. We apply the model to the cases of the Illinois, Michigan and Williston basins. the computed and tectonic subsidence curves agree well for these cases. the proposed model is compatible with the seismic structure of the crust and upper mantle below these basins. the model is also consistent with gravity data. the approach is applicable to other intracratonic basins.     

Stratigraphic signatures of translation of thrust-sheet top basins over low-angle detachment faults

Abstract Low‐angle detachment faults and thrust‐sheet top basins are common features in foreland basins. However, in stratigraphic analysis their influence on sequence architecture is commonly neglected. Usually, only eustatic sea level and changing flexural subsidence are accounted for, and when deformation is considered, the emphasis is on the generation of local thrust‐flank unconformities. This study analyses the effects of detachment angle and repetitive detachment activation on stratigraphic stacking patterns in a large thrust‐sheet top basin by applying a three‐dimensional numerical model. Model experiments show that displacement over low‐angle faults (2–6°) at moderate rates (～5.0 m kyr^?1) results in a vertical uplift component sufficient to counteract the background flexural subsidence rate. Consequently, the basin‐wide accommodation space is reduced, fluvio‐deltaic systems carried by the thrust‐sheet prograde and part of the sediment supply is spilled over towards adjacent basins. The intensity of the forced regression and the interconnectedness of fluvial sheet sandstones increases with the dip angle of the detachment fault or rate of displacement. In addition, the delta plain is susceptible to the formation of incised valleys during eustatic falls because these events are less compensated by regional flexural subsidence, than they would be in the absence of fault displacement.     

Models of the rapid post‐rift subsidence in the eastern Qiongdongnan Basin,South China Sea: implications for the development of the deep thermal anomaly

The Qiongdongnan Basin is one of the largest Cenozoic rifted basins on the northern passive margin of the South China Sea. It is well known that since the Late Miocene, approximately 10 Ma after the end of the syn‐rift phase, this basin has exhibited rapid thermal subsidence. However, detailed analysis reveals a two‐stage anomalous subsidence feature of the syn‐rift subsidence deficit and the well‐known rapid post‐rift subsidence after 10.5 Ma. Heat‐flow data show that heat flow in the central depression zone is 70–105 mW m^?2, considerably higher than the heat flow (<70 mW m^?2) on the northern shelf. In particular, there is a NE‐trending high heat‐flow zone of >85 mW m^?2 in the eastern basin. We used a numerical model of coupled geothermal processes, lithosphere thinning and depositional processes to analyse the origin of the anomalous subsidence pattern. Numerical analysis of different cases shows that the stretching factor β_s based on syn‐rift sequences is less than the observed crustal stretching factor β_c, and if the lithosphere is thinned with β_c during the syn‐rift phase (before 21 Ma), the present basement depth can be predicted fairly accurately. Further analysis does not support crustal thinning after 21 Ma, which indicates that the syn‐rift subsidence is in deficit compared with the predicted subsidence with the crustal stretching factor β_c. The observed high heat flow in the central depression zone is caused by the heating of magmatic injection equivalently at approximately 3–5 Ma, which affected the eastern basin more than the western basin, and the Neogene magmatism might be fed by the deep thermal anomaly. Our results suggest that the causes of the syn‐rift subsidence deficit and rapid post‐rift subsidence might be related. The syn‐rift subsidence deficit might be caused by the dynamic support of the influx of warmer asthenosphere material and a small‐scale thermal upwelling beneath the study area, which might have been persisting for about 10 Ma during the early post‐rift phase, and the post‐rift rapid subsidence might be the result of losing the dynamic support with the decaying or moving away of the deep thermal source, and the rapid cooling of the asthenosphere. We concluded that the excess post‐rift subsidence occurs to compensate for the syn‐rift subsidence deficit, and the deep thermal anomaly might have affected the eastern Qiongdongnan Basin since the Late Oligocene.     

Formation and fragmentation of a late Miocene supradetachment basin in central Crete: implications for exhumation mechanisms of high‐pressure rocks in the Aegean forearc

We present a new lithostratigraphy and chronology for the Miocene on central Crete, in the Aegean forearc. Continuous sedimentation started at ～10.8 Ma in the E–W trending fluvio‐lacustrine Viannos Basin, formed on the hangingwall of the Cretan detachment, which separates high‐pressure (HP) metamorphic rocks from very low‐grade rocks in its hangingwall. Olistostromes including olistoliths deposited shortly before the Viannos Basin submerged into the marine Skinias Basin between 10.4 and 10.3 Ma testifies to significant nearby uplift. Uplift of the Skinias Basin between 9.7 and 9.6 Ma, followed by fragmentation along N–S and E–W striking normal faults, marks the onset of E–W arc‐parallel stretching superimposed on N–S regional Aegean extension. This process continued between 9.6 and 7.36 Ma, as manifested by tilting and subsidence of fault blocks with subsidence events centred at 9.6, 8.8, and 8.2 Ma. Wholesale subsidence of Crete occurred from 7.36 Ma until ～5 Ma, followed by Pliocene uplift and emergence. Subsidence of the Viannos Basin from 10.8 to 10.4 Ma was governed by motion along the Cretan detachment. Regional uplift at ～10.4 Ma, followed by the first reworking of HP rocks (10.4–10.3 Ma) is related to the opening and subsequent isostatic uplift of extensional windows exposing HP rocks. Activity of the Cretan detachment ceased sometime between formation of extensional windows around 10.4 Ma, and high‐angle normal faulting cross‐cutting the detachment at 9.6 Ma. The bulk of exhumation of the Cretan HP‐LT metamorphic rocks occurred between 24 and 12 Ma, before basin subsidence, and was associated with extreme thinning of the hangingwall (by factor ～10), in line with earlier inferences that the Cretan detachment can only explain a minor part of total exhumation. Previously proposed models of buyoant rise of the Cretan HP rocks along the subducting African slab provide an explanation for extension without basin subsidence.     

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.     

Unravelling the multi-stage burial history of the Swiss Molasse Basin: integration of apatite fission track, vitrinite reflectance and biomarker isomerisation analysis

A complex basin evolution was studied using various methods, including thermal constraints based on apatite fission‐track (AFT) analysis, vitrinite reflectance (VR) and biomarker isomerisation, in addition to a detailed analysis of the regional stratigraphic record and of the lithological properties. The study indicates that (1) given the substantial amount of data, the distinction and characterisation of successive stages of heating and burial in the same area are feasible, and (2) the three thermal indicators (AFT, VR and biomarkers) yield internally consistent thermal histories, which supports the validity of the underlying kinetic algorithms and their applicability to natural basins. All data pertaining to burial and thermal evolution were integrated in a basin model, which provides constraints on the thickness of eroded sections and on heat flow over geologic time. Three stages of basin evolution occurred in northern Switzerland. The Permo‐Carboniferous strike–slip basin was characterised by high geothermal gradients (80–100°C km^?1) and maximum temperature up to 160°C. After the erosion of a few hundreds of metres in the Permian, the post‐orogenic, epicontinental Mesozoic basin developed in Central Europe, with subsidence triggered by several stages of rifting. Geothermal gradients in northern Switzerland during Cretaceous burial were relatively high (35–40°C km^?1), and maximum temperature typically reached 75°C (top middle Jurassic) to 100°C (base Mesozoic). At least in the early Cretaceous, a stage of increased heat flow is needed to explain the observed maturity level. After erosion of 600–700 m of Cretaceous and late Jurassic strata during the Paleocene, the wedge‐shaped Molasse Foreland Basin developed. Geothermal gradients were low at this time (≤20°C km^?1). Maximum temperature of Miocene burial exceeded that of Cretaceous burial in proximal parts (<35 km from the Alpine front), but was lower in more distal parts (>45 km). Thus, maximum temperature as well as maximum burial depth ever reached in Mesozoic strata occurred at different times in different regions. Since the Miocene, 750–1050 m were eroded, a process that still continues in the proximal parts of the basin. Current average geothermal gradients in the uppermost 2500 m are elevated (32–47°C km^?1). They are due to a Quaternary increase of heat flow, most probably triggered by limited advective heat transport along Paleozoic faults in the crystalline basement.     

Basin modelling of the Limón Back-arc Basin (Costa Rica): burial history and temperature evolution of an island arc-related basin-system

The Limón back‐arc basin belongs to the southern Central American arc‐trench system and is situated at the east coast of Costa Rica. The basin‐fill consists of Late Cretaceous to Pleistocene sedimentary rocks. A northern and a southern sub‐basin can be defined, separated by the E–W‐trending Trans Isthmic Fault System. The North Limón Basin is nearly undeformed, whereas the South Limón Basin is characterized by a fold‐and‐thrust belt. Both sub‐basins have a very similar sedimentary fill and can act as a natural laboratory for reconstructing controlling factors of arc‐related sedimentary basins as well as the influence of deformation on a basin system. Modelling focused on burial history and temperature evolution. Two‐dimensional simulations were carried out with the software PetroMod^®. The geohistory curve of the North Limón Basin is overall linear, indicating continuous subsidence. The South Limón Basin is also characterized by continuous subsidence, but rates strongly increased at the beginning of the Neogene. Despite a rapid Plio‐Pleistocene deformation of the fold‐and‐thrust belt, the present‐day temperature field is not disturbed in that area. The modelling results indicate a mean heat flow of 60 mW m^?2 for the North Limón Basin and 41 mW m^?2 for the South Limón Basin. These values are low compared with other back‐arc basins. The lower values are attributed to the following effects: (1) underlying basaltic crust, (2) the lack of an initial rift phase, (3) the low extension rates, (4) absence of volcanic activity and (5) insulation effects of a thick sediment pile. The reasons for the locally lower heat flow in the southern sub‐basin can be found in the low‐angle subduction of the Cocos Ridge. Owing to the low subduction angle, the cool fore‐arc mantle‐wedge below the island‐arc is pushed backwards increasing the cooled area.     

Half-graben basin filling models: new constraints on continental extensional basin development

Three end-member models of half-graben development (detachment fault, domino-style, and fault growth) evolve differently through time and produce different basin-filling patterns. The detachment fault model incorporates a basin-bounding fault that soles into a subhorizontal detachment fault; the change in the rate of increase in the volume of the basin during uniform fault displacement is zero. Younger strata consistently pinch out against older synrift strata rather than pre-rift rocks. Both basin-bounding faults and the intervening fault blocks rotate during extension in the domino fault block model; a consequence of this rotation is that the change in the rate of increase of the volume of the basin is negative during uniform extension. Basin fill commonly forms a fanning wedge during fluvial sedimentation, whereas lacustrine strata tend to pinch out against older synrift strata. In the fault growth models, basins grow both wider and longer through time as the basin-bounding faults lengthen and displacement accumulates; the change in the rate of increase in basin volume is positive. Fluvial strata progressively onlap pre-rift rocks of the hanging wall block, whereas lacustrine strata pinch out against older fluvial strata at the centre of the basin but onlap pre-rift rocks along the lateral edges. These fundamental differences may be useful in discriminating among the three end-member models. The transition from fluvial to lacustrine deposition and hanging wall onlap relationships observed in numerous continental extensional basins are best explained by the fault growth models.     

Extension and subsidence of the Pearl River Mouth Basin, northern South China Sea

Abstract The uniform stretching model has been applied to seismic reflection profiles and well-log information from the Pearl River Mouth Basin on the northern flank of the South China Sea. Stretching factors were calculated from subsidence curves determined from the stratigraphy by using the backstripping technique to remove the effects of compaction and sediment loading. Variations in rift topography, palaeobathymetry and global sea-level v/ere taken into account. We argue that the Pearl River Mouth Basin formed by lithospheric extension by a factor of about 1.8, lasting from Late Cretaceous to late Oligocene times. Stretching factors calculated from subsidence agree with those determined from the geometry of normal faulting and from crustal thinning. Thus there is no indication of a significant discrepancy between the different estimates of stretching. The geometry of faulting suggests that considerable amounts of local footwall uplift occurred during the rifting period. Small differences between the observed and calculated subsidence curves (∽ 400 m in the middle Miocene) are best explained by minor amounts of extension ( β ∽ 1.1). The time-temperature history of sediments within the basin has also been calculated so that expected vitrinite reflectance and oil abundance could be determined. The results are consistent with each other and are in reasonable agreement with observations from wells.     

Mudstone compaction curves in basin modelling: a study of Mesozoic and Cenozoic Sediments in the northern North Sea

Basin modelling studies are carried out in order to understand the basin evolution and palaeotemperature history of sedimentary basins. The results of basin modelling are sensitive to changes in the physical properties of the rocks in the sedimentary sequences. The rate of basin subsidence depends, to a large extent, on the density of the sedimentary column, which is largely dependent on the porosity and therefore on the rate of compaction. This study has tested the sensitivity of varying porosity/depth curves and related thermal conductivities for the Cenozoic succession along a cross‐section in the northern North Sea basin, offshore Norway. End‐member porosity/depth curves, assuming clay with smectite and kaolinite properties, are compared with a standard compaction curve for shale normally applied to the North Sea. Using these alternate relationships, basin geometries of the Cenozoic succession may vary up to 15% from those predicted using the standard compaction curve. Isostatic subsidence along the cross‐section varies 2.3–4.6% between the two end‐member cases. This leads to a 3–8% difference in tectonic subsidence, with maximum values in the basin centre. Owing to this, the estimated stretching factors vary up to 7.8%, which further gives rise to a maximum difference in heat flow of more than 8.5% in the basin centre. The modelled temperatures for an Upper Jurassic source rock show a deviation of more than 20 °C at present dependent on the thermal conductivity properties in the post‐rift succession. This will influence the modelled hydrocarbon generation history of the basin, which is an essential output from basin modelling analysis. Results from the northern North Sea have shown that varying compaction trends in sediments with varying thermal properties are important parameters to constrain when analysing sedimentary basins.     

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.     

Effects of active crustal movements on thermal structure in subduction zones

In young suduction zones we observe steady uplift of island arcs. The steady uplift of island arcs is always accompanied by surface erosion. The long duration of uplift and erosion effectively transports heat at depth to shallower parts by advection. If the rates of uplift and erosion are sufficiently large, such a process of heat transportation will strongly affect thermal structure in subduction zones. First, we quantitatively examine the effects of uplift and erosion on thermal structure by using a simple 1-D heat conduction model, based on the assumption that the initial thermal state is in equilibrium. The results show that temperature increase, Δ T  , due to uplift and erosion can be approximately evaluated by Δ T  = ν _e tβ at depth, where ν _e is the rate of uplift (erosion), t is the duration of uplift (erosion), and β is the gradient of the geotherm in the initial state. Next, considering the effects of vertical crustal movements such as uplift and erosion in island arcs and subsidence and sedimentation in ocean trenches, in addition to the effects of radioactive heat generation in the crust, frictional heating at plate boundaries and accretion of oceanic sediments to overriding continental plates, we numerically simulate the evolution process of the thermal structure in subduction zones. The result shows that the temperature beneath the island arc gradually increases as a result of uplift and erosion as plate subduction progresses. Near the ocean trench, on the other hand, the low-temperature region gradually expands as a result of sedimentation and accretion in addition to direct cooling by the cold descending slab. The surface heat flow expected from this model is low in fore-arc basins, high in island arcs and moderately high in back-arc regions.