Our results are surprising and raise questions about fundamental relationships between soil pH and biological processes. While we expected strong positive relationships between soil pH and microbial biomass and activity, we instead observed significant curvilinear relationships with peaks at unexpectedly low values, between pH ~ 3.5 and ~ 4.5. These relationships are based on 20 years of data that span reference and Ca-treated watersheds and organic and mineral soil horizons.
Microbial biomass C is a key functional pool for C and N cycling in organic detritus since it represents the active biological component of soil organic matter. Microbial respiration is a more dynamic, short-term indicator of soil biological activity. The fact that there were many significant relationships between these variables and pH in the Oie horizon, which is the most biologically active zone of the soil profile, suggests that soil pH is a key controller of biological activity at this northern forest site, but the nature of the controls is unexpected. The curvilinear relationships observed may suggest that the soil communities are adapted to the inherently low pH at the site. Addition of calcium to Watershed 1 experimentally increased the pH but clearly did not shift this adaptation to a higher pH. The curvilinear relationships contrast with our expectations and other results from Watershed 1 (Sridevi et al. 2012) showing that decreased soil calcium was correlated with reduced microbial biomass and C sequestration of forest trees.
Our results also show potential adaptation of N cycling processes to low pH in Hubbard Brook soils. We observed significant curvilinear relationships for potential net N mineralization in the Oa/A horizon and for potential net nitrification in the mineral and Oa/A horizons. The most consistent variable was the soil nitrate pool, which was the sole variable that was significantly correlated with pH in all three horizons.
An alternative explanation for the curvilinear relationships observed centers on plant: microbial interactions. Increases in pH could stimulate microbial biomass and activity to a certain level at which plants begin to respond and limit the ability of microbes to respond further. Similar dynamics have been observed in other studies (discussed below) at HBEF. Thus, the complexities of plant: microbial interactions may give the appearance that microbes are adapted to/optimized for the low values while in reality, their ability to respond to increases in pH is being inhibited by plants.
Complex effects of soil pH at a naturally acidic watershed
Forest soils at HBEF are naturally acidic, with low base cation supply due to low rates of mineral weathering and low exchangeable base cation concentrations in the O horizon (Lawrence et al. 1995; Likens et al. 1996). The low levels of exchangeable Ca in mineral soils are due to low rates of mineral weathering and are associated with elevated concentrations of Al in the O horizon and streams at the site (Likens et al. 1998; Palmer et al. 2004).
The addition of calcium silicate to watershed 1 in 1999 rapidly (by 2000) increased exchangeable Ca, soil pH, cation exchange capacity, and base saturation and decreased exchangeable acidity in organic soil horizons (Cho et al. 2010; Cho et al. 2012; Johnson et al. 2014). In the mineral soil layer, short-term increases of exchangeable Ca, soil pH, cation exchange capacity, and base saturation were much smaller and exchangeable acidity exhibited small decreases to no change.
Over the last 20 years, pH has declined in the Oa/A and Oie horizons of the treated watershed, consistent with processing and loss of the added calcium (Johnson et al. 2014). Soil pH has increased in the mineral and Oa/A horizons of the reference watershed, consistent with deacidification and reduction of acid deposition (Likens et al. 2021). Our sampling sites, which range across elevation gradients in both reference and treated watersheds, thus encompass a wide range of pH conditions for the northern hardwood forests represented by the HBEF.
Several studies at HBEF have found that soil pH controls on soil biological processes are more complex that might be expected from consideration of the chemical dynamics of base cations, acidity and aluminum. Fiorentino et al. (Fiorentino et al. 2003) investigated the effects of the calcium addition to Watershed 1 on P pools and cycling and N cycle processes. While there were no significant differences in soil solution or fine root P concentration, microbial biomass P, microbial C to P ratios, available organic and inorganic P fractions were lower in the Oe horizon of the calcium treated watershed than the reference watershed. They concluded that the soil pH changes associated with the Ca addition had increased the rates of P cycling in the forest floor horizons, with enhanced flow of P to plants at the expense of microbes (Fiorentino et al. 2003).
Several studies of the watershed 1 experiment have reported complex effects of calcium additions on N cycling processes (Battles et al. 2014; Juice et al. 2006; Rosi-Marshall et al. 2016). Groffman et al. (Groffman et al. 2006) observed declines in the N content of the microbial biomass, potential net and gross N mineralization rates, and soil inorganic N pools in the Oie horizon of the treated watershed. These results clearly suggested that the calcium addition did not stimulate microbial N cycling, and the authors suggested that plants had responded more dynamically to the calcium addition than microbes, inhibiting the ability of the microbes to respond to the increases in calcium and soil pH. This idea was supported by soil mesocosm studies without plants where calcium additions resulted in marked increases in soil microbial N cycling process rates (Groffman & Fisk 2011a; Groffman & Fisk 2011b).
These results for P and N response to the calcium addition to Watershed 1 lead to the idea that plant:microbial interactions underlie the surprising curvilinear relationships between microbial biomass and activity and soil pH observed in this study. It is logical that increases in pH caused either by natural variation or calcium addition should stimulate microbial biomass and activity. However, if the increases in pH reach a certain level where plants begin to respond, the flow of C, N, and P to microbes could be reduced by plant uptake, resulting in apparent declines in microbial biomass and activity at higher pH. Thus, the complexities of plant: microbial interactions may make it appear that the microbes are adapted to/optimized for low pH while in reality their ability to respond to increases in pH is inhibited by plants.
The previous research at the HBEF highlight that there are important temporal trends in microbial biomass and activity that are not addressed in our analysis of pH controls here. Key trends include marked declines in acid and N deposition, which combined with a lengthening growing season and increases in atmosphere CO2 concentrations have led to a marked reduction in N availability (oligotrophication) at the site (Groffman et al. 2018). These temporal trends are being monitored and analyzed in other studies.
The importance of microbial community composition
Much of the effects of soil pH and acid deposition on ecosystem processes are mediated via effects on microbial community composition. Despite advances in technology and understanding, soil biota have proven to be very difficult to study and characterize. A meta-analysis of over 1300 published datapoints, found that variation in microbial biomass is predictable across biomes with microbial biomass C representing 0.6–1.1% of soil organic C and 1–20% of total plant biomass C (Fierer et al. 2009). This analysis also found that bacterial community composition and fungal: bacterial gene ratios were strongly driven by soil pH and soil C:N ratios (Fierer et al. 2009). A continental scale analysis also found that soil pH was a strong predictor of bacterial community structure (Lauber et al. 2009). This analysis found significant correlations between soil pH and the abundance of acidobacteria, actinobacteria, and bacterioidetes. A study that utilized marker gene sequencing to characterize soil bacterial communities found that 53% of microbes had predictable habitat preferences and strong links between microbial community composition and plant productivity (Delgado-Baquerizo et al. 2018). These studies suggest that there are predictable patterns in soil microbial communities that should be useful for assessing and predicting the effects of environmental changes on these communities.
Previous studies (Sridevi et al. 2012) have found significant differences in bacterial populations in the reference and calcium-treated watersheds at HBEF. Acidobacteriaceae, Comamonadaceae, and Pseudomonadaceae were lower in the Ca-amended soils, while Flavobacteriaceae and Geobacteraceae showed an opposite pattern. Ammonia-oxidizing Nitrosomonadaceae were lower in organic horizon soils of the Ca-amended watershed and an opposite pattern was observed in the mineral soil. The decreases in Acidobacteriacease and Nitrosomonadaceae are consistent with the concept that microbes respond to increases in pH as expected, but that plant domination of the N cycle constrains this response. The presence of Acidobacteria in HBEF soils is important because these bacterial taxa play a major role in the function of the forest soil ecosystem and are viable in acidic, nutrient-limited environments such as forests at HBEF. The presence of specific bacteria can be indicative of how pH impacts the bio-availability of soil nutrients and can potentially be a limiting agent in nutrient pathways.
A metagenomics survey in the reference watershed at the Hubbard Brook focused on the genes involved in denitrification, a process that responds to N enrichment, found evidence that acidity affected the communities that carry out this process (Roco et al. 2019). It has long been recognized that low soil pH inhibits expression of the genes responsible for the last step of the denitrification process (Payne 1981; Zumft 1997). The more recent study showed a series of adaptations to acidity, including removal of toxic intermediates in the denitrification process (nitric oxide) that suggest significant capacity for adaptation to highly acid conditions in HBEF soils (Roco et al. 2019). Further evidence for adaptive capacity was found under experimental soil warming and increased soil freeze/thaw cycles which suppressed bacterial taxa with the genetic potential for the final steps of denitrification (Garcia et al. 2020). This analysis suggested that a key to adaptation to new conditions in HBEF soils is the emergence of taxa that trade-off growth for stress tolerance traits. Another study found that winter snowpack decline associated with global warming increased bacterial richness and phylogenetic diversity (2016). However, any effects of changing soil climate conditions may have been mediated by effects on root growth, which is impaired by soil freezing associated with the snowpack decline (Sorensen et al. 2019).
Much of the microbial community response to changing environmental conditions may be focused in “microzones” or “microbial hotspots” of higher pH in the generally acidic soil matrix at the HBEF (Kuzyakov & Blagodatskaya 2015). These microzones could be associated with Ca-containing mineral particles, or plant roots and could serve as refugia for pH sensitive microbes. Detailed analysis of hotspot characteristics would be required to evaluate this idea.
Regulation of microbial biomass and activity in forest soils
It is undeniable that pH has a strong effect on soil microbial biomass, activity and community composition. However, our results suggest that these effects are strongly mediated by plants and that evaluations of responses of the soil microbial community to environmental change must consider plant: microbial interactions. These interactions greatly complicate assessment and prediction of how ecosystem C and N cycling will respond to environmental change, which will depend on the varying resistance or resilience of taxa within groups or “functional guilds” of microbes that cycle C and N. At the same time, plant responses to environmental change can be more dynamic than environmental responses, as they have access to both aboveground and belowground resources and the ability to shift allocation of these resources in response to environmental stress. Thus, our initial hypothesis that it is useful to view soil pH as a “master variable” controlling belowground biogeochemical processes may be only partially correct. Rather, the effects of pH, and other environmental variables must be evaluated in the context of whole-ecosystem, multi-component interactions and dynamics.