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 at each stage in its development. First formalized by Sclater and Christie in 1980 using borehole data from the North Sea, the method separates the total observed subsidence at a given location into two independent components: the sediment-loading component, which results from the weight of deposited sediment deflecting the lithosphere downward, and the tectonic subsidence component, which reflects genuine crustal extension, thermal cooling, or other tectonic processes unrelated to sediment load. Because the sediment-loading component can be calculated from the known density and thickness of the sedimentary column using Airy or flexural isostasy, subtracting it from the total subsidence leaves the tectonic subsidence curve, which is the geological signal of greatest interest for petroleum system analysis. McKenzie's 1978 lithospheric stretching model relates tectonic subsidence directly to the stretching factor beta (the ratio of post-rift to pre-rift crustal thickness), allowing back-stripping results to constrain when and how severely the crust was extended, which in turn controls the depth and timing of burial heating for source rocks, the generation of petroleum, and the geometry of accommodation space for reservoir deposition. In the Western Canada Sedimentary Basin, back-stripping applied to Devonian and Mississippian strata has elucidated the rift and thermal-subsidence history of the Peace River Arch, the Elk Point Group evaporite basin, and the ancestral Rocky Mountain foredeep, providing calibration constraints for basin models used to predict source rock maturity windows and to locate paleo-structural traps for conventional oil and gas exploration in the Deep Basin and the Foothills.
Key Takeaways
- Sequential layer removal and decompaction correction: The back-stripping procedure begins with the present-day stratigraphic column at a well or seismic location and works backward through time, one formation at a time. The youngest layer is removed first, and the remaining column is restored to the thickness it would have had before being compacted by the overlying layer's weight. This decompaction step uses a porosity-depth relationship of the form phi = phi0 multiplied by exp(-z/c), where phi0 is the depositional surface porosity (typically 0.50 to 0.60 for shales and 0.40 to 0.50 for sands), z is depth in metres, and c is a lithology-specific compaction constant (approximately 1,500 to 3,000 m for shales, 4,000 to 6,000 m for sands). Decompaction restores the true original thickness of each layer as it was at the time of deposition, which is greater than the present compacted thickness by a factor that depends on how deeply the formation was buried. Failure to apply decompaction correction causes systematic underestimation of paleo-water depths and paleo-subsidence rates, particularly in thick shale-dominated sequences where compaction factors can reach 2.5 or greater from surface to maximum burial depth.
- Paleobathymetry correction: To isolate the subsidence of the basin floor itself, each back-stripped horizon must be corrected for the water depth at the time of its deposition, which is known as paleobathymetry. If a carbonate reef complex was deposited in 50 m of water, the apparent subsidence of that horizon below present sea level is overstated by 50 m unless the paleobathymetry correction is applied. Paleobathymetry is reconstructed from multiple proxies: benthic foraminiferal assemblages whose depth ranges are constrained by modern analogues (forams are the most precise for Mesozoic and Cenozoic strata), ichnofacies (trace fossil assemblages that track energy and oxygenation conditions at the seafloor and serve as water-depth indicators in Paleozoic sequences), and seismic-facies geometries such as reef margins, delta clinoforms, and turbidite fan architectures whose water-depth implications are constrained by stratigraphic modelling. In the Devonian Elk Point Basin of the WCSB, paleobathymetry corrections are critical because the basin evolved from a shallow evaporitic platform (0 to 30 m) to an open-marine carbonate ramp (50 to 150 m) within a few million years, and applying incorrect paleobathymetry would create apparent subsidence pulses that do not reflect actual crustal tectonics.
- Isostatic correction: Airy versus flexural isostasy: Isostasy is the principle that the lithosphere floats in gravitational equilibrium on the denser asthenosphere, and that adding or removing sediment load causes the lithosphere to deflect downward or rebound upward. In Airy isostasy, the lithosphere is assumed to respond locally and instantaneously to load changes, so the deflection at any point depends only on the sediment column directly above that point. The Airy correction for sediment loading is W = (rho_s / (rho_m - rho_s)) multiplied by S, where rho_s is the average sediment density, rho_m is the mantle density (approximately 3,300 kg/m3), and S is the sediment thickness. For average sediment densities of 2,300 kg/m3, this gives a subsidence amplification factor of approximately 0.30: for every 1,000 m of sediment deposited, the basin floor subsides an additional 300 m due to sediment-induced isostatic loading alone. Flexural isostasy uses a more realistic elastic-plate model of the lithosphere, in which the lithospheric rigidity (characterised by the elastic thickness Te) spreads the isostatic response laterally over a flexural wavelength of tens to hundreds of kilometres; this approach is required in foreland basins such as the WCSB foredeep where loads from the rising Rocky Mountain thrust belt cause regional depression of the basin floor well beyond the load footprint.
- Tectonic subsidence and the McKenzie stretching model: After removing the sediment loading and paleobathymetry contributions, the residual subsidence curve at a given well location represents pure tectonic subsidence: the component driven by crustal extension, lithospheric thinning, or thermal contraction. McKenzie's 1978 model predicts that tectonic subsidence has two phases following a stretching event: a rapid syn-rift subsidence phase driven by crustal thinning (which thins the buoyant continental crust and allows the denser mantle to take its place, lowering the surface), followed by a prolonged thermal subsidence phase driven by the cooling and contraction of the thermally elevated lithosphere back toward its pre-rift geotherm. The shape of the tectonic subsidence curve, plotted against time, distinguishes these two phases and constrains the stretching factor beta. A beta of 1.5 (crust thinned by one-third) produces a syn-rift subsidence pulse of several hundred metres over 20 to 30 million years, followed by a slow exponential thermal subsidence decay over the next 50 to 150 million years. Back-stripping results that match a McKenzie curve for a known beta value provide independent confirmation of the extensional history inferred from seismic reflection geometries and well log correlations.
- Petroleum system applications: source rock maturity and trap timing: The primary commercial application of back-stripping in petroleum exploration is calibrating basin models used to compute source rock burial depth and temperature through geological time. Because petroleum generation requires sustained heating above the oil-window threshold (approximately 100 to 120 degrees Celsius at geological timescales), knowing when a source rock first reached that temperature and for how long it remained there determines whether a given prospect is oil-charged, gas-charged, or post-mature. Back-stripping provides the burial history that, combined with a paleo-geothermal gradient model calibrated to present-day heat flow measurements and apatite fission track data, gives the temperature-time integral required for maturity calculation. In the Devonian source rocks of northern Alberta, back-stripping studies have shown that the Duvernay Formation reached oil-window maturity (Ro approximately 0.6 to 0.9 percent) during the Mississippian, peaked at gas-condensate maturity (Ro approximately 1.2 to 1.5 percent) during maximum Cretaceous burial, and has been cooling since Laramide erosion removed 1 to 3 km of overburden, a history that explains the present-day condensate-rich nature of the Kaybob South Duvernay play.
Computational Workflow and Data Requirements
A back-stripping analysis begins with the construction of an age-depth table at one or more well locations, assigning a depositional age (top and base in millions of years) and a lithology class to each mapped stratigraphic unit. Accurate chronostratigraphy is the single most important input: errors in formation ages propagate directly into subsidence-rate calculations, and a 5-million-year error on a Devonian horizon can shift the inferred timing of a rift event by one-third of its duration. In the WCSB, Devonian and Mississippian chronostratigraphy is constrained by conodont biozonation supplemented by strontium isotope stratigraphy in carbonate formations and ash bed radiometric dating where Cretaceous tephra layers intersect target stratigraphy. The lithology classification determines the porosity-depth parameters used in decompaction, and using shale parameters for a clean sand or carbonate interval will produce significant errors in restored thickness.
Back-stripping software tools including BasinMod (Platte River Associates), PetroMod (Schlumberger/SLB), Trinity (Ikon Science), and open-source Python implementations (BackPy, pyBackStrip) automate the layer-by-layer decompaction and isostasy calculations, produce tectonic subsidence curves with uncertainty envelopes, and interface with the heat-flow models needed for maturity simulation. Most commercial implementations allow the user to choose between Airy and flexural isostasy and to specify an elastic thickness Te for flexural calculations; for the WCSB foredeep, Te values of 60 to 100 km are appropriate based on the observed wavelength of flexural subsidence documented in seismic cross-sections. Sensitivity analysis on the key inputs (paleobathymetry range, compaction constants, isostasy model) is standard practice and produces a family of tectonic subsidence curves that bracket the uncertainty range, which is then propagated through the maturity model to generate probabilistic maturity estimates rather than single-value predictions.
1D back-stripping at a single well location is the simplest application and is appropriate for preliminary evaluation of a new play area. 2D back-stripping along seismic transects adds the ability to map the lateral variation in subsidence history and to identify basement fault offsets, tilted fault blocks, or differential compaction features that create or destroy structural closure through time. 3D back-stripping, applied to a dense well grid or to horizon picks from a 3D seismic survey, generates paleogeographic maps of basin floor depth and sediment accumulation rate at successive time slices, which are used to reconstruct the sediment provenance systems, identify paleo-highs that focused carbonate reef growth, and pinpoint the fairways where source rocks accumulated the thickest organic-rich intervals. In the Montney play of northeastern British Columbia, 3D back-stripping has been used to reconstruct the Triassic paleogeography that controlled Montney sub-member thickness distribution and thereby the distribution of the highest-porosity Montney C and D intervals that are the primary target for horizontal drilling.
Integration of back-stripping results with present-day observations provides a critical calibration check. Apatite fission track thermochronology (AFTT) and vitrinite reflectance (Ro) measurements in wells provide independent temperature-time constraints that must be consistent with the burial history predicted by back-stripping; if the back-stripped burial depth implies a maximum temperature of 140 degrees Celsius but Ro data indicate only 110 degrees Celsius was reached, either the paleo-geothermal gradient was lower than assumed or the back-stripped burial depth is overestimated. This calibration loop, iterating between the back-stripping burial history, the heat-flow model, and the geochemical maturity data, is the core workflow of quantitative petroleum systems analysis and is the foundation for the resource estimates and play-fairway maps that justify exploration capital allocation in the WCSB and globally.