Variable aging and storage of dissolved black carbon in the ocean.

During wildfires and fossil fuel combustion, biomass is converted to black carbon (BC) via incomplete combustion. BC enters the ocean by rivers and atmospheric deposition contributing to the marine dissolved organic carbon (DOC) pool. The fate of BC is considered to reside in the marine DOC pool, where the oldest BC 14C ages have been measured (>20,000 14C y), implying long-term storage. DOC is the largest exchangeable pool of organic carbon in the oceans, yet most DOC (>80%) remains molecularly uncharacterized. Here, we report 14C measurements on size-fractionated dissolved BC (DBC) obtained using benzene polycarboxylic acids as molecular tracers to constrain the sources and cycling of DBC and its contributions to refractory DOC (RDOC) in a site in the North Pacific Ocean. Our results reveal that the cycling of DBC is more dynamic and heterogeneous than previously believed though it does not comprise a single, uniformly "old" 14C age. Instead, both semilabile and refractory DBC components are distributed among size fractions of DOC. We report that DBC cycles within DOC as a component of RDOC, exhibiting turnover in the ocean on millennia timescales. DBC within the low-molecular-weight DOC pool is large, environmentally persistent and constitutes the size fraction that is responsible for long-term DBC storage. We speculate that sea surface processes, including bacterial remineralization (via the coupling of photooxidation of surface DBC and bacterial co-metabolism), sorption onto sinking particles and surface photochemical oxidation, modify DBC composition and turnover, ultimately controlling the fate of DBC and RDOC in the ocean.

Millennial-scale hydroclimate control of tropical soil carbon storage.

The storage of organic carbon in the terrestrial biosphere directly affects atmospheric concentrations of carbon dioxide over a wide range of timescales. Within the terrestrial biosphere, the magnitude of carbon storage can vary in response to environmental perturbations such as changing temperature or hydroclimate1, potentially generating feedback on the atmospheric inventory of carbon dioxide. Although temperature controls the storage of soil organic carbon at mid and high latitudes2,3, hydroclimate may be the dominant driver of soil carbon persistence in the tropics4,5; however, the sensitivity of tropical soil carbon turnover to large-scale hydroclimate variability remains poorly understood. Here we show that changes in Indian Summer Monsoon rainfall have controlled the residence time of soil carbon in the Ganges-Brahmaputra basin over the past 18,000 years. Comparison of radiocarbon ages of bulk organic carbon and terrestrial higher-plant biomarkers with co-located palaeohydrological records6 reveals a negative relationship between monsoon rainfall and soil organic carbon stocks on a millennial timescale. Across the deglaciation period, a depletion of basin-wide soil carbon stocks was triggered by increasing rainfall and associated enhanced soil respiration rates. Our results suggest that future hydroclimate changes in tropical regions are likely to accelerate soil carbon destabilization, further increasing atmospheric carbon dioxide concentrations.

Aquatic biomass is a major source to particulate organic matter export in large Arctic rivers.

Arctic rivers provide an integrated signature of the changing landscape and transmit signals of change to the ocean. Here, we use a decade of particulate organic matter (POM) compositional data to deconvolute multiple allochthonous and autochthonous pan-Arctic and watershed-specific sources. Constraints from carbon-to-nitrogen ratios (C:N), δ13C, and Δ14C signatures reveal a large, hitherto overlooked contribution from aquatic biomass. Separation in Δ14C age is enhanced by splitting soil sources into shallow and deep pools (mean ± SD: -228 ± 211 vs. -492 ± 173‰) rather than traditional active layer and permafrost pools (-300 ± 236 vs. -441 ± 215‰) that do not represent permafrost-free Arctic regions. We estimate that 39 to 60% (5 to 95% credible interval) of the annual pan-Arctic POM flux (averaging 4,391 Gg/y particulate organic carbon from 2012 to 2019) comes from aquatic biomass. The remainder is sourced from yedoma, deep soils, shallow soils, petrogenic inputs, and fresh terrestrial production. Climate change-induced warming and increasing CO2 concentrations may enhance both soil destabilization and Arctic river aquatic biomass production, increasing fluxes of POM to the ocean. Younger, autochthonous, and older soil-derived POM likely have different destinies (preferential microbial uptake and processing vs. significant sediment burial, respectively). A small (~7%) increase in aquatic biomass POM flux with warming would be equivalent to a ~30% increase in deep soil POM flux. There is a clear need to better quantify how the balance of endmember fluxes may shift with different ramifications for different endmembers and how this will impact the Arctic system.

Climate control on terrestrial biospheric carbon turnover.

Terrestrial vegetation and soils hold three times more carbon than the atmosphere. Much debate concerns how anthropogenic activity will perturb these surface reservoirs, potentially exacerbating ongoing changes to the climate system. Uncertainties specifically persist in extrapolating point-source observations to ecosystem-scale budgets and fluxes, which require consideration of vertical and lateral processes on multiple temporal and spatial scales. To explore controls on organic carbon (OC) turnover at the river basin scale, we present radiocarbon (14C) ages on two groups of molecular tracers of plant-derived carbon-leaf-wax lipids and lignin phenols-from a globally distributed suite of rivers. We find significant negative relationships between the 14C age of these biomarkers and mean annual temperature and precipitation. Moreover, riverine biospheric-carbon ages scale proportionally with basin-wide soil carbon turnover times and soil 14C ages, implicating OC cycling within soils as a primary control on exported biomarker ages and revealing a broad distribution of soil OC reactivities. The ubiquitous occurrence of a long-lived soil OC pool suggests soil OC is globally vulnerable to perturbations by future temperature and precipitation increase. Scaling of riverine biospheric-carbon ages with soil OC turnover shows the former can constrain the sensitivity of carbon dynamics to environmental controls on broad spatial scales. Extracting this information from fluvially dominated sedimentary sequences may inform past variations in soil OC turnover in response to anthropogenic and/or climate perturbations. In turn, monitoring riverine OC composition may help detect future climate-change-induced perturbations of soil OC turnover and stocks.

Molecular 14 C evidence for contrasting turnover and temperature sensitivity of soil organic matter components.

Climate projection requires an accurate understanding for soil organic carbon (SOC) decomposition and its response to warming. An emergent view considers that environmental constraints rather than chemical structure alone control SOC turnover and its temperature sensitivity (i.e., Q10 ), but direct long-term evidence is lacking. Here, using compound-specific radiocarbon analysis of soil profiles along a 3300-km grassland transect, we provide direct evidence for the rapid turnover of lignin-derived phenols compared with slower-cycling molecular components of SOC (i.e., long-chain lipids and black carbon). Furthermore, in contrast to the slow-cycling components whose turnover is strongly modulated by mineral association and exhibits low Q10 , lignin turnover is mainly regulated by temperature and has a high Q10 . Such contrasts resemble those between fast-cycling (i.e., light) and mineral-associated slow-cycling fractions from globally distributed soils. Collectively, our results suggest that warming may greatly accelerate the decomposition of lignin, especially in soils with relatively weak mineral associations.

Multi-molecular 14C evidence for mineral control on terrestrial carbon storage and export.

Compound- and compound class-specific radiocarbon analysis of source-diagnostic 'biomarker' molecules has emerged as a powerful tool to gain insights into terrestrial carbon cycling. While most studies thus far have focused on higher plant biomarkers (i.e. plant leaf-wax n-alkanoic acids and n-alkanes, lignin-derived phenols), tracing paedogenic carbon is crucial given the pivotal role of soils in modulating ecosystem carbon turnover and organic carbon (OC) export. Here, we determine the radiocarbon (14C) ages of glycerol dialkyl glycerol tetraethers (GDGTs) in riverine sediments and compare them to those of higher plant biomarkers as well as markers of pyrogenic (fire-derived) carbon (benzene polycarboxylic acids, BPCAs) to assess their potential as tracers of soil turnover and export. GDGT Δ14C follows similar relationships with basin properties as vegetation-derived lignin phenols and leaf-wax n-alkanoic acids, suggesting that the radiocarbon ages of these compounds are significantly impacted by intermittent soil storage. Systematic radiocarbon age offsets are observable between the studied biomarkers, which are likely caused by different mobilization pathways and/or stabilization by mineral association. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Radiocarbon signatures of carbon phases exported by Swiss rivers in the Anthropocene.

Lateral carbon transport through the land-to-ocean-aquatic-continuum (LOAC) represents a key component of the global carbon cycle. This LOAC involves complex processes, many of which are prone to anthropogenic perturbation, yet the influence of natural and human-induced drivers remains poorly constrained. This study examines the radiocarbon (14C) signatures of particulate and dissolved organic carbon (POC, DOC) and dissolved inorganic carbon (DIC) transported by Swiss rivers to assess controls on sources and cycling of carbon within their watersheds. Twenty-one rivers were selected and sampled during high-flow conditions in summer 2021, a year of exceptionally high rainfall. Δ14C values of POC range from -446‰ to -158‰, while corresponding ranges of Δ14C values for DOC and DIC are -377‰ to -43‰ and -301‰ to -40‰, respectively, indicating the prevalence of pre-aged carbon. Region-specific agricultural practices seem to have an influential effect on all three carbon phases in rivers draining the Swiss Plateau. Based on Multivariate Regression Analysis, mean basin elevation correlated negatively with Δ14C values of all three carbon phases. These contrasts between alpine terrain and the lowlands reflect the importance of overriding ecoregional controls on riverine carbon dynamics within Switzerland, despite high spatial variability in catchment properties. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Degradation and Aging of Terrestrial Organic Carbon within Estuaries: Biogeochemical and Environmental Implications.

Estuaries are action zones for organic carbon (OC) degradation and aging. These processes influence the nature of terrestrial OC (OCterr) export and the magnitude of OCterr burial in marginal seas, with important environmental implications such as CO2 release and hypoxia. In this study, we determined the contents and carbon isotopic compositions (13C and 14C) of bulk OC and fatty acids (FAs) as well as the sedimentological characteristics of suspended particulate matter (SPM) samples collected from two sites over four seasons and of surface sediment samples from three sites in the Pearl River estuary (PRE) to evaluate processes controlling OCterr degradation and aging along an estuarine gradient. We found that the abundance-weighted average C24-32FA 14C ages increased by an average of ∼1170 years for SPM and by an average of ∼3440 years in PR/PRE sediments, along the ∼60 km PRE transect. These increases in the FA age coincided with an 86% decrease in the corresponding mineral surface area-normalized FA loading along the sediment transport pathway, implying that selective degradation of labile and younger OC resulted in apparent OC aging. These measurements reveal an important shift in the nature of OC, with implications for biogeochemical cycling within estuaries and for regional environmental changes.

Nearshore Zone Dynamics Determine Pathway of Organic Carbon From Eroding Permafrost Coasts.

Collapse of permafrost coasts delivers large quantities of particulate organic carbon (POC) to Arctic coastal areas. With rapidly changing environmental conditions, sediment and organic carbon (OC) mobilization and transport pathways are also changing. Here, we assess the sources and sinks of POC in the highly dynamic nearshore zone of Herschel Island-Qikiqtaruk (Yukon, Canada). Our results show that POC concentrations sharply decrease, from 15.9 to 0.3 mg L-1, within the first 100-300 m offshore. Simultaneously, radiocarbon ages of POC drop from 16,400 to 3,600 14C years, indicating rapid settling of old permafrost POC to underlying sediments. This suggests that permafrost OC is, apart from a very narrow resuspension zone (<5 m water depth), predominantly deposited in nearshore sediments. While long-term storage of permafrost OC in marine sediments potentially limits biodegradation and its subsequent release as greenhouse gas, resuspension of fine-grained, OC-rich sediments in the nearshore zone potentially enhances OC turnover.

An Abrupt Aging of Dissolved Organic Carbon in Large Arctic Rivers.

Permafrost thaw in Arctic watersheds threatens to mobilize hitherto sequestered carbon. We examine the radiocarbon activity (F14C) of dissolved organic carbon (DOC) in the northern Mackenzie River basin. From 2003-2017, DOC-F14C signatures (1.00 ± 0.04; n = 39) tracked atmospheric 14CO2, indicating export of "modern" carbon. This trend was interrupted in June 2018 by the widespread release of aged DOC (0.85 ± 0.16, n = 28) measured across three separate catchment areas. Increased nitrate concentrations in June 2018 lead us to attribute this pulse of 14C-depleted DOC to mobilization of previously frozen soil organic matter. We propose export through lateral perennial thaw zones that occurred at the base of the active layer weakened by preceding warm summer and winter seasons. Although we are not yet able to ascertain the broader significance of this "anomalous" mobilization event, it highlights the potential for rapid and large-scale release of aged carbon from permafrost.

Climate warming alters subsoil but not topsoil carbon dynamics in alpine grassland.

Subsoil contains more than half of soil organic carbon (SOC) globally and is conventionally assumed to be relatively unresponsive to warming compared to the topsoil. Here, we show substantial changes in carbon allocation and dynamics of the subsoil but not topsoil in the Qinghai-Tibetan alpine grasslands over 5 years of warming. Specifically, warming enhanced the accumulation of newly synthesized (14 C-enriched) carbon in the subsoil slow-cycling pool (silt-clay fraction) but promoted the decomposition of plant-derived lignin in the fast-cycling pool (macroaggregates). These changes mirrored an accumulation of lipids and sugars at the expense of lignin in the warmed bulk subsoil, likely associated with shortened soil freezing period and a deepening root system. As warming is accompanied by deepening roots in a wide range of ecosystems, root-driven accrual of slow-cycling pool may represent an important and overlooked mechanism for a potential long-term carbon sink at depth. Moreover, given the contrasting sensitivity of SOC dynamics at varied depths, warming studies focusing only on surface soils may vastly misrepresent shifts in ecosystem carbon storage under climate change.

Impacts of Natural and Human-Induced Hydrological Variability on Particulate Organic Carbon Dynamics in the Yellow River.

Natural and human-induced hydrological changes can influence organic carbon (OC) composition in fluvial systems, with biogeochemical consequences in both terrestrial and marine environments. Here, we use bulk and molecular carbon isotopes (13C and 14C) to examine spatiotemporal variations in particulate OC (POC) composition and age from two locations along the course of the Yellow River during 2015 and 2016. Dual carbon isotopes enable deconvolution of modern, pre-aged (millennial age) soil, and fossil inputs, revealing heterogeneous OC sources at both sites. Pre-aged OC predominated at the upstream site (Huayuankou) throughout the study period, mostly reflecting the upper riverine OC. Strong downstream (Kenli) intra-annual variations in modern and pre-aged OC were caused by increased contributions from modern aquatic OC production under the drier and less turbid conditions during this El Niño year. The month of July, which included the human-induced water and sediment regulation (WSR) event at Kenli, accounted for 82% of annual POC flux, with lower modern OC contribution compared with periods of natural seasonal variability. Both natural and human-induced hydrological events clearly exert strong influence on both fluxes and composition of Yellow River POC which, in turn, affect the balance between OC remineralization and burial for this major fluvial system.

Improved Method for Isolation and Purification of Underivatized Amino Acids for Radiocarbon Analysis.

We have improved a method for isolation and purification of individual amino acids for compound-specific radiocarbon analysis (CSRA). To remove high-performance liquid chromatography (HPLC) eluent blanks from isolated amino acid fractions prior to the radiocarbon (Δ14C) measurement, each fraction was filtered through a membrane filter and then washed with diethyl ether twice. Radiocarbon measurements on standard amino acids processed and purified with the above method using elemental analyzer-accelerator mass spectrometry resulted in Δ14C values that were in strong agreement ( R2 = 0.998) with the original Δ14C value of each amino acid standard. From these measurements, we calculate dead and modern carbon contamination contributions as 1.2 ± 0.2 and 0.3 ± 0.1 µgC, respectively, which are consistent with direct assessments of HPLC procedural blanks of 1.0 ± 0.8 µgC per sample. These contamination constraints allow correction of measured Δ14C values for accurate and precise CSRA and are widely applicable to future archeological and biogeochemical studies.

What on Earth Have We Been Burning? Deciphering Sedimentary Records of Pyrogenic Carbon.

Humans have interacted with fire for thousands of years, yet the utilization of fossil fuels marked the beginning of a new era. Ubiquitous in the environment, pyrogenic carbon (PyC) arises from incomplete combustion of biomass and fossil fuels, forming a continuum of condensed aromatic structures. Here, we develop and evaluate 14C records for two complementary PyC molecular markers, benzene polycarboxylic acids (BPCAs) and polycyclic aromatic hydrocarbons (PAHs), preserved in aquatic sediments from a suburban and a remote catchment in the United States (U.S.) from the mid-1700s to 1998. Results show that the majority of PyC stems from local sources and is transferred to aquatic sedimentary archives on subdecadal to millennial time scales. Whereas a small portion stems from near-contemporaneous production and sedimentation, the majority of PyC (∼90%) experiences delayed transmission due to "preaging" on millennial time scales in catchment soils prior to its ultimate deposition. BPCAs (soot) and PAHs (precursors of soot) trace fossil fuel-derived PyC. Both markers parallel historical records of the consumption of fossil fuels in the U.S., yet never account for more than 19% total PyC. This study demonstrates that isotopic characterization of multiple tracers is necessary to constrain histories and inventories of PyC and that sequestration of PyC can markedly lag its production.

Deep CO2 release and the carbon budget of the central Apennines modulated by geodynamics.

Recent studies increasingly recognize the importance of critical-zone weathering during mountain building for long-term CO2 drawdown and release. However, the focus on near-surface weathering reactions commonly does not account for CO2 emissions from the crust, which could outstrip CO2 drawdown where carbonates melt and decarbonize during subduction and metamorphism. We analyse water chemistry from streams in Italy's central Apennines that cross a gradient in heat flow and crustal thickness with relatively constant climatic conditions. We quantify the balance of inorganic carbon fluxes from near-surface weathering processes, metamorphism and the melting of carbonates. We find that, at the regional scale, carbon emissions from crustal sources outpace near-surface fluxes by two orders of magnitude above a tear in the subducting slab characterized by heat flow greater than 150 mW m-2 and crustal thickness of less than 25 km. By contrast, weathering processes dominate the carbon budget where crustal thickness exceeds 40 km and heat flow is lower than 30 mW m-2. The observed variation in metamorphic fluxes is one to two orders of magnitude larger than that of weathering fluxes. We therefore suggest that geodynamic modulations of metamorphic melting and decarbonation reactions are an efficient process by which tectonics can regulate the inorganic carbon cycle.

Spatiotemporal variation of sterols in sediment as markers of primary production and ocean sewage dumping in the southwestern East Sea (Japan Sea).

The sedimentation of organic carbon in the Ulleung Basin, in the southwestern East Sea (Japan Sea) was investigated using radiocarbon and sterols. The accumulation rates of organic carbon and the contents of brassicasterol and dinosterol were higher on the slope than in the central basin, reflecting the surface water productivity, whereas cholesterol showed similar or higher contents in the central basin. The coprostanol concentration in surface sediments reflected the dispersion of sewage dumped in this region. The vertical distribution showed that the coprostanol concentration was the highest in the top 5-cm layer near the Korea Strait, close to one of the two dumping sites. A high coprostanol concentration was also found near the coast further north, where the content peaked at ∼10 cm depth. The vertical distribution of coprostanol helped to estimate the sediment accumulation rate at sites where radiocarbon gradient was too small or the values were too variable.

Vegetal Undercurrents-Obscured Riverine Dynamics of Plant Debris.

Much attention has been focused on fine-grained sediments carried as suspended load in rivers due to their potential to transport, disperse, and preserve organic carbon (OC), while the transfer and fate of OC associated with coarser-grained sediments in fluvial systems have been less extensively studied. Here, sedimentological, geochemical, and biomolecular characteristics of sediments from river depth profiles reveal distinct hydrodynamic behavior for different pools of OC within the Mackenzie River system. Higher radiocarbon (14C) contents, low N/OC ratios, and elevated plant-derived biomarker loadings suggest a systematic transport of submerged vascular plant debris above the active riverbed in large channels both upstream of and within the delta. Subzero temperatures hinder OC degradation promoting the accumulation and waterlogging of plant detritus within the watershed. Once entrained into a channel, sustained flow strength and buoyancy prevent plant debris from settling and keep it suspended in the water column above the riverbed. Helical flow motions within meandering river segments concentrate lithogenic and organic debris near the inner river bends forming a sediment-laden plume. Moving offshore, we observe a lack of discrete, particulate OC in continental shelf sediments, suggesting preferential trapping of coarse debris within deltaic and neritic environments. The delivery of waterlogged plant detritus transport and high sediment loads during the spring flood may reduce oxygen exposure times and microbial decomposition, leading to enhanced sequestration of biospheric OC. Undercurrents enriched in coarse, relatively fresh plant fragments appear to be reoccurring features, highlighting a poorly understood yet significant mechanism operating within the terrestrial carbon cycle.

Rapid northern hemisphere ice sheet melting during the penultimate deglaciation.

The rate and consequences of future high latitude ice sheet retreat remain a major concern given ongoing anthropogenic warming. Here, new precisely dated stalagmite data from NW Iberia provide the first direct, high-resolution records of periods of rapid melting of Northern Hemisphere ice sheets during the penultimate deglaciation. These records reveal the penultimate deglaciation initiated with rapid century-scale meltwater pulses which subsequently trigger abrupt coolings of air temperature in NW Iberia consistent with freshwater-induced AMOC slowdowns. The first of these AMOC slowdowns, 600-year duration, was shorter than Heinrich 1 of the last deglaciation. Although similar insolation forcing initiated the last two deglaciations, the more rapid and sustained rate of freshening in the eastern North Atlantic penultimate deglaciation likely reflects a larger volume of ice stored in the marine-based Eurasian Ice sheet during the penultimate glacial in contrast to the land-based ice sheet on North America as during the last glacial.

Particulate Organic Matter Dynamics in a Permafrost Headwater Stream and the Kolyma River Mainstem.

Ongoing rapid arctic warming leads to extensive permafrost thaw, which in turn increases the hydrologic connectivity of the landscape by opening up subsurface flow paths. Suspended particulate organic matter (POM) has proven useful to trace permafrost thaw signals in arctic rivers, which may experience higher organic matter loads in the future due to expansion and increasing intensity of thaw processes such as thermokarst and river bank erosion. Here we focus on the Kolyma River watershed in Northeast Siberia, the world's largest watershed entirely underlain by continuous permafrost. To evaluate and characterize the present-day fluvial release of POM from permafrost thaw, we collected water samples every 4-7 days during the 4-month open water season in 2013 and 2015 from the lower Kolyma River mainstem and from a small nearby headwater stream (Y3) draining an area completely underlain by Yedoma permafrost (Pleistocene ice- and organic-rich deposits). Concentrations of particulate organic carbon generally followed the hydrograph with the highest concentrations during the spring flood in late May/early June. For the Kolyma River, concentrations of dissolved organic carbon showed a similar behavior, in contrast to the headwater stream, where dissolved organic carbon values were generally higher and particulate organic carbon concentrations lower than for Kolyma. Carbon isotope analysis (δ13C, Δ14C) suggested Kolyma-POM to stem from both contemporary and older permafrost sources, while Y3-POM was more strongly influenced by in-stream production and recent vegetation. Lipid biomarker concentrations (high-molecular-weight n-alkanoic acids and n-alkanes) did not display clear seasonal patterns, yet implied Y3-POM to be more degraded than Kolyma-POM.

Radiocarbon Age Offsets Between Two Surface Dwelling Planktonic Foraminifera Species During Abrupt Climate Events in the SW Iberian Margin.

This study identifies temporal biases in the radiocarbon ages of the planktonic foraminifera species Globigerina bulloides and Globigerinoides ruber (white) in a sediment core from the SW Iberian margin (so-called Shackleton site). Leaching of the outer shell and measurement of the radiocarbon content of both the leachate and leached sample enabled us to identify surface contamination of the tests and its impact on their 14C ages. Incorporation of younger radiocarbon on the outer shell affected both species and had a larger impact downcore. Interspecies comparison of the 14C ages of the leached samples reveal systematic offsets with 14C ages for G. ruber being younger than G. bulloides ages during the last deglaciation and part of the Early and mid-Holocene. The greatest offsets (up to 1,030 years) were found during Heinrich Stadial 1, the Younger Dryas, and part of the Holocene. The potential factors differentially affecting these two planktonic species were assessed by complementary 14C, oxygen and carbon isotopes, and species abundance determinations. The coupled effect of bioturbation with changes in the abundance of G. ruber is invoked to account for the large age offsets. Our results highlight that 14C ages of planktonic foraminifera might be largely compromised even in settings characterized by high sediment accumulation rates. Thus, a careful assessment of potential temporal biases must be performed prior to using 14C ages for paleoclimate investigations or radiocarbon calibrations (e.g., marine calibration curve Marine13, Reimer et al., 2013, https://doi.org/10.2458/azu_js_rc.55.16947).