Search

Process-based ecosystem models are increasingly important for predicting forest dynamics under future environmental conditions, which may encompass non-analogous climate coupled with unprecedented disturbance regimes. However, challenges persist due to the extensive number of model parameters, scarce calibration data, and trade-offs between the local precision and the applicability of the model over a wide range of environmental conditions. In this paper, we describe a protocol that allows a modeller to collect transferable ecosystem properties based on ecosystem characteristic criteria and to compile the parameters that need to be described in the field. We applied the procedure to develop a new parameterisation for European beech (Fagus sylvatica L.) for the Biome-BGCMuSo model, the most advanced member of the Biome-BGC family. For model calibration and testing, we utilised multiyear forest carbon data from 87 plots distributed across five European countries. The initial values of 48 new ecophysiological parameters were defined based on a literature review. The final values of six calibrated parameters were optimised for single sites as well as for multiple sites using generalised likelihood uncertainty estimation (GLUE) and model output conditioning that ensured plausible simulations based on user-defined ranges of carbon stock output variables (carbon stock in aboveground wood biomass, soil, and litter) and finding the intersections of site-specific plausible parameter hyperspaces. To support the model use, we tested the model performance by simulating aboveground tree wood, soil, and litter carbon across a large geographical gradient of central Europe and evaluated the trade-offs between parameters tailored to single plots and parameters estimated using multiple sites. Our findings indicated that parameter sets derived from single sites provided an improved local accuracy of simulations of aboveground wood, soil, and litter carbon stocks by 35 %, 55 %, and 11 % in comparison to the a priori parameter set. However, their broader applicability was very limited. A multi-site optimised parameter set, on the other hand, performed satisfactorily across the entire geographical domain studied here, including on sites not involved in the parameter estimation, but the errors were, on average, 26 %, 35 % and 9 % greater for the aboveground wood, soil, and litter carbon stocks than those obtained with the site-specific parameter sets. Importantly, model simulations demonstrated plausible responses across large-scale environmental gradients, featuring a clear production optimum of beech that aligns with empirical studies. These findings suggest that the model is capable of accurately simulating the dynamics of European beech across its range and can be used for more comprehensive experimentations.

In situ cosmogenic 36Cl analysis by accelerator mass spectrometry (AMS) is routinely employed to date Quaternary surfaces and assess rates of landscape evolution. However, standard laboratory preparation procedures for 36Cl dating require the addition of large amounts of isotopically enriched chlorine spike solution; these solutions are expensive and increasingly difficult to acquire from commercial sources. In addition, the typical workflow for 36Cl dating involves measuring both 35Cl/37Cl and 36Cl/Cl concurrently on the high-energy (post-accelerator) end of the AMS system, but 35Cl/37Cl determinations using this technique can be complicated by isotope fractionation and system memory during measurement. The traditional workflow also does not provide 36Cl extraction laboratories with the data needed to calculate native Cl concentrations in advance of 36Cl/Cl measurements. In light of these concerns, we present an improved workflow for extracting and measuring chlorine in geologic materials. Our initial step is to characterize 35Cl/37Cl on sample aliquots of up to ∼1 g prepared in Ag(Cl, Br) matrices, which greatly reduces the amount of isotopically enriched spike solution required to measure native Cl content in each sample. To avoid potential issues with isotope fractionation through the accelerator, 35Cl/37Cl is measured on the low-energy, pre-accelerator end of the AMS line. Then, for 36Cl/Cl measurements, we extract Cl as AgCl or Ag(Cl, Br) in analytical batches with a consistent total Cl load across all samples; this step is intended to minimize source memory effects during 36Cl/Cl measurements and allows the preparation of AMS standards that are customized to match known Cl contents in the samples. To assess the efficacy of this extraction and measurement workflow, we compare chlorine isotope ratio measurements on seven geologic samples prepared using standard procedures and the updated workflow. Measurements of 35Cl/37Cl and 36Cl/Cl are consistent between the two workflows, and 35Cl/37Cl values measured using our methods have considerably higher precision than those measured following standard protocols. The chemical preparation and measurement workflow presented here (1) reduces the amount of isotopically enriched chlorine spike used per rock sample by up to 95 %; (2) identifies rocks with high native Cl concentrations, which may be lower priority for 36Cl surface exposure dating, at an early stage of analysis; and (3) allows laboratory users to maintain control over the total chlorine content within and across analytical batches. These methods can be incorporated into existing laboratory and AMS protocols for 36Cl analyses and will increase the accessibility of 36Cl dating for geologic applications.

Cosmogenic beryllium-10 and beryllium-7, and the ratio of the two (10Be/7Be), are powerful atmospheric tracers of stratosphere–troposphere exchange (STE) processes; however, measurements are sparse for altitudes well above the tropopause. We present a novel high-altitude balloon campaign aimed to measure these isotopes in the mid-stratosphere called Beryllium Isotopes for Resolving Dynamics in the Stratosphere (BIRDIES). BIRDIES targeted gravity waves produced by tropopause-overshooting convection to study their propagation and impact on STE dynamics, including the production of turbulence in the stratosphere. Two custom-designed payloads called FiSH and GASP were flown at altitudes approaching 30 km to measure in situ turbulence and beryllium isotopes (on aerosols), respectively. These were flown on nine high-altitude balloon flights over Kansas, USA, in summer 2022. The atmospheric samples were augmented with a ground-based rainfall collection targeting isotopic signatures of deep convection overshooting. Our GASP samples yielded mostly negligible amounts of both 10Be and 7Be collected in the mid-stratosphere but led to design improvements to increase aerosol capture in low-pressure environments. Observations from FiSH and the precipitation collection were more fruitful. FiSH showed the presence of turbulent velocity, temperature, and acoustic fluctuations in the stratosphere, including length scales in the infra-sonic range and inertial subrange that indicated times of elevated turbulence. The precipitation collection, and subsequent statistical analysis, showed that large spatial datasets of 10Be/7Be can be measured in individual rainfall events with minimum terrestrial contamination. While the spatial patterns in rainfall suggested some evidence for overshooting convection, inter-event temporal variability was clearly observed and predicted with good agreement using the 3D chemical transport model GEOS-CHEM.

Measurements of cosmic-ray-produced beryllium-10, neon-21, and helium-3 in quartz in a soil profile from a forested landscape in the Oregon Coast Range show that the cosmogenic noble gases 21Ne and 3He are depleted relative to 10Be in the shallow subsurface. The noble gases are mobile in mineral grains via thermally activated diffusion and 10Be is not, implying that noble gas depletion is the result of surface heating by wildfires and subsequent mixing of partially degassed quartz downward into the soil. Cosmogenic noble gas depletion by wildfire heating of soils is a potential means of estimating wildfire intensity and/or frequency over pre-observational timescales.

Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

Here, we present cosmogenic-10Be and cosmogenic-3He data from Ferrar dolerite pyroxenes in surficial rock samples and a bedrock core from the McMurdo Dry Valleys, Antarctica, with the goal of refining the laboratory methods for extracting beryllium from pyroxene, further estimating the 10Be production rate in pyroxene and demonstrating the applicability of 10Be–3He in mafic rock. The ability to routinely measure cosmogenic 10Be in pyroxene will open new opportunities for quantifying exposure durations and Earth surface processes in mafic rocks. We describe scalable laboratory methods for isolating beryllium from pyroxene, which include a simple hydrofluoric acid leaching procedure for removing meteoric 10Be and the addition of a pH 8 precipitation step to reduce the cation load prior to ion exchange chromatography. 10Be measurements in pyroxene from the surface samples have apparent 3He exposure ages of 1–6 Myr. We estimate a spallation production rate for 10Be in pyroxene, referenced to 3He, of 3.6 ± 0.2 atoms g−1 yr−1. 10Be and 3He measurements in the bedrock core yield initial estimates for parameters associated with 10Be and 3He production by negative-muon capture (f10∗=0.00183 and f3∗fCfD=0.00337). Next, we demonstrate that the 10Be–3He pair in pyroxene can be used to simultaneously resolve erosion rates and exposure ages, finding that the measured cosmogenic-nuclide concentrations in our surface samples are best explained by 2–8 Myr of exposure at erosion rates of 0–35 cm Myr−1. Finally, given the low 10Be in our laboratory blanks (average of 5.7 × 103 atoms), the reported measurement precision, and our estimated production rate, it should be possible to measure 2 g samples with 10Be concentrations of 6 × 104 and 1.5 × 104 atoms g−1 with 5 % and 15 % uncertainty, respectively. With this level of precision, Last Glacial Maximum to Late Holocene surfaces can now be dated with 10Be in pyroxene. Application of 10Be in pyroxene, alone or in combination with 3He, will expand possibilities for investigating glacial histories and landscape change in mafic rock.

Deglaciation of the northwestern Laurentide Ice Sheet in the central Mackenzie Valley opened the northern portion of the deglacial Ice-Free Corridor between the Laurentide and Cordilleran ice sheets and a drainage route to the Arctic Ocean. In addition, ice sheet saddle collapse in this section of the Laurentide Ice Sheet has been implicated as a mechanism for delivering substantial freshwater influx into the Arctic Ocean on centennial timescales. However, there is little empirical data to constrain the deglaciation chronology in the central Mackenzie Valley where the northern slopes of the ice saddle were located. Here, we present 30 new 10Be cosmogenic nuclide exposure dates across six sites, including two elevation transects, which constrain the timing and rate of thinning and retreat of the Laurentide Ice Sheet in the area. Our new 10Be dates indicate that the initial deglaciation of the eastern summits of the central Mackenzie Mountains began at ∼15.8 ka (17.1–14.6 ka), ∼1000 years earlier than in previous reconstructions. The main phase of ice saddle collapse occurred between ∼14.9 and 13.6 ka, consistent with numerical modelling simulations, placing this event within the Bølling–Allerød interval (14.6–12.9 ka). Our new dates require a revision of ice margin retreat dynamics, with ice retreating more easterly rather than southward along the Mackenzie Valley. In addition, we quantify a total sea level rise contribution from the Cordilleran–Laurentide ice saddle region of ∼11.2 m between 16 and 13 ka.

We use 25 new measurements of in situ produced cosmogenic 26Al and 10Be in river sand, paired with estimates of dissolved load flux in river water, to characterize the processes and pace of landscape change in central Cuba. Long-term erosion rates inferred from 10Be concentrations in quartz extracted from central Cuban river sand range from 3.4–189 Mg km−2 yr−1 (mean 59, median 45). Dissolved loads (10–176 Mg km−2 yr−1; mean 92, median 97), calculated from stream solute concentrations and modeled runoff, exceed measured cosmogenic-10Be-derived erosion rates in 18 of 23 basins. This disparity mandates that in this environment landscape-scale mass loss is not fully represented by the cosmogenic nuclide measurements. The 26Al / 10Be ratios are lower than expected for steady-state exposure or erosion in 16 of 24 samples. Depressed 26Al / 10Be ratios occur in many of the basins that have the greatest disparity between dissolved loads (high) and erosion rates inferred from cosmogenic nuclide concentrations (low). Depressed 26Al / 10Be ratios are consistent with the presence of a deep, mixed, regolith layer providing extended storage times on slopes and/or burial and extended storage during fluvial transport. River water chemical analyses indicate that many basins with lower 26Al / 10Be ratios and high 10Be concentrations are underlain at least in part by evaporitic rocks that rapidly dissolve. Our data show that when assessing mass loss in humid tropical landscapes, accounting for the contribution of rock dissolution at depth is particularly important. In such warm, wet climates, mineral dissolution can occur many meters below the surface, beyond the penetration depth of most cosmic rays and thus the production of most cosmogenic nuclides. Our data suggest the importance of estimating solute fluxes and measuring paired cosmogenic nuclides to better understand the processes and rates of mass transfer at a basin scale.

Cosmogenic beryllium-10 and beryllium-7, and the ratio of the two (10Be/7Be), are powerful atmospheric tracers of stratosphere–troposphere exchange (STE) processes; however, measurements are sparse for altitudes well above the tropopause. We present a novel high-altitude balloon campaign aimed to measure these isotopes in the mid-stratosphere called Beryllium Isotopes for Resolving Dynamics in the Stratosphere (BIRDIES). BIRDIES targeted gravity waves produced by tropopause-overshooting convection to study their propagation and impact on STE dynamics, including the production of turbulence in the stratosphere. Two custom-designed payloads called FiSH and GASP were flown at altitudes approaching 30 km to measure in situ turbulence and beryllium isotopes (on aerosols), respectively. These were flown on nine high-altitude balloon flights over Kansas, USA, in summer 2022. The atmospheric samples were augmented with a ground-based rainfall collection targeting isotopic signatures of deep convection overshooting. Our GASP samples yielded mostly negligible amounts of both 10Be and 7Be collected in the mid-stratosphere but led to design improvements to increase aerosol capture in low-pressure environments. Observations from FiSH and the precipitation collection were more fruitful. FiSH showed the presence of turbulent velocity, temperature, and acoustic fluctuations in the stratosphere, including length scales in the infra-sonic range and inertial subrange that indicated times of elevated turbulence. The precipitation collection, and subsequent statistical analysis, showed that large spatial datasets of 10Be/7Be can be measured in individual rainfall events with minimum terrestrial contamination. While the spatial patterns in rainfall suggested some evidence for overshooting convection, inter-event temporal variability was clearly observed and predicted with good agreement using the 3D chemical transport model GEOS-CHEM. [ABSTRACT FROM AUTHOR]

Terrestrial biogeochemical models are essential tools to quantify climate–carbon cycle feedback and plant–soil relations from local to global scale. In this study, a theoretical basis is provided for the latest version of the Biome-BGCMuSo biogeochemical model (version 6.2). Biome-BGCMuSo is a branch of the original Biome-BGC model with a large number of developments and structural changes. Earlier model versions performed poorly in terms of soil water content (SWC) dynamics in different environments. Moreover, lack of detailed nitrogen cycle representation was a major limitation of the model. Since problems associated with these internal drivers might influence the final results and parameter estimation, additional structural improvements were necessary. In this paper the improved soil hydrology as well as the soil carbon and nitrogen cycle calculation methods are described in detail. Capabilities of the Biome-BGCMuSo v6.2 model are demonstrated via case studies focusing on soil hydrology, soil nitrogen cycle, and soil organic carbon content estimation. Soil-hydrology-related results are compared to observation data from an experimental lysimeter station. The results indicate improved performance for Biome-BGCMuSo v6.2 compared to v4.0 (explained variance increased from 0.121 to 0.8 for SWC and from 0.084 to 0.46 for soil evaporation; bias changed from −0.047 to −0.007 m3 m−3 for SWC and from −0.68 to −0.2 mm d−1 for soil evaporation). Simulations related to nitrogen balance and soil CO2 efflux were evaluated based on observations made in a long-term field experiment under crop rotation. The results indicated that the model is able to provide realistic nitrate content estimation for the topsoil. Soil nitrous oxide (N2O) efflux and soil respiration simulations were also realistic, with overall correspondence with the observations (for the N2O efflux simulation bias was between −0.13 and −0.1 mgNm-2d-1, and normalized root mean squared error (NRMSE) was 32.4 %–37.6 %; for CO2 efflux simulations bias was 0.04–0.17 gCm-2d-1, while NRMSE was 34.1 %–40.1 %). Sensitivity analysis and optimization of the decomposition scheme are presented to support practical application of the model. The improved version of Biome-BGCMuSo has the ability to provide more realistic soil hydrology representation as well as nitrification and denitrification process estimation, which represents a major milestone.

Long-term erosion rates in Tasmania, at the southern end of Australia's Great Dividing Range, are poorly known; yet, this knowledge is critical for making informed land-use decisions and improving the ecological health of coastal ecosystems. Here, we present quantitative, geologically relevant estimates of erosion rates for the George River basin, in northeast Tasmania, based on in situ-produced 10Be (10Bei) measured from stream sand at two trunk channel sites and seven tributaries (mean: 24.1±1.4 Mgkm-2yr-1; 1σ). These new 10Bei-based erosion rates are strongly related to elevation, which appears to control mean annual precipitation and temperature, suggesting that elevation-dependent surface processes influence rates of erosion in northeast Tasmania. Erosion rates are not correlated with slope in contrast to erosion rates along the mainland portions of Australia's Great Dividing Range. We also extracted and measured meteoric 10Be (10Bem) from grain coatings of sand-sized stream sediment at each site, which we normalize to measured concentrations of reactive 9Be and use to estimate 10Bem-based denudation rates for the George River. 10Bem/9Bereac denudation rates replicate 10Bei erosion rates within a factor of 3 but are highly sensitive to the value of 9Be that is found in bedrock (9Beparent), which was unmeasured in this study. 10Bem/9Bereac denudation rates seem sensitive to recent mining, forestry, and agricultural land use, all of which resulted in widespread topsoil disturbance. Our findings suggest that 10Bem/9Bereac denudation metrics will be most useful in drainage basins that are geologically homogeneous, where recent disturbances to topsoil profiles are minimal, and where 9Beparent is well constrained.

Glaciers preserve climate variations in their geological and geomorphological records, which makes them prime candidates for climate reconstructions. Investigating the glacier–climate system over the past millennia is particularly relevant first because the amplitude and frequency of natural climate variability during the Holocene provides the climatic context against which modern, human-induced climate change must be assessed. Second, the transition from the last glacial to the current interglacial promises important insights into the climate system during warming, which is of particular interest with respect to ongoing climate change. Evidence of stable ice margin positions that record cooling during the past 12 kyr are preserved in two glaciated valleys of the Silvretta Massif in the eastern European Alps, the Jamtal (JAM) and the Laraintal (LAR). We mapped and dated moraines in these catchments including historical ridges using beryllium-10 surface exposure dating (10Be SED) techniques and correlate resulting moraine formation intervals with climate proxy records to evaluate the spatial and temporal scale of these cold phases. The new geochronologies indicate the formation of moraines during the early Holocene (EH), ca. 11.0 ± 0.7 ka (n = 19). Boulder ages along historical moraines (n = 6) suggest at least two glacier advances during the Little Ice Age (LIA; ca. 1250–1850 CE) around 1300 CE and in the second half of the 18th century. An earlier advance to the same position may have occurred around 500 CE. The Jamtal and Laraintal moraine chronologies provide evidence that millennial-scale EH warming was superimposed by centennial-scale cooling. The timing of EH moraine formation coincides with brief temperature drops identified in local and regional paleoproxy records, most prominently with the Preboreal Oscillation (PBO) and is consistent with moraine deposition in other catchments in the European Alps and in the Arctic region. This consistency points to cooling beyond the local scale and therefore a regional or even hemispheric climate driver. Freshwater input sourced from the Laurentide Ice Sheet (LIS), which changed circulation patterns in the North Atlantic, is a plausible explanation for EH cooling and moraine formation in the Nordic region and in Europe.

Gridded model assessments require at least one climatic and one soil database for carrying out the simulations. There are several parallel soil and climate database development projects that provide sufficient, albeit considerably different, observation based input data for crop model based impact studies. The input database related uncertainty of the Biome-BGCMuSo agro-environmental model outputs was investigated using three and four different gridded climatic and soil databases, respectively covering an area of nearly 100.000 km2 with 1104 grid cells. Spatial, temporal, climate and soil database selection related variances were calculated and compared for four model outputs obtained from 30-year-long simulations. The choice of the input database introduced model output variability that was comparable to the variability the year-to-year change of the weather or the spatial heterogeneity of the soil causes. Input database selection could be a decisive factor in carbon sequestration related studies as the soil carbon stock change estimates may either suggest that the simulated ecosystem is a carbon sink or to the contrary a carbon source on the long run. Careful evaluation of the input database quality seems to be an inevitable and highly relevant step towards more realistic plant production and carbon balance simulations.

Process-based ecosystem models are increasingly important for predicting forest dynamics under future environmental conditions, which may encompass non-analogous climate coupled with unprecedented disturbance regimes. However, challenges persist due to the extensive number of model parameters, scarce calibration data, and trade-offs between the local precision and the applicability of the model over a wide range of environmental conditions. In this paper, we describe a protocol that allows a modeller to collect transferable ecosystem properties based on ecosystem characteristic criteria and to compile the parameters that need to be described in the field. We applied the procedure to develop a new parameterisation for European beech (Fagus sylvatica L.) for the Biome-BGCMuSo model, the most advanced member of the Biome-BGC family. For model calibration and testing, we utilised multiyear forest carbon data from 87 plots distributed across five European countries. The initial values of 48 new ecophysiological parameters were defined based on a literature review. The final values of six calibrated parameters were optimised for single sites as well as for multiple sites using generalised likelihood uncertainty estimation (GLUE) and model output conditioning that ensured plausible simulations based on user-defined ranges of carbon stock output variables (carbon stock in aboveground wood biomass, soil, and litter) and finding the intersections of site-specific plausible parameter hyperspaces. To support the model use, we tested the model performance by simulating aboveground tree wood, soil, and litter carbon across a large geographical gradient of central Europe and evaluated the trade-offs between parameters tailored to single plots and parameters estimated using multiple sites. Our findings indicated that parameter sets derived from single sites provided an improved local accuracy of simulations of aboveground wood, soil, and litter carbon stocks by 35 %, 55 %, and 11 % in comparison to the a priori parameter set. However, their broader applicability was very limited. A multi-site optimised parameter set, on the other hand, performed satisfactorily across the entire geographical domain studied here, including on sites not involved in the parameter estimation, but the errors were, on average, 26 %, 35 % and 9 % greater for the aboveground wood, soil, and litter carbon stocks than those obtained with the site-specific parameter sets. Importantly, model simulations demonstrated plausible responses across large-scale environmental gradients, featuring a clear production optimum of beech that aligns with empirical studies. These findings suggest that the model is capable of accurately simulating the dynamics of European beech across its range and can be used for more comprehensive experimentations. [ABSTRACT FROM AUTHOR]

In situ cosmogenic 36Cl analysis by accelerator mass spectrometry (AMS) is routinely employed to date Quaternary surfaces and assess rates of landscape evolution. However, standard laboratory preparation procedures for 36Cl dating require the addition of large amounts of isotopically enriched chlorine spike solution; these solutions are expensive and increasingly difficult to acquire from commercial sources. In addition, the typical workflow for 36Cl dating involves measuring both 35Cl/37Cl and 36Cl/Cl concurrently on the high-energy (post-accelerator) end of the AMS system, but 35Cl/37Cl determinations using this technique can be complicated by isotope fractionation and system memory during measurement. The traditional workflow also does not provide 36Cl extraction laboratories with the data needed to calculate native Cl concentrations in advance of 36Cl/Cl measurements. In light of these concerns, we present an improved workflow for extracting and measuring chlorine in geologic materials. Our initial step is to characterize 35Cl/37Cl on sample aliquots of up to ∼1 g prepared in Ag(Cl, Br) matrices, which greatly reduces the amount of isotopically enriched spike solution required to measure native Cl content in each sample. To avoid potential issues with isotope fractionation through the accelerator, 35Cl/37Cl is measured on the low-energy, pre-accelerator end of the AMS line. Then, for 36Cl/Cl measurements, we extract Cl as AgCl or Ag(Cl, Br) in analytical batches with a consistent total Cl load across all samples; this step is intended to minimize source memory effects during 36Cl/Cl measurements and allows the preparation of AMS standards that are customized to match known Cl contents in the samples. To assess the efficacy of this extraction and measurement workflow, we compare chlorine isotope ratio measurements on seven geologic samples prepared using standard procedures and the updated workflow. Measurements of 35Cl/37Cl and 36Cl/Cl are consistent between the two workflows, and 35Cl/37Cl values measured using our methods have considerably higher precision than those measured following standard protocols. The chemical preparation and measurement workflow presented here (1) reduces the amount of isotopically enriched chlorine spike used per rock sample by up to 95 %; (2) identifies rocks with high native Cl concentrations, which may be lower priority for 36Cl surface exposure dating, at an early stage of analysis; and (3) allows laboratory users to maintain control over the total chlorine content within and across analytical batches. These methods can be incorporated into existing laboratory and AMS protocols for 36Cl analyses and will increase the accessibility of 36Cl dating for geologic applications. [ABSTRACT FROM AUTHOR]

In this article, I reflect on efforts to revise the instruction and evaluation of an undergraduate writing consultant education course. The revisions are motivated by the desire to adopt practices that reflect the writing center’s commitment to social justice for multilingual/translingual students and by a commitment to provide an effective, flexible, and brave environment for writing consultants to continue their professional development. I argue that grounding understanding of multilingual writers in concepts that explicitly explore linguistic diversity and standardized 1 English ideologies as threshold concepts is essential to reconceptualize writing center practices. I also argue for the necessity of adopting a flexible system for reflection, engagement, and evaluation to support writing consultants’ learning and practice. I share prompts used in the course and some of the responses they generated. The responses suggest that although combining threshold concepts with a portfolio system is successful in supporting inclusive practices, there remains a need to expand more inclusive practices across the university.

In this work, we utilize a novel application of cosmogenic 21Ne measurements in chert to compare exposure times measured in eroding surfaces in the central Jordanian Plateau with exposure times from chert pebbles transported by the Miocene Hazeva River. The Miocene Hazeva River was a large fluvial system (estimated catchment size > 100 000 km2) that drained the Arabian Plateau and Sinai Peninsula into the Mediterranean Sea during the early-to-mid Miocene. It was established after the rifting of the Red Sea uplifted the Arabian Plateau during the Oligocene. Following late-Miocene-to-early-Pliocene subsidence along the Dead Sea rift, the Hazeva drainage system was abandoned and dissected, resulting in new drainage divides on either side of the rift. We find modern erosion rates derived from cosmogenic 21Ne, 26Al, and 10Be in exposed in situ chert nodules to be extremely slow (between 2–4 mm kyr−1). Comparison between modern and paleo-erosion rates, measured in chert pebbles, is not straightforward, as cosmogenic 21Ne was acquired partly during bedrock erosion and partly during transport of these pebbles in the Hazeva River. However, 21Ne exposure times calculated in Miocene cherts are generally shorter (ranging between 0-0+59 and 242±113 kyr) compared to exposure times calculated in the currently eroding chert nodules presented here (269±49 and 378±76 kyr) and other chert surfaces currently eroding in hyperarid environments. Miocene exposure times are shorter even when considering that they account for bedrock erosion in addition to maintained transport along this large river. Shorter exposure times in Miocene cherts correspond to faster paleo-erosion rates, which we attribute to a combination of continuous surface uplift and significantly wetter climatic conditions during the early-to-mid Miocene.

The persistence and size of the Greenland Ice Sheet (GrIS) through the Pleistocene is uncertain. This is important because reconstructing changes in the GrIS determines its contribution to sea level rise during prior warm climate periods and informs future projections. To understand better the history of Greenland's ice, we analyzed glacial till collected in 1993 from below 3 km of ice at Summit, Greenland. The till contains plant fragments, wood, insect parts, fungi, and cosmogenic nuclides showing that the bed of the GrIS at Summit is a long-lived, stable land surface preserving a record of deposition, exposure, and interglacial ecosystems. Knowing that central Greenland was tundra-covered during the Pleistocene informs the understanding of Arctic biosphere response to deglaciation. [ABSTRACT FROM AUTHOR]