Statistical modelling of annealing kinetics of fission tracks in zircon; Reassessment of laboratory experiments

We remodel the annealing kinetics of fission tracks (FTs) in zircon to fit data from 43 laboratory experiments in three earlier publications using spontaneous FT lengths in zircon from the Nisatai Dacite. Heating temperature and time are in the ranges of 350–912 °C and 10^− 3–10⁴ h, respectively. The data sets illustrate a transition in annealing process at a certain temperature zone. In order to express the annealing behavior involving the possible transitional process with improved goodness of fit, two types of models are constructed; the hybrid-linear model that consists of the fanning-linear model at high temperature and the parallel-linear model at low temperature, connected with a transitional temperature zone, and the parallel-curvilinear model that approximates smoothly the transition in the annealing process. The validity of the extrapolation of the model to geological time scale is confirmed by comparison with some deep borehole core data. The length estimates by the hybrid-linear model agree well with the borehole data, whilst the estimates by the parallel-curvilinear model show a significant disagreement with highly annealed borehole samples although this model provides good agreement with the experimental data. For the hybrid-linear model, estimates of the partial annealing zone corresponding to isothermal heating for 10⁶ and 10⁷ years are 281–352 °C and 262–330 °C where the top and bottom of the partial annealing zone are defined by the mean length reduction ratio of 0.4 and 0.8, respectively.     

磷灰石裂变径迹退火影响因素研究进展

磷灰石裂变径迹退火是一个繁杂的化学动力学过程,清楚地了解其退火的影响因素对于该技术的应用十分重要。文章概述了磷灰石裂变径迹退火动力学模型的发展史,并结合以往对其退火影响因素的研究,将磷灰石裂变径迹退火影响因素分为自身和外部环境两方面：自身影响因素包括化学成分、晶体结构、径迹长度与半径、径迹与结晶轴的方位关系,其中化学成分对退火起到主导作用;外部环境影响中,温度是主导因素,压力和蚀刻条件的改变也会影响退火。研究成果有利于完善磷灰石裂变径迹的实验室退火模型,提高其作为热历史记录器的精度。     

Modelling the formation of fission tracks in apatite minerals using molecular dynamics simulations

We introduce a simple method to simulate the “ion explosion spike” mechanism of fission track formation within the framework of classical molecular dynamics. The method is applied to six apatite compositions and the resulting tracks are compared with each other as well as with the damage produced by another mechanism—the “Displacement spike”. In contrast to experimentally observed tracks, the radii of simulated tracks are not dependent on their direction in the crystal. Since the simulations model accurately the elastic response of apatites, this suggests that the experimentally observed difference in track radii for tracks along different crystal directions is not entirely caused by anisotropy in the elasticity of apatite. We suggest that anisotropy in the interactions between the electric fields of fission fragments and the crystal ions is a major factor in the final radii of fission tracks. In fluorapatite, the simulations also reveal the formation of small clusters of fluorite-like material in the core of the fission track, a phenomenon which has yet to be confirmed experimentally.

Thermal annealing characteristics of fission tracks in natural zircons of different ages

Fission track (FT) thermochronometry using zircon has widely been applied to the resolution of a variety of geologic problems, for which the understanding of FT annealing behaviour is essential. Thermal annealing experiments were conducted on FTs in natural zircons having different ages (ranging from ~0.6 to ~70 Ma) and radiation damage levels. We measured horizontal confined track lengths on nine zircon concentrates separated from rapidly cooled volcanic rocks, after 1 hr annealing at 400–700°C. As the annealing temperature increases, the observed tracks show a consistent and systematic length reduction in all samples, and the mean track lengths are hardly distinguishable among the nine samples for the same annealing step. Our results suggest that the thermal annealing characteristics at laboratory time‐scale are concordant among the zircons, regardless of their ages, and that identical annealing kinetics may work for Late Mesozoic to Cenozoic zircons.     

Denudation of the Côte d'Ivoire-Ghana transform continental margin from apatite fission tracks

Apatite fission track analysis of samples from the shoulder (marginal ridge) of the Côte d'Ivoire-Ghana transform continental margin reveal a cooling of the margin between 85 and 65 Ma for the central and eastern parts of the ridge. All samples were heated in situ during sedimentary burial with a temperature >120 °C, except for two samples located in the eastern part which were heated between 105 and 120 °C. For the first time, age/depth diagram along a transform margin shows a shape involving erosion starting at the bottom of the continental slope, then stepping backwards towards the edge of the slope. This retrogressive erosion can result from the deepening of the lithospheric plate sliding along the transform margin, from thick continental crust to thin continental crust, and finally to oceanic crust. This process could be at the origin of the shoulder uplift by flexural response to the important crustal discharge (>2 km).     

A Bayesian approach to calibrating apatite fission track annealing models for laboratory and geological timescales

We present a new approach for modelling annealing of fission tracks in apatite, aiming to address various problems with existing models. We cast the model in a fully Bayesian context, which allows us explicitly to deal with data and parameter uncertainties and correlations, and also to deal with the predictive uncertainties. We focus on a well-known annealing algorithm [Laslett, G.M., Green, P.F., Duddy, I.R., Gleadow. A.J.W., 1987. Thermal annealing of fission tracks in apatite. 2. A quantitative-analysis. Chem. Geol., 65 (1), 1-13], and build a hierachical Bayesian model to incorporate both laboratory and geological timescale data as direct constraints. Relative to the original model calibration, we find a better (in terms of likelihood) model conditioned just on the reported laboratory data. We then include the uncertainty on the temperatures recorded during the laboratory annealing experiments. We again find a better model, but the predictive uncertainty when extrapolated to geological timescales is increased due to the uncertainty on the laboratory temperatures. Finally, we explictly include a data set [Vrolijk, P., Donelick, R.A., Quenq, J., Cloos. M., 1992. Testing models of fission track annealing in apatite in a simple thermal setting: site 800, leg 129. In: Larson, R., Lancelet, Y. (Eds.), Proceedings of the Ocean Drilling Program, Scientific Results, vol. 129, pp. 169-176] which provides low-temperature geological timescale constraints for the model calibration. When combined with the laboratory data, we find a model which satisfies both the low-temperature and high-temperature geological timescale benchmarks, although the fit to the original laboratory data is degraded. However, when extrapolated to geological timescales, this combined model significantly reduces the well-known rapid recent cooling artifact found in many published thermal models for geological samples.     

Confined fission track lengths in apatite: a diagnostic tool for thermal history analysis

Fission-track ages in apatite are generally accepted as giving a measure of the time over which a sample has been exposed to temperatures below approximately 100° C. A compilation of the lengths of confined fission tracks in a wide variety of apatites from different geological environments has shown that the distribution of confined track lengths can provide unique thermal history information in the temperature range below about 150° C over times of the order of 10⁶ to 10⁹ years. The distribution of confined lengths of freshly produced induced tracks is characterised by a narrow, symmetrical distribution with a mean length of around 16.3 m and a standard deviation of the distribution of approximately 0.9 m. In volcanic and related rocks which have cooled very rapidly, and never been reheated above about 50° C, the distribution is also narrow and symmetric, but with a shorter mean of 14.5 to 15 m, and a standard deviation of the distribution of approximately 1.0 m. In granitic basement terrains which are thought never to have been significantly disturbed thermally, since their original post-emplacement cooling, the distribution becomes negatively skewed, with a mean around 12 or 13 m and a standard deviation between 1.2 and 2 m.This distribution is thought to characterise slow continuous cooling from temperatures in excess of 120° C, to ambient surface temperatures. More complex thermal histories produce correspondingly complex distributions of confined tracks. The continuous production of tracks through time, coupled with the fact that the length of each track shrinks to a value characteristic of the maximum temperature it has experienced, gives a final length distribution which directly reflects the nature of the variation of temperature with time. Most distinctive of the myriad possible forms of the final distribution are the bimodal distributions, which give clear evidence of a two-stage history, including high and low temperature phases. The study of confined length distributions therefore offers invaluable evidence on the meaning of any fission-track age, and bears the potential of providing rigorous constraints on thermal history in the temperature regime below about 150° C. The results of this study strongly suggest that any apatite fission-track age determination should be supported by a confined track length distribution.     

Annealing studies of fission tracks in allanite

The annealing characteristics of fission tracks in allanite have been studied. All the tracks in the mineral fade at 720° C for one hour heating. The activation energy of allanite ranges from 1.4 to 3.4 eV.     

The influence of artificial radiation damage and thermal annealing on helium diffusion kinetics in apatite

Recent work [Shuster D. L., Flowers R. M. and Farley K. A. (2006) The influence of natural radiation damage on helium diffusion kinetics in apatite. Earth Planet. Sci. Lett.249(3-4), 148-161] revealing a correlation between radiogenic ⁴He concentration and He diffusivity in natural apatites suggests that helium migration is retarded by radiation-induced damage to the crystal structure. If so, the He diffusion kinetics of an apatite is an evolving function of time and the effective uranium concentration in a cooling sample, a fact which must be considered when interpreting apatite (U-Th)/He ages. Here we report the results of experiments designed to investigate and quantify this phenomenon by determining He diffusivities in apatites after systematically adding or removing radiation damage.Radiation damage was added to a suite of synthetic and natural apatites by exposure to between 1 and 100 h of neutron irradiation in a nuclear reactor. The samples were then irradiated with a 220 MeV proton beam and the resulting spallogenic ³He used as a diffusant in step-heating diffusion experiments. In every sample, irradiation increased the activation energy (E_a) and the frequency factor (D_o/a²) of diffusion and yielded a higher He closure temperature (T_c) than the starting material. For example, 100 h in the reactor caused the He closure temperature to increase by as much as 36 °C. For a given neutron fluence the magnitude of increase in closure temperature scales negatively with the initial closure temperature. This is consistent with a logarithmic response in which the neutron damage is additive to the initial damage present. In detail, the irradiations introduce correlated increases in E_a and ln(D_o/a²) that lie on the same array as found in natural apatites. This strongly suggests that neutron-induced damage mimics the damage produced by U and Th decay in natural apatites.To investigate the potential consequences of annealing of radiation damage, samples of Durango apatite were heated in vacuum to temperatures up to 550 °C for between 1 and 350 h. After this treatment the samples were step-heated using the remaining natural ⁴He as the diffusant. At temperatures above 290 °C a systematic change in T_c was observed, with values becoming lower with increasing temperature and time. For example, reduction of T_c from the starting value of 71 to ∼52 °C occurred in 1 h at 375 °C or 10 h at 330 °C. The observed variations in T_c are strongly correlated with the fission track length reduction predicted from the initial holding time and temperature. Furthermore, like the neutron irradiated apatites, these samples plot on the same E_a − ln(D_o/a²) array as natural samples, suggesting that damage annealing is simply undoing the consequences of damage accumulation in terms of He diffusivity.Taken together these data provide unequivocal evidence that at these levels, radiation damage acts to retard He diffusion in apatite, and that thermal annealing reverses the process. The data provide support for the previously described radiation damage trapping kinetic model of Shuster et al. (2006) and can be used to define a model which fully accommodates damage production and annealing.     

Survey of garnets for fossil fission tracks

Andradites and spessartines contain enough uranium to make these minerals interesting for fission track dating. Almandines and pyropes usually are not favourable.     

The effect of long-term low-temperature exposure on apatite fission track stability: A natural annealing experiment in the deep ocean

Since studies on deep-sea cores were carried out in the early 1990s it has been known that ambient temperature may have a marked affect on apatite fission track annealing. Due to sluggish annealing kinetics, this effect cannot be quantified by laboratory annealing experiments. The unknown amount of low-temperature annealing remains one of the main uncertainties for extracting thermal histories from fission track data, particularly for samples which experienced slow cooling in shallow crustal levels. To further elucidate these uncertainties, we studied volcanogenic sediments from five deep-sea drill cores, that were exposed to maximum temperatures between ∼10° and 70 °C over geological time scales of ∼15-120 Ma. Mean track lengths (MTL) and etch pit diameters (Dpar) of all samples were measured, and the chemical composition of each grain analysed for age and track length measurements was determined by electron microprobe analysis. Thermal histories of the sampled sites were independently reconstructed, based on vitrinite reflectance measurements and/or 1D numerical modelling. These reconstructions were used to test the most widely used annealing models for their ability to predict low-temperature annealing. Our results show that long-term exposure to temperatures below the temperature range of the nominal apatite fission track partial annealing zone results in track shortening ranging between 4 and 11%. Both chlorine content and Dpar values explain the downhole annealing patterns equally well. Low chlorine apatite from one drill core revealed a systematic relation between Si-content and Dpar value. The question whether Si-substitution in apatite has direct and systematic effects on annealing properties however, cannot be addressed by our data. For samples, which remained at temperatures <30 °C, and which are low in chlorine, the Laslett et al. [Laslett G., Green P., Duddy I. and Gleadow A. (1987) Thermal annealing of fission tracks in apatite. Chem. Geol. 65, 1-13] annealing model predicts MTL up to 0.6 μm longer than those actually measured, whereas for apatites with intermediate to high chlorine content, which experienced temperatures >30 °C, the predictions of the Laslett et al. (1987) model agree with the measured MTL data within error levels. With few exceptions, predictions by the Ketcham et al. [Ketcham R., Donelick R. and Carlson W. (1999) Variability of apatite fission-track annealing kinetics. III: Extrapolation to geological time scales. Am. Mineral. 84/9, 1235-1255] annealing model are consistent with the measured data for samples which remained at temperatures below ∼30 °C. For samples which experienced maximum temperatures between ∼30 and 70 °C, and which are rich in chlorine, the Ketcham et al. (1999) model overestimates track annealing.     

班公湖- 怒江缝合带东段八宿地区花岗岩体差异性隆升：来自锆石和磷灰石裂变径迹证据

班公湖- 怒江缝合带为青藏高原内部分隔羌塘和拉萨两地块的构造边界,是研究青藏高原构造演化的重要窗口之一。该缝合带自西向东分为西段（班公湖至改则）、中段（安多至东巧）和东段（丁青至怒江）,其中东段的研究程度较低。本次以东段八宿县郭庆乡一条花岗岩高程剖面为研究对象,采用激光剥蚀电感耦合等离子体质谱仪（LA- ICPMS）法对锆石和磷灰石开展裂变径迹测试。花岗岩锆石U- Pb年龄为~180 Ma,指示其结晶时代为早侏罗世。锆石和磷灰石裂变径迹年龄分别为180～130 Ma、86～61 Ma,对应的年龄- 海拔曲线分别为负斜率和正斜率。QTQt模拟显示花岗岩高程剖面顶部在130～60 Ma时剥蚀冷却速率快,中部在130～40 Ma时剥蚀冷却速率居中,而底部在~130 Ma之后一直保持最低的剥蚀冷却速率。这种差异性隆升源自班公湖- 怒江缝合带东段的南向俯冲板片断离早于北向俯冲板片断离。     

Etching of fission tracks in monazite: An experimental study

This study reports a range of etching and annealing experiments to establish the optimum conditions for the etching of fission tracks in monazite. The previously reported concentrated (12 M) HCl etchant at 90°C was found to cause grain loss from epoxy mounts and high degrees of grain corrosion, as did much longer etching times at lower temperatures. Using implanted ²⁵²Cf semi‐tracks, a series of experiments were performed on internal prismatic faces of monazite‐(Ce) crystals from the Palaeozoic Harcourt Granodiorite (Victoria, Australia) using an alternative 6 M HCl etchant, also at 90°C. Step‐etch results show optimal etching at 60–90 min. Further, an isothermal annealing experiment illustrated that the degree of annealing that can be expected during etching at 90°C under laboratory time scales is negligible. The etching rate between grains is not uniform, with a correlation demonstrated between over‐etched grains and high U and Th concentrations.     

Systematics in the fission track annealing of minerals

A linear relationship exists between annealing temperatures and activation energies of fission tracks for different minerals. The temperature and activation energy needed to erase 50% of the tracks in minerals depend on the lattice energy per mole. These relationships make possible predictions about the stability of tracks in silicates, provided their chemical composition is known or a few isochronal annealing data are available.     

Offshore granulites from the Bay of Biscay margins: fission tracks constrain a Proterozoic to Tertiary thermal history

Zircon and apatite fission track ages were determined on granulites dredged along the Bay of Biscay margins. A sample from Ortegal Spur (Iberia margin) yielded 725 ± 67 Ma (zircon). A sample from Le Danois Bank (Iberia margin) yielded 284 ± 58 Ma (zircon), indicating post‐Variscan cooling. Apatite from this sample gave 52 ± 2 Ma, interpreted as final cooling after ‘Pyrenean’ thrust imbrication. Two other samples from Le Danois Bank have Early Cretaceous apatite ages (138 ± 7 and 120 ± 8 Ma), interpreted to result from exhumation during rifting. Finally, a granulite from Goban Spur (Armorican margin) gave 212 ± 10 Ma (apatite), coinciding with a precursory rifting phase. Together with published radiometric results, these data indicate a Precambrian high‐grade terrane at the site of the current margins. The distribution of the granulites on the seafloor reflects tectonic and erosional processes related to (a) Mesozoic rifting and (b) Early Tertiary incipient subduction of the Bay of Biscay beneath Iberia.     

Forward kinematic modelling of a regional transect in the Northern Emirates using geological and apatite fission track age constraints on paleo-burial history

This paper aims to simulate the kinematic evolution of a regional transect crossing the Northern Emirates in the northernmost part of the Semail Ophiolite and the Dibba zone, just south of the Musandam Platform exposures. The studied section comprises, from top to bottom and from inner to outer zones, (1) the erosional remnants of the Semail Ophiolite, mainly made up of serpentinized ultramafics in the west and gabbros in the east, (2) high-grade metamorphic rocks which are currently exposed in the core of a nappe anticline near Masafi, (3) far-travelled Hawasina basinal units and Sumeini paleo-slope units of the Dibba Zone, (4) parautochthonous platform carbonates, which are currently well exposed in the Musandam area, and (5) a flexural basin filled with uppermost Cretaceous to Neogene sediments. Two main compressional episodes are generally identified, resulting first in the obduction of the Semail Ophiolite and then in the stacking of underlying platform carbonate units of the former Arabian passive margin, thus accounting for the present architecture of this transect: (1) first, deformation at the plate boundary initiated in the Late Cretaceous, resulting in the obduction of the Semail Ophiolite and the progressive accretion of the Hawasina and Sumeini tectonic wedge on top of the Arabian foreland, leading to a progressive bending of its lithosphere and development of a wide flexural basin; (2) compression resumed during the Neogene, leading to the tectonic stacking of the parautochthonous platform duplexes of Musandam and Margham trends, the development of out-of-sequence thrusts and triangle zones, refolding of the sole thrust of the former Late Cretaceous accretionary wedge and coeval normal (?) high-angle faulting along the contact between the Musandam and Dibba zones. However, seismic profiles and paleo-thermometers also help in identifying another erosional event at the boundary between the Paleogene Pabdeh and the Neogene Fars series. Evidenced by the local erosional truncation of the Pabdeh series in the vicinity of the frontal triangle zone (i.e. the inner part of the former Late Cretaceous foredeep), this Paleogene uplift/unroofing episode is tentatively interpreted here as an evidence for a continuum of compressional deformation lasting from the Late Cretaceous to the Middle Miocene although one may alternatively speculate that it was related to the detachment of the subducted slab. Although carbonate facies are usually not suitable for apatite fission track (AFT) studies, we were able to extract detrital apatites from quartz-bearing Triassic dolomites in the Musandam area. However, the yield and the quality were both poor and too few fission track lengths could be measured, making it difficult to interpret the meaning of the FT ages. The FT dates obtained in this study are therefore compared with those existing in the literature. Fortunately enough, for each sample, at least ten apatite crystals could be used for fission track dating, except for site 6 with only five datable apatite grains. The obtained apatite fission track dates between 28 and 13 Ma, much younger than the Triassic age of the series, are taken to represent reset fission track ages, implying erosion of an up-to-3-km-thick pile of Jurassic–Cretaceous carbonates and Hawasina allochthon during the Neogene. Apatite fission track dates from the ~95 M-old plagiogranites of the Semail complex (Searle and Cox, Geol Mag 139(3):241–255, 2002) obtained in this study and compared with those recently published provide evidences for more than one cooling event. An early unroofing of the ophiolite during the Late Cretaceous is revealed in fission track dates of 72–76 Ma at the top of the ophiolite in the east, which are coeval and also consistent with the occurrence of paleo-soils, rudists and paleo-reefs on top of serpentinized ultramafics in the west. High-pressure rocks at As Sifah in the southeast near Muscat revealed apatite fission track data ranging from ~46 to 63 Ma (Gray et al. 2006). The leucocratic part of the ophiolite (sample UAE 180) yielded comparable young apatite (40.6?±?3.9 Ma) and zircon (46.6?±?4.3 Ma) FT dates. A Cenozoic (~20–21 Ma) exhumation has been determined for the Bani Hamid metamorphic sole in northern Oman, applying low temperature geochronology and combining apatite FT and apatite (U–Th)/He analyses (Gray et al. 2006). In this study, young apatite fission track dates of 20 Ma have also been found but at the base of the ophiolite near Masafi, in the core of the nappe anticline, thus indicating a Neogene age for the refolding of the allochthon and stacking of underlying parautochthonous platform carbonate units. During the subsequent 2D forward Thrustpack kinematic modelling of the regional transect, these AFT data-set has been used, together with available subsurface information, to reconstruct the past architecture of the structural sections through time, accounting for incremental deformation along the various decollement levels, synorogenic sedimentation and erosion, as well as for successive bending and unbending episodes of the Arabian lithosphere.     

南海北部扩张的地表过程响应：来自海南岛五指山的磷灰石裂变径迹年龄证据

海南岛是南海地块的重要组成部分,其构造演化很大程度上能够反映整个南海的构造活动特征。新生代,海南岛在新构造运动作用下,断块差异升降比较明显,形成了大致以王五 文教断裂为界,北为沉降区,南为断隆区的构造格局。对断隆区隆升过程的研究能够帮助揭示海南岛新构造活动历史,但至今为止,琼中南山地隆升的原因和时限仍存在争议。为了厘清海南岛中部的剥露隆升事件,本文选择琼中南地区海拔最高、高差最大的五指山为研究区,采集8组高程岩石样品,高程范围为203. 55~1153. 52 m。对采集的样品进行了磷灰石裂变径迹测试和热历史模拟分析,结果表明海南岛五指山地区新生代主要经历了两期快速隆升剥露。第一期为渐新世—中新世（32~17 Ma)：隆升速率较快,此时期太平洋板块向欧亚板块俯冲后撤,南海正经历第二次扩张,使得海南岛拉张,活动强烈,造成琼中山地区快速隆升,直到中中新世转为缓慢隆升。第二期为中新世末期（5 Ma)的快速剥蚀隆升阶段：南海扩张已经结束,随着菲律宾板块俯冲亚洲板块,南海北部陆缘整体处于加速热沉降阶段,且全球气候变化加快,造成了海南岛广泛的隆起和加速剥蚀。     

怒江河砂岩屑磷灰石裂变径迹结果与流域地貌演化

河流搬运沉积的河砂作为流域内地质体的平均产物,可以有效地揭示整个流域内区域性的地质体热史演化。对同一河流进行分段采样能够揭示更为详实的热史演化差异。本文对滇西境内怒江上游至下游采集了7个河砂样品进行磷灰石裂变径迹定年,主要年龄峰值依次为: 12.2 Ma和12.8 Ma, 7.7 Ma, 5.3 Ma、4.4 Ma和4.9 Ma, 7.3 Ma,总体上呈现出上游老、中游年轻、下游老的年龄格局。怒江中游河段相对年轻的流域热史经历说明中游河段较其他区域经受了构造地貌的快速演化。由于这一位置对应三江并流,而且气候作用与现代西南季风的水汽通道吻合,降雨量与上、下游流域截然不同。可以推断,气候因素导致了其快而新的构造地貌演化,促成代表年轻冷却历史信息的地质体被迅速剥露造就了怒江河道不同位置流域热史演化的差异。测得的7个样品所有年龄峰值可分为5个区段: 5.3～4.4 Ma、7.7～7.3 Ma、12.8～10.7 Ma、26.8～22.2.3 Ma、48.7～30.1 Ma,基本反映了流域范围新生代以来主要的热史演化阶段。此外,怒江不同河段河砂样品的分段性特征,证实河砂岩屑磷灰石裂变径迹分析方法在解释流域区域热史演化方面具有独特的优点。     

Denudation history of the Snowy Mountains: Constraints from apatite fission track thermochronology

Apatite fission track thermochronology from Early Palaeozoic granitoids centred around the Kosciuszko massif of the Snowy Mountains, records a denudation history that was episodic and highly variable. The form of the apatite fission track age profile assembled from vertical sections and hydroelectric tunnels traversing the mountains, together with numerical forward modelling, provide strong evidence for two episodes of accelerated denudation, commencing in Late Permian—Early Triassic (ca 270–250 Ma) and mid‐Cretaceous (ca 110–100 Ma) times, and a possible third episode in the Cenozoic. Denudation commencing in the Late Permian—Early Triassic was widespread in the eastern and central Snowy Mountains area, continued through much of the Triassic, and amounted to at least ～2.0–2.4 km. This episode was probably the geomorphic response to the Hunter‐Bowen Orogeny. Post‐Triassic denudation to the present in these areas amounted to ～2.0–2.2 km. Unambiguous evidence for mid‐Cretaceous cooling and possible later cooling is confined to a north‐south‐trending sinuous belt, up to ～15 km wide by at least 35 km long, of major reactivated Palaeozoic faults on the western side of the mountains. This zone is the most deeply exposed area of the Kosciuszko block. Denudation accompanying these later events totalled up to ～1.8–2.0 km and ～2.0–2.25 km respectively. Mid‐Cretaceous denudation marks the onset of renewed tectonic activity in the southeastern highlands following a period of relative quiescence since the Late Triassic, and establishes a temporal link with the onset of extension related to the opening of the Tasman Sea. Much of the present day relief of the mountains resulted from surface uplift which disrupted the post‐mid‐Cretaceous apatite fission track profile by variable offsets on faults.