Fission track length histograms from sedimentary samples may contain tracks generated before the time of deposition. They may originate from various source areas with their own thermal histories. After deposition they share thermal history. Inherited tracks, mainly represented by the shortest tracks, are mixed with the young post–sedimentary tracks due to anisotropic annealing, various grain chemistry and uncertainties of measure. Inherited tracks complicate calculating the thermal history. Deconvolution from time sequence analysis [1] is used to separate pre–and post–depositional tracks. Deconvolution reduces the length spread caused by the uncertainties. Hereby tracks are organized according to age, the short tracks being older than the long tracks.

Caprock assemblages associated with salt bodies typically consist of a vertically zoned sequence in ascending order: anhydrite directly above the salt body, a transitional gypsum zone, and occasionally a complex zone of limestone and/or dolomite. Caprock forms when the upper part of a rising diapir is exposed to a crossflow of NaCl-undersaturated water, causing halite to dissolve and the less soluble components, largely anhydrite (CaSO4) and to a lesser extent gypsum (CaSO4•2H2O), to accrete via underplating to the base of the previously formed caprock. If hydrocarbons are present, the CaSO4 minerals are replaced by carbonate minerals in a process mediated by sulfate-reducing bacteria. Previous descriptions of carbonate caprock recognize only two general fabric types: an upper variegated limestone and a lower banded zone that comprises carbonate and sulfate lithologies. Utilizing new facies mapping and petrographic analysis of outcropping caprock from three different salt basins, Paradox Basin and Gulf Coast Region, USA and Flinders Ranges, South Australia, we recognize a wider variety of fabrics and mineralogies. This variety is owed to the location of caprock at the salt-sediment interface, where it is highly prone to substantial chemical and mechanical alteration, resulting in diagenetic textural overprints of precursor caprock fabrics. We propose a new classification based on fabric types in order to facilitate a discussion and interpretation of caprock lithologies in an organized and effective manner. The development of a comprehensive classification is the first step toward deciphering the complex diagenetic processes involved in caprock formation. Understanding the genetic history of caprock fabrics will allow for better identification and prediction of the distribution of caprock mineralogies and fabrics. Our proposed new classification scheme is based on the recognition of four distinct megascopic fabrics: 1) Massive: homogeneous, with micro-to-coarsely crystalline subdivisions; 2) Layered: subdivisions based on thickness of laminae include micro-laminated, laminated, and banded; 3) Brecciated: subdivided into mosaic and disorganized, dictated by clast orientation and volume variations in inter-clast matrix or cement; and 4) Porphyritic: comprising two distinct crystal sizes. These fabrics are not restricted to any one mineralogy (i.e. are found within anhydrite, gypsum and carbonate caprock) and commonly comprise more than one fabric type (e.g. the brecciated clasts of a carbonate caprock can display a layered or massive fabric). To address this issue, the dominant mineralogy and subordinate fabric types is attached to the overall fabric name as prefix-type modifiers (e.g. massive brecciated dolomite caprock). Silicification, carbonate and silica pseudomorphs after gypsum and anhydrite, localized isoclinal folds, secondary dissolution and recrystallization porosity, as well as dead oil are common features found within all fabric types.

Insights from the RAPID–MOCHA observation network in the North Atlantic have motivated a recent focus on the South Atlantic, where water masses are exchanged with the neighboring Indian and Pacific ocean basins. Moreover, the South Atlantic meridional overturning circulation basin-wide array (SAMBA) was recently launched to monitor variability in the South Atlantic MOC (SAMOC) at 34.5ºS. In this study, we are interested in understanding the processes which generate volume transport variability that would be observed at this latitude band. To perform this attribution, we compute sensitivities of the SAMOC at 34ºS to atmospheric state variables (e.g. wind stress, precipitation) using the adjoint of a global ocean model which is fit to a vast number of ocean observations over the past 20 years. These sensitivities isolate the impact from each atmospheric variable, and highlight the oceanic mechanisms, such as Kelvin and Rossby waves, which carry atmospheric forcing perturbations to the SAMOC. The domain of influence for the SAMOC is shown to be quite broad, covering neighboring ocean basins even on short time scales. This result differs from what has previously been shown in the North Atlantic, where Atlantic meridional overturning circulation (AMOC) variability is largely governed by dynamics confined to that basin. We convolve historical forcing variability from ERA-Interim with the computed sensitivities in order to attribute seasonal to interannual SAMOC variability to each atmospheric component. The seasonal cycle of the SAMOC is therefore shown to be largely driven by local zonal wind forcing. Interannual variability, however, is shown to have originated from remote locations across the globe, including a nontrivial component originating from the tropical Pacific. We conclude with preliminary results which employ both modeling results and an analysis of modern altimetry observations to show how El Niño Southern Oscillation variability might influence the South Atlantic.

Hurricane Michael in 2018 was one of the strongest storms to impact the coastal U.S. and was unusual given that it occurred in October (i.e., late in the hurricane season) and intensified over the continental shelf. A potential contributor to this extreme event was thought to be anomalously high heat content on the shelf of the Mississippi Bight, a shelf region significantly impacted by freshwater discharge. Using available long-term time series of regional meteorological and oceanic measurements, water column conditions during the run-up to the rapid intensification of Hurricane Michael were compared to historical conditions in the region. Data for the water column heat content in the western Mississippi Bight were available at a mooring site on the 20 m isobath (Site CP) during August-October of 2018. Unusually high heat content was observed, relative to the typical summer conditions of previous years (N = 13), which resulted from the compounding effects of atmospheric events in the preceding months (August and September): a series of smaller mixing events (e.g., passage of Tropical Storm Gordon) in conjunction with a regional heatwave during most of September. Tropical Storm Gordon traveled across the shelf of the Mississippi Bight, disrupting the stratified water column as observed at Site CP. The stratification breakdown mixed the upper water column heat content deeper into the water column and subsequently allowed for the anomalously warm atmospheric conditions in September to effectively transfer heat deeper into the water column. While the mooring site was significantly distant (250 km) from the center of Hurricane Michael, these processes observed in the western Mississippi Bight likely occurred in the eastern portion of the basin as well. As a result of these compounding atmospheric effects, the shelf water column was primed to support the intensification of following tropical storms, which highlights the need for coupled oceanic-atmospheric forecast models to capture the interaction between the ocean and atmosphere and its effect on water column conditions on the shelf.

Disturbances in the magnetosphere-ionosphere system cause changes in the geomagnetic field that result in ground induced currents (GICs) that are potentially hazardous to electrical systems on Earth. Harmful GICs are driven by magnetic field fluctuations with timescales generally falling in the range of 1-10 minutes; much less attention has been placed on geomagnetic field fluctuations with short timescales (< 60 seconds) because they cause transient induced currents (TICs) that have not been considered to pose a legitimate threat to electrical systems since they are similar to electrical transients due to lightning. On the contrary, short-timescale magnetic field fluctuations have been found to be capable of coupling directly to power grids and electrical systems, inducing substantial voltages without first flowing in the ground. This ionospheric current coupling poses a potential threat to any of these systems, especially electronic equipment with low operating voltage or that does not have surge protection. Transmitting devices that are at risk by such currents are becoming increasingly more prevalent in society with the growth of the Internet of Things (IoT) network. Our characterization of transient magnetic field perturbations provides detail on short-timescale changes of the magnetosphere-ionosphere coupled system and supports the assessment of possible hazards to technological infrastructure on Earth. This research is enabled by modern magnetometers, both ground- and space-based, with high sampling rate capabilities that allow for the measurement of transient surface magnetic field fluctuations with short-timescales. We present the characteristics of transient magnetic field changes observed by the MACCS array in Arctic Canada by selecting events recorded on the ground and analyzing the behavior of the electromagnetic fluctuations within the ionosphere and magnetosphere during such events.

During the last century, descriptions of sediment transport on the surface of Earth have been mostly deterministic and strongly influenced by concepts from continuum mechanics. The assumption that particle motions on hillslopes and in rivers satisfy the continuum hypothesis has provided an important foundation for this topic. Recent studies, however, have recognized that bed load and hillslope sediment transport conditions often are rarefied and do not satisfy continuum assumptions, therein pointing to the need for new ways of describing particle motions and transport. The problem of rarefied sediment transport is probabilistic in nature, and emerging methods for describing particle motions hark back to the pioneering work of Einstein (1938), who conceptualized bed load transport as a probabilistic problem. Here we provide a data set of particle travel distances and supplemental high-speed videos of particle-surface collisions collected during laboratory experiments to assess a theoretical formulation of the probabilistic physics of rarefied particle motions and deposition on rough hillslope surfaces. The formulation is based on a description of the kinetic energy balance of a cohort of particles treated as a rarefied granular gas, and a description of particle deposition that depends on the energy state of the particles. Both laboratory and field-based measurements are consistent with a generalized Pareto distribution of travel distances and predicted variations in behavior associated with the balance between gravitational heating and frictional cooling by particle-surface collisions. These behaviors vary from a truncated distribution associated with rapid thermal collapse to an exponential distribution representing approximately isothermal conditions to a heavy-tailed distribution associated with net heating of particles. The transition to a heavy-tailed distribution likely involves an increasing conversion of translational to rotational kinetic energy leading to larger travel distances with decreasing effectiveness of collisional friction. The analysis points to the need for further clarity concerning how particle size and shape in concert with surface roughness influence the extraction of particle energy and the likelihood of deposition.

Variability in precipitation frequency, intensity and duration deeply affect vegetation and erosion yields. Under climate change scenarios, alterations in precipitation are expected but it is not well understood how it could affect erosion rates and what is the role of vegetation on landscape geomorphic response. Traditional erosion models normally include basic representations of vegetation that do not account for the dynamic character of biomass variability and feedbacks with hydrological and erosion/deposition processes. Hence, using a new modelling frameworks that account for the effect of varying vegetation cover on erosion, and that includes climate change scenarios is needed. Here we use a new model: COPLAS, a tool that couples a Landform Evolution Model with dynamic vegetation and carbon pools modules to investigate the response of landscapes to climate change. The vegetation module includes a coupled photosynthesis-stomatal conductance representation that responds to climatic data inputs as temperature, CO2 concentration and water availability. We use the model to simulate the erosional and geomorphic responses of dynamic vegetation in Howard Springs (Australia) to predicted changes in daily precipitation under future CO2 concentrations of about 940 ppm. The model was calibrated using Ozflux site historical data. Catchment scale simulations were run for a period of 100 years for three scenarios (bare soil, constant and dynamic vegetation) using a daily time step. We found that, for our study case, bare soil produces on average 139% more erosion than the constant vegetation case and 124% more than the dynamic vegetation case. Moreover, an increase in precipitation of around 23% induces an increase of 25%, 35% and 43% in the erosion rates for dynamic vegetation, constant vegetation and bare soil respectively, while a decrease in 26% reduces it in 36%, 60% and 59%. This could be explained by the nonlinear relation between erosion and vegetation (higher rainfall induces higher erosion potential which can be counteracted by an increase in vegetation cover leading to a decrease in soil erodibilty). This finding highlights the importance of considering the dynamic character of vegetation in order to understand the nonlinear relations between fluvial erosion and vegetation cover.

It is now known that many spiral density waves in Saturn’s C ring are forced by resonant gravitational interactions with the planet’s nonradial oscillations. Observations of their periods have the power to reveal new information concerning subtle aspects of the planet’s internal structure. The properties of a certain class of nonradial oscillations, consisting of the inertial and rotational normal modes (i-modes and r-modes), are particularly sensitive to small departures of the internal density stratification from adiabatic, and may be especially useful in revealing such departures and their implications for any compositional gradients residing inside Saturn. We calculate the frequencies and eigenfunctions (shapes) of the i-modes and r-modes for a given suite of internal structure models and determine the Lindblad resonance locations and periods of the density waves they excite in the C ring. We find that for reasonable models of Saturn’s internal structure, the i-mode and r-mode frequencies overlap the range of frequencies of the slowest density waves observed in the C ring, that is, those having pattern speeds similar to Saturn’s rotation rate. In addition, the presence of observable density waves forced by r-modes requires Saturn to possess (at least) two statically stable (non-convective) regions, one in the deep interior and the other in the outer molecular hydrogen envelope. The inner zone is most likely stabilized by a significant compositional gradient, the outer zone by a compositional gradient, by baroclinic processes, or both. Unresolved questions remain. Currently, the model predicts that if energy is equipartitioned among the i-modes, then some slow density waves should be detectable at radii inside ~84500 km in the C ring. None have been observed inside this radius. In addition, the source of excitation of the normal modes remains to be identified. We are currently exploring whether they can draw their energy from baroclinic instability in the outer stable zone. These questions will be addressed at length at the meeting. This work is supported by the NASA SSW program.

CO2 enhanced oil recovery (EOR) is being increasingly deployed for the exploitation of depleted conventional oil and gas fields, due to the twofold benefit of improved recovery efficiency and reduction in carbon dioxide emissions. Geomechanics starts playing an increasingly important role with the injection of CO2 in the subsurface and the subsequent pore pressure buildup. However, the release of pressure through oil production enables a higher amount of CO2 to be introduced into the ground without causing undesirable effects. Understanding the stress perturbation due to fluid injection and withdrawal aids in comprehending the fault activation mechanism and the risks of induced seismicity. The current study evaluates the geomechanical influence of CO2 injection for enhanced oil recovery in a depleted oil field, and assesses the risk of injection and production induced seismicity as a constraint in geologic containment of CO2. The Ankleshwar field in Cambay basin in India is chosen for the study. The field is a potential CO2 EOR site due to its excellent permeability and recovery efficiency. Coupled multiphase fluid flow modeling and geomechanical analysis were carried out to study the hydro-mechanical characteristics of CO2 injection and fluid production. The rise in pore pressure near the injection site induces stress changes even far away from the site of injection due to poroelastic coupling in rocks. The risk of induced seismicity was then analyzed through simulation of fault slip potential in the field. FSP analysis suggests orientation of faults along with proximity to injection site are key parameters influencing fault stability in the Ankleshwar field.

Mantle transition zone (MTZ) is the layer between two seismic upper mantle discontinuities, which are generally believed to be caused by the transformation of high-pressure polymorphs of olivine occurred around depth of 410 km (d410) and 660 km (d660). The fine structure around the discontinuities, especially the undulation of them can be used to infer the variation of temperature and/or composition in the deep mantle. Seismic images have revealed a subhorizontally deflected slab beneath Northeast China, and thus the seismic mapping of the interaction between the slab and the deep mantle is important in understanding the mantle circulation as well as the upper mantle composition. Recently, a temporary NECsaids Array has been installed, and it consists mainly of two nearly orthogonal dense profiles with station spacing of 15–20 km, which provides us an unprecedented opportunity to obtain a high-resolution 3D topography of the MTZ. Considering the dipping and stagnant structure of the slab, we used a novel Ps scatter kernel migration technique to obtain the detailed image of the d410 and d660, which has been proved to be efficient and better for more complex structures. Our results show that the boundaries of the subducting Northwest Pacific slab can be clearly identified in the image, with a dip angle of about 30° in the MTZ. The depth of the d660 depressed 20 - 40 km around the stagnant part of the subducting slab, which could be explained with the lower temperature caused by the cold Pacific slab. Moreover, the d660 is uplifted about 5 - 15 km beneath northwest of Changbaishan and Longgang Volcano, suggesting that the upwelling feeding the volcanoes might originate from the MTZ. Apart from the prominent depression of the d410 observed in the volcanic area, we also observed a wide-distributed depression of the d410 with value ~10 km over most studied areas, seeming inconsistent with the thermal anomaly associated with the subduction environment. We argue that subtle variations of water and/or chemical composition contribute to the observed depth variation of the d410, and the apparent depressed “410” can be related to the wet/dry wadsleyite boundary, which might be caused by the high water content at the top of the MTZ.

The GEOSS Water Strategy–From Observations to Decisions (Lawford et al., 2014) identifies several key water variables as Essential Water Variables (EWVs). This was based on broad meta-surveys (Unninayar et al., 2010) of water-related observational needs for GEO Societal Benefit Aeas (SBAs) that included energy and agriculture among others. This paper summarizes currently defined EWVs required by key research and applications sectors involved in decision support within the Water-Energy-Food (WEF) Nexus. There is a notable overlap between the EWVs and a preliminary assessment of Essential WEF Variables. We also highlight EWVs that are relevant to the indicator monitoring objectives of the UN Sustainable Development Goals (SDGs), and GEOGLOWS priority thematic communities—especially those that overlap with the WEF Nexus. Besides primary EWVs such as precipitation, soil moisture, evapotranspiration and water levels/storage, supplemental EWVs are identified that support the integrated multi-sectoral information needed by WEF decision support applications. Examples of supplemental EWVs include surface meteorology and winds, solar radiation, land use/land cover and vegetation. Both Remote Sensing (RS) platforms and In-Situ observing networks are required to address the broad range of space/time resolutions, accuracies, and data latencies that end-users need. The AGU-2021 is invited to comment on, endorse and/or recommend additional EWVs that should be considered for adoption by GEOGLOWS and the GEO (Group on Earth Observations).

The instrumental data shows 7%-20% increase in annual discharge of the major Eurasian rivers like Ob, Yenisei and Lena between 1936 and 2018 (Wang et al. 2021). The trend has been attributed to increased precipitation and permafrost thawing due to the temperature warming (Walvoord, Kurylyk 2016). However, the instrumental data does not provide the longterm scope of the trend. We modeled seasonal discharge from tree rings for the Yenisei River upstream at Kyzyl gauge and found a remarkable 80% upsurge in winter ‘ow (Nov-Apr) over the last 25 years, which is unprecedented in the last 214 years since 1784 (Panyushkina et al. 2021). In contrast, the annual discharge (Oct-Sept) has only a 7% increase over the last 25 years and shows normal range of variability since 1700 (Fig. 1). Water balance modeling with CRU data at the Yenisei upstream indicates a significant discrepancy between decadal variability of the gauged ‘ow and climate data after 1960 (Fig. 1). The long-term buckhound of the changes in regional hydrology is successfully assessed with tree-ring methods. The tree-ring networks in Eurasian cold climates have a great potential to reconstruct spatial pattern of the seasonal runoff and quantify the long-term impact of permafrost degradation on the hydrological regimes in Siberia. We discuss the impact of melting permafrost on the base ‘ow and enrichment of the surface and groundwater interaction at the Yenisei River basin coupled with the warming temperature and, more importantly, forest fires. Recent increase in the frequency, size, and intensity of boreal fires scale s of its impact on hydrology and permafrost in Siberia. This study demonstrates the complexity of hydrological feedback in Siberia to the Arctic Ampli:cation (AA). The adverse impacts of AA have been and will remain the greatest for the health and socioeconomic of people living in the Pan-Arctic and the geopolitics and macroeconomics of the global society.

Polar regions are highly sensitive to climate change. In recent decades, the Arctic has warmed twice as fast as the world average, which has led to a significant loss of ice cover in the Arctic Ocean. The positive feedback from continental hydrology to Arctic warming amplifies perturbations in the climate system in response to changes in heat and freshwater fluxes. As the surface warms and storm paths change, precipitation and glacier melt have increased pan-Arctic runoff. Melting permafrost adds even more water to river systems. Our research aims to understand 1) whether this feedback impacts the river flow only in the warm season, 2) whether the runoff change in the cold (winter) season can amplify the initial warming, and 3) what is the spatiotemporal pattern of seasonal variations in streamflow across pan-Arctic watersheds. Tree rings serve as critical proxies for quantifying hydrological responses to climate change. We are exploring the potential of modeling seasonal flows, especially winter versus spring and summer flows, and water temperature using intra-seasonal aggregation of climate signals in subsets of ring chronologies derived from quantified tracheal cell parameters: cell area, tangential cell diameter, radial cell diameter, cell wall thickness. The preliminary results demonstrate significant correlation of summer water temperature with cell area and radial cell diameter, whereas the discharge correlates stronger with cell wall thickness and tangential cell diameter. We emphasize the importance of modeling intra-seasonal hydrological parameters versa the hydrological year average to analyze the contribution of continental hydrology to Arctic warming.

In the near-future, new space telescopes like JWST, LUVOIR, and HabEx will begin attempting to explore the properties of atmospheres of potentially habitable planets. This will require a significant amount of time and resources for even a single planet, which makes it essential to prioritize observations by those most-likely to have detectable life. Here we present a statistical method to estimate the probabilities that specific exoplanets have been continuously in the habitable zone of their host stars for more than 2 billion years, the approximate time it took life on Earth to significantly increase the oxygen content of the atmosphere. We introduce the use of statistics of an ensemble of 3D planetary general circulation models to estimate these probabilities, replacing prior 1D model estimates.

The thermochronological method was applied to metamorphic rocks distributed in eastern Nepal to elucidate the denudation process of the upper-crust of the continental collision zone. New results of systematic fission-track (FT) age dating and FT length measurements of zircon and apatite were utilized in the thermochronological inverse analysis to reconstruct the time-temperature (t-T) paths in the temperature range of 60–350°C. Eight t-T paths obtained along the across-strike section in eastern Nepal showed that the cooling process of the metamorphic rocks are characterized by 1) gradual cooling (<30°C/Myr) followed by rapid cooling (~150°C/Myr) and subsequent gradual cooling (gradual-rapid-gradual cooling: GRG cooling), 2) northward-younging of the timing of the rapid cooling since ca. 9 Ma. The observed FT ages and t-T paths were then compared with the FT ages and t-T paths obtained by forward calculations using 3-D thermokinematic models to test the following four tectono-thermal models which have been proposed in the Central and Eastern Himalayas: 1) The denudation of the upper-crust is associated with the movement of the plate boundary fault (Main Himalayan Thrust: MHT) showing flat geometry (the Flat MHT model) and 2) flat-ramp-flat geometry (the Flat-Ramp-Flat MHT model), 3) the denudation of the upper-crust is mainly controlled by the focused uplift associated with the growth of the Lesser Himalayan Duplex (the Duplex 01-03 model) or 4) slip of the splay fault of the MHT (the Splay Fault model). As a result, only the Flat-Ramp-Flat MHT model reconstructed similar t-T paths and age distribution patterns obtained from eastern Nepal. This suggests that the observed FT ages and t-T paths reflect a denudation process driven by the movement of the MHT showing a flat-ramp-flat geometry. The GRG cooling and the northward-younging of the timing of rapid cooling indicate that the flat-ramp-flat geometry of the MHT was established by ca. 9 Ma and has been stable thereafter. The result of the thermokinematic inverse analysis also indicates that the denudation rate and its spatial distribution have been stable since ca. 9 Ma.

We report a new series of step-etch experiments to reveal the influence of microscopy technique on track selection bias. Two different aliquots of induced tracks in Durango apatite were etched for 10-15-20-25-30s and for 20-25-30s in 5.5M HNO3 at 21°C. Three different track selection criteria were applied after the initial etch step of the etching procedure: (1) all tracks with reasonable measurability under transmitted and reflected light switch with 100x objective and 2.5X optovar magnification; (2) ”fully etched” tracks under transmitted light with 100x objective; and (3) ”fully etched” tracks under reflected light with 100x objective. Approach 1 was applied to both aliquots and the approaches 2 and 3 to the latter aliquot. Comparing the mean track lengths, approaches 2 and 3 result in similar values over all experiments, while approach 1 provides a ~0.4 μm higher in the first aliquot and ~0.3 μm lower in the second aliquot than approach 2 and 3 due to higher and lower variations of effective etch times. Comparing the c-axis angles, in the 0-30° range approach 3 provides a severely reduced fraction of tracks due to their weak appearance under reflected light. Furthermore, approach 2 provides ~%14 lower track densities comparing with approaches 1 and 3. We recommend using both transmitted and reflected light during entire track selection and measurement procedures. We are working to develop a new 2+-step etching procedure, where tracks are located after 10s of etching but measured after a second 10s etch step, resulting in better-controlled etching times while reducing the bias associated with analyst choice. Furthermore, this two-step etching procedure can be iterated for more etch-steps, by identifying newly-appeared tracks after each etch step and etching them for 10s more, to increase the number of measured tracks while maintaining consistent selection criteria.

The zircon (U-Th)/He (ZHe) system with a typical closure temperature of ~160-200°C*1, but lower for higher radiation damaged grains*2, offers the potential for evaluating thermal histories in the uppermost ~10 km of the crust. ZHe thermochronometry has been applied to different geological settings in order to estimate tectonics, uplift and denudation, basin evolution, etc.*3, which can also contribute to evaluating long-term tectonic stabilities for the geologic disposal project. So far, the effectivity of ZHe thermochronology has been verified, however improved age standards for the method are required. To date, the method has conventionally employed zircon fission-track age standards such as the Fish Canyon Tuff (FCT) zircon*4. ZHe grain ages are sometimes over-dispersed owing to factors such as zoning of parent nuclei, radiation damage, grain size and He-bearing inclusions*2,5. Considerable parent isotope zonation was reported in some FCT crystals*6, inviting a search for alternative potential ZHe standards*7,8,9. These works reported robust ZHe data with little age dispersion because of homogeneous U-Th distribution in zircon megacrysts, making them possible reference material candidates. However, a practical issue remains because ZHe analyses of unknown samples are carried out grain-by-grain as opposed to analyzing large pieces of a single grain. We have attempted to assess suitable zircon samples as ZHe age standards by using rapid cooling rock samples of relatively young (<100 Ma) age. This is because such rock samples are expected to empirically exhibit simple thermal histories and little radiation damage. Therefore, age dispersion caused by radiation damage can be relatively small. In order to reassess previous data obtained by Tagami et al. (2003)*10, ZHe analyses of the Pliocene Utaosa rhyolite (TRG-04 and -07) and the Miocene Buluk Tuff have been carried out. In addition, OD-3 zircon*11, a zircon U-Pb age standard, was also analyzed. In this presentation, preliminary ZHe age data from these samples will be presented and compared to evaluate their suitability as ZHe reference materials e.g., FCT. References 1: Reiners et al. (2004), Geochim. Cosmochim. Acta, 68, p. 1857–1887 2: Guenthner et al. (2013), Am. J. Sci., 313, p. 145–198 3: Ault et al. (2019), Tectonics, 38, p. 3705–3739 4: Gleadow et al. (2015), Earth Planet Sci. Lett., 424, p. 95–108 5: Danišík et al. (2017), Sci. Adv., 3, p. 1–9 6: Dobson et al. (2008), Geochim. Cosmochim. Acta, 72, p4745–4755 7: Li et al. (2017), Geostand. Geoanal. Res., 41, p. 359–365 8: Yu et al. (2020), Geostand. Geoanal. Res., 44, p. 763–783 9: Kirkland et al. (2020), Geochim. Cosmochim. Acta, 274, p. 1–19 10:Tagami et al. (2003), Earth Planet Sci. Lett., 207, p. 57–67 11:Iwano et al. (2013), Isl. Arc, 22, p. 382–394 Acknowledgements This study was supported by the Ministry of Economy, Trade, and Industry (METI) of Japan.

Optically stimulated luminescence (OSL) thermochronometry is a tool for constraining cooling histories in low-temperature domains (several tens of degree Celsius) during the past 104–105 years [1][2][3]. This method is currently applied only to rapidly denuded regions (~5mm/yr when a general geothermal gradient in is assumed to be ~0.03℃/m) because luminescence signals in slowly denuded regions saturate before the rocks are exhumated to the surface. However, cooling histories in slowly denuded regions may be constrained if unsaturated samples are obtained from deep boreholes. In addition, using deep borehole core enable to compare the results between samples at multiple depths, which is useful to isolate the denudation history from other events, such as faulting or hydrothermal activity. We applied multi-OSL-thermochronometry [2] to the deep borehole core drilled at the Rokko Mountains, Japan, where slow denudation rates (0.1-1.0 mm/yr) are expected from previous studies [4][5][6]. We used the Kabutoyama core collected by National Research Institute for Earth Science and Disaster Resilience [5][7]. The total length of Kabutoyama core is 1,313 m and we collected the samples at 408, 642, 818 and 1048 m for OSL-thermochronometry. Our results showed that useful thermal information can be extracted from the infrared stimulated luminescence signals of samples collected at depths ≥408 m. We found that the sample temperatures remained around the present ambient temperature at each depth for the last 0.1 Myr, indicating that the Rokko Mountains is topographically stable, which was consistent with previous findings. Thus, the thermal denudation history of slowly denuded regions may be constrained by multi-OSL-thermochronometry using samples from deep borehole cores. However, the denudation rates in the Rokko Mountains were too low and could not be determined by this method. Further research is required to quantify the denudation rate. This study was funded by the Ministry of Economy, Trade and Industry (METI), Japan as part of its R&D supporting program titled “Establishment of Advanced Technology for Evaluating the Long-term Geosphere Stability on Geological Disposal Project of Radioactive Waste (Fiscal Years 2019-2021)”. References: [1] Herman et al. (2010). Earth and Planetary Science Letters, 297, 183-189; [2] King et al. (2016). Quaternary Geochronology, 33, 76-87; [3] Herman and King (2018). Elements, 14, 33-38; [4] Huzita (1968). The Quaternary Research, 7, 248-260; [5] Sueoka et al. (2010). Journal of Geography, 119, 84-101; [6] Matsuhi et al. (2014). Transactions, Japanese Geomorphological Union, 35, 165-185; [7] Yamada et al. (2012). Technical Note of the National Research Institute for Earth Science and Disaster Prevention, 371, 27p.