Back-Stripping: Definition, Basin Analysis, and Subsidence History

Back-stripping is a quantitative geological technique that reconstructs the burial and subsidence history of a sedimentary basin by sequentially removing the youngest sedimentary layer, correcting for compaction, paleobathymetry, and isostasy, and iterating backward through geologic time to reveal how and why the basin floor subsided. First formalized by Sclater and Christie (1980) using borehole data from the North Sea, the method separates the total observed subsidence into two independent components: a sediment-loading component that results from the weight of deposited sediment deflecting the lithosphere, and a tectonic subsidence component that reflects genuine crustal extension or thermal cooling. By isolating tectonic subsidence, geoscientists can reconstruct the stretching history of the crust, constrain the timing of syn-rift and post-rift phases, and calibrate basin models used for petroleum system analysis. Back-stripping is today a standard tool in sequence stratigraphy, reservoir characterization modeling, and source rock maturity analysis, with applications ranging from frontier exploration in the Arctic to production optimization in mature basins.

The Geological Foundation: Why Basins Subside

Sedimentary basins form wherever the crust subsides below the surrounding level and accumulates sediment. The mechanisms that drive subsidence are varied: rifting and crustal thinning (extensional basins like the North Sea, Gulf of Mexico, and East African Rift), thermal contraction of the lithosphere after rifting (post-rift sag basins), flexural loading by thrust sheets (foreland basins like the Alberta Basin and Appalachian Basin), and sediment loading in passive margin settings. Understanding which mechanism operated, when it operated, and how intensely it operated is central to predicting the petroleum system: whether a source rock was buried deep enough to generate hydrocarbons, whether a structural trap existed at the time of peak generation, and whether migration pathways were open or blocked. Back-stripping provides the quantitative framework to answer these questions by generating a time-calibrated subsidence curve from rock data that is available today in every exploration or development well.

The fundamental challenge is that the rock record as observed today has been modified by two processes that obscure the original depositional geometry: compaction, which progressively reduces porosity and thickness as sediment is buried under increasing overburden stress; and sea-level change, which alters the apparent water depth in marine sections. Both must be removed before the tectonic signal can be read. Decompaction reverses the porosity loss using an empirical relationship between porosity and depth that varies by lithology. Paleobathymetric correction restores the water depth at the time of deposition using microfossil assemblages as paleodepth proxies. Isostatic correction accounts for the deflection of the lithosphere by the weight of the sediment column itself. Only after all three corrections are applied does the residual subsidence represent the genuine tectonic driving force.

The concept of an accumulation of hydrocarbons in a structural or stratigraphic trap is inseparable from the subsidence history of the basin that produced both the source rock and the reservoir. If the trap formed after the main pulse of hydrocarbon generation and expulsion, the basin may be petroliferous but untrapped: hydrocarbons migrated through and escaped before the structural closure existed. Back-stripping allows geoscientists to date the formation of structural closures and compare that timing to the maturity history of the source rock, directly assessing trap-timing risk in frontier and emerging basins.

The 1D Back-Stripping Method: Step by Step

One-dimensional back-stripping proceeds through a well section from youngest to oldest, removing one stratigraphic layer at a time. For each removal step, the analyst performs four corrections in sequence. First, decompaction restores the removed layer and all underlying layers to their original (pre-burial) thickness using the exponential porosity-depth relationship formulated by Athy (1930) and extended by Sclater and Christie (1980): porosity equals the initial depositional porosity multiplied by e raised to the power of negative c times depth, where c is the compaction coefficient in units of inverse length, which is lithology-dependent. Typical values for c range from 0.27 per kilometer for sandstone to 0.71 per kilometer for shale. Decompacting removes the porosity that was lost during burial, restoring the layer to its original depositional volume.

Second, a paleobathymetric correction is applied to restore the water depth at the time the layer was deposited. This is determined from benthic foraminifera assemblages, trace fossil tiering depth, sedimentary facies analysis, and, where available, oxygen isotope and Mg/Ca paleothermometry data. For continental deposits, paleobathymetry is zero by definition. For deep-water turbidite systems, it may reach 2,000 to 4,000 meters (6,562 to 13,123 feet), introducing a large and uncertain correction that propagates directly into the tectonic subsidence estimate. Third, a paleo-sea level correction is applied using global sea level curves (Haq et al. 1987; Miller et al. 2011) adjusted for local evidence where available. Fourth, an isostatic correction accounts for the deflection of the lithosphere under the sediment load. Two end-member models exist: Airy isostasy, which assumes local isostatic equilibrium with no lateral stress transmission, and flexural isostasy, which treats the lithosphere as a elastic plate with a characteristic flexural rigidity (expressed as the effective elastic thickness Te). For most intracratonic basins and rifts, Airy isostasy is adequate; for foreland basins and passive margins with stiff lithosphere, flexural isostasy produces significantly more accurate reconstructions.

The output of these corrections for each time step is the tectonic subsidence at that moment in the basin's history. Plotting tectonic subsidence against time generates the tectonic subsidence curve, the fundamental diagnostic product of the back-stripping workflow. The shape of this curve reveals the dominant subsidence mechanism. An initial rapid subsidence followed by an exponential decay (thermal subsidence) is diagnostic of rifting followed by thermal cooling, as described by the McKenzie (1978) stretching model. A uniform rate of subsidence sustained over tens of millions of years suggests flexural loading by an advancing thrust sheet. A two-phase pattern with a rapid early phase and a later reactivation may indicate polyphase rifting or inversion tectonics superimposed on an earlier extensional basin.

The McKenzie Stretching Model and the Beta Factor

The McKenzie (1978) pure-shear stretching model provides the theoretical linkage between the tectonic subsidence curve derived from back-stripping and the physical mechanism of crustal extension. The model defines the stretching factor beta as the ratio of the original crustal thickness to the thinned crustal thickness after rifting. A beta of 1.0 represents no extension; a beta of 2.0 represents a two-fold extension of the crust, with crustal thickness halved. The model predicts two phases of subsidence: an instantaneous syn-rift phase caused by the isostatic response to crustal thinning and densification of the lithospheric mantle, and a long-duration post-rift thermal subsidence phase caused by the gradual cooling and thermal contraction of the upwelled asthenosphere. The duration and magnitude of the thermal subsidence phase depend on beta; higher beta values produce faster and deeper thermal subsidence over time scales of 50 to 200 million years.

By fitting the McKenzie model to the tectonic subsidence curve extracted through back-stripping, analysts can determine the best-fit beta value for a given well location or for a map of wells across a basin. This beta map directly describes the pattern of crustal extension, and because extension correlates with the depth to basement and the structural relief on rift-bounding faults, it is a first-order predictor of where the deepest kitchen areas are located (high beta = deepest burial = highest maturity) and where structural traps are likely to have formed on the shoulders of half-grabens (moderate beta). In the North Sea, where back-stripping has been applied to thousands of exploration wells since the 1980s, the spatial variation of beta from approximately 1.3 in the Central Graben flanks to greater than 2.0 in the Viking Graben axis controls the distribution of oil and gas fields as predicted by the McKenzie model and as confirmed by decades of drilling.

The relationship between the asthenosphere and lithospheric thermal state is central to the post-rift subsidence history. Where the asthenosphere is upwelled during rifting, high heat flow drives accelerated maturation of organic matter in source rocks, particularly type II marine kerogens (e.g., Kimmeridge Clay in the North Sea, Niobrara in the Denver Basin, Eagle Ford in Texas). Back-stripping combined with basin thermal modeling allows petroleum system analysts to estimate paleo-heat flow at any point in the burial history, converting the depth-temperature-time path into a vitrinite reflectance equivalent (%Ro) or transformation ratio (TR) for the source rock. This maturity timeline, calibrated against the structural timing derived from back-stripping, is the core input to risk assessments for charge adequacy in frontier exploration wells.

2D and 3D Back-Stripping: Structural Restoration

Two-dimensional back-stripping extends the 1D method along a seismic cross-section by restoring the geometry of each horizon to its interpreted depositional shape, removing compaction effects, and balancing the cross-sectional area to verify that no material has been created or destroyed by the restoration process. Balanced cross-section techniques, originally developed for thrust belt interpretation in the Canadian Rockies and the Appalachians, are combined with the decompaction algorithms of 1D back-stripping to produce chronostratigraphic sections (Wheeler diagrams) that show the temporal and spatial evolution of the basin in a single visualization. These sections reveal accommodation space creation patterns, unconformity development, and the lateral migration of depocenters over time, all of which are critical inputs to sequence stratigraphy analysis and to understanding the distribution of source rocks, reservoir sands, and sealing units in a reservoir characterization model.

Three-dimensional back-stripping using basin-wide grids of well and seismic data is computationally intensive but has become standard practice at major oil companies since the 2000s. Commercial software packages including Petromod (SLB), BasinMod (Zetaware), and Temis (IFP/Beicip-Franlab) implement 3D back-stripping as the first workflow step in any basin model. The 3D approach captures the flexural response of the lithosphere, which varies laterally with effective elastic thickness, and resolves the three-dimensional migration pathways of hydrocarbons from kitchen to trap by tracking the geometry of migration surfaces (top-seal surfaces, unconformities, permeable carrier beds) through time. This capability is particularly important in complex passive margin settings like the Brazilian pre-salt, the West African margins, and the Norwegian Barents Sea, where multiple phases of rifting, salt tectonics, and margin inversion have created structurally complex petroleum systems that cannot be understood with 1D methods alone.