Some limitations of the Rf/φ technique of strain analysis

Finite strain determined from R_f/φ analysis of deformed ellipse distributions cannot be separated into compactional and tectonic parts on the basis of population asymmetry about the orientational mean, as was originally suggested. Furthermore, where there is a ductility contrast between objects and matrix, R_f/φ analysis yields an invalid estimate of the strain suffered by the objects.     

Some limitations of the Rf/φ technique of strain analysis — Discussion

A recent note by De Paor (1980) suggests that finite strain cannot be decomposed into tectonic and compaction components on the basis of asymmetry of the R_fφ data about the mean φ trace. This is true if only the shape of the R_fφ data field is considered, but shown to be incorrect when the distribution of data within this field is considered.     

The strain history of a dalradian slide: Using pebbles with low fluctuations in axis orientation

The strain associated with the Horn Head Slide, a major tectonic break in the Dalradian rocks of NW Ireland, is recorded by pebbles in an adjacent quartzite horizon. Mean X/Y ratios of the deformed pebble shapes in excess of 8.0 are seen closest to the slide and the field of three-dimensional shapes lies along the K = 1 line. The usual methods of separating initial shape ratio (R_i) and tectonic strain ratio (R_s) from the deformed shape ratio (R_f) of ellipsoidal markers are based on measurements of variation in fluctuation (e.g., the (R_f/φ technique). However, due to the high X/Y strains in this situation and since the pebbles initially lay parallel to bedding and to a principal plane of the finite strain ellipsoid, fluctuation is generally very low. Thus, except for the least deformed X/Y data, the R_f/φ technique is inapplicable and other methods are used. For X/Y data with mean (R_f > 4.0: R_s is calculated as the harmonic mean of R_f; maximum R_i values only are obtained from the range of R_f data. For all Y/Z and X/Z data: R_i is calculated from the logarithmic range (ω_log) of R_f; R_s is simply obtained from the geometric mean of R_f modified by R_i. It is concluded from this that a varying prolate tectonic strain (K - 1.5) reaching X/Y values in excess of 8.00 was coaxially superimposed on an initial oblate shape fabric to produce the present field of deformed pebbles in the quartzite near the slide.     

Magnetic susceptibility anisotropy in the Cambrian slate belt of North Wales and correlation with strain

The magnetic susceptibility anisotropy of 275 specimens comprising 38 sites from the Cambrian slate belt in North Wales was measured to determine the magnetic fabric of the slates. The susceptibility ellipsoid is oblate for all sites, and the maximum/intermediate susceptibility plane always coincides with the cleavage plane of the slates which has a Caledonian strike and is nearly vertical. The maximum axes align sub-vertically and the intermediate axes sub-horizontally, trending NE-SW. The minimum susceptibility axes are normal to this foliation plane and coincide with the poles to the slaty cleavage. The orientations of the principal susceptibility axes are found to be in excellent agreement with the orientations of the principal strain directions, determined by X-ray goniometry on one of the samples from almost all of the sites. Correlation of the magnetic susceptibility anisotropy with predicted March strains (March, 1932) shows that the principal magnitudes of susceptibility can be related to those of the strain by: (for i = 1, 2, 3. The orthogonal principal axes), where χ_f and χ₀ are the final and initial susceptibilities along a given axis i and l_f and l_i are final and initial axial dimensions in the same direction i of a principal strain axis. The exponent a for the North Wales slates was found to be 0.145 ± 0.005. Knowledge of such a relationship may permit rapid approximate determinations of a petrofabric in similar rocks from their magnetic fabrics. However, the exponent a will probably have to be recalibrated for each rock type.     

A method for the determination of the initial shapes of deformed xenoliths in granitoids

Measuring strain from deformed xenoliths is problematic due mainly to the large initial shape variations of these markers. A method is described which allows mean initial shape to be determined for a number of xenolith populations by displaying their logarithmic ranges (log R_fmax — log R_fmax) on a Range diagram. The diagram contains a check on validity and allows bad samples to be recognised. Xenolith data from a deformed granite in northwest Ireland is analysed using the method.     

Grain size as an indication of shear strain in brittle fault zones

Shear strain γ in brittle fault zones is related to final and initial grain size parameters, d_f and d_i, respectively, by an expression of the form: where μ_k the coefficient of kinetic friction of crushed rock and σ_n the normal stress across the zone at the time of faulting.Technological literature suggests that ƒ(d_f, d_i) may be given by 10 , where W_i is a material constant. The resulting relationship between shear strain and grain size seems to be compatible with existing experimental data.     

Triaxial stress state deep in orogenic belts: an example from Turkey

A rare metachert pebble containing amphibole grains with microboudin structures in a wide range of orientations provides an opportunity to perform stress analysis in two orthogonal orientations on the foliation surface. The sample was analysed by the microboudin method to infer the triaxial stress state during microboudinage. Stress parameters proportional to the far-field differential stress were determined for sodic amphibole grains in the two orientations. The ratio of the stresses in the two orthogonal orientations (σ₁−σ₂)/(σ₁−σ₃) was calculated to be 0.64, indicating that σ₂ lies closer to the midpoint between σ₁ and σ₃ than to σ₃.     

The importance of the effective intermediate principal stress (σ′2) to fault slip patterns

In normal faulting regimes, the magnitudes and orientations of the maximum and minimum principal compressive stresses may be known with some confidence. However, the magnitude of the intermediate principal compressive stress is generally much more difficult to constrain and is often not considered to be an important factor. In this paper, we show that the slip characteristics of faults and fractures with complex or nonoptimal geometry are highly sensitive to variation or uncertainty in the ambient effective intermediate principal stress (σ′₂). Optimally oriented faults and fractures may be less sensitive to such variations or uncertainties. Slip tendency (T_s) analysis provides a basis for quantifying the effects of uncertainty in the magnitudes and orientations of all principal stresses and in any stress regime, thereby focusing efforts on the most important components of the system. We also show, for a normal faulting stress regime, that the proportion of potential surfaces experiencing high slip tendency (e.g., T_s ≥ 0.6) decreases from a maximum of about 38% where σ′₂ = σ′₃, to a minimum of approximately 14% where σ′₂ is halfway between σ′₃ and σ′₁, and increases to another high of approximately 29% where σ′₂ = σ′₁. This analysis illustrates the influence of the magnitude of σ′₂ on rock mass strength, an observation previously documented by experimental rock deformation studies. Because of the link between fault and fracture slip characteristics and transmissivity in critically stressed rock, this analysis can provide new insights into stress-controlled fault transmissivity.     

Solid bitumen reflectance and Rock-Eval Tmax as maturation indices: an example from the “Nordegg Member”, Western Canada Sedimentary Basin

Marine, organic-rich rock units commonly contain little for vitrinite reflectance (VR₀) measurement, the most commoly used method of assessing thermal maturity. This is true of the Lower Jurassic “Nordegg Member”, a type I/II, sulphur-rich source rock from the Western Canada Sedimentary Basin. This study examines the advantages and pitfalls associated with the use of Rock-Eval T_max and solid bitumen reflectance (BR₀) to determined maturity in the “Nordegg”. Vitrinite reflectance data from Cretaceous coals and known coalification gradients in the study area are used to extrapolate VR₀ values for the “Nordegg”.T_max increases non-linearly with respect to both BR₀ and extrapolated VR₀ values. A sharp increase in the reflectaance of both solid bitumen and vitrinite occurs between T_max 440–450°C, and is coincident with a pronounced decrease in Hydrogen Index values and the loss of solid bitumen and telalginite fluorescence over the same narrow T_max interval. This T_max range is interpreted as the main zone of hydrocarbon generation in the “Nordegg”, and corresponds to extrapolated VR₀ values of 0.55–0.85%. The moderate to high sulphur contents in the kerogen played a significant role in determining the boundaries of the “Nordegg” oil window.A linear relationship between BR₀ and extrapolated VR₀, as proposed elsewhere, is not true for the “Nordegg”. BR₀ increases with respect to extrapolated VR₀ according to Jacob's (1985) formula (VR₀=0.618×(BR₀)+0.40) up to VR₀≈0.72% (BR₀≈0.52%). Beyond this point, BR₀ increases sharply relative to extrapolated VR₀, according to the relatioship VR₀ = 0.277 × (BR₀) + 0.57 (R² = 0.91). The break in the BR₀−VR₀ curve at 0.72%VR₀ is thought to signifiy the peak of hydrocarbon generation and represents a previously unrecognized coalification jump in the solid bitumen analogous to the first coalification jump of liptinites.     

Geometry of the fault zone of the 2003 Ms = 7.5 Chuya earthquake and associated stress fields, Gorny Altai

The co-seismic deformations produced during the September 27, 2003 Chuya earthquake (Ms = 7.5) that affected the Gorny Altai, Russia, are described and discussed along a 30 km long segment. The co-seismic deformations have manifested themselves both in unconsolidated sediments as R- and R′-shears, extension fractures and contraction structures, and in bedrock as the reactivation of preexisting schistosity zones and individual fractures, as well as development of new ruptures and coarse crushing zones. It has been established that the pattern of earthquake ruptures represents a typical fault zone trending NW–SE with a width reaching 4–5 km and a dextral strike–slip kinematics. The initial stress field that produced the whole structural pattern of co-seismic deformations during the Chuya earthquake, is associated with a transcurrent regime with a NNW–SSE, almost N–S, trending of compressional stress axis (σ₁), and a ENE–WSW, almost E–W, trending of tensional stress axis (σ₃). The state of stress in the newly-formed fault zone is relatively uniform. The local stress variations are expressed in insignificant deviation of σ₁ from N–S to NW–SE or NE–SW, in short-term fluctuations of relative stress values in keeping their spatial orientations, or in a local increase of the plunge angle of the σ₁. The geometry of the fault zone associated with the Chuya earthquake has been compared with the mechanical model of fracturing in large continental fault zones with dextral strike–slip kinematics. It is apparent that the observed fracture pattern corresponds to the late disjunctive stage of faulting when the master fault is not fully developed but its segments are already clearly defined. It has been shown that fracturing in widely different rocks follows the common laws of the deformation of solid bodies, even close to the Earth surface, and with high rates of movements.     

Magnitude of σ1, σ2, and σ3 at mid-crustal levels in an orogenic belt: Microboudin method applied to an impure metachert from Turkey

We revised an equation for estimating palaeostress magnitude using the microboudin technique by incorporating the influence of time on the fracture strength of minerals. The equation was used to estimate triaxial palaeostresses from a rare sample of metachert from Turkey that contains microboudinaged, columnar tourmaline grains in a wide range of orientations within the foliation plane. The estimated principal palaeostresses are σ₁ = 605 MPa, σ₂ = 598 MPa, and σ₃ = 597 MPa. As the microboudinage is considered to have occurred immediately before the rock encountered the brittle-plastic transition during exhumation, these stress values correspond to conditions at approximately 18 km depth and 300 °C within a Cretaceous orogenic belt.     

Porphyroclast systems as kinematic indicators

Porphyroclasts of relatively strong minerals in mylonites commonly have an internal monoclinic shape symmetry defined by tails of dynamically recrystallized material. The geometry of a porphyroclast and its tails, called a ‘porhyroclast system’, can serve as a valuable indicator of the sense of vorticity. Porphyroclast systems have been divided into σ- and δ-types on the basis of the geometry of the tails. σ-Types have wedge-shaped recrystallized tails whose median lines lie on opposite sides of a reference plane parallel to the tails and containing the symmetry axis for the system. σ-Types are further subdivided into a σ_a-types, in which the porphyroclast is isolated in a relatively homogeneous matrix, and σ_b-types, in which the porphyroclast system is associated with a shear band foliation in the matrix. δ-Types typically have narrow recrystallized tails whose median lines cross the reference plane adjacent to the porphyroclast. Consequently, embayments of matrix material occur adjacent to the porphyroclasts and the tails display characteristic bends.A porphyroclast system in a mylonite develops when the relatively weak dynamically recrystallized grain aggregate in the porphyroclast mantle changes its shape due to non-coaxial flow in the adjacent matrix. This behaviour has been modelled in shear box experiments. Passive marker lines around rigid cylinders embedded in silicone putty were subjected to simple shear. The experiments were modified to simulate a change in recrystallization rate (R) with respect to rate of deformation (γ) by decreasing the diameter of the rigid cylinder during deformation at variable rates. The ratio R/γ appears to be one of the most important factors in determining which porphyroclast system will develop. At high R/γ values, flow of recrystallized material away from the porphyroclast is continuously appended by the production of new grains and wedge-shaped σ_a-type tails develop. At low R/γ values, relatively few new grains are added to the tails which become thinned and deflected by drag due to the spinning motion of the porphyroclast. In addition, most porphyroclast systems at low shear strains are of σ_a-type or lack monoclinic symmetry, whereas δ-types are only developed at high shear strain values. Complex porphyroclast systems, characterized by two generations of tails, are observed in many of the natural and model shear zones studied and may form due to fluctuating R/γ. Conditions that allow isolated σ_a- and δ-type porphyroclast systems to be used as sense of vorticity indicators are: the systems should have a monoclinic shape symmetry; matrix grain size should be small with respect to porphyroclast size; matrix fabric should be homogeneous; deformation history should be simple, and observations should be made on sections normal to the inferred bulk vorticity vector for the mylonite.     

Changes in bulk properties and molecular compositions within New Zealand Coal Band solvent extracts from early diagenetic to catagenetic maturity levels

The bulk properties and bitumen molecular compositions of a rank-series of 38 humic coals from the New Zealand Coal Band (Cretaceous–Cenozoic) have been analysed to investigate early maturation processes affecting coaly organic matter through diagenesis to moderate catagenesis (Rank(S_r) 0.0–11.8, R_o 0.23–0.81%). The samples comprise a relatively restricted range of vitrinite rich coal types formed largely from higher land plant material under relatively oxic conditions, but with a significant contribution from microbial biomass. With increasing rank, total organic carbon contents show a general increase, whereas moisture and asphaltene contents decrease. Bitumen yields also decrease through the stages of diagenesis and early catagenesis (Rank(S_r) < 9, R_o < 0.55%), indicating partial loss of initial bitumen during early maturation. Thermal generation of hydrocarbons begins slowly at Rank(S_r)  5–6 (R_o  0.40%) as indicated by the constant occurrence and gradual increase of isoprenoids (e.g., pristane and phytane) and hopanoids in their more mature αβ configuration. This early phase of catagenesis, not previously recognised in New Zealand coals, is followed at Rank(S_r)  9 (R_o  0.55%) by the main catagenesis phase characterised by a more rapid increase in the generation of hydrocarbons, including total n-alkanes, isoprenoids and αβ-hopanes. Changes in the maturity of New Zealand coals can be traced by the Carbon Preference Index and several hopane maturity parameters, including 22S/(22S + 22R), αβ/(αβ + βα) and ββ/(αβ + βα + ββ).     

Investigations on the mechanical significance of crenulation cleavage

Measurements of total, incremental and progressive strains associated with the development of small scale crenulation cleavage in some low-grade metamorphic rocks from Australia and Switzerland are applied to a discussion of the mechanical significance of the cleavage.Limits are placed on the amount of incremental and total slip or simple shear possible along the cleavage by the observation that the XY principal plane trace of bulk total crenulation strain coincides within 4° of the crenulation cleavage trace in all cases where this strain has been measured or estimated. The measurements are made on eight specimens using deformed porphyroblasts, crystal fibres in pressure-shadows around pyrite and flattened folds and include deformations with coaxial and non-coaxial histories.Further measurements derived from pressure-shadow fibres (eight specimens) show that the style and orientation of incremental deformation are essentially independent of the crenulation cleavage, except for a limit (43°) to the obliquity of the principal incremental extension axis during a given cleavage episode. The only special deformation related to the cleavage is the coaxial one. An indication of passive cleavage behaviour at high strain is shown by the progressive strain history of one specimen. Evidence for passive rotation of a transected axial plane is shown by another. A model is proposed to account for these observations, especially the conditions necessary for initiation and continued development of a new cleavage fabric.Some further applications of existing strain measurement techniques are described: of the R_f/Ø_f method to heterogeneously superposed tectonic strains and of an improved procedure of t′_α/α flattening analysis.     

The representation and calculation of the deviatoric component of the geological stress tensor

A new diagram for the representation of stress states is proposed and compared with Nadai's stress diagram. The diagram is a graph whose axes are labelled as the differences of the principal stresses (σ₂-σ₃ as ordinate axis, σ₁-σ₂ as the abscissa; where σ₁ > σ₂ > σ₃ are the principal stresses). The design of the plot has been deliberately modelled on that of the ‘log Flinn’ diagram which is used to represent finite strain ellipsoids. The position of the plotted stress state on this diagram depends on the nature of the deviatoric (non-hydrostatic) component of the stress tensor. The distance of the plotted stress from the origin corresponds broadly to the departure of the stress from a hydrostatic state and the parameter R, defined as the gradient of the line joining the plotted state to the origin, expresses the type of symmetry possessed by the stress tensor.It is explained how the diagram can be used to represent calculated palaeostresses and, in particular, how the parameter R can be found directly from some existing methods of stress analysis currently in use.Besides its proposed function to represent the results of such analyses it is felt that the use of the diagram may make clear the essential elements of the definitions of well-known terms used to describe particular stress states (e.g. plane stress, triaxial stress, axial compression, etc.).     

Late Cenozoic stress evolution along the Karasu Valley, SE Turkey

This study defines the Mio-Pliocene to present-day stress regime acting at the northeastern corner of the eastern Mediterranean region along the Karasu Valley (i.e., the Amanos Range), taking in the Antakya, Osmaniye and Kahramanmaras provinces. The inversion slip vectors measured on fault planes and chronologies between striations indicate that the stress regime varied from transpressional initially to transtensional, having consistent NW- and NE-trending σ_Hmax (σ₁) and σ_Hmin (σ₃) axes, respectively; there are significantly different mean stress-ratio (R_m) values however. The older mean stress state is characterized by N151±11°E-trending σ₁ and N59±12°E-trending σ₃ axes, and by a mean arithmetic R_m value of 0.76, indicating that the regional stress regime is transpressional. The younger stress regime is characterized by N154±8°E-trending σ₁ and N243±8°E-trending σ₃ axes, and by a mean arithmetic R_m value of 0.17, indicating a transtensional character for this regional stress regime. The low R values of the stress deviators related to the recent stress state reflect normal-component slips. The earthquake focal mechanism inversions confirm that the younger stress regime continues into the Recent. The inversion identifies a transtensional stress regime representing strike-slip and an extensional stress state with a consistent NE-trending σ_Hmin (σ₃) axis. These stress states are characterized by N66°E and N249°E-trending σ₃ axes, respectively. Both significant regional stress regimes induce left-lateral displacement along the southern part of the East Anatolian Fault (EAF, or Amanos Fault). The temporal change, probably in Quaternary time, within the regional stress regime—from transpression to transtension—resulted from the coeval influences of subduction processes in the west–southwest (i.e., along the Cyprus arc), continental collision in the east, and westward escape of the Anatolian block.     

Determination of stress directions from plagioclase fabrics in high grade deformed rocks (Além Paraíba shear zone, Ribeira fold belt, southeastern Brazil)

A new method to determine stress directions using the preferential orientation of plagioclase mechanical twins has been applied to high-temperature mylonitic rocks from the Além Paraíba shear zone, Ribeira fold belt, southeastern Brazil. We have measured the lattice-preferred orientation of plagioclase grains and calculated the orientation of the stress axes possible for the observed twin orientations. The maximum compressive stress direction (σ₁), determined for all studied samples, is a function of the mechanical twin orientations of a number of distinct plagioclase populations. The σ₁ direction is generally subperpendicular to the (010) plane. The statistical treatment for most of the plagioclase grains examined for each sample shows that σ₁ is almost perpendicular to the foliation plane, suggesting a significant coaxial component in the deformation process of these rocks.     

3-D Mohr circle construction using vein orientation data from Gadag (southern India) – Implications to recognize fluid pressure fluctuation

In this paper orientations of quartz veins from the Archaean age lode-gold bearing region of Gadag (southern India) are used to determine the relative stress and fluid pressure (P_f) conditions by constructing 3-D Mohr circle. Anisotropy of magnetic susceptibility (AMS) analysis of the host massive metabasalt reveals that the magnetic foliation is NW–SE striking, which is related to early NE–SW compression (D1/D2 deformation) that affected the region. The quartz veins have a wide range of orientations, with NW–SE striking veins (steep northeasterly dips) being the most prominent. Vein emplacement is inferred to have taken place under NW–SE compression that is known to have caused late deformation (D3) in the region. It is argued that the NW–SE fabric defined the pre-existing anisotropy and channelized fluid flow during D3. The permeability was initially low, which resulted in high P_f (>σ₂). 3-D Mohr circle analysis indicates that the driving pressure ratio (R′) was 0.94, a condition that favoured fracturing and reactivation of fabric elements (foliations and fractures) having a wide range of orientations. This led to an increase in permeability and fluid flowed (burped) into the fractures. Resulting vein emplacement and sealing of fractures led to a reduction of P_f (<σ₂). It is argued that at this low P_f, NW–SE oriented structures continued to remain susceptible to reactivation and vein emplacement, while fractures of all other orientations were inactive and remained sealed. As a consequence, the study area has a cluster of NW–SE oriented veins. R′ is calculated to be 0.07 from 3-D Mohr circle analysis at low P_f, when fractures with NW–SE orientation only were susceptible to dilation. However, it is envisaged that any emplacement of veins in these fractures would have sealed them, thus reducing the permeability and initiating the next cycle of rise in P_f (>σ₂). Thus, it is concluded that the quartz veins in the Gadag region are a consequence of an interplay between conditions that fluctuated from P_f > σ₂ to P_f < σ₂.     

The role of shear strain in the graphitization of a high-volatile bituminous and an anthracitic coal

A ‘soft’ carbon-based high-volatile bituminous (R_{o max}=0.68%) coal and a ‘hard’ carbon-based Pennsylvania anthracite (R_{o max}=5.27%) were deformed in the steady state at high temperatures and pressures in a series of coaxial and simple shear deformation experiments designed to constrain the role of shear strain and strain energy in the graphitization process. Tests were carried out in a Griggs-t type solid (NaCl) medium apparatus at T=400–900°C, constant displacement rates of 10^-5−10^-6 s⁻¹, at confining pressures of 0.6 GPa (coaxial) or 0.8 and 1.0 GPa (simple shear). Coaxial samples were shortened up to 50%, whereas shear strains up to 4.9 were attained in simple shear tests. Experiments lasted up to 118 h. Deformed, high-volatile bituminous coal was extensively coked and no correlation between strain and R_{o max}, bireflectance or coal texture was observed in any samples. With increasing temperature, R_{o max} and bireflectance increase in highly anisotropic, coarse mosaic units, but remain essentially constant in the fine granular mosaic, which becomes more abundant at higher temperatures. Graphite-like reflectances are observed locally only in highly reactive macerals and in pyrolytic carbon veins. The degree of molecular ordering attained in deformed bituminous coal samples appears to be determined by the heating-pressurization path rather than by subsequent deformation.Graphitization did not occur in coaxially deformed anthracite. Nonetheless, dramatic molecular ordering occurs at T>700°C, with average bireflectance values increasing from 1.68% at 700°C to 6.36% at 900°C. Anisotropy is greatest in zones of high strain at all temperatures. In anthracite samples deformed in simple shear over the 600–900°C range at 1.0 GPa, the average R_{o max} values increase up to 11.9%, whereas average bireflectance values increase up to 10.7%. Bireflectance increases with progressive bedding rotation and, thus, with increasing shear strain. Graphitization occurs in several anthracite samples deformed in simple shear at 900°C. X-ray diffraction and transmission electron microscopy of highly anisotropic material in one sample confirms the presence of graphite with d₀₀₂=0.3363 nm. These data strongly suggest that shear strain, through its tendency to align basic structural units, is the factor responsible for the natural transformation of anthracite to graphite at temperatures far below the 2200°C required in hydrostatic heating experiments at ambient pressure.     

P–T–fO2 controls on a partly inverse chromite-bearing ultramafic intrusive: an evaluation from the Sukinda Massif, India

The chromites from the alpine type ultramafic intrusive of Sukinda, India, display a typical partly inverse spinel form and occur in two distinct zones: Brown Ore Zone (BOZ) and Grey Ore Zone (GOZ). The host ultramafites are mostly altered and are represented by the serpentinite, tremolite-talc(chlorite) schist, talc-serpentine schist and chlorite rock. The less altered variants are dunite, harzburgite and websterite. A dyke of orthopyroxenite runs through the main ultramafic body.The composition of olivine (Fo₉₂), orthopyroxene (En_92–89) and Al₂O₃ contents of the parental liquid (10.40–11.45%) determined from chromites, suggest that the parent melt is of boninitic affinity. The chemical plot of TiO₂ content against cr# of chromites corroborates a boninitic parental melt. The Fe–Mg partitioning in olivine and chromite depicts the temperature for chromitites as 1200 °C. A compositional plot of mg# and cr# suggests crystallization at high pressure conditions, corresponding to the kimberlite xenolith field. From the P–T diagram of pyrolite melting and mineral assemblage, the pressure of crystallization is stipulated to be ≥1.2 GPa. The f_O2 values estimated from Fe³⁺/Cr+Al+Fe³⁺ ratios range from 10^−8.3 to 10^−9.3 for the GOZ and 10^−7.1 to 10^−7.3 for the BOZ. The f_O2 values together with the pressure range suggest crystallization at upper mantle conditions. The heterogeneity in chemical composition and f_O2 conditions for the GOZ and BOZ could be linked to heterogeneity in the upper mantle.