Factors influencing organic-horizon carbon pools in mixed-species stands of central Maine, USA

•A function for predicting O-horizon C content from O-horizon depth was developed.•Using the function, O-horizon C content was calculated for a long-term study site.•Average O-horizon C content did not differ among forest management treatments.•However, variation in O-horizon C content among and within stands was significant.•Most of the variation in O-horizon C content was due to microsite-level factors. The overall goal of this study was to evaluate the correlation of multiple abiotic and biotic factors with organic-horizon (O-horizon) carbon (C) content on the Penobscot Experimental Forest in central Maine, USA. O-horizon samples were collected and their associated depths were recorded from stands managed with a range of silvicultural and harvesting treatments (i.e., selection, shelterwood, and commercial clearcut) and an unmanaged control. The overall mean for O-horizon C content from all samples was 25.6±16.1Mgha−1 (mean±SD). The samples were used to develop a pedotransfer function for predicting O-horizon C content from O-horizon depth (R2=0.47, RMSE=1.6Mgha−1) so that an average of O-horizon C content could be calculated for permanent sample plots on which abiotic and biotic factors were quantified. O-horizon depth measurements recorded along transects on permanent sample plots were used to calculate plot average O-horizon C content. There were no significant differences in average predicted O-horizon C content among selection, shelterwood, and commercial clearcut treatments. However, variation in predicted O-horizon C content between stands where the same treatment was applied was statistically significant and was likely due to the timing of harvests and abundance of dead wood buried within O horizons. Depth to redoximorphic features, cartographic depth to water table or saturated zone, drainage class, fine woody debris mass, downed woody debris volume, tree basal area, and the relative basal area of conifer species were not significant predictors of predicted O-horizon C content at the plot level. When the individual predicted values of O-horizon C content were modeled, within-plot variation accounted for 83.8% of the variance. Bryophyte mass, which was predicted from bryophyte cover, only explained 1.2% of the variation in O-horizon C content at the microsite level. These results highlight the sizeable variability in O-horizon C content within and among these mixed-species stands with various forest management and natural disturbance histories.