Abstract Raman spectroscopy allows thermal maturation of carbonaceous sediments to be determined. The technique has been employed on metamorphic samples exceeding temperatures of 270°C, but recently has been shown to be effective at lower temperatures. Thermal maturation techniques commonly depend on sample size, have varying efficacies at different temperatures and in different conditions. The underlying processes are not well understood, thus data interpretation can be ambiguous. Here we show the efficacy of Raman as a low-temperature thermal marker in a thrust belt. The Bornes region, in a French Subalpine chain, provides an opportunity to test the technique against published vitrinite reflectance data and thermal modelling for the first time. In doing so we show that Raman is an effective thermal marker to temperatures as low as 75°C, has a small error and is consistent with previous work. The Raman data allow us to postulate the relative thickness of the sedimentary succession across the chain, the timing of thermal maturation and the timing and thickness variations of the over-thrust Prealpine nappe. The work establishes Raman as a low-temperature thermal marker for correlation with other techniques to ensure effective and robust interpretation, that can readily be applied to fold–thrust belts in hydrocarbon provinces.

Abstract John G. Ramsay's book Folding and Fracturing of Rocks was first published in 1967. It set the research agenda for, and the approaches to, modern structural geology. It is recognized as one of the most influential texts in the subject area and was underpinned by significant research publications by the author; its weight added to by later texts and articles. John was also an enthusiastic teacher of structural geology, and many of his tutees have gone on to forge their own careers in structural geology, inspired to do so by John, his teaching and approaches. This Special Publication serves as a golden 50 year celebration of the book, and also of John's life-long work and career as an innovative scientist, inspiring mentor, excellent teacher, and great field structural geologist.

Abstract In 1967, the publication of John Ramsay's book, Folding and Fracturing of Rocks , was a landmark event that affected both of us, inspiring us into careers in structural geology. At that time, one of us was a postgraduate student at Imperial College, London, the other a second-year undergraduate at Manchester University. It is hard to convey, 50 years on, how precious this book was then. The true and lasting value is how this book has changed individuals, and influenced the teaching and research in structural geology, in Britain and around the world. In our view, Folding and Fracturing of Rocks marks the birth of modern structural geology.

Abstract On the 50th anniversary of the publication of the book Folding and Fracturing of Rocks by John G. Ramsay we use an analysis of citations to assess the book's impact on research in the field of structural geology. Five topics dealt with in the book have had special success in terms of the number of citations: (1) superimposed folding; (2) fold classification based on layer thickness variation; (3) simple shear deformation; (4) the R f / φ method; and (5) the flexural-slip mechanism of folding. The great impact of this book can be seen in several examples where the book has inspired new lines of research, such as: progressive deformation, strain analysis, shear zones and folding analysis. Ramsay's great merit in writing his book was to have wisely known to apply the quantitative methodology of continuum mechanics to the analysis of the geological structures. We conclude that Folding and Fracturing of Rocks is the most influential textbook in the field of structural geology research.

Abstract The margins to evolving orogenic belts experience near layer-parallel contraction that can evolve into fold–thrust belts. Developing cross-section-scale understanding of these systems necessitates structural interpretation. However, over the past several decades a false distinction has arisen between some forms of so-called fault-related folding and buckle folding. We investigate the origins of this confusion and seek to develop unified approaches for interpreting fold–thrust belts that incorporate deformation arising both from the amplification of buckling instabilities and from localized shear failures (thrust faults). Discussions are illustrated using short case studies from the Bolivian Subandean chain (Incahuasi anticline), the Canadian Cordillera (Livingstone anticlinorium) and Subalpine chains of France and Switzerland. Only fault–bend folding is purely fault-related and other forms, such as fault-propagation and detachment folds, all involve components of buckling. Better integration of understanding of buckling processes, the geometries and structural evolutions that they generate may help to understand how deformation is distributed within fold–thrust belts. It may also reduce the current biases engendered by adopting a narrow range of idealized geometries when constructing cross-sections and evaluating structural evolution in these systems.

Abstract Viscous folding in mechanically heterogeneous layers is modelled numerically in two dimensions for linear and power-law viscous fluids. Viscosity heterogeneities are expressed as circular-shaped variations of the effective viscosity inside and outside the layers. The layers are initially perfectly flat and are shortened in the layer-parallel direction. The viscosity heterogeneities cause a perturbation of the velocity field from the applied bulk pure shear, which perturb geometrically the initially flat-layer interfaces from the first numerical time step. This geometrical perturbation triggers interfacial instabilities, resulting in high-amplitude folding. We compare simulations with heterogeneities with corresponding simulations in which the heterogeneities are removed after the first time step, and, hence, only the initial small geometrical perturbations control wavelength selection and high-amplitude folding. Results for folding in heterogeneous and homogeneous layers are similar, showing that viscosity heterogeneities have a minor to moderate impact on fold wavelength selection and high-amplitude folding. Our results indicate that the interfacial instability is the controlling process for the generation of buckle folds in heterogeneous rock layers. Therefore, existing analytical and numerical solutions for folding in homogeneous layers, in which folding was triggered by geometrical perturbations, are useful and applicable to study folding in natural, heterogeneous rock layers.

Abstract Geometrical heterogeneities along layer interfaces play a key role in determining the geometries of folds developed during shortening of competent layers. We present a series of numerical simulations to investigate the influence of initial sinusoidal perturbations on the folding of single layers. Models consist of a competent viscous single layer embedded in a softer matrix, with the layer orientated parallel to the shortening direction. We first generalize the wide spectrum of sinusoidal perturbations accounting for asymmetries along and across a competent single layer, using two parameters: transversal asymmetry ( A ′) and longitudinal asymmetry ( φ ). These two parameters allow the transition between classical fold shapes and pinch and swell geometries to be studied. The parameter A ′ describes the development of fold hinges with different geometries between the upper and lower layer interfaces, and abnormal curvatures between the outer and inner arcs of fold hinges. The parameter φ induces a strong polarity on the folds, with a systematic preferred orientation of the pinch and swell regions of the layer, even if there is no shear component parallel to the layer. Our results demonstrate the importance of structural inheritance on the resulting fold geometries, and suggest that caution must be taken when using certain types of asymmetrical folds as strain markers and kinematic indicators.

Abstract The detailed analysis of folding in rocks was in part pioneered by John Ramsay, and resulted in a range of techniques and criteria to define folds. Although folding of unlithified or ‘soft’ sediments is typically assumed to produce similar geometries to those in ‘hard rocks’, there has to date been little detailed analysis of such folds. The aim of this paper is therefore to investigate folds developed during soft-sediment deformation (SSD) by applying techniques established for the analysis of tectonic folds during hard-rock deformation (HRD). We use the Late Pleistocene Lisan Formation exposed around the Dead Sea as our case study, as the laminated lake sediments record intricacies of fold detail generated during seismically triggered slumping of mass transport deposits (MTDs) towards the depocentre of the basin. While it is frequently assumed that folds created during SSD are chaotic and form disharmonic structures, we provide analyses that show harmonic fold trains may form during slumping, although larger upright folds cannot be traced for significant distances and are more typically disharmonic. Our analysis also reveals a range of fold styles, with more competent detrital-rich layers displaying buckles (Class 1B), as well as upright Class 1A folds marked by thickened limbs. Class 1A buckle folds are generally considered to be created by flattening that overprints folds with an original Class 1B geometry. As thickened fold limbs are truncated by overlying erosive surfaces, the vertical flattening is considered to have occurred during the slump event. Different fold shapes may partially reflect variable flattening, depending on the original orientation of upright or recumbent folds, together with continued downslope-directed simple-shear deformation that modifies the fold geometry. Analysis of fold wavelength, amplitude and bed thickness allows us to plot strain contour maps, and indicates that beds defining slump folds display viscosity contrasts in the range of 50–250, which are similar to values estimated from folds created during HRD in metamorphic rocks. A range of refold patterns, similar to those established by John Ramsay in metamorphic rocks, are observed within slumps, and are truncated by the overlying sediments, indicating that they formed during a single progressive slump event rather than distinct ‘episodes’ of superimposed deformation. This study confirms that techniques developed for the analysis of folds created during HRD are equally applicable to those formed during SSD, and that resulting folds are generally indistinguishable from one another. Extreme caution should therefore be exercised when interpreting the origin of folds in the rock record where the palaeogeographical and tectonic contexts become increasingly uncertain, thereby leading to potential misidentification of folds created during SSD.

Abstract Classic fold-thrust structures within Carboniferous-age strata at Broadhaven, SW Wales are well-known for their excellent preservation of Variscan deformation. These sites have been important for conceptual model generation of the link between faulting and folding, and are often cited as exemplars of fault-propagation folds following work by Williams & Chapman. Here we employ the virtual outcrop method to digitally map and measure, in detail, the classic Den’s Door outcrop. 3D reconstruction of the site by digital photogrammetry allows us to extract high-density structural measurements, reassess the existing model of structural development for the outcrop, and re-evaluate the link between faulting and folding. We find that digital mapping highlights greater variability in fault displacement and bed thicknesses than previously documented. Fracture analysis shows that fracture intensity is closely linked to structural position and bed-thickness variability, and fracture orientations record the existence of discrete mechanical boundaries through the structure. These results record complex patterns of strain distribution and multi-phase deformation. Evidence for temporal and spatial variability in strain distribution suggests that multiple kinematic and non-kinematic models of deformation are required to faithfully describe even this apparently simple structure. This calls into question the applicability of end-member models of fault-related folding, particularly for multilayered stratigraphy.

Abstract Two models to explain the progressive deformation of syntectonic quartz veins are derived from conventional theories for simple and pure shears. The simple-shear model is based on reorientation and changes in length of linear vein elements, and predicts initial orientations of veins for imposed shear strains, elongations and strain ratios. The pure-shear model considers changes in length of lines variably orientated relative to the maximum compression direction, and yields estimates of elongation strains and strain ratios. Expectations of both models are different, as illustrated by analysis of quartz veins from the Rhoscolyn Anticline, Anglesey, NW Wales. The simple-shear model recognizes three distinct initial orientations, which predict different strains across the fold; the pure-shear model suggests veins were initially sub-parallel to the principal compression direction and predicts effectively constant strains across the fold. In addition, both models predict different patterns of fold vergence: for simple shear, vergence depends on magnitude and direction of shearing and may exhibit complex patterns; for pure shear, vergence patterns are predicted to be essentially constant. In general, the predictions of either model are critically dependent on the origin of the veins, particularly relative to the formation of the Rhoscolyn Anticline.

Abstract Cleavage refraction angles are used to estimate effective viscosity contrasts between interlayered quartzites and phyllites within the Baraboo Syncline, Wisconsin, USA. Both types of layers contain two major phases, quartz and pyrophyllite, with minor amounts of hematite. Quartz (with minor hematite) behaves as the strong phase and pyrophyllite acts as the weak phase. Cleavage refraction directly relates to mineralogy with a linear relationship between bedding/cleavage angle and strong-phase concentration. Mineralogy exerts first-order control over effective viscosity contrasts, which are generally small, in most cases <10. Effective viscosity contrasts are consistent across the fold, so are likely not to be highly strain dependent and indicate approximate linear viscous rheology. Microstructures suggest deformation was dominated by dislocation creep in layers with high quartz concentrations and diffusive mass transfer in layers with lower quartz concentrations, and that the transition of the deformation mechanism is gradual. Thus, the rheological flow laws at the small scale may not reflect the bulk flow law at the large scale over the span of the deformation. Effective viscosity contrasts allow an evaluation of samples compared to theoretical two-phase mixtures. The analysed samples most closely resemble the Reuss bound of two-phase mixtures, regardless of the mineralogy.

Abstract One of the main aims in the field of structural geology is the identification and quantification of deformation or strain. This pursuit has occupied geologists since the 1800s, but has evolved dramatically since those early studies. The quantification of strain in sedimentary lithologies was initially restricted to lithologies of known initial shape, such as fossils or reduction spots. In 1967, Ramsay presented a series of methods and calculations, which allowed populations of clasts to be used as strain markers. These methods acted as a foundation for modern strain analysis, and have influenced thousands of studies. This review highlights the significance of Ramsay's contribution to modern strain analysis. We outline the advances in the field over the 50 years since publication of Folding and Fracturing of Rocks , review the existing limitations of strain analysis methods and look to future developments.

Abstract Geological strain analysis of sedimentary rocks is commonly carried out using clast-based techniques. In the absence of valid strain markers, it can be difficult to identify the presence of an early tectonic fabric development and resulting layer parallel shortening (LPS). In order to identify early LPS, we carried out anisotropy of magnetic susceptibility (AMS) analyses on Mississippian limestones from the Sawtooth Range of Montana. The Sawtooth Range is an arcuate zone of north-trending, closely spaced, west-dipping, imbricate thrust sheets that place Mississippian Madison Group carbonates above Cretaceous shales and sandstones. This structural regime is part of the cordilleran mountain belt of North America, which resulted from accretion of allochthonous terrains to the western edge of the North American continent. Although the region has a general east–west increase in thrust displacement and related brittle deformation, a similar trend in penetrative deformation or the distribution of tectonic fabrics is not observed in the field or in the AMS results. The range of magnetic fabrics identified in each thrust sheet ranges from bedding controlled depositional fabrics to tectonic fabrics at a high angle to bedding.

Abstract Where primary porosity and permeability of a rock are unfavourable for hydrocarbon production, fractures can improve reservoir potential by enhancing permeability. Higher fracture intensity may create a better-connected fracture network, improving fractured-reservoir quality. Investigations into the controls on fracture intensity commonly conclude that either structural or lithological factors have the greatest influence on fracture abundance. We use the Swift Reservoir Anticline in northwestern Montana to investigate how fracture intensity varies throughout the structure and determine that although structural factors do influence fracture intensity, lithology is the main control at outcrop. The Swift Reservoir Anticline exposes bedding surfaces of the Mississippian Castle Reef Formation dolomite. Field data indicates that fracture intensity is highest in the fold forelimb, decreasing into the backlimb except in outcrops of coarse dolomite where fracture intensity is low, regardless of structural position. Field fracture intensity correlates with whole-rock quartz, kaolinite and porosity percentages. We suggest porosity and composition influence bulk-rock mechanical properties, which, in turn, control the fracture intensity at outcrop. Fracture intensity has a stronger relationship with lithological than structural factors, therefore we suggest that the key to predicting fracture intensity in the subsurface here is understanding how lithology varies spatially.

Abstract The Teton anticline and adjacent structures, in the Sawtooth Range, Montana, USA, are fractured in such a way that may be taken as a model for fractures propagating during buckle folding. However, advances in understanding both the process of folding in forelands and the evolution of fracture patterns found within these folds suggest that it is time to reinterpret the nexus between fracturing and folding within these classic structures. With the benefit of seismic lines, the Teton anticline is best described as a fault-propagation fold. Joint propagation initiated with the formation of two major sets whose orientation is controlled by pre-folding, regional stresses. Two more joint sets propagated in local stress fields, developed in response to anticline growth. Some early joints were reactivated as wrench faults during amplification and tightening of the anticlines. The fracture sets identified are consistent with: (a) propagation in a regional stress field, which may be related to stretching in the Sawtooth Range orocline; and (b) tangential longitudinal strain of the backlimb and forcing or trishear of the forelimb during anticline development. Thus, we suggest that fracture networks across folded structures should be interpreted and characterized in the light of the geological history of the entire system.

Abstract A series of analogue models were run to investigate oblique inversion of pre-existing grabens when overprinted by later shortening and the effect of these grabens on development of contractional structures. Obliquity angle ( α ) defining the initial trend of pre-existing grabens relative to the shortening direction, was systematically changed from 0°, 10°, 20°, 30°, 40°, 50°, 65° and 90°. Different structural styles are shown in different models and also in sections cutting across different parts of the models. Model results show that existence of multi-grabens enhances lateral discontinuity of overprinted thrusts in map view. With increasing the obliquity angle, more and longer lateral ramps developed sub-parallel to the graben trends. The pre-existing grabens were apparently rotated from their initial trends during shortening. Some of the normal faults bounding the grabens were partially inverted and resulted in bulging of the syn- and post-rift graben fill sediments. Most normal faults were displaced and rotated by thrusting, and provided relatively weak zones for propagation of thrusts. By comparing with observations from Qingxi graben in western China and from the SW Taiwan fold-and-thrust belt, where oblique inversion occurred, model results can be used to interpret unclear relationships between thrusts and pre-existing extensional structures during superimposed deformation.

Abstract We conducted 39 Ar diffusion experiments using potassium feldspar from the South Cyclades Shear Zone on Ios, in the Cyclades, Aegean Sea, Greece. Irradiated samples were step-heated in an ultra-high-vacuum resistance furnace attached to a mass spectrometer, thus also allowing 40 Ar/ 39 Ar geochronology. Conjoint inversion of these datasets allowed estimation of the relevant diffusion parameters, which were then used to forward model the effect of arbitrary temperature–time histories. Simulations used Monte Carlo methods to improve approximations to the observed age spectra. Two periods of rapid cooling could be inferred. The South Cyclades Shear Zone commenced operation during or shortly after the Eocene–Oligocene transition. Episodes of south-directed movement continued into Early Miocene time, however, with the footwall still hot enough to cause biotite ± garnet metamorphic mineral growth at the base of the overlying, already substantially exhumed, eclogite–blueschist unit. Since its footwall continued to cool, the South Cyclades Shear Zone was an extensional shear zone during both episodes of its operation.

Abstract Global geophysical observations constrain all theories of terrestrial dynamics. We jointly interpret EGM2008 gravity, RET2014 topography and the Global Centroid Moment Tensor database from a structural point of view. We hypothesize that lateral variations of gravity and topography reflect the scale-dependent competence of rocks. We compare the spectral and spatial characteristics of the observed fields with structural predictions from the mechanics of differential grade-2 (DG-2) materials. The results indicate that these viscoelastic materials are a powerful tool for exploring dynamic processes in the Earth. We demonstrate that the known spectral range of Earth's gravity and topography can be explained by the folding, shear banding, faulting and differentiation of the crust, lithosphere and mantle. We show that the low-amplitude long-wavelength bias apparent in the disturbance field can be explained by perturbations to Earth's overall ellipsoidal shape, induced by internal slab loading of the mantle. We find by examining the directional isotropy of the data that the zonal energy in Earth's gravity disturbance is maximized about an axis coincident with the shape-perturbation minimum. The symmetry of tectonic features about this axis, extending from eastern Borneo to Brazil, and its coincidence with the equator suggest the coupling of current plate motions to true polar wander.