1. Field Sites
1.1 Field Site Selection
Four serpentinite bodies in three different climate zones were chosen for investigation to examine climatic influences on chemical weathering of serpentinite bedrock and serpentine soil. These sites include the Trinity Ultramafic Body and Rattlesnake Creek Terrane within the Klamath Mountains of northern California, the Tablelands of Gros Morne National Park in Newfoundland, Canada, and a serpentinite body near the ghost town of Pickhandle Gulch in western Nevada (Table S1, Figure S1).
The modern climate in the Klamath Mountains is mediterranean, featuring hot and dry summers and cold and wet winters. NOAA 1991-2020 climate normals for the Weaverville, CA weather station (Station ID: USC00049490), near the String Bean Creek sampling site in the Rattlesnake Creek Terrane, record mean annual temperatures of 12.8℃, with a mean of 23.2℃ in July and 3.7℃ in December, and an annual average precipitation of ~100 cm1,2. NOAA 1981-2010 climate normals for the Sawyers Bar Ranger Station (Station ID: USC00048025), north of Trinity Lake and within the Trinity Alps, show a mean annual temperature of 12.8℃, with a mean of 22.9℃ in July and 4.4℃ in December, and an annual average precipitation of ~118 cm1,2, 1991-2020 normals were not available for that site. The exact climate at the Eunice Bluff and Deadfall lake sites is likely slightly colder due to the high elevation (~2100 m) compared to ~600-650 m elevation for the monitoring stations, similar high-elevation regions of the Klamath Mountains have mean summertime temperatures (July-October) of ~12.7℃ and mean wintertime (November-May) temperatures hovering around freezing (~0.9℃) though with substantial periods of below-freezing conditions and abundant snowfall3. A dispersed coniferous forest covers much of the Klamath Mountains2,4 and is present at all sampling locations examined within this study. The modern climate of the Tablelands is subarctic; daily mean temperatures at the nearby Cow Head weather station (Station ID: 8401335) maintained by Environment Canada range from -7.0 ℃ in January to 15.7 ℃ in August with a mean yearly temperature of 3.9 ℃, with a yearly average precipitation of 120.4 cm5 from 1980-2010 climate normals. Pickhandle Gulch possesses a desert climate, with the nearest NOAA climate station to Pickhandle Gulch, in Mina, Nevada (Station ID: USC00265168), recording 1991-2020 climate normals of a 14.4 ℃ mean annual temperature, ranging from 2.7 ℃ in December to 27.9 ℃ in July, with an average precipitation of 14.2 cm annually1.
1.2 Klamath Mountains
The Trinity Ultramafic Body covers ~2100 km2 in the eastern Klamath Mountains and is composed of a partially serpentinized mélange of harzburgite and lherzolite (60-70%), dunite (15-20%), and plagioclase lherzolite (10-15%) units cut by clinopyroxene-rich dikes in some areas6,7. Formation of the peridotite sections of the Trinity Ultramafic Body likely occurred around 472 ± 32 Ma with emplacement of the body occurring during the early Devonian8,9. The Rattlesnake Creek Terrane lies southwest of the Trinity Ultramafic Body and is composed of a mixed sequence of serpentinite basement unconformably overlain by Upper Triassic and Lower Jurassic volcanic, hemipelagic, and clastic units with intruded gabbroic and quartz diorites throughout that trend into amphibolite facies and locally higher metamorphic grade materials8,10,11. The serpentinized basement consists of peridotite, greenstone, amphibolite, and pillow basalts that have undergone varying degrees of serpentinization10. The genesis of the Rattlesnake Creek terrane ophiolite and associated facies and their emplacement timing remain poorly understood8.
Samples were collected from five sites within the Klamath Mountains, four sites within the Trinity Ultramafic Body (two sites in high-altitude cirques near Mt. Eddy dubbed Eunice Bluff and Deadfall Lake and two sites within the Swift Creek Valley), and one site within the Rattlesnake Creek Terrane dubbed String Bean Creek (Figure S1). The Trinity Alps, a formerly glaciated region of the Northeast Klamath Mountains where several formerly glaciated valleys feature serpentinite bedrock12, were most recently extensively glaciated coterminous with the late Wisconsinan glacial period of the Laurentide Ice Sheet12. Recent work using cosmogenic techniques has confirmed Younger Dryas deglaciation in high-altitude cirques within the Trinity Alps by ~12.1 ka13,14. The deglaciation date of ~12.1 ka is used as the initiation date for pedogenesis in the high-altitude cirque soils at the Deadfall Lake and Eunice Bluff sites, while within the Swift Creek valley, samples were collected from sites determined by Sharp (1960) to be of Late and Middle Wisconsinan ages. The Late Wisconsinan Cordilleran Ice sheet reached its maximum extent sometime between about 15-23 ka, with initiation of retreat around ~15.6 ka15, with Middle Wisconsinan dates ranging roughly between 25-50 ka16,17. As moraines within Swift Creek lack cosmogenic dates, we use these age ranges as approximate estimates for the ages of the Late and Middle Wisconsinan features described by Sharp (1960)12. In contrast to the Trinity Ultramafic Body, previous glaciation has not been established within the Rattlesnake Creek Terrane, and the soil sampled there is undated but likely substantially older than the soils in the Trinity Ultramafic Body.
1.3 Tablelands
The Gros Morne Tablelands are one of four highland plateaus comprising the Bay of Islands ophiolite, part of the Humber Arm Allochthon located on the western flank of the Island of Newfoundland, Canada18–20. The Bay of Islands ophiolite was emplaced onto the passive continental margin of eastern North America during the Ordovician age Taconian Orogeny19,21. The plateau edges are defined by sharp u-shaped formerly glaciated valleys with associated glacial landforms including moraines, rock glacier deposits, pro-talus lobes, and diamict18,22. Vegetation on the Tablelands is sparse, composed primarily of small ground hugging shrubs with stunted pine trees at lower elevations.
Soil samples were collected from three formerly glaciated valleys on the Tablelands (Figure S1). One sampling site was located just up-slope from and within the protruding lobe of a terminal moraine within the Devil’s Punchbowl cirque. The terminal moraine was dated to 17.6 ± 0.3 ka using chlorine-36 isotope analyses22. One sampling site was located on a terrace overlooking highway 471 in Trout River Gulch. Debris benches in the Trout River Gulch valley were cosmogenically dated using chlorine-36 isotopes to between 15 and 20 ka 22, necessitating the establishment of an ice-free valley prior to that point. Quartzite erratics sourced from over 20 km east of Bonne Bay necessitate westward ice flow at the time of deposition18. As Late Wisconsinan ice flow in the area was directed northward18,23 this westward flow through Trout River Gulch ceased prior to the most recent glacial maximum. The combination of 36Cl ages and glacial erratics indicate that Trout River Gulch was likely not glaciated during the Last Glacial Maximum18,22. Finally, two sampling spots were located on diamict deposits in Winterhouse Gulch, a previously glaciated tributary valley that feeds into Trout River Gulch. One sampling site was located at the mouth of Winterhouse Gulch and one halfway up the valley between the mouth and the cirque. Though cosmogenic ages have not been determined for glacial material in Winterhouse Gulch Canyon, the late Wisconsinan age ranges for other glacial deposits on the Tablelands suggest a similar age range is likely for the presence of a valley glacier in Winterhouse Gulch22, and thus for initiation of soil development following glacial retreat. Reconstructions of Late Wisconsinan ice in Atlantic Canada suggest that the icesheet covering Newfoundland reached a most recent maximum extent around 20 ka23,24 followed by deglaciation that resulted in glaciers retreating away from the coast by ~12 ka 23,24. Initiation of soil development in Winterhouse Gulch likely falls within this age range, concurrent with all other cosmogenically dated glacial landforms surrounding the Tablelands22, with the soil at the mouth of Winterhouse Gulch likely older than the soil closer to the cirque from retreat of the glacier that formed the valley.
1.4 Pickhandle Gulch
Pickhandle Gulch is a ghost town located in the now defunct Candelaria mining district roughly 50 km northwest of Tonopah, Nevada. A serpentinite body roughly 6,800 ft by 300 to 1,350 ft is located southeast of the Northern Belle and Mount Diablo open pit mines with smaller serpentinite lenses scattered nearby, with some exposures of serpentinite within the nearby Northern Belle mine workings25. Soil cover overlying the serpentinite bedrock is thin (<10 cm) and is mantled by a desert pavement composed of serpentinite clasts, with a sparse vegetation cover primarily composed of sagebrush. Previous literature examining the serpentinite body at Pickhandle Gulch is extremely limited and to our knowledge only one report from the Nevada Bureau of Mines attests to the serpentinite’s presence25. The serpentinite is bounded by faults, with a larger mass likely to exist unexposed in the subsurface25. The age and provenance of the serpentinite is currently unknown; however, the body contains inclusions of lower Triassic Candelaria shale material25. The exposure age of the serpentinite is unknown. At Pickhandle Gulch the serpentinite body is bounded by sharp contacts with fluvial channels at its base and toeslope locations could not be identified. Soil samples were therefore collected from a ridgeline position and from two footslope locations on the sample hillslope to provide a degree of spatial variation as time differentiation between soils did not exist in a similar manner as in the Klamath Mountains or Tablelands (Figure S1).
2. Sample Collection and Processing
2.1 Sample Collection
Within the Klamath Mountains and Tablelands, soil pits were located on toeslopes at the base of large hillsides as these have previously been shown to accumulate amorphous material and secondary clay minerals26 and soil development along a hillside is typically greatest at footslope/toeslope positions27. In the Klamath Mountains, there is extensive tree cover, and we chose to locate pits not immediately adjacent to trees to limit the need to cut through large root structures and to avoid such overt biological impacts. In each location, soil pits were excavated by hand using picks and shovels to within a C horizon, contact with the bedrock, point of refusal, or contact with standing water. In the Klamath Mountains, the Eunice Bluff and Deadfall Lake soil pits were excavated to the point of refusal, while the Swift Creek Late, Swift Creek Middle, String Bean Creek soils were excavated until contact with a C horizon determined by decreasing soil redness in relation to overlying soil material 28. In the Tablelands, soil pits were excavated until infilling water was observed in the Winterhouse Gulch Canyon, Winterhouse Gulch Mouth, and Trout River Gulch sites and to contact with bedrock at the Devil’s Punchbowl site. As a result of the shallow nature of the Pickhandle Gulch soils (<~10 cm), soil pits in that field location were excavated to contact with the bedrock.
We collected at least one soil sample from within each soil horizon (see Appendix C). In the Eunice Bluff, Deadfall Lake, and String Bean Creek soils Munsell Color variation allowed for identification of a clear organically enriched A horizon at the surface of each soil pit. We elected to sample in 10 cm depth intervals below the A horizon at those three sites to increase the density of sampling. At the Swift Creek sites in the Klamath Mountains (Late and Middle) one sample was collected from each soil horizon identified in the field. In the Tablelands soil pits, relatively uniform Munsell soil color indicated that A and B horizons had not formed. Because the Tablelands soil profiles were visually undifferentiated, we gathered soil material from four depth intervals within each Tablelands soil pit to provide a higher density of sample material and evaluate possible variation in soil properties with depth that were not visually observable. At Pickhandle Gulch, the presence of small pore spaces indicated the development of Av vesicular surface horizons. We sampled from the Av horizons and from thin underlying C horizons that lacked a similar development of porosity. When referring to our samples in the methods, main text, and supplementary information, we sometimes refer to sampled depth intervals and not to soil horizons, as we sampled multiple different intervals from within some soil horizons in three of the Klamath Mountain sites and from all four Tablelands soils (see Appendix C).
A minimum of ~4 kg of combined bulk soil and gravel material was collected from within each sampled depth interval at all soil pit sampling sites. Parent material bedrock samples were collected from buried clasts from within each soil pit. The collected bulk soil and gravel samples from the Klamath Mountains and Pickhandle Gulch were stored on dry ice during transport and immediately placed in a -20℃ freezer upon return to UNLV. Bulk soil and gravel samples from the Tablelands were packaged and flown back to UNLV and placed in a -20℃ freezer upon arrival but were not packaged on dry ice during transport.
2.2 Parent Material Preparation
We processed one parent material sample each from the Tablelands and Pickhandle Gulch locations as soil pits at these sites were geographically close and developing out of the same ultramafic body. We processed a parent material sample from each Klamath Mountain sampling site as these sites were geographically disparate and developing out of varied ultramafic material. Each parent material sample was trimmed to remove weathering rinds using a cutting saw, crushed within a Bico Chipmunk Jaw Crusher Model 241-36 WD to produce chunks about 3 mm in diameter, and then powdered in a Fritsch pulverisette model 02-101 agate mortar grinder mill. The saw blade used to trim samples was cleaned between uses using a 100% ACS reagent grade ethanol wash. During crushing in the chipmunk jaw crusher, and then again when powdering in the pulverisette, two sub-samples of material were run to pre-contaminate the equipment and those sub-samples were then discarded. A third crushing or powdering run was then conducted using the pre-contaminated instruments and then collected for X-ray diffraction and bulk and selective chemical dissolution analyses. The jaw crusher was thoroughly cleaned between samples using a compressed air blower and a 100% ACS reagent grade ethanol wash. The pulverisette was cleaned between samples with a 100% ACS reagent grade ethanol wash. Following powdering, parent material samples were stored in a -20℃ freezer.
2.3 Soil Preparation and Clay-Size Fraction Separation
After removal from the freezer, each combined soil and gravel sample was first air dried and the <2 mm bulk soil fraction was sieved from the overall sampled material. Gravel content was estimated by massing the <2 mm and >2 mm fractions following sieving. Sub-samples of bulk soil material were then powdered in a Fritsch pulverisette model 02-101 agate mortar grinder mill prior to XRD analysis and application of selective chemical dissolutions. Two sub-samples of bulk soil material were powdered first to pre-contaminate the pulverisette and discarded afterword, and a third sub-sample then powdered and collected for XRD, bulk chemical dissolution, and selective chemical dissolution analyses. The pulverisette was thoroughly washed between samples using 100% ACS reagent grade ethanol. Both powdered and unpowdered bulk soil samples were stored at -20⁰C following powdering.
The clay-size fraction (<2 mm) was separated from the bulk soil following a modification of the methods in Edwards and Bremner (1967). Briefly, ~10 g of soil was mixed with 25 mL of 18.2 MΩ DI water, sonicated for 30 minutes, vortexed, and allowed to settle for 1.5 hours. The suspended load was then pipetted off, NaCl (JT Baker ACS Reagent Grade) was added to the supernatant to a concentration of 1M to promote flocculation, and the samples were allowed to settle for ~30 minutes. Following this, the samples were centrifuged for 3 minutes at 8,000 rpm (8,228 x g) and the clear supernatant decanted from the clay-size fraction pellet. The clay-size fraction pellets then were washed with 18.2 MΩ DI water and centrifuged 3 more times to remove residual salt, frozen at -20℃ for at least 12 hours, and then freeze-dried (LabConco Freezone 4.5 Freeze Dry System).
3. Mineralogical Analyses and Quantification
XRD analyses of randomly oriented powdered bulk soil samples were performed at the UNLV Geosciences Shockwave laboratory on a Proto-AXD Bragg-Brentano type X-ray diffractometer with a Cu K-alpha tube (1.541 Å wavelength), dectris hybrid pixel array detector, 0.4 mm soller slits, and a Ni filter between 4 and 70 ⁰2θ. A knife-blade attachment was added to limit over-saturation of the detector at low angles. Randomly oriented mounts were prepared by top loading the sample into brass or zero background (MTI corp. ZeroSi24D10C1cavity10D) sample stages. Excess material was gently scraped off the top using a razor blade to limit preferred orientation of phyllosilicates. Mineral identification within randomly oriented bulk soil XRD patterns was determined using the program X’Pert HighScore30 and through comparison with d-spacings of crystal lattice planes reported in the American Mineralogist RRUFF database31 and the crystallography open database32.
Randomly oriented clay-size fraction samples extracted as described above were spiked to contain 20 wt.% 0.8 to 1-µm diameter α-Al2O3 (Beantown Chemical 127075-25G) and analyzed from 3 to 80 ⁰2θ at the X-ray Diffraction Laboratory at NASA Johnson Space Center in Houston, Texas on a Malvern PANalytical X’Pert Pro MPD X-ray Diffractometer equipped with a Co K-alpha source (45 kV, 40 mA, ½° anti-scatter slit, ¼° fixed divergence slit, 0.01067° step size, and 0.5 s/step). Spiked samples were prepared by combining pre-weighed α-Al2O3 and powdered parent material or clay-size fraction material and mixing the two materials together within an agate mortar and pestle by a combination of pressing material in toward the center and a circular motion. Mineral identification within randomly oriented clay-size fraction XRD patterns was determined by using the programs X’Pert HighScore30 and Profex 33 and through comparison with d-spacings of crystal lattice planes reported in the American Mineralogist RRUFF database31 and the crystallography open database32.
Oriented mount XRD measurements were conducted on a Proto-AXD Bragg-Brentano type X-ray diffractometer in the UNLV Geosciences Shockwave laboratory with a Cu K-alpha tube (1.541 Å wavelength), dectris hybrid pixel array detector, 0.4 mm soller slits, a Ni filter, and a knife-blade attachment between 2 and 14 ⁰2θ. XRD analyses were conducted on oriented clay-size fraction material following cation saturation and various chemical and heat treatments to determine the presence or absence of expandable 2:1 clay minerals in each soil. Phyllosilicate mineralogy determinations from oriented mount analysis as described below assist with determinations of phyllosilicate mineral presence, and particularly the presence of expandable 2:1 clay minerals, within XRD patterns acquired from randomly oriented mounts.
To identify whether expandable 2:1 clay minerals were present in the clay-size fraction, oriented mounts of the clay-size fraction was analyzed by XRD after 3 different treatments: 1) Saturation with Mg2+ followed by air-drying, 2) ethylene glycol vapor solvation of the air-dried Mg2+ saturated mounts, and 3) K+ saturation followed by heating to 550 °C (see Appendix C). Mg2+ or K+ saturations were prepared34–38 by suspending 300 mg of clay-size fraction material in solutions of either 1N MgCl2 (VWR ACS Reagent grade CAS 7791-18-6) or 0.1N KCl (OmniPur EMG Chemical Argentometric grade CAS 7447-40-7) in 18.2 MΩ water, shaking intermittently by hand for ~1 hour followed by centrifugation for 3 minutes at 8,000 rpm (8,228 x g), after which the supernatant was decanted, and the samples then washed with a 50:50 18.2 MΩ DI water-to-ethanol mixture and centrifuged for 3 minutes at 8,000 rpm 3 times to remove excess salts.
Three mL of 70% ethanol 30% water solution (Aldon Corp. EE0069 Denatured Ethyl Alcohol) were then added to each sample, and the container vortexed to suspend the clay-size particles. ~1 mL of the resulting slurry was then pipetted onto fused silica slides (Alfa Aesar Quartz microscope slides, Part #42295) and allowed to air dry in a covered container overnight in a fume hood.
Previous methods have called for 3-5 overnight wash cycles using cation saturated solutions followed by rinsing35,39, or direct application of cation solution to already oriented clays through a suction device onto ceramic slides followed by a DI-water rinse38. To test whether a 1-hour treatment produced similar saturation to the 3+ day treatment, we compared XRD patterns from a 3-day and 1-hour Mg-saturation procedure. We observed similar peak locations between air dried and glycolated XRD patterns from the two saturation durations indicating that a 1-hour saturation procedure produces cation saturation like a 72-hour saturation procedure (Figure S2).
Mg-saturated mounts were analyzed after air drying and again after solvation with ethylene glycol (VWR ethylene glycol ≥99%, CAS: 107-21-1) vapor in a closed container maintained at 60-70 ℃ for at least 12 hours35,40. K-saturated oriented mounts were analyzed after heating to 550 ℃ for at least 30 minutes40. Table S2 gives reference d-spacings used to identify the presence of various phyllosilicates within the oriented mount XRD patterns. Smectite (001) planes expand upon solvation with ethylene glycol35 with the degree of expansion largely controlled by layer charge34,35,37,41,42, with typically observed d(001) expansion of low-charge smectites following solvation with ethylene glycol to between ~16 to ~18 Å35,37,43. High charge smectites and vermiculite typically do not exhibit (001) expansion or exhibit extremely limited expansion to between ~14 to ~14.5 Å following ethylene glycol vapor solvation34,35,37,44. Both high- and low-charge smectite (001) planes collapse to ~10 Å following heating to 550℃34–36.
Serpentine mineral (lizardite, antigorite, chrysotile) XRD peaks commonly overlap, and distinguishing whether one or multiple serpentine polymorphs is present from powder XRD patterns can be difficult. Lizardite, antigorite, and chrysotile presence in XRD patterns was interpreted primarily through the methodology of Peacock (1987)7. Identification of antigorite was based on the lack of a asymmetric (020) peaks at ~4.62 Ǻ to 4.57 Ǻ and the presence of a (16 0 1) peak at ~2.53 Å. Although transmission electron microscopy (TEM) analyses detected chrysotile fibers in the clay-size fraction of four soil samples (Figure S), chrysotile identification is difficult and generally not definitive by XRD in when multiple serpentine polymorphs are present and in various states of disorder or incorporate variable amounts of Mg, Fe, Al, Ni, or other elements. Determination from XRD alone may not be definitive when mixtures of these phases are present within a sample7. Additional complications in distinguishing individual serpentine phases arise from the common co-occurrence of the serpentine polymorphs45–47 and from overlap of chlorite and talc peaks in the range of the lizardite and chrysotile (020) peaks and from olivine, which has a major (131) peak at ~2.51 Ǻ7, and from a major goethite peak (111) at ~2.44Ǻ48. While TEM imagery suggests widespread chrysotile presence, it’s presence is not definitive based on XRD alone.
Clay-size fraction material from at least one depth interval within A and B (when present) and C soil horizons from all soil pits was analyzed through Rietveld refinement (Appendix C) of randomly oriented XRD patterns49,50 using Profex 5.133, based on the BGMN refinement architecture51, with structure files from the Crystallography Open Database 32 (see Table S3 for crystal structure files used during refinement). Rietveld refinements were conducted on the clay-size fraction samples spiked with 20 wt.% 0.8 to 1-µm diameter α-Al2O3 (Beantown Chemical 127075-25G). The Crystallography Open Database crystallographic information files used for Rietveld refinement, as well as the origin of files used for smectites, are given in Table S3. For the Rietveld refinements, XRD patterns were analyzed between 5 and 80 ⁰2θ for samples that lacked smectite content, or when samples possessed smectites from 10 to 80 ⁰2θ to exclude the highly variable smectite (001) planes as in Smith et al. (2018)52. Though chrysotile was present in TEM imagery (Figure S49), it was not included in Rietveld refinement due to unavoidable significant overlap between major peaks between different serpentine polymorphs. During fitting, we refined unit cell parameters, scale factor, crystallite size, microstrain, and preferred orientation. Halite resulting from the Na-saturation utilized during the clay-size fraction extraction procedure was observed in several clay-size samples utilized for Rietveld refinement. XRD patterns that exhibited the presence of peaks attributable to halite were cut between 31.6 to 32.2 ⁰2θ, 36.4 to 37.4 ⁰2θ, 52.6 to 54.0 ⁰2θ, and 66.2 to 67.2 ⁰2θ as these represent peaks attributed to halite from comparison with a halite crystallographic information file from the Crystallography Open Database (COD 9006374)53.
Rietveld refinements were utilized to fit a background to each clay-size fraction XRD pattern for the purpose of determining relative amorphous material contributions to each XRD pattern. Rietveld refinement was not utilized to directly calculate X-ray amorphous material abundances via the internal standard method as we could not accurately estimate the density of the various potential amorphous materials present. The prevalence of phyllosilicate phases with highly overlapping peaks and a high potential degree of disorder and chemical variation in octahedral and tetrahedral sites similarly would induce a high degree of uncertainty in any calculated phase abundances. Instead, relative amorphous contributions were determined by calculating a potential amorphous contribution to the fitted Rietveld background as amorphous material is typically part of the background fit of XRD Rietveld refinements52,54. Following fitting of a background during Rietveld refinement, a baseline was calculated for each fitted background utilizing a SNIP function with 100 iterations in the Profex program33. This baseline was then subtracted from the fitted Rietveld background to produce an amorphous hump pattern. The area under the central peak of each amorphous hump was then integrated in OriginPro software. As each XRD samples was analyzed using identical instrument parameters and processing methods, the relative area underneath each amorphous hump contribution to the fitted Rietveld backgrounds provides a relative indication of the abundance of amorphous material within each clay-size fraction.
To ensure the presence of X-ray amorphous material in the samples, we analyzed select samples of clay-size fraction material from the Klamath Mountains and the Tablelands via high-resolution transmission electron microscopy (HRTEM). We analyzed material from the Eunice Bluff soil BC horizon (23-33 cm depth interval) from the Klamath Mountains and the Devil’s Punchbowl soil C horizon (20-35 cm depth interval) from the Tablelands. HRTEM sample preparation was performed as follows. A small amount (~15 mg) of clay-size fraction material was added to ~15 mL of 70% ethanol solution. The solution was then sonicated for 5 minutes to disperse aggregates of clay-size particles. A carbon-coated copper TEM grid was then swirled in the mixture several times and allowed to air dry. Prepared samples were then analyzed using a Titan 300/80 (FEI) aberration-corrected TEM at 300 kV accelerating voltage at Arizona State University. Electron-dispersive spectroscopy results were normalized without Cu or C due to contributions to those elements from the grid material. During our analytical run, the selected area diffraction (SAED) detector exhibited problems and could not be utilized. The presence or absence of amorphous material was determined using fast Fourier transform (FFT) imagery. FFT images produce similar 2D patterns to SAED imagery and can be analyzed in a similar manner. Amorphous material was identified based on the absence of lattice fringes in FFT diffractograms and HRTEM images. Nanocrystalline material, still falling under the X-ray amorphous component of XRD data55,56, was identified based on the presence of crystalline packets less than 100 nm across in HRTEM images with diffuse spots in FFT diffractograms.
4. Chemical Analyses
Bulk soil pH was measured via a 1:1 mass-to-volume ratio of unpowdered <2 mm bulk soil fraction material to 18.2 MΩ water with a Mettler Toledo InLab Expert Pro pH probe57. All samples were dried at 105 ℃ overnight for at least 12 hours to dehydrate the samples using ceramic crucibles that had also previously been dried at 105 ℃ overnight for at least 12 hours. Loss on ignition at 550℃ (LOI550), a proxy for organic carbon content, was applied to unpowdered bulk soil material following the method of Heiri et al. (2001)58. Estimated error for the method is ±2 wt. %58. Briefly, 5 grams of unpowdered air-dried bulk soil material was loaded into oven dried ceramic crucibles, heated to 105℃ overnight for at least 12 hours to dehydrate the samples, cooled and weighed, heated to 550℃ for 4 hours to burn off organic carbon, then cooled and weighed again. Loss on ignition at 950℃ (LOI950) as in Jones et al. (2000)59 was applied to all clay-size fraction and powdered parent material samples to examine total volatile content. 1 gram of freeze-dried clay-size fraction material or powdered parent material oven dried at 105℃ overnight for at least 12 hours to ensure dehydration was heated to 950 ℃ for 2 hours, allowed to cool overnight, then weighed.
To examine trends in crystallinity of Fe-containing materials within each soil site as well as the bulk chemical composition of the parent bedrock and clay-size fraction material, a suite of total and selective chemical dissolutions was applied to material from all sampled depth intervals in all soil pits (depth intervals and soil horizon designations are given in Appendix C).
ICP-MS analyses for bulk major and trace element chemistry of the total digestions of the powdered parent material and clay-size fraction samples were measured at the Nevada Plasma Facility Lab (NPFL) at UNLV using a ThermoFisher ScientificTM iCAP Qc ICP-MS after the following digestion— 50 mg of sample material was dissolved in acid cleaned Savillex beakers in a 1:1 mixture of sub-boiled double distilled HNO3 and Optima grade HF. Samples were analyzed concurrently with the USGS standard BHVO-1, AGV-1, and BCR-1, and total procedural blanks. Reproducibility and detailed methods of NPFL data is presented in DeFelice et al. (2019)60. Off gassing of SiF4 during dissolution of silicates in HF precludes accurate direct measurements of SiO2 in our samples61. SiO2 content was instead calculated by subtracting the total major element oxides measured by ICP-MS from 100% as in Udry et al. (2019)62 after incorporating the volatile content determined by loss on ignition at 950 ℃ from 100% as in Jones et al. (2000)59.
Na2O was excluded from ICP-MS measurements of the clay-size fraction chemistry as the clay-size fraction extraction procedure involved the addition of NaCl to promote flocculation. Na2O content within parent material samples was measured (see Appendix B) to estimate whether exclusion of Na2O would substantially alter calculated SiO2 values, and minimal Na2O content is present in the parent material samples for Pickhandle Gulch (0.00 wt.%), the Tablelands (0.02 wt.%), and for the Eunice Bluff (0.06 wt.%), Swift Creek (0.01 wt.%), and String Bean Creek (0.02 wt.%) sampling sites in the Klamath Mountains. Na2O content was slightly higher within the Deadfall Lake parent material (1.40 wt.%) because of the presence of plagioclase, likely including albite, not observed in the bedrock at other sites and likely not contributing to this extent to the soil formed at Deadfall lake based on the decreased Al present in the clay-size fraction from Deadfall lake relative to the parent material, and the presence of olivine in the clay size fraction despite the absence of olivine in the parent material (Figure S3 and Table 1). As SiO2 content is calculated based on the measured major element oxide concentrations not including Na2O it should be read as a potential maximum value and is likely very accurate for the soil pits with minimal parent material Na2O values, but potentially slightly over-estimates the SiO2 content for the Deadfall Lake soil. Additionally, if plagioclase found within some soils included a significant albite component this would also lead to overestimates of the SiO2 content.
Selective chemical dissolution techniques were applied to powdered bulk soil samples, oven dried at 105⁰C for 12 hours to ensure dehydration, from all sampled depth intervals from all soil pits (Appendix C) to examine the distribution of secondary Fe within different soil reservoirs. Hydroxylamine hydrochloride, citrate dithionite, and sodium pyrophosphate extractants were applied to study Fe within amorphous material63,64, amorphous and crystalline (oxyhydr)oxides57,64,65, and organically bound64,66 Fe-reservoirs respectively. The hydroxylamine hydrochloride supernatant was also analyzed for Si content.
A sodium pyrophosphate extraction was applied to analyze organically complexed Fe (FeP) following the method delineated within Shang and Zelazny (2008)62 a modification of the method used by McKeague (1967)64. Whether the pyrophosphate technique results in dissolution of purely organic compounds or whether the supernatant potentially includes some Fe from colloidal Fe-(oxyhydr)oxide particles has been called into question65. Fe measured by the pyrophosphate technique thus represents a potential maximum amount of organically complexed Fe. Briefly, 0.1M Na4P2O7∙10H2O (Alfa Aesar ACS grade #33385-22) was brought to pH 10 with Aristar Plus trace metal grade HCl (VWR chemicals #87003-259). 30 mL of the pyrophosphate extractant was added to 300 mg of 105⁰C oven dried powdered bulk soil material and the mixture was shaken overnight for 16 hours at 180 rpm (Eberbach Mod. 6010 shaker table). The samples were then centrifuged at 15,550 x g (10,980 rpm) for 10 minutes and the supernatant decanted from overtop the soil pellet and filtered through a 0.2 µm filter (VWR #28145-499). Samples were analyzed for Fe within five days of the extraction procedure by atomic absorption spectroscopy on a Thermo iCE 3000 Flame Atomic Absorption Spectrometer. The sodium pyrophosphate results are presented in the supplementary online material (Appendix B).
The hydroxylamine hydrochloride extraction to analyze Fe and Si within amorphous material was applied following the methodology of Shang and Zelazny (2008)64, a modification of the methods used by Ross et al. (1985)63 and Chao and Zhou (1983)67. The hydroxylamine hydrochloride technique acts as a combined reduction and acid protonation technique and is selective for Fe and Si containing amorphous material63,64,67–69. The hydroxylamine hydrochloride technique was chosen rather than the ammonium oxalate in darkness (AOD) method as it extracts Fe in comparable amounts from amorphous material63,64, because the hydroxylamine technique does not attack Fe within magnetite or chlorites and exhibits minimal potential for dissolving crystalline Fe-(oxyhydr)oxides such as lepidocrocite63,64,67, and because the hydroxylamine supernatant produces minimal interference during absorbance measurements compared to the oxalate supernatant64. 25 mL of 0.25 M ACS grade NH2OH∙HCl (Beantown Chemical Catalogue #144935) + 0.25 M Aristar Plus trace metals grade HCl (VWR chemicals #87003-253) solution was added to 100 mg of 105⁰C oven dried powdered bulk soil material from each sampled depth interval (Appendix C) within acid washed 50 mL centrifuge tubes. Samples were shaken for 16 hours at 180 rpm at room temperature (Eberbach Mod. 6010 shaker table) and the supernatant was then filtered through a sterile 0.2 µm filter (VWR #28145-499). The filtered supernatant samples were acidified to 2% v/v with trace metals grade HNO3 (VWR chemicals Aristar Plus #87003-259). The samples were then diluted 10x using a 2% v/v trace metals grade HNO3 in 18.2 MΩ DI-water solution and analyzed for Fe and Si content via atomic absorption spectroscopy on a Thermo iCE 3000 Flame Atomic Absorption Spectrometer.
The citrate dithionite extraction was applied to examine Fe within both amorphous and crystalline (oxyhydr)oxides64,65 based on the procedure of the 2004 Soil Survey Laboratory Manual57. The dithionite acts as a reducing agent for Fe3+ in Fe3+-containing amorphous and crystalline materials at near neutral pH conditions64,65,70–72. Some work has suggested that in addition to reduction of Fe from within oxides and oxyhydroxides that reduction with dithionite extracts minor amounts of Fe3+ from vermiculite interlayers73–75 and limited but measurable quantities of structural Fe from within smectites64,70–72. Briefly, 25 mL of 0.57 M Na3C6H507·2H2O, 400 mg of Na2S2O4, and 750 mg of 105⁰C oven-dried powdered bulk soil material from each sampled soil depth interval (Appendix C) were combined within acid-washed 50 mL centrifuge tubes and shaken overnight for 12-16 hours at 180 rpm (Eberbach Mod. 6010 shaker table), then agitated vigorously by hand for approximately 30 seconds to ensure suspension of soil particles in solution and allowed to incubate without shaking overnight at room temperature for 12-16 hours. Samples were then centrifuged at 4,000 rpm (2,057 x g) for 15 minutes and the supernatant was decanted from the soil pellet and filtered through a 0.2 µm filter (VWR #28145-499). Samples were acidified to 2% v/v with trace metal grade Aristar Plus HNO3 (VWR chemicals #87003-259), diluted 100x using a solution of 2% v/v trace metal grade HNO3 and 18.2 MΩ DI water, and then analyzed for Fe content via atomic absorption spectroscopy on a Thermo iCE 3000 Flame Atomic Absorption Spectrometer.
References
- Arguez, A. et al. NOAA’s U.S. Climate Normals (1981-2020). NOAA National Centers for Environmental Information vol. 10 Preprint at https://doi.org/10.7289/V5PN93JP (2010).
- Skinner, C. N., Taylor, A. H. & Agee, J. K. Klamath Mountains bioregion. Fire in California’s Ecosystems 170–194 (2006).
- Garwood, J. M. & Welsh, H. H. Ecology of the Cascades Frog (Rana cascadae) and Interactions with Garter Snakes and Nonnative Trout in the Trinity Alps Wilderness, California. (2007).
- Mohr, J. A., Whitlock, C. & Skinner, C. N. Postglacial vegetation and fire history, eastern Klamath Mountains, California, USA. Holocene 10, 587–601 (2000).
- Environment Canada. Canadian Climate Normals 1981-2010 Station Data. Cow Head. Climate ID: 8401335. Preprint at (2010).
- Quick, J. E. The origin and significance of large, tabular dunite bodies in the Trinity peridotite, northern California. Contributions to Mineralogy and Petrology 78, 413–422 (1982).
- Peacock, S. M. Serpentinization and infiltration metasomatism in the Trinity peridotite, Klamath province, northern California: implications for subduction zones. Contributions to Mineralogy and Petrology 95, 55–70 (1987).
- Snoke, A. W. & Barnes, C. G. The development of tectonic concepts for the Klamath Mountains province, California and oregon. Special Paper of the Geological Society of America 410, 1–29 (2006).
- Wallin, E. T. & Metcalf, R. V. Supra-subduction zone ophiolite formed in an extensional forearc: Trinity terrane, Klamath Mountains, California. Geology 106, 591–608 (1998).
- Wright, J. E. & Wyld, S. J. The Rattlesnake Creek Terrane, Klamath Mountains, California: an early Mesozoic volcanic arc and its basement of tectonically disrupted oceanic crust. Geol Soc Am Bull 106, 1033–1056 (1994).
- Irwin, W. P. Terranes of the western Paleozoic and Triassic belt in the southern Klamath Mountains, California. in U.S. Geological Survey Professional Paper 800-C C103–C111 (1972).
- Sharp, R. P. Pleistocene Glaciation in the Trinity Alps of Northern California. Am J Sci 258, 305–340 (1960).
- Daniels, M. L., Anderson, S. & Whitlock, C. Vegetation and fire history since the Late Pleistocene from the Trinity Mountains, northwestern California, USA. Holocene (2005) doi:10.1191/0959683605hl878ra.
- Dickey, N. W. Chronology and Paleoclimate of Late Pleistocene Glaciation in The Klamath Mountains, CA. (Colorado State University, Northridge, 2016).
- Blaise, B., Clague, J. J. & Mathewes, R. W. Time of maximum Late Wisconsin glaciation, West Coast of Canada. Quat Res 34, 282–295 (1990).
- Anderson, P. M. & Lozhkin, A. The stage 3 interstadial complex (Karginskii/middle Wisconian interval) of Beringia: Variations in paleoenvironments and implications for paleoclimatic interpretations. Quat Sci Rev 20, 93–125 (2001).
- Dredge, L. A. & Thorleifson, H. L. The Middle Wisconsinan History of the Laurentide Ice Sheet. Quat Sci Rev 41, 215–235 (1987).
- Brookes, I. A. Table Mountain, Gros Morne National Park, Newfoundland. The Canadian Geographer 37, 69–75 (1993).
- Williams, H., Cawood, P. A. & Dohar, V. Geology, Humer Arm Allochthon, Newfoundland. (Geological Survey of Canada, 1989).
- Suen, C. J., Frey, F. A. & Malpas, J. Bay of Islands Ophiolite Suite, Newfoundland Petrologic and Geochemical Characteristics with Emphasis on Rare Earth Element Geochemistry. Earth Planet Sci Lett 45, 337–348 (1979).
- Suhr, G. & Cawood, P. A. Structural history of ophiolite obduction, Bay of Islands, Newfoundland. Geol Soc Am Bull 105, 399–410 (1993).
- Osborn, G., Spooner, I., Gosse, J. & Clark, D. Alpine glacial geology of the Tablelands, Gros Morne National Park, Newfoundland. Can J Earth Sci 44, 819–834 (2007).
- Shaw, J. Submarine moraines in Newfoundland coastal waters: Implications for the deglaciation of Newfoundland and adjacent areas. Quaternary International 99–100, 115–134 (2003).
- Stea, R. R., Piper, D. J. W., Fader, G. B. J. & Boyd, R. Wisconsinan glacial and sea-level history of Maritime Canada and the adjacent continental shelf: a correlation of land and sea events. Bulletin of the Geological Society of America 110, 821–845 (1998).
- Page, B. M. Geology of the Candelaria Mining District, Mineral County, Nevada. Nevada Bureau of Mines Bulletin 56 (1959).
- Istok, J. D. D. & Harward, M. E. E. Influence of Soil Moisture on Smectite Formation in Soils Derived from Serpentinite. Soil Science Society of America Journal 46, 1106–1108 (1982).
- Birkeland, P. W. Soils and Geomorphology. (Oxford University Press, 1984).
- Soil Science Division Staff. Soil Survey Manual. Government Printing Office (2017). doi:10.1097/00010694-195112000-00022.
- Edwards, A. P. & Bremner, J. M. Dispersion of Soil Particles By Sonic Vibration. Journal of Soil Science 18, 47–63 (1967).
- Degen, T., Sadki, M., Bron, E., König, U. & Nénert, G. The high score suite. in Powder Diffraction vol. 29 (2014).
- Lafuente, B., Downs, R. T., Yang, H. & Stone, N. The power of databases: The RRUFF project. Highlights in Mineralogical Crystallography (2016). doi:10.1515/9783110417104-003.
- Graulis, S. et al. Crystallography Open Database - An open-access collection of crystal structures. J Appl Crystallogr 42, (2009).
- Doebelin, N. & Kleeberg, R. Profex: A graphical user interface for the Rietveld refinement program BGMN. J Appl Crystallogr 48, 1573–1580 (2015).
- Harward, M. E., Carstea, D. D. & Sayegh, A. H. Properties of vermiculites and smectites: Expansion and collapse. Clays Clay Miner 16, 437–447 (1969).
- Moore, D. M. & Reynolds, R. C. X-ray Diffraction and the Identification and Analysis of Clay Minerals. (Oxford University Press, 1997).
- Lee, B. D., Sears, S. K., Graham, R. C., Amrhein, C. & Vali, H. Secondary Mineral Genesis from Chlorite and Serpentine in an Ultramafic Soil Toposequence. Soil Science Society of America Journal 67, 1309–1317 (2003).
- Mosser-Ruck, R., Devineau, K., Charpentier, D. & Cathelineau, M. Effects of ethylene glycol saturation protocols on XRD patterns: A critical review and discussion. Clays Clay Miner 53, 631–638 (2005).
- Harris, W. & White, G. N. X-ray diffraction techniques for soil mineral identification. in Methods of Soil Analysis, Part 5: Mineralogical Methods vol. 5 (2015).
- Whittig, L. D. & Allardice, W. R. X-ray diffraction techniques. in X-ray Diffraction Techniques 331–362 (wiley, 1986). doi:10.2136/sssabookser5.1.2ed.c12.
- Poppe, L. J., Paskevich, V. F., Hathaway, J. C. & Blackwood, D. S. A Laboratory Manual For X-Ray Powder Diffraction. US Geological Survey Open-File Report 1.041 1–88 (2001).
- Szczerba, M., Kłapyta, Z. & Kalinichev, A. Ethylene glycol intercalation in smectites. Molecular dynamics simulation studies. Appl Clay Sci 91–92, 87–97 (2014).
- Sato, T., Watanabe, T. & Otsuka, R. Effects of layer charge, charge location, and energy change on expansion properties of dioctahedral smectites. Clays Clay Miner 40, 103–113 (1992).
- Christidis, G. E. & Eberl, D. D. Determination of layer-charge characteristics of smectites. Clays Clay Miner 51, 644–655 (2003).
- Vali, H. & Hesse, R. Identification of vermiculite by transmission electron microscopy and X-ray diffraction. Clay Miner 27, 185–192 (1992).
- Evans, B. W. The Serpentinite Multisystem Revisited: Chrysotile Is Metastable. Int Geol Rev 46, (2004).
- Evans, B. W., Hattori, K. & Baronnet, A. Serpentinite: What, why, where? Elements 9, 99–106 (2013).
- Baumeister, J. L. et al. Biogeochemical weathering of serpentinites: An examination of incipient dissolution affecting serpentine soil formation. Applied Geochemistry 54, 74–84 (2015).
- Nagai, T., Kagi, H. & Yamanaka, T. Variation of hydrogen bonded O...O distances in goethite at high pressure. American Mineralogist 88, 1423–1427 (2003).
- Rietveld, H. M. A profile refinement method for nuclear and magnetic structures. J Appl Crystallogr 2, 65–71 (1969).
- Bish, D. L. & Post, J. E. Quantitative mineralogical analysis using the Rietveld full-pattern fitting method. American Mineralogist 78, 932–940 (1993).
- Bergmann, J., Friedel, P. & Kleeberg, R. BGMN—a new fundamental parameters based Rietveld program for laboratory X-ray sources, its use in quantitative analysis and structure investigations. IUCr Commision on Powder Diffraction Newsletter 20, 5–8 (1998).
- Smith, R. J., Rampe, E. B., Horgan, B. H. N. & Dehouck, E. Deriving amorphous component abundance and composition of rocks and sediments on Earth and Mars. J Geophys Res Planets 123, 2485–2505 (2018).
- Wang, K. & Reeber, R. R. Thermal expansion of alkali halides at high pressure: NaCI as an example. Phys Chem Minerals vol. 23 (1996).
- Klug, H. P. & Alexander, L. E. X-ray Diffraction Procedures: for Polycrystalline and Amorphous Materials. (Wiley, 1974).
- Smith, R. J. & Horgan, B. H. N. Nanoscale Variations in Natural Amorphous and Nanocrystalline Weathering Products in Mafic to Intermediate Volcanic Terrains on Earth: Implications for Amorphous Detections on Mars. J Geophys Res Planets 126, 1–30 (2021).
- Holder, C. F. & Schaak, R. E. Tutorial on Powder X-ray Diffraction for Characterizing Nanoscale Materials. ACS Nano vol. 13 7359–7365 Preprint at https://doi.org/10.1021/acsnano.9b05157 (2019).
- Soil Survey Staff. Soil Survey Laboratory Methods Manual. Soil Survey Investigations Report 42, 700 (2004).
- Heiri, O., Lotter, A. F. & Lemcke, G. Loss on ignition as a method for estimating organic and carbonate content in sediments: reproducibility and comparability of results Oliver. J Paleolimnol 25, 101–110 (2001).
- Jones, R. C., Babcock, C. J. & Knowlton, W. B. Estimation of the Total Amorphous Content of Hawai’i Soils by the Rietveld Method. Soil Science Society of America Journal 64, 1100–1108 (2000).
- DeFelice, C., Mallick, S., Saal, A. E. & Huang, S. An isotopically depleted lower mantle component is intrinsic to the Hawaiian mantle plume. Nat Geosci 12, 487–492 (2019).
- Sparks, D. L., Page, A. L., Helmke, P. A. & Loeppert, R. H. Methods of Soil Analysis, part 3: Chemical methods (vol 14). in (2020).
- Udry, A. et al. Reclassification of four aubrites as enstatite chondrite impact melts: Potential geochemical analogs for Mercury. Meteorit Planet Sci 54, 785–810 (2019).
- Ross, G. J., Wang, C. & Schuppli, P. A. Hydroxylamine and Ammonium Oxalate Solutions As Extractants For Iron and Aluminum From Soils. Soil Science Society of America Journal 49, 783–785 (1985).
- Shang, C. & Zelazny, L. W. Selective Dissolution Techniques for Mineral Analysis of Soils and Sediments. in Methods of Soil Analysis Part 5 - Mineralogical Methods (eds. ULERY, A. L. & RICHARD DREES, L.) 33–80 (Soil Science Society of America, 2008). doi:10.2136/sssabookser5.5.c3.
- Holmgren, G. G. S. A Rapid Citrate-Dithionite Extractable Iron Procedure. Soil Science Society of America Journal 31, 210–211 (1967).
- McKeague, J. A. An Evaluation of 0.1 M Pyrophosphate and Pyrophosphate-Dithionite in Comparison With Oxalate as Extractancs of the Accumulation Products in Podzols and Some Other Soils. Can J Soil Sci 47, 95–99 (1967).
- Chao, T. T. & Zhou, L. Extraction Techniques for Selective Dissolution of Amorphous Iron Oxides from Soils and Sediments. Soil Science Society of America Journal 47, 225–232 (1983).
- Chester, R. & Hughes, M. J. A CHEMICAL TECHNIQUE FOR THE SEPARATION OF FERRO-MANGANESE MINERALS, CARBONATE MINERALS AND ADSORBED TRACE ELEMENTS FROM PELAGIC SEDIMENTS. Chem Geol 2, 249–262 (1967).
- Chao, T. T. Selective Dissolution of Manganese Oxides from Soils and Sediments with Acidified Hydroxylamine Hydrochloride. Soil Science Society of America Journal 36, 764–768 (1972).
- Stucki, J. W., Golden, D. C. & Roth, C. B. Preparation and Handling of Dithionite-Reduced Smectite Suspensions. Clays Clay Miner 32, 191–197 (1984).
- Russell, J. D. Infrared and Mössbauer Studies of Reduced Nontronites. Clays Clay Miner 27, 63–71 (1979).
- Ericsson, T., Linares, J. & Lotse, E. A Mossbauer Study of the Effect of Dithionite/Citrate/Bicarbonate Treatment on a Vermiculite, A Smectite, and a Soil. Clay Miner 19, 85–91 (1984).
- Roth, C. B., Jackson, M. L. & Syers, J. K. Deferration Effect on Structural Ferrous-Ferric Iron Ratio and CEC of Vermiculites and Soils. Clays and Clays Minerals 17, 253–264 (1969).
- Roth, C. B., Jackson, M. L., Lotse, E. G. & Syers, J. K. Ferrous-Ferric Ratio and CEC Changes on Deferration of Weathered Micaceous Vermiculite. Isr J Chem 6, 261–273 (1968).
- Veith, J. A. & Jackson, M. L. Iron oxidation and reduction effects on structural hydroxyl and layer charge in aqueous suspensions of micaceous vermiculites. Clays Clay Miner 22, 345–353 (1974).