Nonparametric approaches for estimating regional lake nutrient thresholds

Nonparametric approaches including a classification and regression tree (CART), a nonparametric changepoint analysis (nCPA) and a Bayesian hierarchical modeling (BHM) method were developed to determine ecoregional nutrient response thresholds. A CART analysis revealed that hierarchical structure was important for predicting Chl a concentrations from total nitrogen (TN) and total phosphorus (TP). The nCPA and BHM methods confirmed the CART results for each node in the tree, and the 90% confidence interval for each threshold was calculated to quantify uncertainty. The CART, nCPA, and BHM methods suggested that the nutrient criteria differed significantly within certain nutrient ecoregions and that numerical nutrient criteria of 0.0150-0.222 mg/L TP and 0.300-1.766 mg/L TN may control Chl a concentrations in the various lake ecoregions. The results of this analysis suggest that the integration of CART, nCPA and BHM might be useful for determining nutrient thresholds.