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 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.     

Tectonic and thermal evolution of Queen Charlotte Basin: lithospheric deformation and subsidence models

A two-layer lithospheric stretching model that includes the effects of decompression melting was used to estimate the deformation and thermal evolution of the Queen Charlotte Basin, British Columbia. The basin contains up to 6 km of Tertiary fill and is postulated to have been formed during a transtensional stage of Cenozoic plate motion between the Pacific and North American plates. Several models of basin formation have been proposed to explain the sediment distribution, contemporaneous volcanism and high present-day heat flow. We used bathymetry, Tertiary sediment thickness and crustal thickness to calculate the amount of stretching in the crust and lower lithosphere, and the volume of melt generated during advection of mantle rocks. A second set of calculations traced the thermal evolution of the sediments and lithosphere, and we show maps of estimated present-day heat flow and sediment maturity. This study differs significantly from previous work in the use of gridded data that provide coverage over a large region and permit lateral variations in lithospheric deformation and thermal properties to be clearly defined, a difficult quest in studies based on single-point or profile data. In addition, the use of crustal thickness, derived from a regional interpretation of gravity data and constrained by seismic refraction results, as an input allows reliable estimates of extension to be made despite recent deformation of sedimentary strata in Hecate Strait. We present results for a model which used a prerift crustal thickness of ≈34 km and a short rifting period from 25 to 20 Ma. This model infers that significant thinning occurred beneath south-western Hecate Strait and southern Queen Charlotte Sound, and several kilometres of igneous crust were added at these sites, without requiring elevated asthenospheric temperatures prior to extension. Net lithospheric extension is surprisingly uniform within the basin and averages 76%, or ≈50 km, across the margin. This amount is consistent with other estimates of extension and may provide information useful in refining models of plate motion along this margin.     

Interpretation of migrated seismic reflection profiles across the northern Appalachians in Maine

Summary. A simple line-migration technique has been developed and applied to deep crustal reflection data collected in the northern Appalachians of Quebec and central Maine. Preliminary interpretation of these data, combined with the results of gravity and refraction studies, shows thinning of the crust in two distinct steps beneath the Silurian-Devonian Merrimack synclinorium and shows evidence for Mesozoic crustal stretching and dike intrusion.     

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.     

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.     

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.     

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.     

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.     

Older crust underlies Iceland

The oldest rocks outcropping in northwest Iceland are ∼16 Myr old and in east Iceland ∼13 Myr. The full plate spreading rate in this region during the Cenozoic has been ∼2 cm a⁻¹, and thus these rocks are expected to be separated by ∼290 km. They are, however, ∼500 km apart. The conclusion is inescapable that an expanse of older crust ∼210 km wide underlies Iceland, submerged beneath younger lavas. This conclusion is independent of any considerations regarding spreading ridge migrations, jumps, the simultaneous existence of multiple active ridges, three-dimensionality, or subsidence of the lava pile. Such complexities bear on the distribution and age of the older crust, but not on its existence or its width. If it is entirely oceanic its maximum age is most likely 26–37 Ma. It is at least 150 km in north–south extent, but may taper and extend beneath south Iceland. Part of it might be continental—a southerly extension of the Jan Mayen microcontinent. This older crust contributes significantly to crustal thickness beneath Iceland and the ∼40 km local thickness measured seismically is thus probably an overestimate of present-day steady-state crustal production at Iceland.     

Anomalous passive subsidence of deep‐water sedimentary basins: a prearc basin example,southern New Caledonia Trough and Taranaki Basin,New Zealand

Stratigraphic data from petroleum wells and seismic reflection analysis reveal two distinct episodes of subsidence in the southern New Caledonia Trough and deep‐water Taranaki Basin. Tectonic subsidence of ~2.5 km was related to Cretaceous rift faulting and post‐rift thermal subsidence, and ~1.5 km of anomalous passive tectonic subsidence occurred during Cenozoic time. Pure‐shear stretching by factors of up to 2 is estimated for the first phase of subsidence from the exponential decay of post‐rift subsidence. The second subsidence event occured ~40 Ma after rifting ceased, and was not associated with faulting in the upper crust. Eocene subsidence patterns indicate northward tilting of the basin, followed by rapid regional subsidence during the Oligocene and Early Miocene. The resulting basin is 300–500 km wide and over 2000 km long, includes part of Taranaki Basin, and is not easily explained by any classic model of lithosphere deformation or cooling. The spatial scale of the basin, paucity of Cenozoic crustal faulting, and magnitudes of subsidence suggest a regional process that acted from below, probably originating within the upper mantle. This process was likely associated with inception of nearby Australia‐Pacific plate convergence, which ultimately formed the Tonga‐Kermadec subduction zone. Our study demonstrates that shallow‐water environments persisted for longer and their associated sedimentary sequences are hence thicker than would be predicted by any rift basin model that produces such large values of subsidence and an equivalent water depth. We suggest that convective processes within the upper mantle can influence the sedimentary facies distribution and thermal architecture of deep‐water basins, and that not all deep‐water basins are simply the evolved products of the same processes that produce shallow‐water sedimentary basins. This may be particularly true during the inception of subduction zones, and we suggest the term ‘prearc’ basin to describe this tectonic setting.     

Crustal structure and riftogenesis of the Valencia Trough (north-western Mediterranean Sea)

The Valencia Trough is a rift formed during the late Oligocene – early Miocene opening of the western Mediterranean Sea. In this paper, we focus on the crustal structure and on the deep structure of the basin which is hard to delineate because of the widespread volcanism that conceals part of the basement. This work is the result of the study of a dense network of seismic profiling surveys and exploratory wells made in the region. The structure of the deep basement reveals the importance of transfer fracture zones which represent steps in the deepening of the basin. The thinning of the crust follows the basement deepening and we find the same partitioning of structural blocks at the crustal level. Transfer faults also represent limits in the thinning of the crust and each compartment thus delineated has a different thinning and different extensional ratios. Such a discrepancy between the thinning of the upper crust and the thinning of the lower crust may be common in many other rift zones, but is seldom as well imaged as in this study of the Valencia Trough. The transfer zones are related to extensional processes but a simple shear opening is envisaged to explain the discrepancies between thinning and extension and the asymmetry of the margins. The more efficient thinning in the lower crust can be explained by a thermal anomaly in accordance with the recent evolution of the trough. The steady thinning of the margins is discussed in terms of a marginal basin in a compressional context.     

Teleseismic delay-time tomography of the upper mantle beneath southeastern India: imprint of Indo–Antarctica rifting

A 3-D P -velocity map of the crust and upper mantle beneath the southeastern part of India has been reconstructed through the inversion of teleseismic traveltimes. Salient geological features in the study region include the Archean Dharwar Craton and Eastern Ghat metamorphic belt (EGMB), and the Proterozoic Cuddapah and Godavari basins. The Krishna–Godavari basin, on the eastern coastal margin, evolved in response to the Indo–Antarctica breakup. A 24-station temporary network provided 1161 traveltimes, which were used to model 3-D P -velocity variation. The velocity model accounts of 80 per cent of the observed data variance. The velocity picture to a depth of 120 km shows two patterns: a high velocity beneath the interior domain (Dharwar craton and Cuddapah basin), and a lower velocity beneath the eastern margin region (EGMB and coastal basin). Across the array velocity variations of 7–10 per cent in the crust (0–40 km) and 3–5 per cent in the uppermost mantle (40–120 km) are observed. At deeper levels (120–210 km) the upper-mantle velocity differences are insignificant among different geological units. The presence of such a low velocity along the eastern margin suggests significantly thin lithosphere (<100 km) beneath it compared to a thick lithosphere (>200 km) beneath the eastern Dharwar craton. Such lithospheric thinning could be a consequence of Indo–Antarctica break-up.     

Anomalous seismic crustal structure of oceanic fracture zones

Summary. The seismic structure of crust found within fracture zones falls outside the range of velocity structures observed for normal oceanic crust in the North Atlantic. The crust in fracture zones is frequently very thin and is characterized by low crustal velocities and by the conspicuous absence of a refractor with a velocity typical of oceanic layer 3. Anomalous crust is present in both large- and small-offset fracture zones. Since they are among the most common tectonic features in the ocean basins, and are particularly closely spaced on slow-spreading ridges, fracture zones represent a major source of seismic crustal heterogeneity. We interpret the anomalous crust as a thin, intensely fractured, faulted and hydrothermally altered basaltic and gabbroic section overlying ultramafics that, in places, are extensively serpentinized. The unusually thin crust found within fracture zones and the gradual crustal thinning over a distance of several tens of kilometres on either side of the fracture zones can be explained by two main processes; firstly the cold lithosphere edge opposite the spreading centre at the ridgetransform intersection modifies the normal intrusive and extrusive processes of the spreading centre leading to the accretion of an anomalous and thin igneous section; and secondly each spreading ridge segment is fed from a separate subcrustal magma supply point, so as the magma flows laterally down the spreading centre it generates a crustal section of decreasing thickness, culminating in the very thin crust of the fracture zones at either end of the ridge segment.     

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.     

The seismic structure of thinned continental crust in the northern Bay of Biscay

Summary. The stretching and thinning of the continental crust, which occurs during the formation of passive continental margins, may cause important changes in the velocity structure of such crust. Further, crust attenuated to a few kilometres' thickness, can be found underlying 'oceanic' water depths. This paper poses the question of whether thinned continental crust can be distinguished seismically from normal oceanic crust of about the same thickness. A single seismic refraction line shot over thinned continental crust as part of the North Biscay margin transect in 1979 was studied in detail. Tau— p inversion suggested that there are differences between oceanic and continental crust in the lower crustal structure. This was confirmed when synthetic seismograms were calculated. The thinned continental crust (β± 7.0) exhibits a two-gradient structure in the non-sedimentary crust with velocities between 5.9 and 7.4 km s⁻¹; an upper 0.8 s⁻¹ layer overlies a 0.4 s⁻¹ layer. No layer comparable to oceanic layer 3 was detected. The uppermost mantle also contains a low-velocity zone.     

Tectonics of the late Mesozoic wide extensional basin system in the China–Mongolia border region

The China–Mongolia border region contains many late Mesozoic extensional basins that together constitute a regionally extensive basin system. Individual basins within the system are internally composed of a family of sub‐basins filled with relatively thin sedimentary piles mostly less than 5 km in thickness. There are two types of sub‐basins within the basins, failed and combined, respectively. The failed sub‐basins are those that failed to continue developing with time. In contrast, the combined ones are those that succeeded in growing by coalescing adjacent previously isolated sub‐basins. Thus, a combined sub‐basin is bounded by a linked through‐going normal fault that usually displays a corrugated trace on map view and a shallower dip on cross‐section. Along‐strike existence of discrete depocenters and alternation of sedimentary wedges of different types validate the linkage origin of combined sub‐basins. Localized high‐strain extension resulted in large‐amount displacement on linked faults, but contemporaneously brought about the cessation of some isolated fault segments and the formation of corresponding failed sub‐basins in intervening areas between active linked faults. Some combined sub‐basins might have evolved into supradetachment basins through time, concurrent with rapid denudation of footwall rocks and formation of metamorphic core complexes in places. A tectonic scenario of the broad basin system can be envisioned as an evolution from early‐stage distributed isolated sub‐basins to late‐stage focused combined or/and supradetachment sub‐basins bounded by linked faults, accompanied by synchronous cessation of some early‐formed sub‐basins. Initiation of the late Mesozoic extension is believed to result from gravitational collapse of the crust that had been overthickened shortly prior to the extension. Compression, arising from collision of Siberia and the amalgamated North China–Mongolia block along the Mongol–Okhotsk suture in the time interval from the Middle to Late Jurassic, led to significant shortening and thickening over a broad area and subsequent extensional collapse. Pre‐ and syn‐extensional voluminous magmatism must have considerably reduced the viscosity of the overthickened crust, thereby not only facilitating the gravitational collapse but enabling the lower‐middle crust to flow as well. Flow of a thicker crustal layer is assumed to have occurred coevally with upper‐crustal stretching so as to diminish the potential contrast of crustal thickness by repositioning materials from less extended to highly extending regions. Lateral middle‐ and lower‐crustal flow and its resultant upward push upon the upper crust provide a satisfying explanation for a number of unusual phenomena, such as supracrustal activity of the extension, absence or negligibleness of postrift subsidence of the basin system, less reduction of crustal thickness after extension, and non‐compression‐induced basin inversion, all of which have been paradoxical in the previous study of the late Mesozoic basin tectonics in the China–Mongolia border region.     

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.     

Crustal P- and S-velocity structure of the Serbomacedonian Massif (Northern Greece) obtained by non-linear inversion of traveltimes

In the present study, the P - and S -velocity structure of the crust and uppermost mantle in the area of central Macedonia (northern Greece) is presented, as derived from the inversion of traveltimes of local events. An appropriate preconditioning of the final linearized system is used in order to reduce ray density effects on the results. The study focuses mainly on the structure of the broader area of the Serbomacedonian Massif. Interesting features and details of the crustal structure can be recognized in the final tomographic images. The crustal thickness shows strong variations. Under the Serbomacedonian and western Rhodope massifs the crust has a thickness that exceeds 30  km. On the other hand, the North Aegean Trough exhibits a fairly thin crust (25–27  km). Moreover, the Serbomacedonian Massif is bounded by two regions that trend parallel to the Axios river–Thermaikos gulf and the Strymon river–Orfanou gulf, respectively, which show significant crustal thinning (25–28  km). The observed match between the direction of this crustal thinning and the basins' axes indicates that they have been generated by the same extensional deformation episode.     

Development of sedimentary basins: differential stretching,phase transitions,shear heating and tectonic pressure

Classical models of lithosphere thinning predict deep synrift basins covered by wider and thinner post‐rift deposits. However, synextensional uplift and/or erosion of the crust are widely documented in nature (e.g. the Base Cretaceous unconformity of the NE Atlantic), and generally the post‐rift deposits dominate basins fills. Accordingly, several basin models focus on this discrepancy between observations and the classical approach. These models either involve differential thinning, where the mantle thins more than the crust thereby increasing average temperature of the lithosphere, or focus on the effect of metamorphic reactions, showing that such reactions decrease the density of lithospheric rocks. Both approaches result in less synrift subsidence and increased post‐rift subsidence. The synextensional uplift in these two approaches happens only for special cases, that is for a case of initially thin crust, specific mineral assemblage of the lithospheric mantle or extensive differential thinning of the lithosphere. Here, we analyse the effects of shear heating and tectonic underpressure on the evolution of sedimentary basins. In simple 1D models, we test the implications of various mechanisms in regard to uplift, subsidence, density variations and thermal history. Our numerical experiments show that tectonic underpressure during lithospheric thinning combined with pressure‐dependent density is a widely applicable mechanism for synextensional uplift. Mineral phase transitions in the subcrustal lithosphere amplify the effect of underpressure and may result in more than 1 km of synextensional erosion. Additional heat from shear heating, especially combined with mineral phase transitions and differential thinning of the lithosphere, greatly decreases the amount of synrift deposits.