Soil Phosphorus Speciation and Availability in Meadows and Forests in Alpine Lake Watersheds With Different Parent Materials

Heron, Thomas; Strawn, Daniel G.; Dobre, Mariana; Cade-Menun, Barbara J.; Deval, Chinmay; Brooks, Erin S.; Piaskowski, Julia; Gasch, Caley; Crump, Alex

doi:10.3389/ffgc.2020.604200

ORIGINAL RESEARCH article

Front. For. Glob. Change, 23 February 2021

Sec. Forest Soils

Volume 3 - 2020 | https://doi.org/10.3389/ffgc.2020.604200

This article is part of the Research Topic Changes in Forest Ecosystem Nutrition View all 25 articles

Soil Phosphorus Speciation and Availability in Meadows and Forests in Alpine Lake Watersheds With Different Parent Materials

$\nThomas Heron$ Thomas Heron¹

Daniel G. Strawn¹^*

Mariana Dobre¹

Barbara J. Cade-Menun²

Chinmay Deval¹

Erin S. Brooks¹

Julia Piaskowski¹

Caley Gasch³

Alex Crump¹

¹College of Agricultural and Life Sciences, University of Idaho, Moscow, ID, United States
²Agriculture and Agri-Food Canada, Swift Current Research and Development Centre, Swift Current, SK, Canada
³Department of Soil Science, North Dakota State University, Fargo, ND, United States

In the Lake Tahoe Basin in California and Nevada (USA), managing nutrient export from watersheds into streams and the lake is a significant challenge that needs to be addressed to improve water quality. Leaching and runoff of phosphorus (P) from soils is a major nutrient source to the lake, and P loading potential from different watersheds varies as a function of landscape and ecosystem properties, and how the watershed is managed. In this research, P availability and speciation in forest and meadow soils in the Lake Tahoe Basin were measured at two watersheds with different parent material types. Soils developed on andesitic parent materials had approximately twice as much total P compared to those developed on granitic parent materials. Regardless of parent material, organic P was 79–92% of the total P in the meadow soils, and only 13–47% in the forest soils. Most of the soil organic P consisted of monoester P compounds, but a significant amount, especially in meadow soils, was diester P compounds (up to 30% of total extracted P). Water extractable P (WEP) concentrations were ~10 times greater in the granitic forest soils compared to the andesitic forest soils, which had more poorly crystalline aluminosilicates and iron oxides that retain P and thus restrict WEP export. In the meadow soils, microbial biomass P was approximately seven times greater than the forest soils, which may be an important sink for P leached from upland forests. Results show that ecosystem and parent material are important attributes that control P speciation and availability in the Lake Tahoe Basin, and that organic P compounds are a major component of the soil P and are available for leaching from the soils. These factors can be used to develop accurate predictions of P availability and more precise forest management practices to reduce P export into Lake Tahoe.

Introduction

Lake Tahoe, located in the Sierra Nevada Mountain range in California and Nevada, is the sixth largest lake by volume in the United States. It is classified as an ultra-oligotrophic lake, meaning that it has naturally low nutrient concentrations and low primary production, and it is renowned for the clarity of its water (Hatch et al., 2001; Goldberg et al., 2015). In recent years, however, water clarity in Lake Tahoe has declined, with Secchi depth readings decreasing from ~31 m in 1968 to 21.6 m in 2018 (Schladow, 2019). As a result of non-point-source nutrient loading, primary production in Lake Tahoe has increased by ~6% per year (Jassby et al., 1999; Roberts and Reuter, 2010).

Historically, algae growth in Lake Tahoe has been co-limited by nitrogen (N) and phosphorus (P) (Hatch et al., 1999). However, as has been observed for lakes worldwide (Elser et al., 2009), increased atmospheric N loading and N deposition have altered plankton species and N:P stoichiometry, shifting nutrient limitation in Lake Tahoe to P (Hatch et al., 1999; Goldberg et al., 2015). A recent lake-clarity model demonstrated that a return to the historical Secchi depth reading in Lake Tahoe would be possible within 20 years if P loading were reduced by at least 2.75% per year (Sahoo et al., 2010). However, to control P sources and subsequent loading into surface waters, a full understanding of P cycling and species in soils in the Lake Tahoe Basin is required.

The physical forms of P that can enter and cycle in lakes are defined as particulate P that is >0.45 μm, and solution P that can pass through a 0.45-μm filter that consists of dissolved and colloidal P (Bol et al., 2016). Colloidal P particles are 1 to 1,000 nm in diameter and the colloids less than 450 nm can pass through a 0.45 μm filter (Jiang et al., 2017), and can thus be mobilized through soils and remain suspended in surface waters. Dissolved, colloidal, or particulate P species can be organic [bound to a carbon (C) group] or inorganic (singular or multiple phosphate groups). Short-term changes in Lake Tahoe primary productivity are well-explained by dissolved inorganic and organic P loads from Lake Tahoe Basin streams, which contribute up to 1,000 kg of dissolved P annually (Hatch et al., 1999). Different forms of P, both physically and chemically, differ in their mobility, environmental reactivity, and bioavailability. Dissolved molybdate-reactive phosphate (MRP), also known as soluble reactive phosphate (SRP), is the most readily bioavailable form (Hatch et al., 1999; Sahoo et al., 2010). Thus, to manage Lake Tahoe Basin landscapes for P-load reduction requires knowledge of soil P species and pools, and their potential for release and transport into surface waters.

About 6% of the Lake Tahoe Basin is considered urban and has been developed for residential and commercial use. Phosphorus inputs within the urban portion of the Basin account for 18% of total P inputs into the lake (Sahoo et al., 2013). In the non-urban regions of the Basin, 50% are covered by forests with yellow pine associations (TRPA, 2015) containing Jeffrey pine (Pinus jeffreyi Balf.), white fir [Abies concolor (Gord. and Glend.) Lindl. Ex Hildebr.], incense cedar [Calocedrus decurrens (Torr.) Florin], and sugar pine (Pinus labertiana Douglas); and another 17% of the Tahoe Basin landscape are red fir associations containing red fir (Abies magnifica A. Murray bis), Jeffrey pine and Lodgepole pine (Pinus contorta Douglas ex Loudon). Wet and dry meadows comprise 2 and 1% of non-urban land, respectively, and contain grasses, sedges, and rushes (TRPA, 2015). Soils in the Lake Tahoe Basin are developed on andesitic, granitic, or mixed parent materials (Coats et al., 2016). Studies in the Eastern Sierra Nevada have shown that forest soils on granitic parent materials can have substantially higher extractable-P concentrations than those developed on andesitic parent materials (Johnson et al., 1997; Coats et al., 2016). However, it is unclear the mechanistic processes that are responsible for these differences.

In soils, inorganic and organic P molecular species have distinct potential for uptake by vegetation or mobilization out of the soil profile. Uselman et al. (2012) suggested that the amount of dissolved organic P in soil solution is largely dependent on the type and amount of above- and below-ground organic matter. Some forest-soil P is exported as particulate P (Prairie and Kalff, 1988). The amount of eroded particulate P that is exported from a site depends on three factors: (1) site geography (slope, climate, and geology); (2) site management (harvest, thinning, and development); and (3) wildfire history (Miller et al., 2006). Particulate P that enters streams and lakes is not directly available for uptake by aquatic organisms, although it can be released as dissolved P from the particles and then is bioavailable to support aquatic algae growth (Young et al., 1985; Reid et al., 2018).

A significant fraction of forest P exists in the plant litter and O horizons that can be illuviated into lower depths in the soil profile or be lost in runoff (Miller et al., 2010; Bol et al., 2016). Miller et al. (2005, 2006) observed that organic horizons on forest floors in Lake Tahoe Basin have high levels of water-soluble P that may be a source of P loading to streams via overland or subsurface flow, the latter of which moves through and reacts with soils. Phosphorus leached from O horizons can be transported into the soils through several mechanisms, depending on the soil physical properties that facilitate preferential vs. matrix flow (Julich et al., 2017; Luo et al., 2019). In alpine environments, spring snowmelt runoff is an important mechanism of P loading because it transports dissolved, colloidal, and particulate P from decomposed forest litter and soils, which can then emerge as subsurface P loading to streams and lakes (Backnäs et al., 2012).

Estimates suggest that groundwater sources make up 15% by mass of total P loading to Lake Tahoe (Roberts and Reuter, 2010). Furthermore, 61% (3,700 kg) of the annual total dissolved P that is found in Lake Tahoe Basin groundwater is believed to be derived from natural sources from unimpacted non-urban areas, predominantly from overlying forest litter P pools and P released from the adsorbed and mineral-bound soil P pool (U.S. Army Corp of Engineers, 2003). Sohrt et al. (2019) used an end-member mixing model that included soil water input to predict that up to 92% of the stream P in a mixed deciduous/evergreen forest in Europe was leached from the mineral soil horizons. Considering the hydrologic interfaces in the soil, it follows that in the Lake Tahoe Basin forests, soil P biogeochemistry is an important factor that controls P discharge into surface waters, and to reduce P loads released from forests it is imperative to understand how site and management factors impact P solubility and mobility.

Although meadows comprise a small area of the Lake Tahoe Basin watershed, they are important controllers of P entering streams because they are transitional zones connecting terrestrial and aquatic ecosystems and are commonly located adjacent to forests (Roby et al., 2015). Some meadows in the Lake Tahoe Basin are categorized as stream environment zones (SEZ), which is a designation used by the Lake Tahoe Basin Management Unit for an area of high value and management priority based on ecosystem services, including the filtering and storage of nutrients in runoff (Roby et al., 2015). Forest-derived P is commonly hydrologically transported through meadow ecosystems, which can act as either sinks that intercept P or sources that release P to streams and lakes. Several groups have studied the capacity of riparian systems to perform these functions. For example, Casey and Klaine (2001) studied P adsorption behavior in meadow soils, including Cumulic Humaquepts (similar taxa are found in some Lake Tahoe Basin meadows), and demonstrated the importance of sorption capacity as a mechanism of nutrient attenuation. They determined that soil P concentration was 100 times below concentrations that would cause soil solution P levels to exceed U.S. Environmental Protection Agency recommendations for lentic waters. In contrast, Hoffmann et al. (2006) found a net loss of P via leaching from soils in riparian meadows during two of three sampling years. In a later study, Hoffmann et al. (2009) concluded that, although sedimentation in riparian buffers is an important mechanism of P retention, these buffers may eventually become significant sources of dissolved reactive P release to surface or groundwater. Gergans et al. (2011) studied nutrient flow through a Lake Tahoe watershed that included a riparian meadow ecosystem and observed that the meadow soils were sources of phosphate into a nearby stream and that the release varied with season. The contrasting reports of nutrient retention and release from meadows highlight the complex nature of meadow biogeochemical processes that can make them either sources or sinks of P into surface waters.

Organic P species in forest soils have been shown to be a dominant loading factor to surface waters (Condron et al., 2005; Sohrt et al., 2017). Backnäs et al. (2012) observed higher soluble organic P (labile monoester and diester P species) in surface horizons of Podzol soils in a mixed-coniferous forest in Finland compared to deeper soils. Anderson and Magdoff (2005) observed higher levels of labile organic P than inorganic P in leachate from packed soil columns leached with DNA (diester P) and orthophosphate solutions. Missong et al. (2016) separated bulk soil extractions from forest soils into colloidal and electrolytic fractions and found most of the extractable P was organic bound P (diesters) on colloids. Brödlin et al. (2019b) studied P forms in soils from three different parent materials in deciduous forests and observed a tendency for organic P to dominate mobilized dissolved P. Bol et al. (2016) reviewed organic P in forested soils and concluded that, although it is a significant component of P cycling, the lack of knowledge of organic P species creates a “blind spot in ecosystem research.” Therefore, the dynamics and vulnerability of P leached from both forest and adjacent meadow soils needs to be investigated to understand the potential impact on water quality. This is especially true in watersheds like in the Lake Tahoe Basin, where nutrients leached through forest soils are major inputs into the lake.

In this paper, we investigated the influence of parent material and ecosystem type on soil P species and solubility in the Lake Tahoe Basin. We hypothesized that there would be distinct P biogeochemistry in forest and meadow ecosystems, and that granitic and andesitic parent materials would influence total and available soil P, as well as the amount and type of organic P. Speciation of P in the soils was determined by extraction and ³¹P nuclear magnetic resonance spectroscopy (P-NMR) to elucidate organic and inorganic P species. Water-extractable, exchangeable, and microbial-biomass bound soil P were measured to determine soil P fractions that are potentially soluble and labile. These extractions are good predictors of P immobilization and potential runoff from soils (Pote et al., 1996; Campo et al., 1998; Vadas et al., 2005; Wang et al., 2010; Pistocchi et al., 2018).

Methods

Study Sites and Sample Collection

Soils were sampled from two subalpine meadow systems (Paige Meadow and Meeks Meadow) and their adjacent forests. The research sites are located on the west shore of the Lake Tahoe Basin (Figure 1). Paige Meadow is an alluvial floodplain surrounded by forested hillslopes of terminal moraines at elevation ~2,115 m. Meeks Meadow is situated in an elongated glacial valley trough floodplain (elevation ~1,905 m), confined on both sides by steep forested hillslopes of lateral moraines. Separate lobes of the Sierran Ice Cap extended over the present-day locations of both meadows, carving out their current floodplain topographic environments (Ehlers and Gibbard, 2003). At Paige Meadow and its surrounding forest, soils developed on glacial deposits of eroded basaltic and andesitic rocks from Miocene- through Pleistocene-age volcanic activity (Kortemeier et al., 2018). The geologic substrate of the Meeks Meadow watershed is primarily granodiorite eroded from a glacial drift of till and outwash (Saucedo, 2005). Both meadows contain perennial grasses mixed with sedges, rushes, and forbs (Soil Survey Staff, 2007). The forest surrounding Paige Meadow is a red fir forest association, while the forest next to Meeks Meadow consists of a yellow pine association (Soil Survey Staff, 2007; TRPA, 2015). Climatic data from Tahoe City and Rubicon SNOTEL stations show approximate cumulative precipitation of 900 mm and a mean annual temperature of 7.5°C (USDA-NRCS, 2019). Soils at nearby SNOTEL stations fall under a xeric soil moisture and frigid soil temperature regimes. At each location, 3–5 soil profiles were viewed to 20–40 cm depth and characterized using either a shovel or corer. The descriptions of the soils were done following USDA NRCS soil description methods (Schoeneberger et al., 2012).

FIGURE 1

Figure 1. Andesitic watershed at Paige Meadow (a) and granitic watershed at Meeks Meadow (b). Solid lines are watershed boundaries. Dashed lines separate forest soils from meadow soils, based on map units from the SSURGO database. Circles are meadow sampling sites and triangles are forest sampling sites.

Eight locations from each ecosystem-parent material type were sampled in June, August, and October of 2018 (Figure 1). Samples were taken from the top 15 cm of the soil below the O horizon using a 10-cm diameter soil auger. At each of the eight replicate locations, a composite sample was collected by mixing three sub-samples from 1 m radius. After sampling, the soils were stored on ice while in transport to the lab. A portion of each sample was sieved (<2 mm) and oven-dried at 50°C, and the remainder was stored undried in re-sealable plastic bags at 4°C and sieved (<2 mm) immediately prior to analysis. At each site seven to fourteen 7-cm diameter by 15-cm depth intact cores were sampled for bulk density measurement. O-horizon samples were taken by compositing material from several locations at each forest site into a bag that was thoroughly mixed. For P-NMR analysis and P availability extractions, field-moist samples were used, and P concentrations were adjusted using the percent moisture content determined by the difference in mass of water between the field-moist and oven-dried samples. Other analyses used 50°C oven-dried soils.

Laboratory Analyses

Soil Characterization

Replicate samples from the June 2018 sampling were analyzed for pH, percent sand, and concentrations of total organic C (TOC), total N (TN), and oxalate-extractable iron (Fe), aluminum (Al), silicon (Si), and P. Soils pH was measured on soils at 1:1 soil to 18-megaohm deionized water mass ratio. Percent sand by mass was measured by sieving the <2 mm soil fraction through a 63 μm sieve. Bulk density was measured in cores dried at ~50°C and corrected for rocks using granite and andesite density of 2.65 and 2.60 g cm⁻³, respectively (Soil Survey Staff, 2014). Concentrations of TOC and TN were measured using a CNS dry combustion analyzer (Shimadzu Corporation, Oregon). Soils were extracted for poorly crystalline iron and aluminum oxides in a 1:50 solid-solution ratio of 0.2 M ammonium oxalate solution in darkness (Soil Survey Staff, 2014), shaken for 4 h, allowed to settle overnight, centrifuged (1,500 × g for 30 min), filtered (0.22 μm diameter PES membrane filter), and analyzed by inductively coupled plasma-atomic emission spectrometry (ICP-AES, Thermo Scientific, Waltham, Massachusetts) that was calibrated using ISO traceable standards.

Soil Total P

A subset of samples from granitic meadow (n = 6), granitic forest (n = 6), andesitic meadow (n = 5), and andesitic forest (n = 5) soils were analyzed for total P (TP) concentration by an analytical laboratory (Bureau Veritas, Inc.; Vancouver, BC; ISO/IEC 17025 and ISO 9001) using a two-step multi-acid (HNO₃-HClO₄-HF, and HCl) heated-digestion and analysis by ICP-mass spectrometry.

Total Organic Soil P by Ignition

Total soil organic P concentrations of the same subset of samples used for P-NMR analysis were measured using the ignition method (Saunders and Williams, 1955; Cade-Menun and Lavkulich, 1997). Duplicate 0.5 g subsamples of oven-dried soil were weighed. One replicate was incinerated at 550°C over a 2 h ramp-up period and maintained at this temperature for an additional 1 h followed by a 2-hr cool down. Both samples were then extracted in 1:60 solid-solution ratio of 1 N H₂SO₄, shaken for ~16 h, centrifuged at 1,500 × g for 15 min, and the supernatant was decanted and analyzed colorimetrically (Murphy and Riley, 1962). Total organic P was calculated as the difference between incinerated and non-incinerated samples. The P concentration in the incinerated sample is an estimate of soil total P (TP_inc).

Soil P Speciation by P-NMR Analysis

A subset of samples that included at least two replicates from each soil/ecosystem type were selected for P NMR analysis to identify concentrations and speciation of organic P in the soils. Following standard extraction procedures for P NMR (Cade-Menun and Preston, 1996; Cade-Menun and Liu, 2014), 2 g dry-mass equivalent undried soil subsamples were suspended in 25 ml of 0.5 M NaOH and 0.1 M Na₂-EDTA solution, shaken for 4 h, centrifuged at 1,500 × g for 20 min, and the supernatant was decanted and freeze-dried. A 1 ml aliquot of extract was taken from each sample, diluted 1:10 with deionized water, and analyzed by ICP-AES for total P, Fe, and manganese (Mn) concentrations. The P-NMR spectroscopy was conducted at the University of Idaho's Department of Chemistry. Approximately 0.24 g of freeze-dried extract powder from each sample was dissolved in 0.9 ml of NaOH-EDTA solution and 0.1 ml and D₂O and 0.5 ml of this solution was placed in a 5-mm NMR tube. The NMR spectra were obtained at 202.48 MHz on a 500 MHz Bruker Avance III spectrometer equipped with a 5-mm broadband probe. The 1D ³¹P spectra were acquired with 67.5° pulses, at 30°C, with proton decoupling, and a total recycle delay (pre-scan delay plus acquisition time) of 4 s, for 3,000–8,000 scans, determined by signal-to-noise ratios. This delay time will be sufficient for relaxation based on the ratio of P/Fe+Mn in these samples (McDowell et al., 2006; Cade-Menun and Liu, 2014). Spectra were plotted with 7 Hz line-broadening for the main spectra and 2 Hz line-broadening to assess finer details. Peak areas were computed by integration and visual inspection using NUTS software (Acorn NMR, Livermore CA, 2000 edition), with correction for the degradation of orthophosphate diesters (Cade-Menun and Liu, 2014; Schneider et al., 2016). Peak assignments were made from the literature and confirmed using phytate and β-glycerophosphate spikes (Cade-Menun, 2015).

Extractable Soil P

Concentrations of labile soil P were measured using water-extractable P (WEP), Bray-1 P (B1P), and microbial biomass P (MBP) methods. Field-moist soil samples were extracted for WEP in a 1:10 solid-solution ratio of 18 megaohm deionized water, shaken for 1 h, centrifuged (1,500 x g) for 10 min, and filtered through 0.45-μm diameter PES membrane filters (Kuo, 1996; Self-Davis et al., 2009). An aliquot was subsampled from the filtered extract for molybdate colorimetry (Murphy and Riley, 1962). Colorimetry measures phosphate that reacts with molybdate (MRP), which is used as an estimation of inorganic P in solution. However, some organic P compounds may hydrolyze during the colorimetric reaction and are included in the MRP measurement, while complex inorganic P compounds such as polyphosphates will not react with molybdate (Haygarth and Sharpley, 2000; Worsfold et al., 2016). Therefore, we hereafter refer to WEP MRP as WEP_MR. The total P in the WEP was analyzed by ICP-AES. The difference between the total WEP and WEP_MR concentrations is operationally defined as molybdate-unreactive (WEP_MU), which primarily consists of P associated with organic, non-hydrolysable, and colloidal forms (Haygarth et al., 1997; Haygarth and Sharpley, 2000). In addition to soil extraction, five subsamples from composite O-horizon samples from each forest were ground, passed through a 2-mm sieve, and extracted at 1:50 solid solution ratio for WEP and WEP_MR.

Field-moist soils were extracted for Bray-1 P (B1P) as described in Sims (2009). An aliquot of the Bray-1 extract was filtered through a 0.45-μm PES membrane filter and measured colorimetrically (B1P_MR) and by ICP-AES (B1P). The difference between B1P and B1P_MR is the B1P molybdate unreactive (B1P_MU).

Microbial biomass P (MBP) was measured by treating a 1 g dry-mass equivalent sample of undried soil using 1 ml of chloroform, placing it under a vacuum with a beaker of ~30 ml of chloroform, allowing it to evaporate for 24 h, and then extracting with the Bray-1 P extractant (Voroney et al., 2008; Reddy et al., 2013). In acidic soils, Bray-1 is a better extract for microbial biomass P than Na-bicarbonate extract (Oberson et al., 1997; Wu et al., 2000). Microbial biomass P was calculated as the difference between chloroform-fumigated and unfumigated samples, without an efficiency correction factor.

Statistical Analyses

The three seasonal samples of WEP, B1P, and MBP were pooled in the statistical analysis using a mixed model to estimate the random and fixed effects. Extract concentrations that were below the method detection limit (MDL) of the ICP (0.05 mg kg⁻¹ for WEP and B1P and 0.1 mg kg⁻¹ for MBP) were assigned values ½ MDL. The extract data were analyzed with a generalized linear mixed model using a log-normal distribution. Landscape type, parent material, and their interaction were evaluated as fixed effects, and sample point was evaluated as a random effect. Repeated measurements on the concentrations from the sample points were modeled using a compound symmetry covariance structure. For all variables except total WEP and WEP_MU, sample identification effects were estimated at each time point. Model fit was assessed by examining the log-likelihoods and inspecting residual plots. All analyses were performed in R version 3.6 (R Core Team, 2019) using the packages “nlme” (Pinheiro, 2019) for model building and ANOVA and “emmeans” (Lenth, 2019) for finding the estimated marginal means and conducting comparisons. Tukey honest significance difference (HSD) test (p < 0.05) was used to test significance for the following paired comparisons: andesitic-meadow vs. andesitic forest, granitic-meadow vs. granitic forest, andesitic forest vs. granitic forest, and andesitic meadow vs. granitic meadow.

All other soil data were tested for significance by fitting the data to the analysis of variance (ANOVA) linear models, and Tukey HSD test was used for assessing statistical differences (p < 0.05) between treatment means. Pearson's correlation coefficients were used to evaluate the strength of relationships between soil properties (Origin Lab, Northampton, MA).

Results

Soil Characterization

Upper soil profile descriptions for each watershed are listed in Supplementary Table 1. Both meadow soils are mapped as Inceptisols, and typically have aquic conditions in spring and early summer (Soil Survey Staff, 2007). Meadow soils had a darker chroma of 1 compared to chroma between 2 and 3 in the forest soils. The taxonomic descriptions of the meadow soils include subgroups Cumulic Humaquept at Paige Meadow and Cumulic Humaquept and Aquic/Oxyaquic Dystroxerept at Meeks Meadow. The difference between these subgroups is a higher seasonal water table and an epipedon thick enough to qualify as either mollic or umbric in Cumulic Humaquepts. The Aquic and Oxyaquic Dystroxerepts at Meeks Meadow have slightly deeper water tables and dark ochric epipedons (~ 15 cm) that verge on meeting the thickness requirement of a mollic or umbric epipedon. The forest soils surrounding Paige Meadow are mapped as Humic Vitrixerands. The forest soils at Meeks are mapped as Humic Dystroxerepts.

Organic horizons were ~6 cm thick at the granitic forest sites (Supplementary Table 1), and 1.5–2 cm thick at the andesitic sites. At the granitic forest sites, the decomposed litter could be separated into Oi and Oe horizons. In contrast, at sampling time (October 2018), only an Oe horizon was present in the litter at the andesitic forest site, suggesting a greater litter decomposition rate. Forest canopy coverage at the sites are similar: 49% at the andesitic site and 41% at the granitic site (Landfire, 2020). In the meadows, O horizons were 0.5–3 cm thick.

The pH of the meadow and forest soils ranged from pH 5.3 to 6.0, with meadows slightly lower than forests (Table 1). Average sand content was similar in the granitic meadow, granitic forest, and andesitic meadow (84, 87, and 88%, respectively), but it was significantly lower (75%) in the andesitic forest soils (Table 1). Bulk density of the forest and meadow soils ranged from 0.78 to 1.46 g cm⁻³. The andesitic soils had significantly lower bulk densities than the granitic soils (Table 1). The andesitic forest soils contain the most poorly crystalline iron and aluminum oxides (measured by oxalate extraction), which is consistent with Andisol classification by the USDA NRCS (Soil Survey Staff, 1999). Oxalate-extractable Fe and Al were not significantly different among the other three ecosystem-parent material types (Table 1). Oxalate-extractable Si concentrations followed the same patterns as Fe and Al concentrations. Oxalate-extractable P concentrations were significantly higher in andesitic forest soil than the other soils, and significantly lower in granitic meadow soils compared to granitic forest soils. Soil TOC concentrations were significantly different between the two forests (Table 1), but not the two meadows, and were not significantly different between meadows and forests within each watershed. Granitic forest soils contained the lowest average TOC concentration of all four ecosystem-parent material types, while andesitic forest soil had the highest average TOC concentration. The average TN concentration was approximately three times higher in the andesitic meadow soils compared to the andesitic forest soils (Table 1) and was higher in the andesitic forest soils than the granitic forest. Total N concentrations in the granitic forest and granitic meadow soils were not significantly different.

TABLE 1

Table 1. Soil physicochemical properties of replicates sampled in June 2018.

Soil Total P

The mean total soil P concentrations (TP) in the soils developed on andesitic parent materials were significantly higher than those for the soils developed on the granitic parent materials (Figure 2). Differences between forest and meadow soils within either watershed were not significant. The estimated total P concentrations via incineration and H₂SO₄ extraction for the soil samples analyzed by P-NMR were similar to the total P measured from three acid digestion (slope = 0.95, r² = 0.95; Supplementary Table 2). The total P stocks for the 0–15 cm mineral soils calculated using the mean bulk densities were 0.69, 1.10, 1.22, and 1.81 Mg Ha⁻¹ for the granitic meadow, granitic forest, andesitic meadow, and andesitic forest, respectively.

FIGURE 2

Figure 2. Mean total soil phosphorus. Error bars are standard errors. Values with the same letter are not significantly different (α =0.05), n = 6 for granitic parent material sites and n = 5 for andesitic parent material sites.

Speciation of Soil P

In the meadow soils, total organic P concentrations determined from incineration and H₂SO₄ extraction were 79–92% of the TP_inc (Figure 3, Supplementary Table 2). Organic P concentrations in the forest soils were much lower (13–47% of the TP_inc) than in the meadow soils.

FIGURE 3

Figure 3. Organic and inorganic P in forest and meadow soils from andesitic and granitic parent materials determined by difference of H₂SO₄ extraction of non-incinerated and incinerated soils. Numbers within bars are percent organic P.

The NaOH-EDTA extraction efficiency ranged from 33 to 75% of total soil P (Supplementary Table 2). The P not extracted by NaOH-EDTA is considered to be predominantly mineral-bound inorganic P and not readily available to the soil solution or for biological cycling (Cade-Menun et al., 2015).

Example P-NMR spectra are shown in Figure 4, the concentrations (and percentage of extracted P) are shown in Supplementary Table 3, the grouping of these P species into pools (total organic and inorganic P) and compound classes (total polyphosphates, etc.) are shown in Supplementary Table 4, and the chemical shifts of the identified P compounds are shown in Supplementary Table 5. The concentrations of the main P compound classes within each ecosystem and parent material type are shown in Figure 5.

FIGURE 4

Figure 4. NMR spectra from NaOH-EDTA extracts of forest and meadow soils from the granitic site.

FIGURE 5

Figure 5. Mean P species in soils measured using 31-NMR analysis of NaOH-EDTA extracts of soils. The monoester and diester values were corrected for degradation (denoted with “c” prefix). Inorganic P is the sum of the inorganic orthophosphate and polyphosphate compounds.

Inorganic P compounds identified in the NaOH-EDTA extracts by NMR include orthophosphate, pyrophosphate, and polyphosphates. Pyrophosphate and polyphosphates were grouped together as total polyphosphates (Supplementary Table 4), and all three were summed together as inorganic P (Supplementary Table 4, Figure 5). For all soils, orthophosphate was the dominant inorganic P form, and for forest soils developed on both parent materials it comprised the majority of P in the NaOH-EDTA extracts (71.2–84.7% for granitic forests; 51.7–69.3% for andesitic forests). In contrast, the percentages and concentrations of all inorganic P compounds were much lower in meadow soils developed on both parent materials, averaging 21% of extracted P (Figure 5). There were no clear trends among the ecosystem and parent material soil types for pyrophosphate or polyphosphates, which were present in all samples, ranging from at 2.1–4.9% of NaOH-EDTA extracted P.

The percentage organic P determined by P-NMR on the soil extracts was directly correlated with the percentage determined using incineration and H₂SO₄ extraction (r² = 0.95, Supplementary Figure 1). For all soil types, all the major organic P compound classes were identified: phosphonates, orthophosphate monoesters (hereafter called monoesters), and orthophosphate diesters (hereafter called diesters). The phosphonates included several different peaks (Supplementary Table 5), indicating that a number of different compounds were present, but these were not specifically identified. Concentrations of phosphonates ranged from 1.4 to 10.3 mg kg⁻¹ (0.7–4.9% of extracted P) and were generally higher in meadows than forests.

Monoesters identified in the P-NMR spectra included four stereoisomers of inositol hexakisphosphate (IHP): myo-IHP (phytate), scyllo-IHP, neo-IHP, and D-chiro-IHP. Of these, myo-IHP was the predominant P form and was generally more abundant in meadows than forests. For most of the soil samples, myo-IHP exceeded the sum of the other three stereoisomers. Other specifically identified monoesters were glucose 6-phosphate (0.6–2.1% of extracted P), choline phosphate (0.3–1.3%), α-glycerophosphate (0.3–2.8%), β-glycerophosphate (0.7–5.7%), nucleotides (1.3–12.8%), and an unidentified peak at ~5 ppm, which was present in all samples at 0.7–9.2% of the extracted P. Although peaks for α-glycerophosphate, β-glycerophosphate, and nucleotides are present in the monoester region of spectra, they originate during NaOH-EDTA extraction and P-NMR analysis as a result of degradation of diesters in the original soil samples (Cade-Menun, 2015; Schneider et al., 2016). Thus, the peak areas from these compounds were subtracted from the monoester peak areas and included with the diesters.

Peaks representing diester compounds were separated into DNA (0.5–5.9%), Diester 1 (2.33 to −0.27 ppm, 0.6–8.7%), and Diester 2 (−0.9 to −3.72 ppm, 0.2–3.4%). The Diester 1 region included phospholipids and lipoteichoic acids, while the compounds in the Diester 2 region have not been specifically identified. The proportions and concentrations of P in these three diester regions were generally greater in meadows than forests for both parent materials. Total diesters (cDiesters), calculated by including the degradation compounds from the monoesters, confirmed that the percentages of cDiesters were greater in meadows than forests, and concentrations were greater in andesitic forests and meadows than the granitic forests and meadows, respectively (Supplementary Table 4, Figure 5).

The average of the replicate NMR results (Figure 5) shows that in the meadow soils the three major organic P compound classes were greater in concentration than in the forest soils. In meadow soils, cMonester was 50.3% (107 mg kg⁻¹), cDiester was 25.8% (52.1 mg kg⁻¹), and phosphonates were 3.0% (5.6 mg kg⁻¹). In forest soils, cMonester was 15.7% (42.5 mg kg⁻¹), cDiester was 8.5% (22.2 mg kg⁻¹), and phosphonates were 1.3% (3.4 mg kg⁻¹). The ratio of cMonesters to cDiesters was over 1 for all soils, indicating that monoesters were the dominant P compound class in both ecosystems. In meadow soils, total IHP concentrations comprised about one third of the cMonoesters 17.3% (33.5 mg kg⁻¹) but were half of the cMonoesters in forest soils 7.7% (20.9 mg kg⁻¹).

Effects of Parent Material and Ecosystem on Extractable P

In evaluating the extractable P concentrations, main effects were calculated as well as interactions (Supplementary Table 6). For all the WEP and B1P extractions, there were crossover interactions between parent material and ecosystem type (Supplementary Figure 3). Thus, the effect of ecosystem should be interpreted in context of parent material, and vice versa.

The mean total WEP (WEP_Total) concentrations from all the paired interactions (parent material within ecosystem type and ecosystem type within parent material type) were significantly different from each other (Supplementary Table 6). The granitic forest soils had the greatest mean WEP_Total concentration of the four ecosystem-parent material types (Figure 6). The andesitic forest soils had the lowest mean WEP_Total concentration. WEP_MU comprised the largest fraction of WEP (Supplementary Table 6) in all ecosystem-parent material types. This suggests that most of the WEP exists as soluble organic P compounds or as inorganic P complexed to colloids instead of as dissolved phosphate.

FIGURE 6

Figure 6. Estimated marginal means for water extractable P (WEP) (bottom), Bray 1 extractable P (B1P) (middle), and microbial biomass P (MBP) (top) from soil samples from the different ecosystem and parent material sites. WEP and B1P extract subscripts indicate total solution P, molybdate reactive P (MR), and molybdate unreactive P (MU). Error bars are the standard errors.

The mean total B1P concentration from the granitic forest soils was more than 10 times greater than the mean B1P concentration from the other soils (Figure 6). The B1P concentrations at the two forest soils were significantly different, as were the meadow to forest comparisons (Supplementary Table 6). However, the B1P from the andesitic meadow and granitic meadow soils were not significantly different. Most of the B1P was molybdate reactive P (Supplementary Table 6), suggesting B1P was predominantly inorganic P extracted from the soil.

The mean MBP concentrations for the meadow soils were approximately seven times more than the mean MBP concentrations from the forest soils (17.4 mg kg⁻¹ compared to 2.55 mg kg⁻¹, Supplementary Table 6). All paired comparisons for the interactions were significantly different. The andesitic meadow soils had the greatest MBP concentrations, followed by the granitic meadow soils (Figure 6).

Water Soluble P of Soil O Horizons

The WEP and WEP_MR concentrations from composite O-horizon samples are shown in Table 2. The less decomposed Oi composite sample from the granitic site had greater WEP concentrations than the Oe sample.

TABLE 2

Table 2. Mean^a WEP and WEP_MR extract concentrations from composite samples of the O horizons from the forested sites.

Discussion

Parent Material Effects on Soil P

Total P concentrations in unfertilized soils are typically linked to P content of the parent material, which is the source of P and is also an important control of soil mineralogy that affects P retention (Gardner, 1990; Porder and Ramachandran, 2013; Hahm et al., 2014). The granitic Meeks Meadow soils developed on granodiorite eroded from glacial till and outwash, while the andesitic soils at Paige Meadow developed on glacial deposits of eroded basaltic and andesitic volcanic rock. In a database of P concentrations of common rock types, mean total P concentrations of granodiorite and granite are 810 and 568 mg kg⁻¹, respectively, while mean total P concentrations of basalt, andesite, and basaltic-andesite are 1,304, 1,150, and 1,551 mg kg⁻¹ (Porder and Ramachandran, 2013). Total P concentrations in Tahoe soils developed in andesitic parent materials were over two times greater than in soils developed in the granitic parent materials (Figure 2), which is consistent with expected differences based on parent material.

In addition to P inputs from parent materials, other P inputs and losses from the soil also affect soil TP concentrations. Phosphorus released from parent material by weathering is taken up by plants and microbes and converted to other P forms and can also be leached out of the soil profile. It has been noted that total P concentrations in forest soils decrease over time (Yang and Post, 2011; Deiss et al., 2018; Nelson et al., 2020). There are no clear trends for WEP, B1P, or MBP concentrations between the granitic or andesitic watersheds (Figure 6, Supplementary Table 6). However, there were significant differences for WEP, B1P, and MBP when the interaction of parent material type and ecosystem were considered (Supplementary Figure 3). The granitic forest soils had much greater WEP and B1P concentrations than the andesitic forest soils, despite having lower TP concentrations. This suggests greater P availability, or more labile P, in the granitic forest soils than in the andesitic forest soils. The labile P is available to either efflux out of the soil or be taken up by plants, thus causing a decrease in TP concentration in the soil. Johnson et al. (1997) reported significantly higher concentrations for a B1P-type extraction from andesitic forest soils than granitic forest soils at other locations in the Lake Tahoe Basin, but did not report concentrations for other P pools such as TP. Brödlin et al. (2019a) observed that in sandy soils derived from glacial-till parent materials in a European deciduous forest, greater P leaching occurred than in soils derived from volcanic rock. Thus, for the granitic forest soils in Lake Tahoe, the decreased TP concentration is both due to the decreased parent material rock-P inputs and export of available P out of the soil. Since the meadow soils are depositional environments, the TP concentrations in these soils reflect the parent material trends as well: andesitic meadows have more TP than granitic meadows. But differences in labile P are not as great as the differences between the forest soils. Thus in the meadows, other ecosystem processes are more important in controlling the availability of the soil P.

In the soils of the current study, soil properties were not clearly delineated by parent material alone. Soil pH was more closely related to ecosystem type than parent materials (Table 1), and concentrations of oxalate-extractable Al, Fe, Si and P, and sand content of the andesitic forest were significantly different from the andesitic meadow and the granitic forest and meadow (Table 1). The andesitic forest soil properties are due to the andesitic-parent material contributions, while the andesitic meadow is in a depositional environment collecting eroded soils transported from surrounding landscapes. In addition, the hydraulic conditions and vegetative community dynamics of the meadow create pedogenic processes that can overshadow the andic soil properties. Thus, the meadow soil has less andesitic parent material influences, even though it is within a predominantly andesitic parent material landscape. The A horizon in the andesitic forest soils is only 6 cm thick, and thus some of the soil properties in the 15 cm cores sampled from this site were impacted by Bw horizon material (Supplementary Table 1). At nearly all of the other sites, except one andesitic meadow profile, A horizons are more than 15 cm thick, and thus the cores represent the pedogenic horizons. The inclusion of Bw material in the andesitic forest soil sample may have contributed to its distinct properties as compared to the other soil cores; however, the mineral properties of this sample as measured by the oxalate extractable Fe, Al, and Si (Table 1) are indicative of andic soils, which have high concentrations of poorly crystalline aluminum and iron oxides in both the surface and subsurface horizons. Thus, the andesitic soil samples, while composed of A and B horizon soils, have mineral properties that are indicative of the andesitic parent materials from which the soils developed and the soil P properties from this sample are indicative of the top 15 cm of the soil.

Soil P sorption capacity is related to clay minerals and iron and aluminum oxides (Gerard, 2016). The andesitic forest soils had the lowest sand content and the highest concentrations of amorphous iron hydroxides and aluminosilicates (Table 1), suggesting that P adsorption on mineral surfaces will be strong in these soils (Khare et al., 2005; Redel et al., 2008). Phosphorus complexed as Al or Fe-organic matter ternary complexes are also common in andesitic soils (Gerke and Hermann, 1992; Gerke, 2010). Oxalate-extractable P was positively correlated to oxalate extractable Fe, Al, and Si (r = 0.82, 0.85, 0.82, respectively; Supplementary Figure 2) and the highest concentrations of oxalate-extractable P were in the andesitic forest soils (Table 1), confirming the relation of the soil P to the andesitic soil minerals.

Sorption capacity will influence both inorganic and organic P forms. Total organic P in the Lake Tahoe soils does not appear to be related to the parent material, beyond the differences in TP already discussed, and is clearly more influenced by ecosystem (Figure 3; Supplementary Table 2). For specific organic P forms, correlations of concentrations of myo-IHP and other IHP stereoisomers with concentrations of oxalate-extractable Fe and Al have been widely reported for soils and are thought to demonstrate the sorption of these compounds to amorphous iron hydroxides and aluminosilicates (Jorgensen et al., 2015; Nelson et al., 2020). The average of the percentage of total extracted IHP (Supplementary Table 4) for the eight granitic soils was similar to that for the four andesitic soils (13.1 and 11.3% respectively), while the average concentrations of total IHP in andesitic soils were double those of the granitic soils (41.2 vs. 20.2 mg kg⁻¹). The same trend was also observed for DNA and the general Diester 1 category (Supplementary Table 3). In acidic soils, adsorption of DNA occurs, and its NaOH-EDTA-extracted concentrations have been correlated with those of oxalate-extracted Fe and Al (Condron et al., 2005; Nelson et al., 2020). However, the Diester 1 category includes phospholipids and lipoteichoic acids, which do not sorb to soil minerals (Condron et al., 2005). The differences in concentrations of these compounds and compound classes are consistent with differences in total concentrations of P in the NaOH-EDTA extracts (Supplementary Table 2), which are consistent with total concentrations of soil P (Figure 2), so may simply reflect trends in total P rather than selective binding of P compounds.

Ecosystem Type and Soil P

Although TP concentration did not differ between ecosystem types at the two different parent material watersheds, total organic P measured by the incineration method was substantially greater in meadow soils than forest soils (85.5 vs. 36.2%, averaging data from Figure 3). The same trend was observed for total organic P determined by P-NMR, even though the recovery of total P was lower in NaOH-EDTA extracts from meadow soils (40.8%) than forest (61.9%, averaging data in Supplementary Table 2), which could underestimate inorganic P. Chiu et al. (2005) measured P-NMR spectra in NaOH-EDTA extracts from subalpine grassland and forest soils and observed a similar fractionation of inorganic and organic P forms between the two ecosystems. In addition to differences in total organic P, P-NMR revealed differences in P forms and compound classes between the two ecosystem types (Figure 3, Supplementary Table 4).

The forest and meadow ecosystems differ in vegetation, elevation, and slope position, all of which influence P cycling. The plant species in the meadows have above and below-ground vegetation that can readily decompose compared to forests (Margalef et al., 2017). In forests, litter is deposited onto the forest floor and gets incorporated into an Oi horizon that decomposes to an Oe horizon. An Oi horizon with identifiable pine needles was observed in the granitic soil, but not the andesitic forest. This may be due to andesitic soil properties that increase soil moisture retention, which facilitates greater microbial decomposition rates and thus quicker breakdown of forest litter (Sun et al., 2017). Different species of plants and even the same species of plants growing under different soil fertility conditions will contribute different P species (Noack et al., 2014) that can change with depth in the soil profile (Nelson et al., 2020).

The vegetation from these ecosystems was not analyzed by P-NMR, so we cannot say with certainty the P forms input from plants. However, myo-IHP is widely recognized as a plant P compound (Condron et al., 2005). Other compounds may originate from plants or microbes or can be produced by alteration of plant-P compounds (Condron et al., 2005). In the soils of this study, microbial P (MBP) concentrations were greater in meadow soils than they were in forest soils. However, in coniferous forests, the majority of microbial activity and P cycling occurs in the O horizon, associated with the hyphal mat of ectomycorrhizal fungi at the soil-organic matter interface (Plassard et al., 2011; Nelson et al., 2020). This may account for difference in P forms and MBP concentrations between these ecosystems and would be consistent with the substantially higher WEP concentrations in O horizons (Table 2) than mineral soils in the forest.

Availability of soil P controls P immobilization into microbial biomass (Olander and Vitousek, 2004; Yang and Post, 2011; Spohn and Widdig, 2017; Pistocchi et al., 2018). For example, Pistocchi et al. (2018) observed that during incubation of a deciduous forest soil with low available P, P cycling between soil and microbial biomass was conservative, while in soil with higher available P there was more exchange between microbial-bound P and inorganic soil P pools (i.e., mineral-bound P). Thus, in forest soils, when P availability exceeds biological demands, geochemical processes (adsorption and precipitation) predominate over immobilization by microbes (Olander and Vitousek, 2004); this implies that geochemical processes control P availability for leaching or root uptake in the mineral soil horizons. In the andesitic soils, the decreased labile P concentrations (Figure 6) inhibit microbial P fixation, causing low MBP. Aluminum toxicity is another cause of decreased MBP in the andic forest soils because high soluble Al concentrations inhibit microbial enzyme production, including phosphatase, thereby limiting P immobilization (Kunito et al., 2016).

Another possible cause of the MBP increase in the meadow soils compared to the forest soils is the increase in N availability in the meadow soils (Table 2). Microbes mineralize organic P for microbial uptake using phosphatase enzymes, which require N for production (Vitousek et al., 2010; Marklein and Houlton, 2012). Total N and MBP were significantly correlated (r = 0.81) (Supplementary Figure 3). Thus, the occurrence of sufficient N availability for phosphatase generation in the meadows facilitates degradation of organic P compounds, and subsequently, the biologically available P can be immobilized by microbes. Soil moisture may also play a factor, with greater moisture in the meadows increasing microbial activity.

Mycorrhizal association may also greatly influence P cycling. The ectomycorrhizal fungi found in forests will produce more phosphatases than endomycorrhizae associated with meadow plants (Plassard et al., 2011; Margalef et al., 2017). They will also produce organic acids such as oxalate (Plassard et al., 2011). These will desorb both inorganic and organic P, and both organic acids and phosphatases may need to be present simultaneously to mineralize organic P (Giles et al., 2018). This could also account for the reduced organic P concentrations in these forests compared to meadows.

In meadows, high seasonal water tables can have a significant impact on soil properties that influences P cycling. This, combined with high organic C concentrations, can mask the influence of mineralogy on labile P (Sah et al., 1989; Johnston et al., 1995). Sah et al. (1989) observed that C availability controlled P availability in wetland soils: when total organic C concentration exceeded 0.8%, it promoted the reduction of ferric (Fe³⁺) oxides, which decreased sorption capacity of the soil for P. Alternatively, during periods of flooding, the precipitation of ferrous iron-phosphate minerals such as vivianite may occur (Zhang et al., 2003; Heiberg et al., 2010; Rothe et al., 2016), which is less soluble upon drainage because the P remains occluded by oxidized iron-hydroxide minerals that form when the vivianite oxidizes (Sah and Mikkelsen, 1986b). The seasonal redox cycles that occur in wetland soils may decrease P leaching compared to unflooded soils, even after soils remained drained over 4 months (Sah and Mikkelsen, 1986a). In contrast, Gergans et al. (2011) proposed that excess sulfate in Tahoe Basin wetland soils facilitates production of iron sulfides during reducing conditions that make Fe unavailable when the soils re-oxidize to adsorb phosphate and organic P compounds, thus making soil P more available for continued leaching in these wetland soils. Based on the varying results of the studies discussed above, there are several factors that influence P speciation in Lake Tahoe Basin meadows, including seasonal flooding, high organic matter, and Fe biogeochemical transformations. These factors may be more important than parent materials for controlling P mobility. The presence of higher concentrations of phosphonates in meadow soils compared to forest soils is also consistent with higher moisture levels (Condron et al., 2005), and higher diester concentrations have also been reported for poorly drained soils compared to well-drained soils (Young et al., 2013).

Labile Soil P and Potential Loss to Lake Tahoe

Soil P buffering capacity is the degree to which soil can adsorb or release P from exchange sites to maintain dissolved P concentrations in the soil solution (Holford, 1997). Soils with larger total P reserves are considered to have greater buffering capacity to replenish P taken up by plants or leached out of the soil (Daly et al., 2015). An estimate of soil P buffering capacity is the P saturation index (PSI), which is calculated as concentrations of oxalated-extracted P divided by the sum of the oxalate-extracted Fe and Al (Schoumans, 2009). For the Lake Tahoe watershed soils, PSI was 0.029 and 0.072 for the granitic meadow and forest soils, and 0.047 and 0.029 for the andesitic meadow and forest soil (Table 1).

Of the extractants used in this study, WEP measures the most labile P, B1P extracts less labile P that is sorbed to the soil, and oxalate or NaOH-EDTA extract both labile P and P that is more tightly held by soil through either adsorption complexes, mineral-bound P, or larger organic P compounds. In all of the Lake Tahoe Basin soils, B1P concentrations were approximately an order of magnitude greater than WEP concentrations (Figure 6), suggesting a large amount of adsorbed P is released by the B1P extractant that is not released by water extraction. Although a stronger P buffering capacity is expected in the andesitic forest soils, where the highest TP was observed, the lower PSI in these soils suggests there is excess P sorption capacity on high-adsorption affinity iron and aluminum oxides, thus causing the lower concentrations of labile P (WEP and B1P) as compared to the granitic forest soils. The B1P extractant apparently did not access the P in the andesitic soils that was either strongly adsorbed, had formed Al or Fe-P mineral phases with low solubility (Negrín et al., 1996), or was complexed as Al or Fe-organic matter ternary complexes, which are common in andesitic soils (Gerke and Hermann, 1992; Gerke, 2010). The NaOH-EDTA extracts removed a much greater amount of the total P than the B1P extracts (Supplementary Table 6) but were similar in concentration to those of oxalate-extractable P (Table 1). Thus, there appears to be a large reserve of P associated with iron oxide and allophane minerals in the andesitic forest soils that is not readily available for release to the soil solution.

Analysis of the NaOH-EDTA soil extracts from the andesitic forests by P-NMR showed that 26–46% of the extracted P was organic P species (Supplementary Table 4). This organic P fraction in the andesitic soils may be slowly available to plants and microbes that release organic acids, which enhance P release through competitive exchange (Harrold and Tabatabai, 2006), and which may be synergistic with hydrolysis of organic P compounds by phosphatases (Giles et al., 2018).

Most of the WEP in forests and meadow soils from Lake Tahoe was not reactive with molybdate blue chemistry (WEP_MU) (Figure 6, Supplementary Table 6). The source of water-extractable organic P compounds are inputs from plants and soil microbes. The predominance of WEP_MU from both forest and meadow soils in the Lake Tahoe Basin is a potential source of mobile P that most likely consists of organic P compounds (Worsfold et al., 2016). The labile organic P compounds can be transported to Lake Tahoe by vertical and lateral transport processes, especially during high intensity events that cause preferential flow through macropores, which are common in coniferous forest soils (Luo et al., 2019). Organic P species have been reported in soil leachate, snowmelt runoff, and samples of river inlet and floodplain waters during flooding events (Toor et al., 2003; Cade-Menun et al., 2006; Wiens et al., 2019), and in water-extractable colloids from grasslands and forests (Missong et al., 2016; Jiang et al., 2017).

Water-soluble P from leaf litter is an important source of labile P return to the soil (Uselman et al., 2012; Sohrt et al., 2019). In a separate study of soluble P from O horizons from Lake Tahoe Basin (unpublished data), nine samples were collected from forest and meadow sites near the Paige Meadow and Meeks Bay watersheds. WEP concentrations in these samples ranged from 54 to 209 mg kg⁻¹ (mean = 122 mg kg⁻¹, standard deviation = 45 mg kg⁻¹), indicating that WEP concentrations in the O-horizon samples are highly variable throughout the two watersheds. Based on the two Oe composite samples (Table 2), molybdate-reactive P (inorganic P) was the predominant phase of WEP (67% in both the granitic and andesitic forests). The concentration of P from the composite Oi horizon sample from the granitic site was more than twice that of the Oe horizon sample, suggesting P is lost from the litter as it decomposes (Table 2).

The Oe WEP concentration was ~58 times greater than soil A horizon WEP concentration in the granitic forest soils, and 690 times greater than in andesitic forest soils. Miller et al. (2005) measured soluble P from O horizons of Lake Tahoe Basin forest soils using laboratory simulated precipitation and snowmelt leaching experiments and observed 46 mg kg⁻¹ of soluble P leached from the less decomposed Oi horizon and 28 mg kg⁻¹ from the more decomposed Oe horizon. The O horizons in the Miller et al. (2005) study came from Jeffrey and Sugar Pine forests in a granitic watershed. Although the Miller et al. (2005) water extraction methods were different than those used in this study, both the Miller et al. (2005) leaching experiment and this study's granitic forest O-horizon samples (Table 2) have more water-soluble P in the minimally decomposed Oi horizons compared to the more decomposed Oe horizons.

Approximately one-third of the WEP from the O horizon samples was WEP_MU, which could be a considerable source of organic P mobilized into the soil and possibly to the surface water, depending on the organic P species and reactivity. Both the granitic and andesitic sites have a similar vegetative density (41 vs. 49%) (Landfire, 2020), thus, vegetative P inputs to the soils should be similar in the two watersheds. However, based on the greater WEP and B1P concentrations in the granitic forest soils, P outputs to streams and groundwater are expected to be much greater from these systems; this is due to the lower sorption capacity of the soils that allows for a significant amount of P release in the extractions. Uhlig and von Blanckenburg (2019) estimated that P inventory of the forest litter in montane, temperate forest ecosystems can only sustain vegetative demand for a few decades, and that continuous release of P from parent rocks must occur to sustain forest growth. The different adsorption capacities of the andesitic and granitic soils in the Lake Tahoe forests can have a major influence on the timescales of P availability and its cycling between the forests, litter, and parent material.

There are several sinks for WEP from forest litter: it can be taken up by plants and microorganisms for internal cycling; leached into the soil where it may adsorb, be immobilized, or further leached into ground water; or be transported off site via surface runoff of dissolved P or eroded P-containing particles. Although concentrations of WEP in the O horizons are much greater than the soils, it is a much smaller total P pool in the ecosystem than soil P, which is large and stores much of the WEP leached into it from O horizons (Yang and Post, 2011). The high concentrations of soluble P in the Lake Tahoe Basin O horizon samples indicate that a large flux of available P can enter the soil. Much of this flux occurs during spring snowmelt. The fate of this O horizon-sourced P is a function of the characteristics of soil biological and physical properties and site hydrology.

Ohara et al. (2011) recorded that more than 90% of field-observed hillslope drainage in a Lake Tahoe watershed occurred as subsurface lateral flow through the soil. Thus, soil reactions are important processes controlling P transport to surface waters, which would be especially high during periods of continuous snowmelt. These processes are impacted by the species of soluble P in the soils, which both the soil extractions and soil P-NMR analyses suggest are both inorganic and organic P species.

When streams near our research sites experience peak discharge, molybdate-unreactive fractions make up 61–67% of filterable (<0.45 μm) P [Supplementary Figure 4 (USGS, 2016)]. Therefore, molybdate-reactive and unreactive fractions in these nearby streams during snowmelt more closely reflect the WEP fractionation of soils (50–74% WEP_MU; Supplementary Table 6) than WEP_MU from O horizons (28–36%; Table 2). A possible explanation of this is that during high-flow periods, the inorganic (molybdate-reactive) P species are attenuated by forest and meadow soils leading to net exports of organic forms. Thus, the flowing solution reflects the WEP_MU fraction leaching from the soil. Bol et al. (2016) conducted an extensive review of P fluxes in forested ecosystems and concluded that P loss as colloidal-organic P that is exported from soil profiles through macropores during high-intensity rainfall events is likely a critical factor in P export. Colloidal organic P would be included in the WEP_MU fraction in this study. Considering that organic P increases in the stream in the Lake Tahoe Basin during high-intensity events (Supplementary Figure 4), preferential-flow path loading is a likely scenario occurring in forest and meadow watersheds in the Lake Tahoe Basin. However, to explain the molybdate unreactive ratio of the stream water, there must be attenuation of the inorganic P as it moves through the preferential flow paths; otherwise the ratio of inorganic to organic P forms in the stream water would be more similar to the ratios in the O-horizon extracts. Alternatively, forest soils may be transporting P-laden water through preferential flow paths where P attenuation is minimal, but as the flow continues toward the streams, it is intercepted by riparian meadows that have fewer preferential flow paths, enhanced groundwater storage, and greater microbial activity that immobilizes orthophosphate, causing the soil water that exfiltrates into the streams to have a greater proportion of WEP_MU than what is leached from the forest O horizons. A third mechanism of inorganic P attenuation that may enrich Lake Tahoe Basin stream waters with molybdate unreactive P is preferential adsorption of inorganic P within the stream on suspended particles eroded from soils. Since the highest total suspended solids occurs during high runoff events, adsorption may be significant enough to alter the dissolved inorganic and organic solution composition during these periods.

Because molybdate unreactive P (organic P) accounted for the majority of WEP from the Lake Tahoe Basin soils (Figure 6), it is likely the most vulnerable for transport as lateral flow during spring snowmelt or exfiltration from meadows, thereby increasing the P load in surface waters. The high concentrations of P released from litter suggests that forest management practices that remove timber and deposit deep layers of chopped fresh organic matter (mastication) to prevent erosion may be creating a potential source of P that can be leached into surface waters—at least in the short-term time it takes for the material to degrade. In watersheds with soils developed on granitic parent materials, this would be especially problematic. A beneficial focus of future research would be an examination of the speciation of P in the forest O horizons and comparison to P forms in both forest and meadow soils, as well as measurement of soil and stream water samples for P-species composition. Additionally, evaluation of the subsoil deeper than 15 cm should be done to account for how P reactions influence leaching through the deeper soil profile.

Conclusion

In soils in the Lake Tahoe Basin, P storage shifts from sorption on minerals in forests, to immobilization in microbial biomass in meadows. In forested hillslopes, adsorbed P may be gradually depleted if it is leached from the soil into ground and surface water. The degree of P depletion depends on the parent materials from which the soils developed (granitic vs. andesitic). In soils developed on andesitic parent materials, forest cycling of P is mediated by the high adsorption capacity of P on andic minerals, while in soils derived from granitic parent materials the increased resistance to weathering creates coarser-textured soils and fewer soil clays, causing a decreased P adsorption capacity. As a result, granitic soils have greater potential P mobilization into groundwater and lateral runoff into surface waters.

Organic P was a predominant water-extractable fraction from all soils. Total organic P concentration was greater in meadow soils than forest soils, and in all soils of this study, orthophosphate monoesters were the main organic P compound class, even after correcting for diester degradation during analysis. The organic P compounds in the soils can be leached into the surface waters. Once in the surface water, mineralization of the organic P compounds can make phosphate available to aquatic organisms, causing surface water quality degradation.

Results from this study provide insights into speciation of P in forest and meadow soils and show the importance of parent materials on P availability. This information can be used to better understand which ecosystems present the most risks for P loading into Lake Tahoe, which will allow for better forest management practices to prevent P export into feeder streams and groundwater that discharge into the lake. Current management strategies use controlled burns and erosion prevention strategies to prevent P loss from Lake Tahoe Basin soils into the lake. Resource managers need to consider the highly variable sources of P in the Basin to decide which watersheds are most vulnerable to P loss, such as the granitic forest watersheds, and match management strategy to site properties to optimize the site management for decreased P loss.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

TH conducted the research and reported it as part of a MS thesis at the University of Idaho. DS is the corresponding author and led all aspects of the research and writing of this manuscript. MD, CD, EB, and CG designed the research and assisted with experiments and data interpretation. BC-M conducted NMR data analyses and interpretation as well as helped interpret all of the experimental data. JP did the mixed-model statistical analyses of the data and consulted in data interpretation. AC assisted with conducting experiments and collecting and analyzing data. All authors contributed to writing and editing of the article and approved the submitted version.

Funding

This work was supported by AFRI program (Grant No. 2016-67020-25320/project accession no. 1009827) from the USDA National Institute of Food and Agriculture.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We appreciate assistance with NMR experiments by Alex Blumenfeld in the University of Idaho Chemistry Department, and assistance with O horizon water extractable P measurements by Tiffany Perez.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/ffgc.2020.604200/full#supplementary-material

References

Anderson, B. H., and Magdoff, F. R. (2005). Relative movement and soil fixation of soluble organic and inorganic phosphorus. J. Env. Qual. 34, 2228–2233. doi: 10.2134/jeq2005.0025