US Watershed Carbon - Land Properities
Katie Fankhauser
2024-06-04
Change log
- Summarized land cover from 2021 NLCD product
- Correlations between seasonal TOC summaries and land cover
- TOC regressed on streamflow and percent cropland
Lab Data (Yampa)
Land Use Land Cover (NLCD)
Downloaded from NLCD 2021 from MRLC Viewer: https://www.mrlc.gov/viewer/
## [1] "Demo for sensor: Phippsburg - Sonde"
## LocationName lulc n_pixels pct
## <char> <fctr> <num> <num>
## 1: Morrison-Yampa - Sonde cultivated 88 0.0039
## 2: Howelsen tunnel 1 - Sonde cultivated 282 0.0113
## 3: Upstream of Stagecoach 1 - Sonde cultivated 733 0.0338
## 4: Downstream Catamount - Sonde cultivated 1007 0.0586
## 5: Upstream of Stagecoach 2 - Sonde cultivated 1355 0.0612
## 6: Yampa - Sonde cultivated 3165 0.0693
## 7: Phippsburg - Sonde cultivated 3548 0.1451
## 8: Chuck Lewis 1 - Sonde cultivated 8498 0.3842Correlations
Include flag for whether sensor is directly downstream of reservoir. Morrison-Yampa is downstream of Stagecoach and Downstream Catamount - Sonde is downstream of Catamount.
Streamflow
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: pred_value ~ streamflow + downstream_of_reservoir + (1 | mw_id)
## Data: toc_yampa_tx
##
## REML criterion at convergence: 341
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.1319 -0.5331 0.0097 0.5221 2.7839
##
## Random effects:
## Groups Name Variance Std.Dev.
## mw_id (Intercept) 1.200 1.096
## Residual 0.239 0.488
## Number of obs: 211, groups: mw_id, 9
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 4.1603 0.4230 7.3137 9.84 0.000018 ***
## streamflow 0.1855 0.0171 204.5623 10.85 < 2e-16 ***
## downstream_of_reservoirY 1.1752 0.8825 6.8518 1.33 0.23
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) strmfl
## streamflow -0.179
## dwnstrm_f_Y -0.467 0.018
Streamflow moderated by percent cultivated land (pasture/hay or cultivated crops)
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: pred_value ~ streamflow + pct_cultivated + (1 | mw_id)
## Data: toc_yampa_tx_lulc
##
## REML criterion at convergence: 339
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.1438 -0.5303 0.0072 0.5312 2.7769
##
## Random effects:
## Groups Name Variance Std.Dev.
## mw_id (Intercept) 1.384 1.176
## Residual 0.239 0.488
## Number of obs: 211, groups: mw_id, 9
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 4.6521 0.4971 7.1184 9.36 3e-05 ***
## streamflow 0.1858 0.0171 204.1989 10.86 <2e-16 ***
## pct_cultivated -2.7134 3.4583 6.8210 -0.78 0.46
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) strmfl
## streamflow -0.138
## pct_cultvtd -0.593 -0.011Streamflow moderated by all land cover classes
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: pred_value ~ streamflow + pct_water + pct_developed + pct_forest +
## pct_shrubland + pct_herbaceous + pct_cultivated + (1 | mw_id)
## Data: toc_yampa_tx_lulc
##
## REML criterion at convergence: 299
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.185 -0.534 0.010 0.529 2.747
##
## Random effects:
## Groups Name Variance Std.Dev.
## mw_id (Intercept) 1.812 1.346
## Residual 0.238 0.488
## Number of obs: 211, groups: mw_id, 9
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 5.0066 1.3511 2.0075 3.71 0.065 .
## streamflow 0.1889 0.0172 201.1958 10.99 <2e-16 ***
## pct_water -11.3204 24.3232 1.9991 -0.47 0.687
## pct_developed -13.9591 9.1256 2.0260 -1.53 0.264
## pct_forest 1.6438 3.7576 1.9998 0.44 0.704
## pct_shrubland 1.5648 5.2032 1.9984 0.30 0.792
## pct_herbaceous -689.2468 1387.9019 1.9974 -0.50 0.669
## pct_cultivated 0.8093 5.5626 2.0020 0.15 0.898
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) strmfl pct_wt pct_dv pct_fr pct_sh pct_hr
## streamflow -0.045
## pct_water -0.058 0.003
## pct_develpd -0.151 -0.048 0.708
## pct_forest -0.364 -0.003 -0.810 -0.664
## pct_shrblnd -0.274 0.023 -0.132 -0.148 0.188
## pct_herbacs 0.004 -0.009 0.661 0.612 -0.664 -0.682
## pct_cultvtd -0.267 0.004 -0.557 -0.604 0.655 0.124 -0.521
## fit warnings:
## Some predictor variables are on very different scales: consider rescaling