1 Set up

1.1 1) Load packages

1.2 2) Import data

1.3 3) Data Wrangling

1.4 4) Depth Wrangling

1.5 5) Theme James

1.6 Comments

2 Biological Graphs

2.1 B-Glucosidase

2.1.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        406.69  203.34     2 19.197  1.4314 0.26336  
## precip         762.23  762.23     1 29.735  5.3656 0.02762 *
## landuse:precip  55.90   27.95     2 19.197  0.1967 0.82304  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2919.0 1459.50     2    27 29.2361 1.753e-07 ***
## precip         2990.0 1495.00     2    27 29.9472 1.403e-07 ***
## landuse:precip 1979.9  494.98     4    27  9.9152 4.539e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean   SE   df lower.CL upper.CL
##  472      18.8 2.04 13.5     14.4     23.2
##  579      40.2 2.04 13.5     35.8     44.5
##  850      23.9 2.04 13.5     19.5     28.3
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  579      40.2 2.04 13.5     35.8     44.5  a    
##  850      23.9 2.04 13.5     19.5     28.3   b   
##  472      18.8 2.04 13.5     14.4     23.2   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE   df lower.CL upper.CL .group
## 1    579 40.15929 2.039632 13.5 35.76946 44.54911     a 
## 3    850 23.88515 2.039632 13.5 19.49532 28.27497      b
## 2    472 18.78914 2.039632 13.5 14.39932 23.17897      b

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        132.975  66.488     2 11.647  1.6235 0.2388
## precip          62.711  62.711     1 29.377  1.5313 0.2257
## landuse:precip 125.403  62.701     2 10.936  1.5311 0.2593

2.1.2 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        77.959  38.980     2 21.012  4.9186 0.017694 * 
## precip         85.348  85.348     1 28.016 10.7696 0.002765 **
## landuse:precip 34.545  17.272     2 21.012  2.1795 0.137998   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        165.906  82.953     2    27 13.9815 6.799e-05 ***
## precip         214.256 107.128     2    27 18.0562 1.052e-05 ***
## landuse:precip  88.463  22.116     4    27  3.7275   0.01534 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        8.65 0.703 24.3     7.20     10.1
##  EA        9.13 0.703 24.3     7.68     10.6
##  NP       13.42 0.703 24.3    11.97     14.9
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       13.42 0.703 24.3    11.97     14.9  a    
##  EA        9.13 0.703 24.3     7.68     10.6   b   
##  AG        8.65 0.703 24.3     7.20     10.1   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse    lsmean        SE   df  lower.CL upper.CL .group
## 1      NP 13.424425 0.7031503 24.3 11.974142 14.87471     a 
## 2      EA  9.125530 0.7031503 24.3  7.675246 10.57581      b
## 3      AG  8.652411 0.7031503 24.3  7.202128 10.10269      b

##  precip lsmean    SE   df lower.CL upper.CL
##  472      8.22 0.703 13.5     6.71     9.73
##  579     13.81 0.703 13.5    12.29    15.32
##  850      9.17 0.703 13.5     7.66    10.69
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean    SE   df lower.CL upper.CL .group
##  579     13.81 0.703 13.5    12.29    15.32  a    
##  850      9.17 0.703 13.5     7.66    10.69   b   
##  472      8.22 0.703 13.5     6.71     9.73   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip    lsmean        SE   df  lower.CL  upper.CL .group
## 1    579 13.806733 0.7031503 13.5 12.293368 15.320097     a 
## 3    850  9.174260 0.7031503 13.5  7.660896 10.687625      b
## 2    472  8.221373 0.7031503 13.5  6.708009  9.734738      b

2.2 Glucosaminidase

2.2.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse         0.174   0.087     2 21.216  0.0105 0.98958  
## precip         48.786  48.786     1 29.692  5.8617 0.02180 *
## landuse:precip 66.237  33.118     2 21.216  3.9793 0.03408 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        1021.33  510.66     2 19.493 67.2371 1.786e-09 ***
## precip           61.41   30.71     2  6.888  4.0429   0.06895 .  
## landuse:precip   96.68   24.17     4 19.493  3.1823   0.03617 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean   SE   df lower.CL upper.CL
##  472      12.0 1.28 7.59     9.07     15.0
##  579      13.3 1.28 7.59    10.35     16.3
##  850      10.1 1.28 7.59     7.15     13.1
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  579      13.3 1.28 7.59    10.35     16.3  a    
##  472      12.0 1.28 7.59     9.07     15.0  ab   
##  850      10.1 1.28 7.59     7.15     13.1   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean      SE       df  lower.CL upper.CL .group
## 3    579 13.31379 1.27514 7.586872 10.345208 16.28237     a 
## 1    472 12.04215 1.27514 7.586872  9.073567 15.01073     ab
## 2    850 10.12164 1.27514 7.586872  7.153057 13.09022      b

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse         0.174   0.087     2 21.216  0.0105 0.98958  
## precip         48.786  48.786     1 29.692  5.8617 0.02180 *
## landuse:precip 66.237  33.118     2 21.216  3.9793 0.03408 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        1021.33  510.66     2 19.493 67.2371 1.786e-09 ***
## precip           61.41   30.71     2  6.888  4.0429   0.06895 .  
## landuse:precip   96.68   24.17     4 19.493  3.1823   0.03617 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean  SE   df lower.CL upper.CL
##  AG        8.53 1.7 18.7     4.97    12.10
##  EA       11.22 1.7 18.7     7.66    14.78
##  NP       16.37 1.7 18.7    12.81    19.94
## 
## precip = 579:
##  landuse lsmean  SE   df lower.CL upper.CL
##  AG        5.69 1.7 18.7     2.13     9.26
##  EA       14.10 1.7 18.7    10.53    17.66
##  NP       20.15 1.7 18.7    16.59    23.71
## 
## precip = 850:
##  landuse lsmean  SE   df lower.CL upper.CL
##  AG        2.11 1.7 18.7    -1.45     5.67
##  EA        9.33 1.7 18.7     5.77    12.89
##  NP       18.92 1.7 18.7    15.36    22.49
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean  SE   df lower.CL upper.CL .group
##  NP       16.37 1.7 18.7    12.81    19.94  a    
##  EA       11.22 1.7 18.7     7.66    14.78   b   
##  AG        8.53 1.7 18.7     4.97    12.10   b   
## 
## precip = 579:
##  landuse lsmean  SE   df lower.CL upper.CL .group
##  NP       20.15 1.7 18.7    16.59    23.71  a    
##  EA       14.10 1.7 18.7    10.53    17.66   b   
##  AG        5.69 1.7 18.7     2.13     9.26    c  
## 
## precip = 850:
##  landuse lsmean  SE   df lower.CL upper.CL .group
##  NP       18.92 1.7 18.7    15.36    22.49  a    
##  EA        9.33 1.7 18.7     5.77    12.89   b   
##  AG        2.11 1.7 18.7    -1.45     5.67    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean       SE       df  lower.CL  upper.CL .group
## 1      NP    472 16.374757 1.700531 18.68356 12.811419 19.938094     a 
## 2      EA    472 11.218369 1.700531 18.68356  7.655032 14.781707      b
## 3      AG    472  8.533324 1.700531 18.68356  4.969986 12.096662      b
## 4      NP    579 20.151593 1.700531 18.68356 16.588256 23.714931    a  
## 5      EA    579 14.095959 1.700531 18.68356 10.532622 17.659297     b 
## 6      AG    579  5.693820 1.700531 18.68356  2.130482  9.257157      c
## 7      NP    850 18.924386 1.700531 18.68356 15.361048 22.487724    a  
## 8      EA    850  9.331009 1.700531 18.68356  5.767671 12.894346     b 
## 9      AG    850  2.109527 1.700531 18.68356 -1.453810  5.672865      c
##   precip landuse    lsmean       SE       df  lower.CL  upper.CL .group
## 1    472      NP 16.374757 1.700531 18.68356 12.811419 19.938094     a 
## 2    472      EA 11.218369 1.700531 18.68356  7.655032 14.781707      b
## 3    472      AG  8.533324 1.700531 18.68356  4.969986 12.096662      b
## 4    579      NP 20.151593 1.700531 18.68356 16.588256 23.714931    a  
## 5    579      EA 14.095959 1.700531 18.68356 10.532622 17.659297     b 
## 6    579      AG  5.693820 1.700531 18.68356  2.130482  9.257157      c
## 7    850      NP 18.924386 1.700531 18.68356 15.361048 22.487724    a  
## 8    850      EA  9.331009 1.700531 18.68356  5.767671 12.894346     b 
## 9    850      AG  2.109527 1.700531 18.68356 -1.453810  5.672865      c
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

2.2.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        28.497  14.248     2 24.271  7.1049  0.003724 ** 
## precip          0.013   0.013     1 28.446  0.0066  0.935626    
## landuse:precip 76.829  38.414     2 23.760 19.1554 1.111e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        28.497  14.248     2 24.271  7.1049  0.003724 ** 
## precip          0.013   0.013     1 28.446  0.0066  0.935626    
## landuse:precip 76.829  38.414     2 23.760 19.1554 1.111e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        4.76 0.813 11.6     2.98     6.54
##  EA        4.69 0.722 10.1     3.08     6.30
##  NP       11.34 0.722 10.1     9.73    12.95
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       11.34 0.722 10.1     9.73    12.95  a    
##  AG        4.76 0.813 11.6     2.98     6.54   b   
##  EA        4.69 0.722 10.1     3.08     6.30   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse    lsmean        SE       df lower.CL  upper.CL .group
## 2      NP 11.340304 0.7218120 10.05769 9.733256 12.947352     a 
## 1      AG  4.761182 0.8125304 11.63811 2.984705  6.537658      b
## 3      EA  4.687993 0.7218120 10.05769 3.080945  6.295041      b

## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF   DenDF F value   Pr(>F)    
## landuse        316.19 158.094     2 19.0440 72.1280  1.3e-09 ***
## precip           0.05   0.025     2  9.1102  0.0114 0.988725    
## landuse:precip  69.00  17.250     4 15.9060  7.8700 0.001066 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        6.97 1.15 11.4    4.455     9.48
##  EA        4.24 1.15 11.4    1.725     6.76
##  NP        9.48 1.15 11.4    6.966    12.00
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        5.40 1.15 11.4    2.881     7.91
##  EA        4.87 1.15 11.4    2.358     7.39
##  NP       10.64 1.15 11.4    8.130    13.16
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        2.12 1.15 11.4   -0.393     4.64
##  EA        5.02 1.15 11.4    2.504     7.53
##  NP       13.96 1.15 11.4   11.448    16.48
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        9.48 1.15 11.4    6.966    12.00  a    
##  AG        6.97 1.15 11.4    4.455     9.48  ab   
##  EA        4.24 1.15 11.4    1.725     6.76   b   
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       10.64 1.15 11.4    8.130    13.16  a    
##  AG        5.40 1.15 11.4    2.881     7.91   b   
##  EA        4.87 1.15 11.4    2.358     7.39   b   
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       13.96 1.15 11.4   11.448    16.48  a    
##  EA        5.02 1.15 11.4    2.504     7.53   b   
##  AG        2.12 1.15 11.4   -0.393     4.64    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean       SE       df   lower.CL  upper.CL .group
## 2      NP    472  9.481236 1.148036 11.44541  6.9663719 11.996100     a 
## 1      AG    472  6.969692 1.148036 11.44541  4.4548283  9.484556     ab
## 3      EA    472  4.240214 1.148036 11.44541  1.7253496  6.755077      b
## 5      NP    579 10.644967 1.148036 11.44541  8.1301032 13.159831     a 
## 4      AG    579  5.395481 1.148036 11.44541  2.8806173  7.910345      b
## 6      EA    579  4.872860 1.148036 11.44541  2.3579962  7.387724      b
## 7      NP    850 13.962659 1.148036 11.44541 11.4477952 16.477523    a  
## 8      EA    850  5.018855 1.148036 11.44541  2.5039916  7.533719     b 
## 9      AG    850  2.122223 1.148036 11.44541 -0.3926413  4.637087      c
##   precip landuse    lsmean       SE       df   lower.CL  upper.CL .group
## 1    472      NP  9.481236 1.148036 11.44541  6.9663719 11.996100     a 
## 2    472      AG  6.969692 1.148036 11.44541  4.4548283  9.484556     ab
## 3    472      EA  4.240214 1.148036 11.44541  1.7253496  6.755077      b
## 4    579      NP 10.644967 1.148036 11.44541  8.1301032 13.159831     a 
## 5    579      AG  5.395481 1.148036 11.44541  2.8806173  7.910345      b
## 6    579      EA  4.872860 1.148036 11.44541  2.3579962  7.387724      b
## 7    850      NP 13.962659 1.148036 11.44541 11.4477952 16.477523    a  
## 8    850      EA  5.018855 1.148036 11.44541  2.5039916  7.533719     b 
## 9    850      AG  2.122223 1.148036 11.44541 -0.3926413  4.637087      c
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

2.2.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         3.751  1.8756     2    30  1.5555   0.22764    
## precip          4.931  4.9309     1    30  4.0892   0.05215 .  
## landuse:precip 33.375 16.6873     2    30 13.8389 5.517e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF  F value    Pr(>F)    
## landuse        265.143 132.571     2    27 127.5786 1.745e-14 ***
## precip          10.113   5.057     2    27   4.8662 0.0156764 *  
## landuse:precip  36.310   9.078     4    27   8.7357 0.0001166 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        3.59 0.51 27     2.55     4.64
##  EA        2.47 0.51 27     1.43     3.52
##  NP        6.50 0.51 27     5.45     7.54
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        3.16 0.51 27     2.12     4.21
##  EA        3.86 0.51 27     2.82     4.91
##  NP        8.96 0.51 27     7.92    10.01
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        2.26 0.51 27     1.21     3.31
##  EA        2.75 0.51 27     1.71     3.80
##  NP       10.87 0.51 27     9.82    11.91
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        6.50 0.51 27     5.45     7.54  a    
##  AG        3.59 0.51 27     2.55     4.64   b   
##  EA        2.47 0.51 27     1.43     3.52   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        8.96 0.51 27     7.92    10.01  a    
##  EA        3.86 0.51 27     2.82     4.91   b   
##  AG        3.16 0.51 27     2.12     4.21   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       10.87 0.51 27     9.82    11.91  a    
##  EA        2.75 0.51 27     1.71     3.80   b   
##  AG        2.26 0.51 27     1.21     3.31   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean        SE df lower.CL  upper.CL .group
## 2      NP    472  6.495736 0.5096899 27 5.449939  7.541534     a 
## 1      AG    472  3.592756 0.5096899 27 2.546958  4.638553      b
## 3      EA    472  2.473458 0.5096899 27 1.427661  3.519255      b
## 4      NP    579  8.962847 0.5096899 27 7.917049 10.008644     a 
## 5      EA    579  3.863588 0.5096899 27 2.817790  4.909385      b
## 6      AG    579  3.163233 0.5096899 27 2.117436  4.209030      b
## 7      NP    850 10.865018 0.5096899 27 9.819221 11.910815     a 
## 8      EA    850  2.752753 0.5096899 27 1.706956  3.798551      b
## 9      AG    850  2.259754 0.5096899 27 1.213957  3.305552      b
##   precip landuse    lsmean        SE df lower.CL  upper.CL .group  location
## 1    472      NP  6.495736 0.5096899 27 5.449939  7.541534     a    Tribune
## 2    472      AG  3.592756 0.5096899 27 2.546958  4.638553      b   Tribune
## 3    472      EA  2.473458 0.5096899 27 1.427661  3.519255      b   Tribune
## 4    579      NP  8.962847 0.5096899 27 7.917049 10.008644     a       Hays
## 5    579      EA  3.863588 0.5096899 27 2.817790  4.909385      b      Hays
## 6    579      AG  3.163233 0.5096899 27 2.117436  4.209030      b      Hays
## 7    850      NP 10.865018 0.5096899 27 9.819221 11.910815     a  Manhattan
## 8    850      EA  2.752753 0.5096899 27 1.706956  3.798551      b Manhattan
## 9    850      AG  2.259754 0.5096899 27 1.213957  3.305552      b Manhattan

2.3 Acid Phosphatase

2.3.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse         9933.8  4966.9     2 19.488  2.8414 0.08263 .
## precip         13022.4 13022.4     1 29.633  7.4498 0.01057 *
## landuse:precip  9553.6  4776.8     2 19.488  2.7327 0.08992 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF   DenDF F value    Pr(>F)    
## landuse        265172  132586     2 20.2421 355.867 < 2.2e-16 ***
## precip          32076   16038     2  8.2925  43.047 4.175e-05 ***
## landuse:precip  39710    9927     4 20.2421  26.646 8.166e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean   SE   df lower.CL upper.CL
##  472       110 6.73 9.68     94.7      125
##  579       182 6.73 9.68    167.2      197
##  850       119 6.73 9.68    104.3      134
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  579       182 6.73 9.68    167.2      197  a    
##  850       119 6.73 9.68    104.3      134   b   
##  472       110 6.73 9.68     94.7      125   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE       df  lower.CL upper.CL .group
## 1    579 182.2288 6.734382 9.676881 167.15546 197.3021     a 
## 3    850 119.3474 6.734382 9.676881 104.27409 134.4208      b
## 2    472 109.7360 6.734382 9.676881  94.66269 124.8094      b

2.3.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        9352.2  4676.1     2    30  6.2452 0.005401 **
## precip         4163.5  4163.5     1    30  5.5605 0.025084 * 
## landuse:precip 6416.6  3208.3     2    30  4.2849 0.023074 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         85297   42649     2    27 97.4347 4.478e-13 ***
## precip           7908    3954     2    27  9.0333 0.0009917 ***
## landuse:precip  13316    3329     4    27  7.6056 0.0003070 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        46.5 6.04 14.8     33.6     59.4
##  EA        60.2 6.04 24.8     47.7     72.6
##  NP       155.9 6.04 24.8    143.5    168.4
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       155.9 6.04 24.8    143.5    168.4  a    
##  EA        60.2 6.04 24.8     47.7     72.6   b   
##  AG        46.5 6.04 14.8     33.6     59.4   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse    lsmean       SE    df  lower.CL  upper.CL .group
## 1      NP 155.92509 6.039549 24.75 143.48004 168.37015     a 
## 2      EA  60.16690 6.039549 24.75  47.72184  72.61196      b
## 3      AG  46.52442 6.039549 14.85  33.64009  59.40874      b

##  precip lsmean   SE   df lower.CL upper.CL
##  472      88.4 6.04 24.8     76.0    100.9
##  579     105.2 6.04 14.8     92.3    118.1
##  850      69.0 6.04 14.8     56.1     81.8
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  579     105.2 6.04 14.8     92.3    118.1  a    
##  472      88.4 6.04 24.8     76.0    100.9  a    
##  850      69.0 6.04 14.8     56.1     81.8   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip    lsmean       SE    df lower.CL  upper.CL .group
## 3    579 105.22521 6.039549 14.85 92.34088 118.10954     a 
## 1    472  88.43687 6.039549 24.75 75.99181 100.88193     a 
## 2    850  68.95433 6.039549 14.85 56.07001  81.83866      b

## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       75.83 10.5 27     54.4     97.3
##  EA       52.64 10.5 27     31.2     74.1
##  NP      136.85 10.5 27    115.4    158.3
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       54.09 10.5 27     32.6     75.6
##  EA       61.56 10.5 27     40.1     83.0
##  NP      200.03 10.5 27    178.6    221.5
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        9.66 10.5 27    -11.8     31.1
##  EA       66.30 10.5 27     44.8     87.8
##  NP      130.90 10.5 27    109.4    152.4
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP      136.85 10.5 27    115.4    158.3  a    
##  AG       75.83 10.5 27     54.4     97.3   b   
##  EA       52.64 10.5 27     31.2     74.1   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP      200.03 10.5 27    178.6    221.5  a    
##  EA       61.56 10.5 27     40.1     83.0   b   
##  AG       54.09 10.5 27     32.6     75.6   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP      130.90 10.5 27    109.4    152.4  a    
##  EA       66.30 10.5 27     44.8     87.8   b   
##  AG        9.66 10.5 27    -11.8     31.1    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip     lsmean       SE df  lower.CL  upper.CL .group
## 2      NP    472 136.845849 10.46081 27 115.38205 158.30965     a 
## 1      AG    472  75.825188 10.46081 27  54.36139  97.28899      b
## 3      EA    472  52.639563 10.46081 27  31.17576  74.10336      b
## 4      NP    579 200.028036 10.46081 27 178.56423 221.49184     a 
## 5      EA    579  61.561062 10.46081 27  40.09726  83.02486      b
## 6      AG    579  54.086524 10.46081 27  32.62272  75.55033      b
## 7      NP    850 130.901396 10.46081 27 109.43759 152.36520    a  
## 8      EA    850  66.300074 10.46081 27  44.83627  87.76388     b 
## 9      AG    850   9.661535 10.46081 27 -11.80227  31.12534      c
##   precip landuse     lsmean       SE df  lower.CL  upper.CL .group  location
## 1    472      NP 136.845849 10.46081 27 115.38205 158.30965     a    Tribune
## 2    472      AG  75.825188 10.46081 27  54.36139  97.28899      b   Tribune
## 3    472      EA  52.639563 10.46081 27  31.17576  74.10336      b   Tribune
## 4    579      NP 200.028036 10.46081 27 178.56423 221.49184     a       Hays
## 5    579      EA  61.561062 10.46081 27  40.09726  83.02486      b      Hays
## 6    579      AG  54.086524 10.46081 27  32.62272  75.55033      b      Hays
## 7    850      NP 130.901396 10.46081 27 109.43759 152.36520    a   Manhattan
## 8    850      EA  66.300074 10.46081 27  44.83627  87.76388     b  Manhattan
## 9    850      AG   9.661535 10.46081 27 -11.80227  31.12534      c Manhattan

2.3.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        2692.8  1346.4     2    30  2.4621 0.102314   
## precip         7199.4  7199.4     1    30 13.1651 0.001048 **
## landuse:precip  877.1   438.5     2    30  0.8019 0.457838   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         33885 16942.6     2    27 59.8329 1.197e-10 ***
## precip           9346  4673.0     2    27 16.5026 2.079e-05 ***
## landuse:precip   7491  1872.6     4    27  6.6132 0.0007607 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean   SE   df lower.CL upper.CL
##  472      75.5 4.86 13.5     65.1     86.0
##  579      84.2 4.86 13.5     73.8     94.7
##  850      46.5 4.86 13.5     36.1     57.0
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  579      84.2 4.86 13.5     73.8     94.7  a    
##  472      75.5 4.86 13.5     65.1     86.0  a    
##  850      46.5 4.86 13.5     36.1     57.0   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE   df lower.CL upper.CL .group
## 3    579 84.21449 4.857691 13.5 73.75947 94.66952     a 
## 1    472 75.53836 4.857691 13.5 65.08334 85.99339     a 
## 2    850 46.53298 4.857691 13.5 36.07795 56.98801      b

2.4 Alkaline Phosphatase

2.4.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        6274.1 3137.06     2 21.902 14.7820 8.645e-05 ***
## precip          758.2  758.22     1 21.925  3.5728   0.07203 .  
## landuse:precip 2161.8 1080.88     2 21.902  5.0932   0.01526 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF  F value    Pr(>F)    
## landuse        26411.8 13205.9     2 21.808 128.5598 8.521e-13 ***
## precip          1347.3   673.7     2 10.606   6.5582     0.014 *  
## landuse:precip  4519.6  1129.9     4 21.808  10.9997 4.880e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        24.2 3.27 18.5     17.3     31.1
##  EA        35.1 3.27 18.5     28.2     41.9
##  NP        86.3 3.27 18.5     79.5     93.2
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        86.3 3.27 18.5     79.5     93.2  a    
##  EA        35.1 3.27 18.5     28.2     41.9   b   
##  AG        24.2 3.27 18.5     17.3     31.1    c  
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE       df lower.CL upper.CL .group
## 1      NP 86.31681 3.273376 18.53298 79.45386 93.17977    a  
## 2      EA 35.08142 3.273376 18.53298 28.21846 41.94438     b 
## 3      AG 24.19335 3.273376 18.53298 17.33039 31.05631      c

## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        11.8 5.44 25.3    0.637     23.0
##  EA        18.2 5.44 25.3    7.008     29.4
##  NP        82.7 5.44 25.3   71.536     93.9
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        12.9 5.44 25.3    1.751     24.1
##  EA        56.1 5.44 25.3   44.885     67.3
##  NP        94.1 5.44 25.3   82.860    105.2
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        47.8 5.44 25.3   36.619     59.0
##  EA        31.0 5.44 25.3   19.778     42.2
##  NP        82.2 5.44 25.3   70.981     93.4
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        82.7 5.44 25.3   71.536     93.9  a    
##  EA        18.2 5.44 25.3    7.008     29.4   b   
##  AG        11.8 5.44 25.3    0.637     23.0   b   
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        94.1 5.44 25.3   82.860    105.2  a    
##  EA        56.1 5.44 25.3   44.885     67.3   b   
##  AG        12.9 5.44 25.3    1.751     24.1    c  
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        82.2 5.44 25.3   70.981     93.4  a    
##  AG        47.8 5.44 25.3   36.619     59.0   b   
##  EA        31.0 5.44 25.3   19.778     42.2    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE      df   lower.CL  upper.CL .group
## 1      NP    472 82.72739 5.436833 25.2876 71.5364806  93.91831     a 
## 2      EA    472 18.19938 5.436833 25.2876  7.0084641  29.39029      b
## 3      AG    472 11.82769 5.436833 25.2876  0.6367806  23.01861      b
## 4      NP    579 94.05087 5.436833 25.2876 82.8599553 105.24178    a  
## 5      EA    579 56.07564 5.436833 25.2876 44.8847300  67.26656     b 
## 6      AG    579 12.94232 5.436833 25.2876  1.7514103  24.13324      c
## 8      NP    850 82.17218 5.436833 25.2876 70.9812688  93.36310    a  
## 7      AG    850 47.81004 5.436833 25.2876 36.6191264  59.00095     b 
## 9      EA    850 30.96925 5.436833 25.2876 19.7783355  42.16016      c
##   precip landuse   lsmean       SE      df   lower.CL  upper.CL .group
## 1    472      NP 82.72739 5.436833 25.2876 71.5364806  93.91831     a 
## 2    472      EA 18.19938 5.436833 25.2876  7.0084641  29.39029      b
## 3    472      AG 11.82769 5.436833 25.2876  0.6367806  23.01861      b
## 4    579      NP 94.05087 5.436833 25.2876 82.8599553 105.24178    a  
## 5    579      EA 56.07564 5.436833 25.2876 44.8847300  67.26656     b 
## 6    579      AG 12.94232 5.436833 25.2876  1.7514103  24.13324      c
## 7    850      NP 82.17218 5.436833 25.2876 70.9812688  93.36310    a  
## 8    850      AG 47.81004 5.436833 25.2876 36.6191264  59.00095     b 
## 9    850      EA 30.96925 5.436833 25.2876 19.7783355  42.16016      c
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

2.4.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        5806.7 2903.36     2    27 48.8669 1.066e-09 ***
## precip          697.6  348.80     2    27  5.8708  0.007639 ** 
## landuse:precip 3257.7  814.43     4    27 13.7078 3.170e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        5806.7 2903.36     2    27 48.8669 1.066e-09 ***
## precip          697.6  348.80     2    27  5.8708  0.007639 ** 
## landuse:precip 3257.7  814.43     4    27 13.7078 3.170e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        22.9 2.23 14.8    18.17     27.7
##  EA        13.0 2.23 24.8     8.38     17.5
##  NP        43.5 2.23 24.8    38.88     48.0
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        43.5 2.23 24.8    38.88     48.0  a    
##  AG        22.9 2.23 14.8    18.17     27.7   b   
##  EA        13.0 2.23 24.8     8.38     17.5    c  
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE    df  lower.CL upper.CL .group
## 2      NP 43.46246 2.225116 24.75 38.877396 48.04752    a  
## 1      AG 22.91538 2.225116 14.85 18.168486 27.66228     b 
## 3      EA 12.96014 2.225116 24.75  8.375082 17.54520      c

##  precip lsmean   SE   df lower.CL upper.CL
##  472      22.2 2.23 24.8     17.7     26.8
##  579      24.6 2.23 14.8     19.8     29.3
##  850      32.5 2.23 14.8     27.8     37.3
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  850      32.5 2.23 14.8     27.8     37.3  a    
##  579      24.6 2.23 14.8     19.8     29.3   b   
##  472      22.2 2.23 24.8     17.7     26.8   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE    df lower.CL upper.CL .group
## 1    850 32.52672 2.225116 14.85 27.77982 37.27361     a 
## 2    579 24.56131 2.225116 14.85 19.81442 29.30821      b
## 3    472 22.24995 2.225116 24.75 17.66489 26.83501      b

## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        8.84 3.85 27    0.929     16.7
##  EA       14.07 3.85 27    6.159     22.0
##  NP       43.85 3.85 27   35.938     51.8
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       11.90 3.85 27    3.989     19.8
##  EA       15.41 3.85 27    7.501     23.3
##  NP       46.38 3.85 27   38.471     54.3
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       48.01 3.85 27   40.104     55.9
##  EA        9.40 3.85 27    1.497     17.3
##  NP       40.16 3.85 27   32.255     48.1
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       43.85 3.85 27   35.938     51.8  a    
##  EA       14.07 3.85 27    6.159     22.0   b   
##  AG        8.84 3.85 27    0.929     16.7   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       46.38 3.85 27   38.471     54.3  a    
##  EA       15.41 3.85 27    7.501     23.3   b   
##  AG       11.90 3.85 27    3.989     19.8   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  AG       48.01 3.85 27   40.104     55.9  a    
##  NP       40.16 3.85 27   32.255     48.1  a    
##  EA        9.40 3.85 27    1.497     17.3   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean       SE df   lower.CL upper.CL .group
## 1      NP    472 43.845743 3.854013 27 35.9379608 51.75353     a 
## 2      EA    472 14.066993 3.854013 27  6.1592110 21.97478      b
## 3      AG    472  8.837121 3.854013 27  0.9293387 16.74490      b
## 4      NP    579 46.378364 3.854013 27 38.4705820 54.28615     a 
## 5      EA    579 15.408482 3.854013 27  7.5006995 23.31626      b
## 6      AG    579 11.897095 3.854013 27  3.9893128 19.80488      b
## 8      AG    850 48.011935 3.854013 27 40.1041524 55.91972     a 
## 9      NP    850 40.163259 3.854013 27 32.2554765 48.07104     a 
## 7      EA    850  9.404951 3.854013 27  1.4971690 17.31273      b
##   precip landuse    lsmean       SE df   lower.CL upper.CL .group  location
## 1    472      NP 43.845743 3.854013 27 35.9379608 51.75353     a    Tribune
## 2    472      EA 14.066993 3.854013 27  6.1592110 21.97478      b   Tribune
## 3    472      AG  8.837121 3.854013 27  0.9293387 16.74490      b   Tribune
## 4    579      NP 46.378364 3.854013 27 38.4705820 54.28615     a       Hays
## 5    579      EA 15.408482 3.854013 27  7.5006995 23.31626      b      Hays
## 6    579      AG 11.897095 3.854013 27  3.9893128 19.80488      b      Hays
## 7    850      AG 48.011935 3.854013 27 40.1041524 55.91972     a  Manhattan
## 8    850      NP 40.163259 3.854013 27 32.2554765 48.07104     a  Manhattan
## 9    850      EA  9.404951 3.854013 27  1.4971690 17.31273      b Manhattan

2.5 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2020.8 1010.39     2    30 10.4876 0.0003519 ***
## precip          424.9  424.90     1    30  4.4103 0.0442347 *  
## landuse:precip 1790.5  895.26     2    30  9.2926 0.0007232 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2738.95 1369.48     2    27 16.6134 1.978e-05 ***
## precip          433.98  216.99     2    27  2.6323 0.0902803 .  
## landuse:precip 2446.00  611.50     4    27  7.4182 0.0003628 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        26.2 2.62 24.3     20.8     31.6
##  EA        18.3 2.62 24.3     12.8     23.7
##  NP        39.4 2.62 24.3     34.0     44.8
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        39.4 2.62 24.3     34.0     44.8  a    
##  AG        26.2 2.62 24.3     20.8     31.6   b   
##  EA        18.3 2.62 24.3     12.8     23.7    c  
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE   df lower.CL upper.CL .group
## 2      NP 39.40352 2.620945 24.3 33.99769 44.80935    a  
## 1      AG 26.21479 2.620945 24.3 20.80896 31.62062     b 
## 3      EA 18.25195 2.620945 24.3 12.84612 23.65779      c

##  precip lsmean   SE   df lower.CL upper.CL
##  472      31.9 2.62 13.5     26.3     37.6
##  579      28.4 2.62 13.5     22.8     34.1
##  850      23.5 2.62 13.5     17.8     29.1
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  472      31.9 2.62 13.5     26.3     37.6  a    
##  579      28.4 2.62 13.5     22.8     34.1  ab   
##  850      23.5 2.62 13.5     17.8     29.1   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE   df lower.CL upper.CL .group
## 3    472 31.94195 2.620945 13.5 26.30099 37.58292     a 
## 2    579 28.44832 2.620945 13.5 22.80736 34.08929     ab
## 1    850 23.47999 2.620945 13.5 17.83902 29.12095      b

## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       26.00 4.54 27     16.7     35.3
##  EA       24.34 4.54 27     15.0     33.7
##  NP       45.49 4.54 27     36.2     54.8
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       15.01 4.54 27      5.7     24.3
##  EA       24.21 4.54 27     14.9     33.5
##  NP       46.13 4.54 27     36.8     55.4
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       37.64 4.54 27     28.3     47.0
##  EA        6.21 4.54 27     -3.1     15.5
##  NP       26.59 4.54 27     17.3     35.9
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       45.49 4.54 27     36.2     54.8  a    
##  AG       26.00 4.54 27     16.7     35.3   b   
##  EA       24.34 4.54 27     15.0     33.7   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       46.13 4.54 27     36.8     55.4  a    
##  EA       24.21 4.54 27     14.9     33.5   b   
##  AG       15.01 4.54 27      5.7     24.3   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  AG       37.64 4.54 27     28.3     47.0  a    
##  NP       26.59 4.54 27     17.3     35.9  a    
##  EA        6.21 4.54 27     -3.1     15.5   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean      SE df  lower.CL upper.CL .group
## 2      NP    472 45.490348 4.53961 27 36.175838 54.80486     a 
## 1      AG    472 25.997263 4.53961 27 16.682753 35.31177      b
## 3      EA    472 24.338253 4.53961 27 15.023743 33.65276      b
## 4      NP    579 46.130202 4.53961 27 36.815692 55.44471     a 
## 5      EA    579 24.205126 4.53961 27 14.890616 33.51964      b
## 6      AG    579 15.009646 4.53961 27  5.695136 24.32416      b
## 8      AG    850 37.637466 4.53961 27 28.322956 46.95198     a 
## 9      NP    850 26.590013 4.53961 27 17.275503 35.90452     a 
## 7      EA    850  6.212484 4.53961 27 -3.102026 15.52699      b
##   precip landuse    lsmean      SE df  lower.CL upper.CL .group  location
## 1    472      NP 45.490348 4.53961 27 36.175838 54.80486     a    Tribune
## 2    472      AG 25.997263 4.53961 27 16.682753 35.31177      b   Tribune
## 3    472      EA 24.338253 4.53961 27 15.023743 33.65276      b   Tribune
## 4    579      NP 46.130202 4.53961 27 36.815692 55.44471     a       Hays
## 5    579      EA 24.205126 4.53961 27 14.890616 33.51964      b      Hays
## 6    579      AG 15.009646 4.53961 27  5.695136 24.32416      b      Hays
## 7    850      AG 37.637466 4.53961 27 28.322956 46.95198     a  Manhattan
## 8    850      NP 26.590013 4.53961 27 17.275503 35.90452     a  Manhattan
## 9    850      EA  6.212484 4.53961 27 -3.102026 15.52699      b Manhattan

2.6 Arylsulfatase

2.6.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        707.48  353.74     2    30 11.2066 0.0002319 ***
## precip         245.16  245.16     1    30  7.7669 0.0091407 ** 
## landuse:precip  17.44    8.72     2    30  0.2762 0.7605832    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF  F value    Pr(>F)    
## landuse        15832.6  7916.3     2 21.560 1021.497 < 2.2e-16 ***
## precip           252.5   126.2     2 10.372   16.289 0.0006306 ***
## landuse:precip   705.4   176.4     4 21.560   22.757 1.774e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        3.21 0.876 19.8     1.38     5.04
##  EA        8.45 0.876 19.8     6.62    10.28
##  NP       50.09 0.876 19.8    48.26    51.92
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       50.09 0.876 19.8    48.26    51.92  a    
##  EA        8.45 0.876 19.8     6.62    10.28   b   
##  AG        3.21 0.876 19.8     1.38     5.04    c  
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse    lsmean       SE       df  lower.CL  upper.CL .group
## 1      NP 50.086544 0.875904 19.76627 48.258055 51.915034    a  
## 2      EA  8.452328 0.875904 19.76627  6.623838 10.280818     b 
## 3      AG  3.211371 0.875904 19.76627  1.382881  5.039861      c

##  precip lsmean   SE   df lower.CL upper.CL
##  472      18.8 0.93 10.3     16.7     20.8
##  579      18.5 0.93 10.3     16.5     20.6
##  850      24.5 0.93 10.3     22.4     26.5
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  850      24.5 0.93 10.3     22.4     26.5  a    
##  472      18.8 0.93 10.3     16.7     20.8   b   
##  579      18.5 0.93 10.3     16.5     20.6   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean        SE       df lower.CL upper.CL .group
## 2    850 24.46099 0.9296914 10.27579 22.39702 26.52496     a 
## 1    472 18.76439 0.9296914 10.27579 16.70042 20.82836      b
## 3    579 18.52486 0.9296914 10.27579 16.46089 20.58883      b

2.6.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         42.16   21.08     2 21.883  2.8845   0.07729 .  
## precip         423.76  423.76     1 29.445 57.9862 1.939e-08 ***
## landuse:precip 547.18  273.59     2 21.142 37.4374 1.129e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        5991.2 2995.61     2 22.867 480.154 < 2.2e-16 ***
## precip          674.0  336.98     2 16.466  54.013 5.846e-08 ***
## landuse:precip  642.1  160.52     4 21.095  25.729 7.502e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean    SE    df lower.CL upper.CL
##  472      8.26 0.820 17.32     6.53     9.99
##  579     14.45 0.891  9.13    12.44    16.46
##  850     19.99 0.891  9.13    17.98    22.00
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean    SE    df lower.CL upper.CL .group
##  850     19.99 0.891  9.13    17.98    22.00  a    
##  579     14.45 0.891  9.13    12.44    16.46   b   
##  472      8.26 0.820 17.32     6.53     9.99    c  
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean        SE       df  lower.CL  upper.CL .group
## 1    850 19.99076 0.8906200  9.12986 17.980397 22.001122    a  
## 2    579 14.45133 0.8906200  9.12986 12.440971 16.461696     b 
## 3    472  8.25981 0.8195958 17.32099  6.533052  9.986568      c

## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        1.65 1.35 24.4   -1.146     4.44
##  EA        3.98 1.35 24.4    1.189     6.77
##  NP       19.15 1.35 24.4   16.361    21.94
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        3.47 1.35 24.4    0.679     6.26
##  EA        6.73 1.35 24.4    3.941     9.52
##  NP       33.15 1.35 24.4   30.359    35.94
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        6.71 1.35 24.4    3.923     9.51
##  EA        7.63 1.35 24.4    4.839    10.42
##  NP       45.63 1.35 24.4   42.836    48.42
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       19.15 1.35 24.4   16.361    21.94  a    
##  EA        3.98 1.35 24.4    1.189     6.77   b   
##  AG        1.65 1.35 24.4   -1.146     4.44   b   
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       33.15 1.35 24.4   30.359    35.94  a    
##  EA        6.73 1.35 24.4    3.941     9.52   b   
##  AG        3.47 1.35 24.4    0.679     6.26   b   
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       45.63 1.35 24.4   42.836    48.42  a    
##  EA        7.63 1.35 24.4    4.839    10.42   b   
##  AG        6.71 1.35 24.4    3.923     9.51   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean       SE       df   lower.CL  upper.CL .group
## 1      NP    472 19.153109 1.353888 24.42584 16.3613982 21.944820     a 
## 2      EA    472  3.980589 1.353888 24.42584  1.1888776  6.772300      b
## 3      AG    472  1.645733 1.353888 24.42584 -1.1459779  4.437444      b
## 4      NP    579 33.150908 1.353888 24.42584 30.3591970 35.942619     a 
## 5      EA    579  6.732545 1.353888 24.42584  3.9408335  9.524255      b
## 6      AG    579  3.470547 1.353888 24.42584  0.6788364  6.262258      b
## 7      NP    850 45.627650 1.353888 24.42584 42.8359394 48.419361     a 
## 8      EA    850  7.630217 1.353888 24.42584  4.8385061 10.421928      b
## 9      AG    850  6.714410 1.353888 24.42584  3.9226987  9.506121      b
##   precip landuse    lsmean       SE       df   lower.CL  upper.CL .group
## 1    472      NP 19.153109 1.353888 24.42584 16.3613982 21.944820     a 
## 2    472      EA  3.980589 1.353888 24.42584  1.1888776  6.772300      b
## 3    472      AG  1.645733 1.353888 24.42584 -1.1459779  4.437444      b
## 4    579      NP 33.150908 1.353888 24.42584 30.3591970 35.942619     a 
## 5    579      EA  6.732545 1.353888 24.42584  3.9408335  9.524255      b
## 6    579      AG  3.470547 1.353888 24.42584  0.6788364  6.262258      b
## 7    850      NP 45.627650 1.353888 24.42584 42.8359394 48.419361     a 
## 8    850      EA  7.630217 1.353888 24.42584  4.8385061 10.421928      b
## 9    850      AG  6.714410 1.353888 24.42584  3.9226987  9.506121      b
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

2.6.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         24.141  12.070     2 24.486  1.4358 0.2572630    
## precip          68.493  68.493     1 23.836  8.1474 0.0087831 ** 
## landuse:precip 178.201  89.101     2 24.486 10.5988 0.0004832 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF  F value    Pr(>F)    
## landuse        4219.7 2109.87     2 21.467 283.6221 3.663e-16 ***
## precip           65.5   32.73     2 10.466   4.3998 0.0410429 *  
## landuse:precip  220.3   55.08     4 21.467   7.4043 0.0006567 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean    SE df lower.CL upper.CL
##  472      12.9 0.871 11     11.0     14.9
##  579      13.5 0.871 11     11.6     15.4
##  850      16.2 0.871 11     14.3     18.1
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean    SE df lower.CL upper.CL .group
##  850      16.2 0.871 11     14.3     18.1  a    
##  579      13.5 0.871 11     11.6     15.4   b   
##  472      12.9 0.871 11     11.0     14.9   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean        SE      df lower.CL upper.CL .group
## 1    850 16.19983 0.8706956 11.0358 14.28420 18.11546     a 
## 2    579 13.52797 0.8706956 11.0358 11.61234 15.44360      b
## 3    472 12.94841 0.8706956 11.0358 11.03278 14.86404      b

## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        7.05 1.41 26.5    4.145     9.95
##  EA        7.87 1.41 26.5    4.968    10.77
##  NP       23.93 1.41 26.5   21.024    26.83
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        3.74 1.41 26.5    0.837     6.64
##  EA        7.37 1.41 26.5    4.469    10.27
##  NP       29.47 1.41 26.5   26.569    32.37
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        7.38 1.41 26.5    4.474    10.28
##  EA        6.05 1.41 26.5    3.145     8.95
##  NP       35.18 1.41 26.5   32.272    38.08
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       23.93 1.41 26.5   21.024    26.83  a    
##  EA        7.87 1.41 26.5    4.968    10.77   b   
##  AG        7.05 1.41 26.5    4.145     9.95   b   
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       29.47 1.41 26.5   26.569    32.37  a    
##  EA        7.37 1.41 26.5    4.469    10.27   b   
##  AG        3.74 1.41 26.5    0.837     6.64   b   
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       35.18 1.41 26.5   32.272    38.08  a    
##  AG        7.38 1.41 26.5    4.474    10.28   b   
##  EA        6.05 1.41 26.5    3.145     8.95   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean       SE       df   lower.CL  upper.CL .group
## 1      NP    472 23.927095 1.413488 26.49604 21.0242733 26.829917     a 
## 2      EA    472  7.870503 1.413488 26.49604  4.9676812 10.773325      b
## 3      AG    472  7.047637 1.413488 26.49604  4.1448147  9.950459      b
## 4      NP    579 29.471810 1.413488 26.49604 26.5689884 32.374632     a 
## 5      EA    579  7.371796 1.413488 26.49604  4.4689742 10.274618      b
## 6      AG    579  3.740303 1.413488 26.49604  0.8374805  6.643124      b
## 8      NP    850 35.175205 1.413488 26.49604 32.2723830 38.078027     a 
## 7      AG    850  7.376330 1.413488 26.49604  4.4735079 10.279152      b
## 9      EA    850  6.047956 1.413488 26.49604  3.1451338  8.950778      b
##   precip landuse    lsmean       SE       df   lower.CL  upper.CL .group
## 1    472      NP 23.927095 1.413488 26.49604 21.0242733 26.829917     a 
## 2    472      EA  7.870503 1.413488 26.49604  4.9676812 10.773325      b
## 3    472      AG  7.047637 1.413488 26.49604  4.1448147  9.950459      b
## 4    579      NP 29.471810 1.413488 26.49604 26.5689884 32.374632     a 
## 5    579      EA  7.371796 1.413488 26.49604  4.4689742 10.274618      b
## 6    579      AG  3.740303 1.413488 26.49604  0.8374805  6.643124      b
## 7    850      NP 35.175205 1.413488 26.49604 32.2723830 38.078027     a 
## 8    850      AG  7.376330 1.413488 26.49604  4.4735079 10.279152      b
## 9    850      EA  6.047956 1.413488 26.49604  3.1451338  8.950778      b
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

2.7 Phosphodiesterase

2.7.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         819.72  409.86     2    30  14.433 4.062e-05 ***
## precip         1582.48 1582.48     1    30  55.727 2.549e-08 ***
## landuse:precip   74.62   37.31     2    30   1.314    0.2838    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF   DenDF  F value    Pr(>F)    
## landuse        14835.3  7417.6     2 20.4399 359.6981 < 2.2e-16 ***
## precip          1377.0   688.5     2  8.9759  33.3867 6.963e-05 ***
## landuse:precip   307.1    76.8     4 20.4399   3.7232   0.01975 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        15.4 1.39 21.3     12.5     18.2
##  EA        19.5 1.39 21.3     16.6     22.4
##  NP        60.3 1.39 21.3     57.4     63.2
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        60.3 1.39 21.3     57.4     63.2  a    
##  EA        19.5 1.39 21.3     16.6     22.4   b   
##  AG        15.4 1.39 21.3     12.5     18.2    c  
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE       df lower.CL upper.CL .group
## 1      NP 60.32737 1.387609 21.26087 57.44383 63.21090    a  
## 2      EA 19.47552 1.387609 21.26087 16.59198 22.35906     b 
## 3      AG 15.35040 1.387609 21.26087 12.46687 18.23394      c

##  precip lsmean   SE   df lower.CL upper.CL
##  472      24.2 1.45 11.1     21.0     27.4
##  579      30.5 1.45 11.1     27.4     33.7
##  850      40.4 1.45 11.1     37.2     43.6
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  850      40.4 1.45 11.1     37.2     43.6  a    
##  579      30.5 1.45 11.1     27.4     33.7   b   
##  472      24.2 1.45 11.1     21.0     27.4    c  
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE      df lower.CL upper.CL .group
## 1    850 40.42103 1.446041 11.0832 37.24122 43.60083    a  
## 2    579 30.54402 1.446041 11.0832 27.36422 33.72383     b 
## 3    472 24.18824 1.446041 11.0832 21.00844 27.36804      c

2.7.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         205.01  102.50     2 25.987  3.7910  0.035922 *  
## precip         1638.23 1638.23     1 29.954 60.5880 1.112e-08 ***
## landuse:precip  420.44  210.22     2 25.483  7.7748  0.002312 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        3484.2 1742.10     2 24.626 107.858 6.563e-13 ***
## precip         1652.1  826.06     2 21.224  51.144 7.640e-09 ***
## landuse:precip  766.1  191.52     4 23.732  11.858 1.953e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        15.6 1.36 10.1    12.58     18.6
##  EA        11.5 1.27 19.2     8.89     14.2
##  NP        34.2 1.27 19.2    31.55     36.9
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        34.2 1.27 19.2    31.55     36.9  a    
##  AG        15.6 1.36 10.1    12.58     18.6   b   
##  EA        11.5 1.27 19.2     8.89     14.2    c  
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE       df  lower.CL upper.CL .group
## 2      NP 34.21412 1.272457 19.16915 31.552430 36.87581    a  
## 1      AG 15.60217 1.355610 10.05137 12.583769 18.62056     b 
## 3      EA 11.54995 1.272457 19.16915  8.888262 14.21165      c

##  precip lsmean   SE   df lower.CL upper.CL
##  472      12.5 1.27 19.2      9.8     15.1
##  579      19.0 1.36 10.1     16.0     22.1
##  850      29.9 1.36 10.1     26.9     32.9
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  850      29.9 1.36 10.1     26.9     32.9  a    
##  579      19.0 1.36 10.1     16.0     22.1   b   
##  472      12.5 1.27 19.2      9.8     15.1    c  
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE       df  lower.CL upper.CL .group
## 1    850 29.87663 1.355610 10.05137 26.858234 32.89503    a  
## 2    579 19.03206 1.355610 10.05137 16.013660 22.05045     b 
## 3    472 12.45756 1.272457 19.16915  9.795863 15.11925      c

## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        9.32 2.13 25.6     4.94     13.7
##  EA        8.57 2.13 25.6     4.19     12.9
##  NP       19.49 2.13 25.6    15.11     23.9
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        8.34 2.13 25.6     3.96     12.7
##  EA       11.50 2.13 25.6     7.12     15.9
##  NP       37.26 2.13 25.6    32.88     41.6
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG       29.15 2.13 25.6    24.77     33.5
##  EA       14.59 2.13 25.6    10.21     19.0
##  NP       45.90 2.13 25.6    41.52     50.3
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       19.49 2.13 25.6    15.11     23.9  a    
##  AG        9.32 2.13 25.6     4.94     13.7   b   
##  EA        8.57 2.13 25.6     4.19     12.9   b   
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       37.26 2.13 25.6    32.88     41.6  a    
##  EA       11.50 2.13 25.6     7.12     15.9   b   
##  AG        8.34 2.13 25.6     3.96     12.7   b   
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP       45.90 2.13 25.6    41.52     50.3  a    
##  AG       29.15 2.13 25.6    24.77     33.5   b   
##  EA       14.59 2.13 25.6    10.21     19.0    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean       SE       df  lower.CL upper.CL .group
## 2      NP    472 19.488273 2.128294 25.55399 15.109782 23.86676     a 
## 1      AG    472  9.316488 2.128294 25.55399  4.937997 13.69498      b
## 3      EA    472  8.567905 2.128294 25.55399  4.189414 12.94640      b
## 4      NP    579 37.257613 2.128294 25.55399 32.879123 41.63610     a 
## 5      EA    579 11.496888 2.128294 25.55399  7.118398 15.87538      b
## 6      AG    579  8.341668 2.128294 25.55399  3.963178 12.72016      b
## 8      NP    850 45.896480 2.128294 25.55399 41.517990 50.27497    a  
## 7      AG    850 29.148341 2.128294 25.55399 24.769850 33.52683     b 
## 9      EA    850 14.585071 2.128294 25.55399 10.206580 18.96356      c
##   precip landuse    lsmean       SE       df  lower.CL upper.CL .group
## 1    472      NP 19.488273 2.128294 25.55399 15.109782 23.86676     a 
## 2    472      AG  9.316488 2.128294 25.55399  4.937997 13.69498      b
## 3    472      EA  8.567905 2.128294 25.55399  4.189414 12.94640      b
## 4    579      NP 37.257613 2.128294 25.55399 32.879123 41.63610     a 
## 5    579      EA 11.496888 2.128294 25.55399  7.118398 15.87538      b
## 6    579      AG  8.341668 2.128294 25.55399  3.963178 12.72016      b
## 7    850      NP 45.896480 2.128294 25.55399 41.517990 50.27497    a  
## 8    850      AG 29.148341 2.128294 25.55399 24.769850 33.52683     b 
## 9    850      EA 14.585071 2.128294 25.55399 10.206580 18.96356      c
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

2.7.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        327.34  163.67     2    30  3.7784 0.03440 *
## precip         254.75  254.75     1    30  5.8810 0.02154 *
## landuse:precip 356.95  178.47     2    30  4.1201 0.02624 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        1749.81  874.90     2    27 31.2774 9.337e-08 ***
## precip          281.36  140.68     2    27  5.0292 0.0139137 *  
## landuse:precip  874.60  218.65     4    27  7.8166 0.0002549 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        15.6 1.53 24.3     12.5     18.8
##  EA        13.7 1.53 24.3     10.6     16.9
##  NP        29.4 1.53 24.3     26.2     32.5
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        29.4 1.53 24.3     26.2     32.5  a    
##  AG        15.6 1.53 24.3     12.5     18.8   b   
##  EA        13.7 1.53 24.3     10.6     16.9   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE   df lower.CL upper.CL .group
## 2      NP 29.37031 1.526772 24.3 26.22126 32.51935     a 
## 1      AG 15.62377 1.526772 24.3 12.47472 18.77281      b
## 3      EA 13.72205 1.526772 24.3 10.57301 16.87110      b

##  precip lsmean   SE   df lower.CL upper.CL
##  472      16.0 1.53 13.5     12.7     19.3
##  579      19.8 1.53 13.5     16.6     23.1
##  850      22.9 1.53 13.5     19.6     26.1
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  850      22.9 1.53 13.5     19.6     26.1  a    
##  579      19.8 1.53 13.5     16.6     23.1  ab   
##  472      16.0 1.53 13.5     12.7     19.3   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE       df  lower.CL upper.CL .group
## 1    850 29.87663 1.355610 10.05137 26.858234 32.89503    a  
## 2    579 19.03206 1.355610 10.05137 16.013660 22.05045     b 
## 3    472 12.45756 1.272457 19.16915  9.795863 15.11925      c

## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       13.82 2.64 27     8.39     19.2
##  EA       12.84 2.64 27     7.42     18.3
##  NP       21.41 2.64 27    15.98     26.8
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        7.55 2.64 27     2.12     13.0
##  EA       17.68 2.64 27    12.25     23.1
##  NP       34.29 2.64 27    28.86     39.7
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       25.51 2.64 27    20.08     30.9
##  EA       10.64 2.64 27     5.22     16.1
##  NP       32.42 2.64 27    26.99     37.8
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       21.41 2.64 27    15.98     26.8  a    
##  AG       13.82 2.64 27     8.39     19.2  ab   
##  EA       12.84 2.64 27     7.42     18.3   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       34.29 2.64 27    28.86     39.7  a    
##  EA       17.68 2.64 27    12.25     23.1   b   
##  AG        7.55 2.64 27     2.12     13.0    c  
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       32.42 2.64 27    26.99     37.8  a    
##  AG       25.51 2.64 27    20.08     30.9  a    
##  EA       10.64 2.64 27     5.22     16.1   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE df  lower.CL upper.CL .group
## 2      NP    472 21.40554 2.644448 27 15.979580 26.83150     a 
## 1      AG    472 13.81916 2.644448 27  8.393198 19.24511     ab
## 3      EA    472 12.84488 2.644448 27  7.418923 18.27084      b
## 4      NP    579 34.28563 2.644448 27 28.859670 39.71159    a  
## 5      EA    579 17.68023 2.644448 27 12.254271 23.10619     b 
## 6      AG    579  7.54616 2.644448 27  2.120202 12.97212      c
## 8      NP    850 32.41976 2.644448 27 26.993800 37.84572     a 
## 7      AG    850 25.50598 2.644448 27 20.080025 30.93194     a 
## 9      EA    850 10.64105 2.644448 27  5.215088 16.06700      b
##   precip landuse   lsmean       SE df  lower.CL upper.CL .group  location
## 1    472      NP 21.40554 2.644448 27 15.979580 26.83150     a    Tribune
## 2    472      AG 13.81916 2.644448 27  8.393198 19.24511     ab   Tribune
## 3    472      EA 12.84488 2.644448 27  7.418923 18.27084      b   Tribune
## 4    579      NP 34.28563 2.644448 27 28.859670 39.71159    a        Hays
## 5    579      EA 17.68023 2.644448 27 12.254271 23.10619     b       Hays
## 6    579      AG  7.54616 2.644448 27  2.120202 12.97212      c      Hays
## 7    850      NP 32.41976 2.644448 27 26.993800 37.84572     a  Manhattan
## 8    850      AG 25.50598 2.644448 27 20.080025 30.93194     a  Manhattan
## 9    850      EA 10.64105 2.644448 27  5.215088 16.06700      b Manhattan

2.8 Protein

2.8.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse         9.860   4.930     2    30  1.0907 0.348949   
## precip         55.397  55.397     1    30 12.2561 0.001473 **
## landuse:precip  7.114   3.557     2    30  0.7869 0.464425   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF   DenDF F value    Pr(>F)    
## landuse        217.444 108.722     2 19.1466  72.152 1.215e-09 ***
## precip          44.081  22.041     2  7.0513  14.627  0.003096 ** 
## landuse:precip  90.177  22.544     4 19.1466  14.961 1.058e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean    SE   df lower.CL upper.CL
##  472      7.53 0.421 9.91     6.59     8.47
##  579      5.68 0.421 9.91     4.74     6.61
##  850      4.32 0.421 9.91     3.39     5.26
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean    SE   df lower.CL upper.CL .group
##  472      7.53 0.421 9.91     6.59     8.47  a    
##  579      5.68 0.421 9.91     4.74     6.61   b   
##  850      4.32 0.421 9.91     3.39     5.26    c  
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean        SE       df lower.CL upper.CL .group
## 3    472 7.530212 0.4206309 9.909329 6.591824 8.468600    a  
## 2    579 5.675769 0.4206309 9.909329 4.737382 6.614157     b 
## 1    850 4.324931 0.4206309 9.909329 3.386543 5.263319      c

2.8.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)    
## landuse         9.2553  4.6276     2 22.952  5.9863 0.008087 ** 
## precip          1.9296  1.9296     1 29.178  2.4961 0.124910    
## landuse:precip 16.7774  8.3887     2 22.286 10.8516 0.000512 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF   DenDF F value    Pr(>F)    
## landuse        28.6485 14.3242     2 21.0171 22.3799 6.209e-06 ***
## precip          8.0651  4.0325     2  9.6116  6.3003   0.01785 *  
## landuse:precip 11.6571  2.9143     4 17.1425  4.5532   0.01094 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean    SE    df lower.CL upper.CL
##  AG        2.54 0.392  6.84     1.61     3.47
##  EA        2.93 0.330 10.75     2.20     3.66
##  NP        4.69 0.330 10.75     3.96     5.41
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean    SE    df lower.CL upper.CL .group
##  NP        4.69 0.330 10.75     3.96     5.41  a    
##  EA        2.93 0.330 10.75     2.20     3.66   b   
##  AG        2.54 0.392  6.84     1.61     3.47   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean        SE        df lower.CL upper.CL .group
## 1      NP 4.686371 0.3299719 10.747280 3.958018 5.414723     a 
## 2      EA 2.929352 0.3299719 10.747280 2.200999 3.657704      b
## 3      AG 2.536165 0.3915567  6.835987 1.605758 3.466572      b

## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        4.37 0.51 16.5   3.2935     5.45
##  EA        3.60 0.51 16.5   2.5216     4.68
##  NP        5.08 0.51 16.5   4.0000     6.16
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        2.22 0.51 16.5   1.1406     3.30
##  EA        2.52 0.51 16.5   1.4426     3.60
##  NP        3.76 0.51 16.5   2.6852     4.84
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        1.02 0.51 16.5  -0.0606     2.10
##  EA        2.67 0.51 16.5   1.5889     3.75
##  NP        5.22 0.51 16.5   4.1389     6.30
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        5.08 0.51 16.5   4.0000     6.16  a    
##  AG        4.37 0.51 16.5   3.2935     5.45  ab   
##  EA        3.60 0.51 16.5   2.5216     4.68   b   
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        3.76 0.51 16.5   2.6852     4.84  a    
##  EA        2.52 0.51 16.5   1.4426     3.60   b   
##  AG        2.22 0.51 16.5   1.1406     3.30   b   
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  NP        5.22 0.51 16.5   4.1389     6.30  a    
##  EA        2.67 0.51 16.5   1.5889     3.75   b   
##  AG        1.02 0.51 16.5  -0.0606     2.10    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean        SE       df    lower.CL upper.CL .group
## 2      NP    472 5.078364 0.5098937 16.49031  4.00004423 6.156684     a 
## 1      AG    472 4.371856 0.5098937 16.49031  3.29353629 5.450176     ab
## 3      EA    472 3.599918 0.5098937 16.49031  2.52159787 4.678238      b
## 4      NP    579 3.763522 0.5098937 16.49031  2.68520156 4.841842     a 
## 5      EA    579 2.520902 0.5098937 16.49031  1.44258197 3.599222      b
## 6      AG    579 2.218873 0.5098937 16.49031  1.14055311 3.297193      b
## 7      NP    850 5.217226 0.5098937 16.49031  4.13890632 6.295546    a  
## 8      EA    850 2.667235 0.5098937 16.49031  1.58891498 3.745555     b 
## 9      AG    850 1.017765 0.5098937 16.49031 -0.06055468 2.096085      c
##   precip landuse   lsmean        SE       df    lower.CL upper.CL .group
## 1    472      NP 5.078364 0.5098937 16.49031  4.00004423 6.156684     a 
## 2    472      AG 4.371856 0.5098937 16.49031  3.29353629 5.450176     ab
## 3    472      EA 3.599918 0.5098937 16.49031  2.52159787 4.678238      b
## 4    579      NP 3.763522 0.5098937 16.49031  2.68520156 4.841842     a 
## 5    579      EA 2.520902 0.5098937 16.49031  1.44258197 3.599222      b
## 6    579      AG 2.218873 0.5098937 16.49031  1.14055311 3.297193      b
## 7    850      NP 5.217226 0.5098937 16.49031  4.13890632 6.295546    a  
## 8    850      EA 2.667235 0.5098937 16.49031  1.58891498 3.745555     b 
## 9    850      AG 1.017765 0.5098937 16.49031 -0.06055468 2.096085      c
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        0.51448 0.257239     2 20.736  0.6189 0.5482
## precip         0.02700 0.026997     1 21.355  0.0650 0.8013
## landuse:precip 0.52772 0.263862     2 20.736  0.6349 0.5400

2.9 Permanganate Oxidizable Organic Carbon

2.9.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        21874.6 10937.3     2 21.585  12.794 0.0002165 ***
## precip          9475.2  9475.2     1 26.513  11.084 0.0025662 ** 
## landuse:precip 30249.6 15124.8     2 21.585  17.693 2.825e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         11652  5825.9     2    27  14.275 5.891e-05 ***
## precip          34811 17405.3     2    27  42.648 4.401e-09 ***
## landuse:precip  39551  9887.7     4    27  24.228 1.347e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        1409 5.83 24.3     1397     1421
##  EA        1453 5.83 24.3     1441     1465
##  NP        1433 5.83 24.3     1421     1445
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  EA        1453 5.83 24.3     1441     1465  a    
##  NP        1433 5.83 24.3     1421     1445   b   
##  AG        1409 5.83 24.3     1397     1421    c  
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE   df lower.CL upper.CL .group
## 1      EA 1453.179 5.831805 24.3 1441.150 1465.207    a  
## 3      NP 1432.897 5.831805 24.3 1420.868 1444.925     b 
## 2      AG 1409.156 5.831805 24.3 1397.128 1421.185      c

##  precip lsmean   SE   df lower.CL upper.CL
##  472      1398 5.83 13.5     1385     1410
##  579      1473 5.83 13.5     1460     1485
##  850      1425 5.83 13.5     1412     1437
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean   SE   df lower.CL upper.CL .group
##  579      1473 5.83 13.5     1460     1485  a    
##  850      1425 5.83 13.5     1412     1437   b   
##  472      1398 5.83 13.5     1385     1410    c  
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE   df lower.CL upper.CL .group
## 1    579 1472.798 5.831805 13.5 1460.247 1485.350    a  
## 3    850 1424.867 5.831805 13.5 1412.316 1437.419     b 
## 2    472 1397.566 5.831805 13.5 1385.014 1410.118      c

## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        1372 10.1 27     1351     1393
##  EA        1390 10.1 27     1369     1410
##  NP        1431 10.1 27     1410     1452
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        1488 10.1 27     1467     1509
##  EA        1463 10.1 27     1442     1484
##  NP        1468 10.1 27     1447     1488
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        1368 10.1 27     1347     1388
##  EA        1507 10.1 27     1486     1528
##  NP        1400 10.1 27     1379     1421
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        1431 10.1 27     1410     1452  a    
##  EA        1390 10.1 27     1369     1410   b   
##  AG        1372 10.1 27     1351     1393   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  AG        1488 10.1 27     1467     1509  a    
##  NP        1468 10.1 27     1447     1488  a    
##  EA        1463 10.1 27     1442     1484  a    
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  EA        1507 10.1 27     1486     1528  a    
##  NP        1400 10.1 27     1379     1421   b   
##  AG        1368 10.1 27     1347     1388    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE df lower.CL upper.CL .group
## 1      NP    472 1430.925 10.10098 27 1410.199 1451.650     a 
## 2      EA    472 1389.748 10.10098 27 1369.022 1410.473      b
## 3      AG    472 1372.025 10.10098 27 1351.300 1392.751      b
## 5      AG    579 1487.880 10.10098 27 1467.154 1508.605      a
## 6      NP    579 1467.616 10.10098 27 1446.890 1488.341      a
## 4      EA    579 1462.900 10.10098 27 1442.174 1483.625      a
## 7      EA    850 1506.888 10.10098 27 1486.163 1527.614    a  
## 9      NP    850 1400.150 10.10098 27 1379.425 1420.876     b 
## 8      AG    850 1367.564 10.10098 27 1346.838 1388.289      c
##   precip landuse   lsmean       SE df lower.CL upper.CL .group  location
## 1    472      NP 1430.925 10.10098 27 1410.199 1451.650     a    Tribune
## 2    472      EA 1389.748 10.10098 27 1369.022 1410.473      b   Tribune
## 3    472      AG 1372.025 10.10098 27 1351.300 1392.751      b   Tribune
## 4    579      AG 1487.880 10.10098 27 1467.154 1508.605      a      Hays
## 5    579      NP 1467.616 10.10098 27 1446.890 1488.341      a      Hays
## 6    579      EA 1462.900 10.10098 27 1442.174 1483.625      a      Hays
## 7    850      EA 1506.888 10.10098 27 1486.163 1527.614    a   Manhattan
## 8    850      NP 1400.150 10.10098 27 1379.425 1420.876     b  Manhattan
## 9    850      AG 1367.564 10.10098 27 1346.838 1388.289      c Manhattan

2.9.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         38797 19398.5     2    30 10.6366 0.0003224 ***
## precip           1674  1674.1     1    30  0.9179 0.3456774    
## landuse:precip  54061 27030.3     2    30 14.8213 3.338e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         14172    7086     2    27  32.198 7.094e-08 ***
## precip          40129   20065     2    27  91.171 9.814e-13 ***
## landuse:precip  64376   16094     4    27  73.129 4.350e-14 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        1407 4.28 14.8     1398     1416
##  EA        1453 4.28 24.8     1444     1462
##  NP        1416 4.28 24.8     1408     1425
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  EA        1453 4.28 24.8     1444     1462  a    
##  NP        1416 4.28 24.8     1408     1425   b   
##  AG        1407 4.28 14.8     1398     1416   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean       SE    df lower.CL upper.CL .group
## 1      EA 1453.039 4.282481 24.75 1444.214 1461.863     a 
## 3      NP 1416.333 4.282481 24.75 1407.509 1425.158      b
## 2      AG 1407.099 4.282481 14.85 1397.963 1416.235      b

## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        1375 7.42 27     1359     1390
##  EA        1375 7.42 27     1360     1391
##  NP        1409 7.42 27     1394     1425
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        1491 7.42 27     1476     1506
##  EA        1453 7.42 27     1438     1468
##  NP        1460 7.42 27     1445     1475
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        1356 7.42 27     1341     1371
##  EA        1531 7.42 27     1516     1546
##  NP        1380 7.42 27     1364     1395
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        1409 7.42 27     1394     1425  a    
##  EA        1375 7.42 27     1360     1391   b   
##  AG        1375 7.42 27     1359     1390   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  AG        1491 7.42 27     1476     1506  a    
##  NP        1460 7.42 27     1445     1475   b   
##  EA        1453 7.42 27     1438     1468   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  EA        1531 7.42 27     1516     1546  a    
##  NP        1380 7.42 27     1364     1395   b   
##  AG        1356 7.42 27     1341     1371    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE df lower.CL upper.CL .group
## 1      NP    472 1409.322 7.417474 27 1394.102 1424.541     a 
## 2      EA    472 1375.343 7.417474 27 1360.124 1390.563      b
## 3      AG    472 1374.584 7.417474 27 1359.365 1389.803      b
## 5      AG    579 1490.951 7.417474 27 1475.732 1506.170     a 
## 6      NP    579 1460.074 7.417474 27 1444.854 1475.293      b
## 4      EA    579 1452.924 7.417474 27 1437.704 1468.143      b
## 7      EA    850 1530.849 7.417474 27 1515.629 1546.068    a  
## 9      NP    850 1379.605 7.417474 27 1364.386 1394.824     b 
## 8      AG    850 1355.762 7.417474 27 1340.542 1370.981      c
##   precip landuse   lsmean       SE df lower.CL upper.CL .group  location
## 1    472      NP 1409.322 7.417474 27 1394.102 1424.541     a    Tribune
## 2    472      EA 1375.343 7.417474 27 1360.124 1390.563      b   Tribune
## 3    472      AG 1374.584 7.417474 27 1359.365 1389.803      b   Tribune
## 4    579      AG 1490.951 7.417474 27 1475.732 1506.170     a       Hays
## 5    579      NP 1460.074 7.417474 27 1444.854 1475.293      b      Hays
## 6    579      EA 1452.924 7.417474 27 1437.704 1468.143      b      Hays
## 7    850      EA 1530.849 7.417474 27 1515.629 1546.068    a   Manhattan
## 8    850      NP 1379.605 7.417474 27 1364.386 1394.824     b  Manhattan
## 9    850      AG 1355.762 7.417474 27 1340.542 1370.981      c Manhattan

2.9.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         23205   11602     2 22.612  23.518 2.994e-06 ***
## precip          19448   19448     1 24.069  39.421 1.702e-06 ***
## landuse:precip  36378   18189     2 22.612  36.869 7.633e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         19507  9753.3     2 21.532  68.081 4.901e-10 ***
## precip          56085 28042.4     2 10.929 195.744 2.769e-09 ***
## landuse:precip  44618 11154.5     4 21.532  77.862 1.736e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean    SE    df lower.CL upper.CL
##  AG        1414 3.557 22.76     1407     1422
##  EA        1459 3.557 22.76     1451     1466
##  NP        1406 3.557 22.76     1398     1413
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean    SE    df lower.CL upper.CL .group
##  EA        1459 3.557 22.76     1451     1466  a    
##  AG        1414 3.557 22.76     1407     1422   b   
##  NP        1406 3.557 22.76     1398     1413   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean      SE       df lower.CL upper.CL .group
## 3      EA 1458.766 3.55651 22.76111 1451.404 1466.127     a 
## 1      AG 1414.266 3.55651 22.76111 1406.905 1421.628      b
## 2      NP 1405.642 3.55651 22.76111 1398.281 1413.004      b

##  precip lsmean    SE   df lower.CL upper.CL
##  472      1380 3.636 12.1     1372     1387
##  579      1480 3.636 12.1     1472     1488
##  850      1419 3.636 12.1     1411     1427
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean    SE   df lower.CL upper.CL .group
##  579      1480 3.636 12.1     1472     1488  a    
##  850      1419 3.636 12.1     1411     1427   b   
##  472      1380 3.636 12.1     1372     1387    c  
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean       SE       df lower.CL upper.CL .group
## 1    579 1479.694 3.635528 12.09504 1471.780 1487.608    a  
## 3    850 1419.405 3.635528 12.09504 1411.491 1427.319     b 
## 2    472 1379.575 3.635528 12.09504 1371.661 1387.490      c

## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        1376 6.09 26.9     1363     1388
##  EA        1384 6.09 26.9     1372     1397
##  NP        1379 6.09 26.9     1366     1391
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        1503 6.09 26.9     1490     1515
##  EA        1472 6.09 26.9     1460     1485
##  NP        1464 6.09 26.9     1452     1477
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL
##  AG        1364 6.09 26.9     1352     1377
##  EA        1520 6.09 26.9     1507     1532
##  NP        1374 6.09 26.9     1362     1387
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  EA        1384 6.09 26.9     1372     1397  a    
##  NP        1379 6.09 26.9     1366     1391  a    
##  AG        1376 6.09 26.9     1363     1388  a    
## 
## precip = 579:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  AG        1503 6.09 26.9     1490     1515  a    
##  EA        1472 6.09 26.9     1460     1485   b   
##  NP        1464 6.09 26.9     1452     1477   b   
## 
## precip = 850:
##  landuse lsmean   SE   df lower.CL upper.CL .group
##  EA        1520 6.09 26.9     1507     1532  a    
##  NP        1374 6.09 26.9     1362     1387   b   
##  AG        1364 6.09 26.9     1352     1377   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE      df lower.CL upper.CL .group
## 1      EA    472 1384.425 6.090471 26.8725 1371.925 1396.924      a
## 3      NP    472 1378.523 6.090471 26.8725 1366.023 1391.022      a
## 2      AG    472 1375.779 6.090471 26.8725 1363.280 1388.278      a
## 6      AG    579 1502.843 6.090471 26.8725 1490.343 1515.342     a 
## 5      EA    579 1472.117 6.090471 26.8725 1459.618 1484.616      b
## 4      NP    579 1464.123 6.090471 26.8725 1451.624 1476.623      b
## 7      EA    850 1519.755 6.090471 26.8725 1507.256 1532.255     a 
## 9      NP    850 1374.281 6.090471 26.8725 1361.782 1386.781      b
## 8      AG    850 1364.178 6.090471 26.8725 1351.678 1376.677      b
##   precip landuse   lsmean       SE      df lower.CL upper.CL .group  location
## 1    472      EA 1384.425 6.090471 26.8725 1371.925 1396.924      a   Tribune
## 2    472      NP 1378.523 6.090471 26.8725 1366.023 1391.022      a   Tribune
## 3    472      AG 1375.779 6.090471 26.8725 1363.280 1388.278      a   Tribune
## 4    579      AG 1502.843 6.090471 26.8725 1490.343 1515.342     a       Hays
## 5    579      EA 1472.117 6.090471 26.8725 1459.618 1484.616      b      Hays
## 6    579      NP 1464.123 6.090471 26.8725 1451.624 1476.623      b      Hays
## 7    850      EA 1519.755 6.090471 26.8725 1507.256 1532.255     a  Manhattan
## 8    850      NP 1374.281 6.090471 26.8725 1361.782 1386.781      b Manhattan
## 9    850      AG 1364.178 6.090471 26.8725 1351.678 1376.677      b Manhattan

2.10 Soil Respiration

2.10.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        3.0520  1.5260     2 19.819  0.6252 0.5454
## precip         4.5394  4.5394     1 29.553  1.8598 0.1829
## landuse:precip 5.4906  2.7453     2 19.819  1.1248 0.3446

2.10.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        6.1095 3.05476     2 22.985 16.4145 3.732e-05 ***
## precip         0.0704 0.07043     1 24.556  0.3785    0.5441    
## landuse:precip 6.2650 3.13250     2 22.754 16.8322 3.261e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        11.8739  5.9369     2 24.564  42.896 9.682e-09 ***
## precip         12.2047  6.1023     2 23.306  44.091 1.199e-08 ***
## landuse:precip  7.5323  1.8831     4 24.093  13.606 6.169e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG       0.686 0.112 13.0    0.443    0.929
##  EA       0.830 0.110 23.2    0.602    1.058
##  NP       1.973 0.110 23.2    1.745    2.201
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       1.973 0.110 23.2    1.745    2.201  a    
##  EA       0.830 0.110 23.2    0.602    1.058   b   
##  AG       0.686 0.112 13.0    0.443    0.929   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse    lsmean        SE       df  lower.CL  upper.CL .group
## 1      NP 1.9727342 0.1101649 23.22087 1.7449606 2.2005079     a 
## 2      EA 0.8298641 0.1101649 23.22087 0.6020905 1.0576377      b
## 3      AG 0.6859999 0.1123319 13.00457 0.4433303 0.9286695      b

## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG       0.229 0.189 26.9  -0.1583    0.617
##  EA       0.423 0.189 26.9   0.0357    0.811
##  NP       0.339 0.189 26.9  -0.0482    0.727
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG       0.715 0.189 26.9   0.3271    1.102
##  EA       1.501 0.189 26.9   1.1133    1.889
##  NP       2.908 0.189 26.9   2.5202    3.296
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG       1.114 0.189 26.9   0.7261    1.501
##  EA       0.565 0.189 26.9   0.1775    0.953
##  NP       2.671 0.189 26.9   2.2831    3.058
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  EA       0.423 0.189 26.9   0.0357    0.811  a    
##  NP       0.339 0.189 26.9  -0.0482    0.727  a    
##  AG       0.229 0.189 26.9  -0.1583    0.617  a    
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       2.908 0.189 26.9   2.5202    3.296  a    
##  EA       1.501 0.189 26.9   1.1133    1.889   b   
##  AG       0.715 0.189 26.9   0.3271    1.102    c  
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       2.671 0.189 26.9   2.2831    3.058  a    
##  AG       1.114 0.189 26.9   0.7261    1.501   b   
##  EA       0.565 0.189 26.9   0.1775    0.953    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean        SE       df    lower.CL  upper.CL .group
## 1      EA    472 0.4233481 0.1889066 26.87877  0.03566189 0.8110343      a
## 3      NP    472 0.3394918 0.1889066 26.87877 -0.04819445 0.7271780      a
## 2      AG    472 0.2294218 0.1889066 26.87877 -0.15826445 0.6171080      a
## 4      NP    579 2.9079069 0.1889066 26.87877  2.52022071 3.2955931    a  
## 5      EA    579 1.5010137 0.1889066 26.87877  1.11332753 1.8886999     b 
## 6      AG    579 0.7147885 0.1889066 26.87877  0.32710232 1.1024747      c
## 8      NP    850 2.6708041 0.1889066 26.87877  2.28311788 3.0584903    a  
## 7      AG    850 1.1137894 0.1889066 26.87877  0.72610319 1.5014756     b 
## 9      EA    850 0.5652306 0.1889066 26.87877  0.17754438 0.9529168      c
##   precip landuse    lsmean        SE       df    lower.CL  upper.CL .group
## 1    472      EA 0.4233481 0.1889066 26.87877  0.03566189 0.8110343      a
## 2    472      NP 0.3394918 0.1889066 26.87877 -0.04819445 0.7271780      a
## 3    472      AG 0.2294218 0.1889066 26.87877 -0.15826445 0.6171080      a
## 4    579      NP 2.9079069 0.1889066 26.87877  2.52022071 3.2955931    a  
## 5    579      EA 1.5010137 0.1889066 26.87877  1.11332753 1.8886999     b 
## 6    579      AG 0.7147885 0.1889066 26.87877  0.32710232 1.1024747      c
## 7    850      NP 2.6708041 0.1889066 26.87877  2.28311788 3.0584903    a  
## 8    850      AG 1.1137894 0.1889066 26.87877  0.72610319 1.5014756     b 
## 9    850      EA 0.5652306 0.1889066 26.87877  0.17754438 0.9529168      c
##    location
## 1   Tribune
## 2   Tribune
## 3   Tribune
## 4      Hays
## 5      Hays
## 6      Hays
## 7 Manhattan
## 8 Manhattan
## 9 Manhattan

2.10.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        1.8994  0.9497     2    30  3.0116 0.064279 . 
## precip         3.8505  3.8505     1    30 12.2102 0.001499 **
## landuse:precip 3.6003  1.8002     2    30  5.7085 0.007927 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        9.2111  4.6055     2    27  75.650 8.567e-12 ***
## precip         7.8223  3.9112     2    27  64.244 5.438e-11 ***
## landuse:precip 7.4452  1.8613     4    27  30.574 1.135e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  precip lsmean     SE   df lower.CL upper.CL
##  472     0.233 0.0712 13.5   0.0794    0.386
##  579     1.235 0.0712 13.5   1.0818    1.388
##  850     1.208 0.0712 13.5   1.0542    1.361
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean     SE   df lower.CL upper.CL .group
##  579     1.235 0.0712 13.5   1.0818    1.388  a    
##  850     1.208 0.0712 13.5   1.0542    1.361  a    
##  472     0.233 0.0712 13.5   0.0794    0.386   b   
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip    lsmean         SE   df   lower.CL  upper.CL .group
## 1    579 1.2350554 0.07122692 13.5 1.08175634 1.3883545     a 
## 3    850 1.2075208 0.07122692 13.5 1.05422174 1.3608199     a 
## 2    472 0.2327438 0.07122692 13.5 0.07944469 0.3860429      b

## precip = 472:
##  landuse lsmean    SE df lower.CL upper.CL
##  AG       0.128 0.123 27  -0.1250    0.381
##  EA       0.428 0.123 27   0.1752    0.681
##  NP       0.142 0.123 27  -0.1113    0.395
## 
## precip = 579:
##  landuse lsmean    SE df lower.CL upper.CL
##  AG       0.270 0.123 27   0.0167    0.523
##  EA       1.041 0.123 27   0.7880    1.294
##  NP       2.394 0.123 27   2.1411    2.647
## 
## precip = 850:
##  landuse lsmean    SE df lower.CL upper.CL
##  AG       0.936 0.123 27   0.6826    1.189
##  EA       0.426 0.123 27   0.1730    0.679
##  NP       2.261 0.123 27   2.0075    2.514
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean    SE df lower.CL upper.CL .group
##  EA       0.428 0.123 27   0.1752    0.681  a    
##  NP       0.142 0.123 27  -0.1113    0.395  a    
##  AG       0.128 0.123 27  -0.1250    0.381  a    
## 
## precip = 579:
##  landuse lsmean    SE df lower.CL upper.CL .group
##  NP       2.394 0.123 27   2.1411    2.647  a    
##  EA       1.041 0.123 27   0.7880    1.294   b   
##  AG       0.270 0.123 27   0.0167    0.523    c  
## 
## precip = 850:
##  landuse lsmean    SE df lower.CL upper.CL .group
##  NP       2.261 0.123 27   2.0075    2.514  a    
##  AG       0.936 0.123 27   0.6826    1.189   b   
##  EA       0.426 0.123 27   0.1730    0.679    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean        SE df    lower.CL  upper.CL .group
## 1      EA    472 0.4282920 0.1233687 27  0.17516042 0.6814235      a
## 3      NP    472 0.1417868 0.1233687 27 -0.11134478 0.3949183      a
## 2      AG    472 0.1281525 0.1233687 27 -0.12497902 0.3812841      a
## 4      NP    579 2.3942570 0.1233687 27  2.14112543 2.6473886    a  
## 5      EA    579 1.0410972 0.1233687 27  0.78796563 1.2942288     b 
## 6      AG    579 0.2698121 0.1233687 27  0.01668053 0.5229437      c
## 8      NP    850 2.2606769 0.1233687 27  2.00754531 2.5138084    a  
## 7      AG    850 0.9357511 0.1233687 27  0.68261950 1.1888826     b 
## 9      EA    850 0.4261345 0.1233687 27  0.17300297 0.6792661      c
##   precip landuse    lsmean        SE df    lower.CL  upper.CL .group  location
## 1    472      EA 0.4282920 0.1233687 27  0.17516042 0.6814235      a   Tribune
## 2    472      NP 0.1417868 0.1233687 27 -0.11134478 0.3949183      a   Tribune
## 3    472      AG 0.1281525 0.1233687 27 -0.12497902 0.3812841      a   Tribune
## 4    579      NP 2.3942570 0.1233687 27  2.14112543 2.6473886    a        Hays
## 5    579      EA 1.0410972 0.1233687 27  0.78796563 1.2942288     b       Hays
## 6    579      AG 0.2698121 0.1233687 27  0.01668053 0.5229437      c      Hays
## 7    850      NP 2.2606769 0.1233687 27  2.00754531 2.5138084    a   Manhattan
## 8    850      AG 0.9357511 0.1233687 27  0.68261950 1.1888826     b  Manhattan
## 9    850      EA 0.4261345 0.1233687 27  0.17300297 0.6792661      c Manhattan

3 Chemical

3.1 NH4+

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        1.2012  0.6006     2    29  1.8387 0.1770834    
## precip         4.3932  4.3932     1    29 13.4495 0.0009791 ***
## landuse:precip 0.8080  0.4040     2    29  1.2368 0.3051679    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        0.278002 0.139001     2 21.341  8.0597 0.002469 **
## precip         0.045177 0.045177     1 29.873  2.6195 0.116071   
## landuse:precip 0.191443 0.095722     2 20.564  5.5502 0.011817 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.055792 0.027896     2    30  3.8259 0.03311 *
## precip         0.016342 0.016342     1    30  2.2413 0.14482  
## landuse:precip 0.033857 0.016929     2    30  2.3217 0.11548  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.27472 0.13736     2    29  3.5490   0.04181 *  
## precip         0.84933 0.84933     1    29 21.9440 6.093e-05 ***
## landuse:precip 0.20333 0.10167     2    29  2.6268   0.08944 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.167393 0.083696     2    30 10.4425 0.0003614 ***
## precip         0.059303 0.059303     1    30  7.3990 0.0107502 *  
## landuse:precip 0.121247 0.060623     2    30  7.5637 0.0021888 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq   Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        0.043897 0.0219486     2 24.084  5.9488 0.007951 **
## precip         0.009368 0.0093679     1 22.921  2.5390 0.124761   
## landuse:precip 0.024428 0.0122141     2 24.084  3.3104 0.053681 . 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.2 NO3-

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        259.52  129.76     2    29  2.4548 0.1035391    
## precip         959.70  959.70     1    29 18.1552 0.0001962 ***
## landuse:precip 149.18   74.59     2    29  1.4111 0.2601302    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)    
## landuse        17.358   8.679     2 18.754  9.1407 0.001694 ** 
## precip         36.235  36.235     1 29.846 38.1636 8.68e-07 ***
## landuse:precip  8.841   4.421     2 17.913  4.6559 0.023536 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         4.7669  2.3834     2    30  4.5836   0.01832 *  
## precip         15.9468 15.9468     1    30 30.6675 5.104e-06 ***
## landuse:precip  1.9554  0.9777     2    30  1.8802   0.17010    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.70450 0.35225     2 19.881  2.5104 0.1066342    
## precip         2.37101 2.37101     1 26.947 16.8976 0.0003311 ***
## landuse:precip 0.09917 0.04958     2 19.869  0.3534 0.7066455    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.50448 0.25224     2 21.840  6.3321  0.006778 ** 
## precip         2.35232 2.35232     1 27.859 59.0506 2.346e-08 ***
## landuse:precip 0.09064 0.04532     2 21.183  1.1377  0.339381    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.1601  0.0801     2    30  1.7019    0.1995    
## precip         3.1848  3.1848     1    30 67.7004 3.483e-09 ***
## landuse:precip 0.0036  0.0018     2    30  0.0381    0.9627    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.3 Total Nitrogen

3.3.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq   Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.032787 0.0163936     2    30  3.2473 0.05289 .
## precip         0.003716 0.0037165     1    30  0.7362 0.39769  
## landuse:precip 0.002981 0.0014904     2    30  0.2952 0.74650  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.3.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                   Sum Sq   Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.0057792 0.0028896     2    30  4.0590 0.02753 *
## precip         0.0026554 0.0026554     1    30  3.7300 0.06294 .
## landuse:precip 0.0002458 0.0001229     2    30  0.1726 0.84228  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean      SE   df lower.CL upper.CL
##  AG       0.158 0.00792 15.4    0.141    0.175
##  EA       0.161 0.00773 26.9    0.145    0.177
##  NP       0.260 0.00773 26.9    0.244    0.276
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean      SE   df lower.CL upper.CL .group
##  NP       0.260 0.00773 26.9    0.244    0.276  a    
##  EA       0.161 0.00773 26.9    0.145    0.177   b   
##  AG       0.158 0.00792 15.4    0.141    0.175   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse    lsmean          SE       df  lower.CL  upper.CL .group
## 1      NP 0.2600000 0.007726970 26.85141 0.2441415 0.2758585     a 
## 2      EA 0.1608333 0.007726970 26.85141 0.1449748 0.1766919      b
## 3      AG 0.1583333 0.007921595 15.40516 0.1414874 0.1751792      b

3.3.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                   Sum Sq   Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.0022613 0.0011306     2    30  2.1009 0.13998  
## precip         0.0017053 0.0017053     1    30  3.1687 0.08519 .
## landuse:precip 0.0045665 0.0022832     2    30  4.2427 0.02384 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq   Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.037006 0.0185028     2    27 41.9811 5.171e-09 ***
## precip         0.001839 0.0009194     2    27  2.0861  0.143727    
## landuse:precip 0.008678 0.0021694     4    27  4.9223  0.004117 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG       0.140 0.0105 27    0.118    0.162
##  EA       0.145 0.0105 27    0.123    0.167
##  NP       0.175 0.0105 27    0.153    0.197
## 
## precip = 579:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG       0.122 0.0105 27    0.101    0.144
##  EA       0.142 0.0105 27    0.121    0.164
##  NP       0.222 0.0105 27    0.201    0.244
## 
## precip = 850:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG       0.163 0.0105 27    0.141    0.184
##  EA       0.125 0.0105 27    0.103    0.147
##  NP       0.225 0.0105 27    0.203    0.247
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  NP       0.175 0.0105 27    0.153    0.197  a    
##  EA       0.145 0.0105 27    0.123    0.167  ab   
##  AG       0.140 0.0105 27    0.118    0.162   b   
## 
## precip = 579:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  NP       0.222 0.0105 27    0.201    0.244  a    
##  EA       0.142 0.0105 27    0.121    0.164   b   
##  AG       0.122 0.0105 27    0.101    0.144   b   
## 
## precip = 850:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  NP       0.225 0.0105 27    0.203    0.247  a    
##  AG       0.163 0.0105 27    0.141    0.184   b   
##  EA       0.125 0.0105 27    0.103    0.147    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean         SE df  lower.CL  upper.CL .group
## 1      NP    472 0.1750 0.01049691 27 0.1534621 0.1965379     a 
## 2      EA    472 0.1450 0.01049691 27 0.1234621 0.1665379     ab
## 3      AG    472 0.1400 0.01049691 27 0.1184621 0.1615379      b
## 4      NP    579 0.2225 0.01049691 27 0.2009621 0.2440379     a 
## 5      EA    579 0.1425 0.01049691 27 0.1209621 0.1640379      b
## 6      AG    579 0.1225 0.01049691 27 0.1009621 0.1440379      b
## 8      NP    850 0.2250 0.01049691 27 0.2034621 0.2465379    a  
## 7      AG    850 0.1625 0.01049691 27 0.1409621 0.1840379     b 
## 9      EA    850 0.1250 0.01049691 27 0.1034621 0.1465379      c
##   precip landuse lsmean         SE df  lower.CL  upper.CL .group  location
## 1    472      NP 0.1750 0.01049691 27 0.1534621 0.1965379     a    Tribune
## 2    472      EA 0.1450 0.01049691 27 0.1234621 0.1665379     ab   Tribune
## 3    472      AG 0.1400 0.01049691 27 0.1184621 0.1615379      b   Tribune
## 4    579      NP 0.2225 0.01049691 27 0.2009621 0.2440379     a       Hays
## 5    579      EA 0.1425 0.01049691 27 0.1209621 0.1640379      b      Hays
## 6    579      AG 0.1225 0.01049691 27 0.1009621 0.1440379      b      Hays
## 7    850      NP 0.2250 0.01049691 27 0.2034621 0.2465379    a   Manhattan
## 8    850      AG 0.1625 0.01049691 27 0.1409621 0.1840379     b  Manhattan
## 9    850      EA 0.1250 0.01049691 27 0.1034621 0.1465379      c Manhattan

3.4 Total Organic Carbon

3.4.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        3.6739 1.83695     2 24.041  2.5457 0.09936 .
## precip         0.0401 0.04007     1 24.088  0.0555 0.81569  
## landuse:precip 1.8953 0.94766     2 24.041  1.3133 0.28754  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.4.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.16035 0.08017     2    30  0.7959 0.4604606    
## precip         1.69163 1.69163     1    30 16.7939 0.0002915 ***
## landuse:precip 0.75036 0.37518     2    30  3.7247 0.0359050 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        14.3807  7.1903     2    27 130.536 1.319e-14 ***
## precip          2.4270  1.2135     2    27  22.030 2.120e-06 ***
## landuse:precip  1.5496  0.3874     4    27   7.033 0.0005146 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean    SE df lower.CL upper.CL
##  AG        1.67 0.117 27    1.427     1.91
##  EA        1.42 0.117 27    1.179     1.66
##  NP        2.52 0.117 27    2.279     2.76
## 
## precip = 579:
##  landuse lsmean    SE df lower.CL upper.CL
##  AG        0.99 0.117 27    0.749     1.23
##  EA        1.45 0.117 27    1.212     1.69
##  NP        2.59 0.117 27    2.352     2.83
## 
## precip = 850:
##  landuse lsmean    SE df lower.CL upper.CL
##  AG        1.96 0.117 27    1.722     2.20
##  EA        1.52 0.117 27    1.282     1.76
##  NP        3.41 0.117 27    3.172     3.65
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean    SE df lower.CL upper.CL .group
##  NP        2.52 0.117 27    2.279     2.76  a    
##  AG        1.67 0.117 27    1.427     1.91   b   
##  EA        1.42 0.117 27    1.179     1.66   b   
## 
## precip = 579:
##  landuse lsmean    SE df lower.CL upper.CL .group
##  NP        2.59 0.117 27    2.352     2.83  a    
##  EA        1.45 0.117 27    1.212     1.69   b   
##  AG        0.99 0.117 27    0.749     1.23    c  
## 
## precip = 850:
##  landuse lsmean    SE df lower.CL upper.CL .group
##  NP        3.41 0.117 27    3.172     3.65  a    
##  AG        1.96 0.117 27    1.722     2.20   b   
##  EA        1.52 0.117 27    1.282     1.76    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean        SE df  lower.CL upper.CL .group
## 2      NP    472 2.5200 0.1173492 27 2.2792193 2.760781     a 
## 1      AG    472 1.6675 0.1173492 27 1.4267193 1.908281      b
## 3      EA    472 1.4200 0.1173492 27 1.1792193 1.660781      b
## 4      NP    579 2.5925 0.1173492 27 2.3517193 2.833281    a  
## 5      EA    579 1.4525 0.1173492 27 1.2117193 1.693281     b 
## 6      AG    579 0.9900 0.1173492 27 0.7492193 1.230781      c
## 8      NP    850 3.4125 0.1173492 27 3.1717193 3.653281    a  
## 7      AG    850 1.9625 0.1173492 27 1.7217193 2.203281     b 
## 9      EA    850 1.5225 0.1173492 27 1.2817193 1.763281      c
##   precip landuse lsmean        SE df  lower.CL upper.CL .group  location
## 1    472      NP 2.5200 0.1173492 27 2.2792193 2.760781     a    Tribune
## 2    472      AG 1.6675 0.1173492 27 1.4267193 1.908281      b   Tribune
## 3    472      EA 1.4200 0.1173492 27 1.1792193 1.660781      b   Tribune
## 4    579      NP 2.5925 0.1173492 27 2.3517193 2.833281    a        Hays
## 5    579      EA 1.4525 0.1173492 27 1.2117193 1.693281     b       Hays
## 6    579      AG 0.9900 0.1173492 27 0.7492193 1.230781      c      Hays
## 7    850      NP 3.4125 0.1173492 27 3.1717193 3.653281    a   Manhattan
## 8    850      AG 1.9625 0.1173492 27 1.7217193 2.203281     b  Manhattan
## 9    850      EA 1.5225 0.1173492 27 1.2817193 1.763281      c Manhattan

3.4.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.3615 0.18075     2 23.085  2.4825  0.105551    
## precip         2.7765 2.77650     1 22.620 38.1335 2.872e-06 ***
## landuse:precip 1.2321 0.61602     2 23.085  8.4607  0.001752 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        6.9614  3.4807     2 22.464 64.6753 4.780e-10 ***
## precip         3.0641  1.5321     2 12.489 28.4674 2.228e-05 ***
## landuse:precip 1.6607  0.4152     4 22.464  7.7142 0.0004555 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG       1.325 0.118 26.8    1.082     1.57
##  EA       1.300 0.118 26.8    1.057     1.54
##  NP       1.760 0.118 26.8    1.517     2.00
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG       0.945 0.118 26.8    0.702     1.19
##  EA       1.230 0.118 26.8    0.987     1.47
##  NP       2.112 0.118 26.8    1.869     2.36
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG       1.880 0.118 26.8    1.637     2.12
##  EA       1.387 0.118 26.8    1.144     1.63
##  NP       2.953 0.118 26.8    2.709     3.20
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       1.760 0.118 26.8    1.517     2.00  a    
##  AG       1.325 0.118 26.8    1.082     1.57   b   
##  EA       1.300 0.118 26.8    1.057     1.54   b   
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       2.112 0.118 26.8    1.869     2.36  a    
##  EA       1.230 0.118 26.8    0.987     1.47   b   
##  AG       0.945 0.118 26.8    0.702     1.19   b   
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP       2.953 0.118 26.8    2.709     3.20  a    
##  AG       1.880 0.118 26.8    1.637     2.12   b   
##  EA       1.387 0.118 26.8    1.144     1.63    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean       SE       df  lower.CL upper.CL .group
## 2      NP    472 1.7600 0.118403 26.82631 1.5169835 2.003016     a 
## 1      AG    472 1.3250 0.118403 26.82631 1.0819835 1.568016      b
## 3      EA    472 1.3000 0.118403 26.82631 1.0569835 1.543016      b
## 4      NP    579 2.1125 0.118403 26.82631 1.8694835 2.355516     a 
## 5      EA    579 1.2300 0.118403 26.82631 0.9869835 1.473016      b
## 6      AG    579 0.9450 0.118403 26.82631 0.7019835 1.188016      b
## 8      NP    850 2.9525 0.118403 26.82631 2.7094835 3.195516    a  
## 7      AG    850 1.8800 0.118403 26.82631 1.6369835 2.123016     b 
## 9      EA    850 1.3875 0.118403 26.82631 1.1444835 1.630516      c
##   precip landuse lsmean       SE       df  lower.CL upper.CL .group  location
## 1    472      NP 1.7600 0.118403 26.82631 1.5169835 2.003016     a    Tribune
## 2    472      AG 1.3250 0.118403 26.82631 1.0819835 1.568016      b   Tribune
## 3    472      EA 1.3000 0.118403 26.82631 1.0569835 1.543016      b   Tribune
## 4    579      NP 2.1125 0.118403 26.82631 1.8694835 2.355516     a       Hays
## 5    579      EA 1.2300 0.118403 26.82631 0.9869835 1.473016      b      Hays
## 6    579      AG 0.9450 0.118403 26.82631 0.7019835 1.188016      b      Hays
## 7    850      NP 2.9525 0.118403 26.82631 2.7094835 3.195516    a   Manhattan
## 8    850      AG 1.8800 0.118403 26.82631 1.6369835 2.123016     b  Manhattan
## 9    850      EA 1.3875 0.118403 26.82631 1.1444835 1.630516      c Manhattan

3.5 Extractable Ca Content (mg/kg)

3.5.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        10526795  5263397     2 22.781  60.384 7.837e-10 ***
## precip         29875092 29875092     1 21.561 342.742 1.014e-14 ***
## landuse:precip 15040461  7520230     2 22.781  86.276 2.347e-11 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         7009233  3504617     2    27  49.290 9.728e-10 ***
## precip         31173574 15586787     2    27 219.216 < 2.2e-16 ***
## landuse:precip 15354522  3838631     4    27  53.987 1.710e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean  SE df lower.CL upper.CL
##  AG        1445 133 27     1172     1719
##  EA        1549 133 27     1275     1822
##  NP        2129 133 27     1855     2402
## 
## precip = 579:
##  landuse lsmean  SE df lower.CL upper.CL
##  AG        2633 133 27     2359     2907
##  EA        2224 133 27     1951     2498
##  NP        2886 133 27     2612     3159
## 
## precip = 850:
##  landuse lsmean  SE df lower.CL upper.CL
##  AG        5797 133 27     5524     6071
##  EA        2860 133 27     2586     3133
##  NP        3246 133 27     2972     3519
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean  SE df lower.CL upper.CL .group
##  NP        2129 133 27     1855     2402  a    
##  EA        1549 133 27     1275     1822   b   
##  AG        1445 133 27     1172     1719   b   
## 
## precip = 579:
##  landuse lsmean  SE df lower.CL upper.CL .group
##  NP        2886 133 27     2612     3159  a    
##  AG        2633 133 27     2359     2907  a    
##  EA        2224 133 27     1951     2498   b   
## 
## precip = 850:
##  landuse lsmean  SE df lower.CL upper.CL .group
##  AG        5797 133 27     5524     6071  a    
##  NP        3246 133 27     2972     3519   b   
##  EA        2860 133 27     2586     3133   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE df lower.CL upper.CL .group
## 1      NP    472 2128.800 133.3251 27 1855.239 2402.361     a 
## 2      EA    472 1548.825 133.3251 27 1275.264 1822.386      b
## 3      AG    472 1445.200 133.3251 27 1171.639 1718.761      b
## 5      NP    579 2885.750 133.3251 27 2612.189 3159.311     a 
## 4      AG    579 2632.950 133.3251 27 2359.389 2906.511     a 
## 6      EA    579 2224.500 133.3251 27 1950.939 2498.061      b
## 8      AG    850 5797.475 133.3251 27 5523.914 6071.036     a 
## 9      NP    850 3245.700 133.3251 27 2972.139 3519.261      b
## 7      EA    850 2859.800 133.3251 27 2586.239 3133.361      b
##   precip landuse   lsmean       SE df lower.CL upper.CL .group  location
## 1    472      NP 2128.800 133.3251 27 1855.239 2402.361     a    Tribune
## 2    472      EA 1548.825 133.3251 27 1275.264 1822.386      b   Tribune
## 3    472      AG 1445.200 133.3251 27 1171.639 1718.761      b   Tribune
## 4    579      NP 2885.750 133.3251 27 2612.189 3159.311     a       Hays
## 5    579      AG 2632.950 133.3251 27 2359.389 2906.511     a       Hays
## 6    579      EA 2224.500 133.3251 27 1950.939 2498.061      b      Hays
## 7    850      AG 5797.475 133.3251 27 5523.914 6071.036     a  Manhattan
## 8    850      NP 3245.700 133.3251 27 2972.139 3519.261      b Manhattan
## 9    850      EA 2859.800 133.3251 27 2586.239 3133.361      b Manhattan

3.5.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        10832552  5416276     2 26.715  62.460 8.467e-11 ***
## precip         21485011 21485011     1 29.995 247.765 4.855e-16 ***
## landuse:precip 17045109  8522555     2 26.262  98.282 6.399e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         9036133  4518066     2 24.347  62.697 2.491e-10 ***
## precip         21699428 10849714     2 22.878 150.561 6.797e-14 ***
## landuse:precip 17326979  4331745     4 23.817  60.111 4.237e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         9036133  4518066     2 24.347  62.697 2.491e-10 ***
## precip         21699428 10849714     2 22.878 150.561 6.797e-14 ***
## landuse:precip 17326979  4331745     4 23.817  60.111 4.237e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean  SE   df lower.CL upper.CL
##  AG        1610 137 26.8     1330     1891
##  EA        1700 137 26.8     1420     1980
##  NP        2225 137 26.8     1945     2506
## 
## precip = 579:
##  landuse lsmean  SE   df lower.CL upper.CL
##  AG        2844 137 26.8     2564     3125
##  EA        2453 137 26.8     2172     2733
##  NP        2436 137 26.8     2155     2716
## 
## precip = 850:
##  landuse lsmean  SE   df lower.CL upper.CL
##  AG        5812 137 26.8     5532     6093
##  EA        2550 137 26.8     2270     2831
##  NP        2934 137 26.8     2654     3215
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean  SE   df lower.CL upper.CL .group
##  NP        2225 137 26.8     1945     2506  a    
##  EA        1700 137 26.8     1420     1980   b   
##  AG        1610 137 26.8     1330     1891   b   
## 
## precip = 579:
##  landuse lsmean  SE   df lower.CL upper.CL .group
##  AG        2844 137 26.8     2564     3125  a    
##  EA        2453 137 26.8     2172     2733  ab   
##  NP        2436 137 26.8     2155     2716   b   
## 
## precip = 850:
##  landuse lsmean  SE   df lower.CL upper.CL .group
##  AG        5812 137 26.8     5532     6093  a    
##  NP        2934 137 26.8     2654     3215   b   
##  EA        2550 137 26.8     2270     2831   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE       df lower.CL upper.CL .group
## 1      NP    472 2225.350 136.5988 26.84412 1944.996 2505.704     a 
## 2      EA    472 1700.125 136.5988 26.84412 1419.771 1980.479      b
## 3      AG    472 1610.475 136.5988 26.84412 1330.121 1890.829      b
## 6      AG    579 2844.325 136.5988 26.84412 2563.971 3124.679     a 
## 5      EA    579 2452.750 136.5988 26.84412 2172.396 2733.104     ab
## 4      NP    579 2435.500 136.5988 26.84412 2155.146 2715.854      b
## 8      AG    850 5812.175 136.5988 26.84412 5531.821 6092.529     a 
## 9      NP    850 2934.275 136.5988 26.84412 2653.921 3214.629      b
## 7      EA    850 2550.275 136.5988 26.84412 2269.921 2830.629      b
##   precip landuse   lsmean       SE       df lower.CL upper.CL .group  location
## 1    472      NP 2225.350 136.5988 26.84412 1944.996 2505.704     a    Tribune
## 2    472      EA 1700.125 136.5988 26.84412 1419.771 1980.479      b   Tribune
## 3    472      AG 1610.475 136.5988 26.84412 1330.121 1890.829      b   Tribune
## 4    579      AG 2844.325 136.5988 26.84412 2563.971 3124.679     a       Hays
## 5    579      EA 2452.750 136.5988 26.84412 2172.396 2733.104     ab      Hays
## 6    579      NP 2435.500 136.5988 26.84412 2155.146 2715.854      b      Hays
## 7    850      AG 5812.175 136.5988 26.84412 5531.821 6092.529     a  Manhattan
## 8    850      NP 2934.275 136.5988 26.84412 2653.921 3214.629      b Manhattan
## 9    850      EA 2550.275 136.5988 26.84412 2269.921 2830.629      b Manhattan

3.5.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2916289 1458145     2    30  4.7086 0.0166551 *  
## precip         5287323 5287323     1    30 17.0735 0.0002653 ***
## landuse:precip 5961045 2980522     2    30  9.6245 0.0005900 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        7358980 3679490     2    27 12.6716 0.0001315 ***
## precip         5663302 2831651     2    27  9.7518 0.0006491 ***
## landuse:precip 7035408 1758852     4    27  6.0572 0.0012981 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean  SE df lower.CL upper.CL
##  AG        2596 269 27     2043     3149
##  EA        2116 269 27     1563     2668
##  NP        2381 269 27     1828     2934
## 
## precip = 579:
##  landuse lsmean  SE df lower.CL upper.CL
##  AG        3058 269 27     2505     3611
##  EA        2944 269 27     2391     3497
##  NP        2586 269 27     2033     3139
## 
## precip = 850:
##  landuse lsmean  SE df lower.CL upper.CL
##  AG        4824 269 27     4271     5376
##  EA        2450 269 27     1897     3003
##  NP        2734 269 27     2181     3286
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean  SE df lower.CL upper.CL .group
##  AG        2596 269 27     2043     3149  a    
##  NP        2381 269 27     1828     2934  a    
##  EA        2116 269 27     1563     2668  a    
## 
## precip = 579:
##  landuse lsmean  SE df lower.CL upper.CL .group
##  AG        3058 269 27     2505     3611  a    
##  EA        2944 269 27     2391     3497  a    
##  NP        2586 269 27     2033     3139  a    
## 
## precip = 850:
##  landuse lsmean  SE df lower.CL upper.CL .group
##  AG        4824 269 27     4271     5376  a    
##  NP        2734 269 27     2181     3286   b   
##  EA        2450 269 27     1897     3003   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean      SE df lower.CL upper.CL .group
## 2      AG    472 2595.875 269.431 27 2043.048 3148.702      a
## 3      NP    472 2381.250 269.431 27 1828.423 2934.077      a
## 1      EA    472 2115.500 269.431 27 1562.673 2668.327      a
## 6      AG    579 3058.025 269.431 27 2505.198 3610.852      a
## 5      EA    579 2944.000 269.431 27 2391.173 3496.827      a
## 4      NP    579 2586.000 269.431 27 2033.173 3138.827      a
## 8      AG    850 4823.625 269.431 27 4270.798 5376.452     a 
## 9      NP    850 2733.550 269.431 27 2180.723 3286.377      b
## 7      EA    850 2449.725 269.431 27 1896.898 3002.552      b
##   precip landuse   lsmean      SE df lower.CL upper.CL .group  location
## 1    472      AG 2595.875 269.431 27 2043.048 3148.702      a   Tribune
## 2    472      NP 2381.250 269.431 27 1828.423 2934.077      a   Tribune
## 3    472      EA 2115.500 269.431 27 1562.673 2668.327      a   Tribune
## 4    579      AG 3058.025 269.431 27 2505.198 3610.852      a      Hays
## 5    579      EA 2944.000 269.431 27 2391.173 3496.827      a      Hays
## 6    579      NP 2586.000 269.431 27 2033.173 3138.827      a      Hays
## 7    850      AG 4823.625 269.431 27 4270.798 5376.452     a  Manhattan
## 8    850      NP 2733.550 269.431 27 2180.723 3286.377      b Manhattan
## 9    850      EA 2449.725 269.431 27 1896.898 3002.552      b Manhattan

3.6 Extractable Cu Content (mg/kg)

3.6.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.80480 0.40240     2 21.599  19.063 1.697e-05 ***
## precip         0.26667 0.26667     1 26.774  12.633  0.001433 ** 
## landuse:precip 0.64562 0.32281     2 21.599  15.293 7.289e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.35389 0.17694     2    27   9.800 0.0006312 ***
## precip         0.75056 0.37528     2    27  20.785 3.433e-06 ***
## landuse:precip 0.74444 0.18611     4    27  10.308 3.361e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG        1.30 0.0672 27    1.162     1.44
##  EA        1.48 0.0672 27    1.337     1.61
##  NP        1.05 0.0672 27    0.912     1.19
## 
## precip = 579:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG        1.55 0.0672 27    1.412     1.69
##  EA        1.90 0.0672 27    1.762     2.04
##  NP        1.40 0.0672 27    1.262     1.54
## 
## precip = 850:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG        1.25 0.0672 27    1.112     1.39
##  EA        1.30 0.0672 27    1.162     1.44
##  NP        1.55 0.0672 27    1.412     1.69
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  EA        1.48 0.0672 27    1.337     1.61  a    
##  AG        1.30 0.0672 27    1.162     1.44  a    
##  NP        1.05 0.0672 27    0.912     1.19   b   
## 
## precip = 579:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  EA        1.90 0.0672 27    1.762     2.04  a    
##  AG        1.55 0.0672 27    1.412     1.69   b   
##  NP        1.40 0.0672 27    1.262     1.54   b   
## 
## precip = 850:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  NP        1.55 0.0672 27    1.412     1.69  a    
##  EA        1.30 0.0672 27    1.162     1.44   b   
##  AG        1.25 0.0672 27    1.112     1.39   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean         SE df  lower.CL upper.CL .group
## 3      EA    472  1.475 0.06718548 27 1.3371468 1.612853     a 
## 1      AG    472  1.300 0.06718548 27 1.1621468 1.437853     a 
## 2      NP    472  1.050 0.06718548 27 0.9121468 1.187853      b
## 6      EA    579  1.900 0.06718548 27 1.7621468 2.037853     a 
## 4      AG    579  1.550 0.06718548 27 1.4121468 1.687853      b
## 5      NP    579  1.400 0.06718548 27 1.2621468 1.537853      b
## 7      NP    850  1.550 0.06718548 27 1.4121468 1.687853     a 
## 8      EA    850  1.300 0.06718548 27 1.1621468 1.437853      b
## 9      AG    850  1.250 0.06718548 27 1.1121468 1.387853      b
##   precip landuse lsmean         SE df  lower.CL upper.CL .group  location
## 1    472      EA  1.475 0.06718548 27 1.3371468 1.612853     a    Tribune
## 2    472      AG  1.300 0.06718548 27 1.1621468 1.437853     a    Tribune
## 3    472      NP  1.050 0.06718548 27 0.9121468 1.187853      b   Tribune
## 4    579      EA  1.900 0.06718548 27 1.7621468 2.037853     a       Hays
## 5    579      AG  1.550 0.06718548 27 1.4121468 1.687853      b      Hays
## 6    579      NP  1.400 0.06718548 27 1.2621468 1.537853      b      Hays
## 7    850      NP  1.550 0.06718548 27 1.4121468 1.687853     a  Manhattan
## 8    850      EA  1.300 0.06718548 27 1.1621468 1.437853      b Manhattan
## 9    850      AG  1.250 0.06718548 27 1.1121468 1.387853      b Manhattan

3.6.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        0.35934 0.179670     2 23.054  9.0448 0.001260 **
## precip         0.04352 0.043516     1 25.197  2.1907 0.151244   
## landuse:precip 0.32113 0.160566     2 22.752  8.0831 0.002228 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.37167 0.18583     2    27  7.8093  0.002108 ** 
## precip         0.78167 0.39083     2    27 16.4241 2.154e-05 ***
## landuse:precip 0.93167 0.23292     4    27  9.7879 5.010e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG       1.400 0.0771 27    1.242     1.56
##  EA       1.325 0.0771 27    1.167     1.48
##  NP       0.975 0.0771 27    0.817     1.13
## 
## precip = 579:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG       1.450 0.0771 27    1.292     1.61
##  EA       1.825 0.0771 27    1.667     1.98
##  NP       1.500 0.0771 27    1.342     1.66
## 
## precip = 850:
##  landuse lsmean     SE df lower.CL upper.CL
##  AG       1.150 0.0771 27    0.992     1.31
##  EA       1.550 0.0771 27    1.392     1.71
##  NP       1.650 0.0771 27    1.492     1.81
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  AG       1.400 0.0771 27    1.242     1.56  a    
##  EA       1.325 0.0771 27    1.167     1.48  a    
##  NP       0.975 0.0771 27    0.817     1.13   b   
## 
## precip = 579:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  EA       1.825 0.0771 27    1.667     1.98  a    
##  NP       1.500 0.0771 27    1.342     1.66   b   
##  AG       1.450 0.0771 27    1.292     1.61   b   
## 
## precip = 850:
##  landuse lsmean     SE df lower.CL upper.CL .group
##  NP       1.650 0.0771 27    1.492     1.81  a    
##  EA       1.550 0.0771 27    1.392     1.71  a    
##  AG       1.150 0.0771 27    0.992     1.31   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean         SE df  lower.CL upper.CL .group
## 3      AG    472  1.400 0.07713024 27 1.2417418 1.558258     a 
## 2      EA    472  1.325 0.07713024 27 1.1667418 1.483258     a 
## 1      NP    472  0.975 0.07713024 27 0.8167418 1.133258      b
## 4      EA    579  1.825 0.07713024 27 1.6667418 1.983258     a 
## 6      NP    579  1.500 0.07713024 27 1.3417418 1.658258      b
## 5      AG    579  1.450 0.07713024 27 1.2917418 1.608258      b
## 7      NP    850  1.650 0.07713024 27 1.4917418 1.808258     a 
## 8      EA    850  1.550 0.07713024 27 1.3917418 1.708258     a 
## 9      AG    850  1.150 0.07713024 27 0.9917418 1.308258      b
##   precip landuse lsmean         SE df  lower.CL upper.CL .group  location
## 1    472      AG  1.400 0.07713024 27 1.2417418 1.558258     a    Tribune
## 2    472      EA  1.325 0.07713024 27 1.1667418 1.483258     a    Tribune
## 3    472      NP  0.975 0.07713024 27 0.8167418 1.133258      b   Tribune
## 4    579      EA  1.825 0.07713024 27 1.6667418 1.983258     a       Hays
## 5    579      NP  1.500 0.07713024 27 1.3417418 1.658258      b      Hays
## 6    579      AG  1.450 0.07713024 27 1.2917418 1.608258      b      Hays
## 7    850      NP  1.650 0.07713024 27 1.4917418 1.808258     a  Manhattan
## 8    850      EA  1.550 0.07713024 27 1.3917418 1.708258     a  Manhattan
## 9    850      AG  1.150 0.07713024 27 0.9917418 1.308258      b Manhattan

3.6.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)  
## landuse        0.23902 0.11951     2 19.288  1.2545 0.3075  
## precip         0.40712 0.40712     1 29.123  4.2733 0.0477 *
## landuse:precip 0.22971 0.11486     2 19.288  1.2056 0.3211  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.7 Extractable Mg Content (mg/kg)

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse         73281   36640     2    30  5.2138 0.011393 * 
## precip           5291    5291     1    30  0.7529 0.392456   
## landuse:precip 112623   56312     2    30  8.0130 0.001628 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        110597   55299     2    30  7.1464 0.0028962 ** 
## precip          16017   16017     1    30  2.0700 0.1605757    
## landuse:precip 157633   78816     2    30 10.1856 0.0004208 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        128756   64378     2 18.657  7.7128 0.003621 **
## precip           8519    8519     1 23.553  1.0206 0.322647   
## landuse:precip 154460   77230     2 18.657  9.2525 0.001616 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.8 Extractable Mn Content (mg/kg)

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2677.4  1338.7     2    30  5.4179  0.009799 ** 
## precip         6811.9  6811.9     1    30 27.5693 1.148e-05 ***
## landuse:precip 2857.1  1428.6     2    30  5.7818  0.007519 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         684.33  342.17     2    30  2.3057 0.1170911    
## precip         2422.07 2422.07     1    30 16.3213 0.0003421 ***
## landuse:precip  937.55  468.77     2    30  3.1589 0.0568907 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        110.197  55.099     2    30  0.5042 0.6090
## precip         261.908 261.908     1    30  2.3967 0.1321
## landuse:precip  88.779  44.390     2    30  0.4062 0.6698
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        1.9658  0.9829     2    30  8.2622 0.0013855 ** 
## precip         5.6622  5.6622     1    30 47.5955 1.167e-07 ***
## landuse:precip 2.3822  1.1911     2    30 10.0123 0.0004667 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.73645 0.36823     2    30  2.9824   0.06586 .  
## precip         2.48313 2.48313     1    30 20.1118 9.938e-05 ***
## landuse:precip 1.09313 0.54657     2    30  4.4268   0.02067 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.17184 0.085919     2    30  0.5615 0.5762
## precip         0.23635 0.236351     1    30  1.5446 0.2236
## landuse:precip 0.13117 0.065583     2    30  0.4286 0.6553

3.9 Extractable Na Content (mg/kg)

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse          66.01   33.01     2 23.510  2.5208   0.10189    
## precip         1058.94 1058.94     1 24.823 80.8761 2.779e-09 ***
## landuse:precip   93.64   46.82     2 23.510  3.5758   0.04407 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         767.56  383.78     2    30  5.8723 0.0070437 ** 
## precip         2029.74 2029.74     1    30 31.0575 4.624e-06 ***
## landuse:precip 1238.34  619.17     2    30  9.4741 0.0006468 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         6868.1  3434.0     2 25.879  6.4274 0.0054168 ** 
## precip          7963.3  7963.3     1 24.994 14.9046 0.0007083 ***
## landuse:precip 10154.4  5077.2     2 25.879  9.5029 0.0008044 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.18675 0.09337     2 17.612  1.6027    0.2294    
## precip         2.96662 2.96662     1 27.113 50.9206 1.103e-07 ***
## landuse:precip 0.23069 0.11535     2 17.612  1.9799    0.1676    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.6502  0.3251     2    30  2.4771    0.1010    
## precip         5.5264  5.5264     1    30 42.1119 3.582e-07 ***
## landuse:precip 1.1598  0.5799     2    30  4.4189    0.0208 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        4.9461  2.4730     2 24.479   9.503 0.0008825 ***
## precip         8.7795  8.7795     1 24.103  33.736 5.381e-06 ***
## landuse:precip 6.4365  3.2182     2 24.479  12.367 0.0001941 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.10 Extractable P Content (mg/kg)

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        4022.1 2011.03     2 22.361  3.7651 0.03897 *
## precip           52.3   52.28     1 24.712  0.0979 0.75702  
## landuse:precip 2820.6 1410.30     2 22.361  2.6404 0.09345 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse          7677    3838     2    30  2.6693   0.08572 .  
## precip          33278   33278     1    30 23.1426 3.975e-05 ***
## landuse:precip   6154    3077     2    30  2.1399   0.13528    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        3435.1 1717.57     2 22.127  4.8386 0.018064 * 
## precip           60.7   60.66     1 25.600  0.1709 0.682762   
## landuse:precip 4844.3 2422.15     2 22.127  6.8234 0.004917 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        2.6248 1.31238     2 22.029  4.0362 0.032101 * 
## precip         0.0258 0.02582     1 25.488  0.0794 0.780367   
## landuse:precip 5.1732 2.58660     2 22.029  7.9551 0.002509 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         6.1486  3.0743     2    30  3.2986   0.05071 .  
## precip         23.6475 23.6475     1    30 25.3730 2.098e-05 ***
## landuse:precip  7.6806  3.8403     2    30  4.1205   0.02624 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        5.0455  2.5227     2 22.030  7.3776 0.0035270 ** 
## precip         0.1243  0.1243     1 25.773  0.3634 0.5518766    
## landuse:precip 7.7727  3.8863     2 22.030 11.3654 0.0004061 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.11 Extractable K Content (mg/kg)

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        243718  121859     2 20.875  11.264 0.0004809 ***
## precip         543302  543302     1 28.478  50.218 9.439e-08 ***
## landuse:precip 235695  117847     2 20.875  10.893 0.0005756 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF  F value    Pr(>F)    
## landuse          52814   26407     2 21.421   2.8658   0.07891 .  
## precip         1169188 1169188     1 29.981 126.8859 2.687e-12 ***
## landuse:precip   49780   24890     2 20.649   2.7012   0.09077 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF  F value    Pr(>F)    
## landuse          37264   18632     2 19.520   2.7769   0.08683 .  
## precip         1217054 1217054     1 18.019 181.3872 7.596e-11 ***
## landuse:precip   31742   15871     2 19.520   2.3654   0.12027    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF  F value    Pr(>F)    
## landuse         5.722   2.861     2    30   9.9069 0.0004973 ***
## precip         39.911  39.911     1    30 138.2090 9.302e-13 ***
## landuse:precip  6.430   3.215     2    30  11.1326 0.0002419 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF  F value    Pr(>F)    
## landuse         2.353   1.177     2    30   5.6981  0.007987 ** 
## precip         45.499  45.499     1    30 220.3300 2.316e-15 ***
## landuse:precip  2.808   1.404     2    30   6.7998  0.003669 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF  F value    Pr(>F)    
## landuse         1.270   0.635     2    30   3.5712   0.04062 *  
## precip         43.675  43.675     1    30 245.7125 5.409e-16 ***
## landuse:precip  1.556   0.778     2    30   4.3767   0.02149 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.12 Extractable Zn Content (mg/kg)

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        6.5639  3.2820     2 23.060  1.8274 0.1834
## precip         3.9597  3.9597     1 25.815  2.2048 0.1497
## landuse:precip 2.1397  1.0698     2 23.060  0.5957 0.5594
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        0.73378 0.36689     2 26.774  1.2090 0.3143
## precip         0.66870 0.66870     1 29.985  2.2034 0.1481
## landuse:precip 0.73000 0.36500     2 26.336  1.2027 0.3163
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        0.36864 0.184322     2 20.070  2.6635 0.09420 .
## precip         0.07967 0.079669     1 21.808  1.1513 0.29502  
## landuse:precip 0.40302 0.201509     2 20.070  2.9119 0.07756 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.28547 0.142736     2    30  2.8620 0.07285 .
## precip         0.08426 0.084260     1    30  1.6895 0.20356  
## landuse:precip 0.11955 0.059778     2    30  1.1986 0.31564  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.095955 0.047977     2    30  1.7297 0.1946
## precip         0.060541 0.060541     1    30  2.1827 0.1500
## landuse:precip 0.086563 0.043281     2    30  1.5604 0.2266
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        0.071159 0.035580     2 20.706  2.8258 0.08221 .
## precip         0.018298 0.018298     1 23.084  1.4533 0.24021  
## landuse:precip 0.070124 0.035062     2 20.706  2.7847 0.08491 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.13 Extractable Fe Content (mg/kg)

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2419.8 1209.92     2    30  9.0976 0.0008162 ***
## precip            0.3    0.27     1    30  0.0020 0.9645424    
## landuse:precip 4399.9 2199.96     2    30 16.5418 1.439e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        6504.4  3252.2     2 24.791  19.432 8.386e-06 ***
## precip         2731.7  2731.7     1 29.957  16.322 0.0003427 ***
## landuse:precip 9592.0  4796.0     2 24.225  28.656 4.126e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        3931.7  1965.8     2 23.996  10.479 0.0005357 ***
## precip         8960.2  8960.2     1 23.997  47.764 3.799e-07 ***
## landuse:precip 5957.0  2978.5     2 23.996  15.877 4.049e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2.5850 1.29252     2    30 16.5234 1.452e-05 ***
## precip         0.2069 0.20687     1    30  2.6446    0.1144    
## landuse:precip 4.6526 2.32631     2    30 29.7392 7.604e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        5.3158  2.6579     2 24.394  25.858 9.392e-07 ***
## precip         0.6818  0.6818     1 29.987   6.633   0.01518 *  
## landuse:precip 7.9518  3.9759     2 23.794  38.681 3.330e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        2.8847  1.4424     2 23.478  8.0877 0.0021280 ** 
## precip         5.8035  5.8035     1 23.083 32.5416 8.164e-06 ***
## landuse:precip 4.3114  2.1557     2 23.478 12.0875 0.0002461 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3.14 pH

3.14.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         4.9834  2.4917     2    30  24.884 4.259e-07 ***
## precip         11.7672 11.7672     1    30 117.519 6.758e-12 ***
## landuse:precip  5.8123  2.9061     2    30  29.024 9.685e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         2.4267  1.2133     2    27  13.537 8.473e-05 ***
## precip         11.7950  5.8975     2    27  65.799 4.163e-11 ***
## landuse:precip  6.3683  1.5921     4    27  17.763 2.961e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         2.4267  1.2133     2    27  13.537 8.473e-05 ***
## precip         11.7950  5.8975     2    27  65.799 4.163e-11 ***
## landuse:precip  6.3683  1.5921     4    27  17.763 2.961e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  landuse lsmean     SE   df lower.CL upper.CL
##  AG        6.32 0.0864 24.3     6.14     6.49
##  EA        5.75 0.0864 24.3     5.57     5.93
##  NP        6.28 0.0864 24.3     6.11     6.46
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  landuse lsmean     SE   df lower.CL upper.CL .group
##  AG        6.32 0.0864 24.3     6.14     6.49  a    
##  NP        6.28 0.0864 24.3     6.11     6.46  a    
##  EA        5.75 0.0864 24.3     5.57     5.93   b   
## 
## Results are averaged over the levels of: precip 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse   lsmean         SE   df lower.CL upper.CL .group
## 2      AG 6.316667 0.08642416 24.3 6.138412 6.494921     a 
## 3      NP 6.283333 0.08642416 24.3 6.105079 6.461588     a 
## 1      EA 5.750000 0.08642416 24.3 5.571746 5.928254      b

##  precip lsmean     SE   df lower.CL upper.CL
##  472      5.51 0.0864 13.5     5.32     5.69
##  579      5.96 0.0864 13.5     5.77     6.14
##  850      6.88 0.0864 13.5     6.70     7.07
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
##  precip lsmean     SE   df lower.CL upper.CL .group
##  850      6.88 0.0864 13.5     6.70     7.07  a    
##  579      5.96 0.0864 13.5     5.77     6.14   b   
##  472      5.51 0.0864 13.5     5.32     5.69    c  
## 
## Results are averaged over the levels of: landuse 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   precip   lsmean         SE   df lower.CL upper.CL .group
## 1    850 6.883333 0.08642416 13.5 6.697326 7.069341    a  
## 2    579 5.958333 0.08642416 13.5 5.772326 6.144341     b 
## 3    472 5.508333 0.08642416 13.5 5.322326 5.694341      c

## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        5.30 0.15 27     4.99     5.61
##  EA        5.17 0.15 27     4.87     5.48
##  NP        6.05 0.15 27     5.74     6.36
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        5.75 0.15 27     5.44     6.06
##  EA        5.85 0.15 27     5.54     6.16
##  NP        6.28 0.15 27     5.97     6.58
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        7.90 0.15 27     7.59     8.21
##  EA        6.22 0.15 27     5.92     6.53
##  NP        6.53 0.15 27     6.22     6.83
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        6.05 0.15 27     5.74     6.36  a    
##  AG        5.30 0.15 27     4.99     5.61   b   
##  EA        5.17 0.15 27     4.87     5.48   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        6.28 0.15 27     5.97     6.58  a    
##  EA        5.85 0.15 27     5.54     6.16  ab   
##  AG        5.75 0.15 27     5.44     6.06   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  AG        7.90 0.15 27     7.59     8.21  a    
##  NP        6.53 0.15 27     6.22     6.83   b   
##  EA        6.22 0.15 27     5.92     6.53   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean       SE df lower.CL upper.CL .group
## 2      NP    472  6.050 0.149691 27 5.742859 6.357141     a 
## 1      AG    472  5.300 0.149691 27 4.992859 5.607141      b
## 3      EA    472  5.175 0.149691 27 4.867859 5.482141      b
## 4      NP    579  6.275 0.149691 27 5.967859 6.582141     a 
## 5      EA    579  5.850 0.149691 27 5.542859 6.157141     ab
## 6      AG    579  5.750 0.149691 27 5.442859 6.057141      b
## 8      AG    850  7.900 0.149691 27 7.592859 8.207141     a 
## 9      NP    850  6.525 0.149691 27 6.217859 6.832141      b
## 7      EA    850  6.225 0.149691 27 5.917859 6.532141      b
##   precip landuse lsmean       SE df lower.CL upper.CL .group  location
## 1    472      NP  6.050 0.149691 27 5.742859 6.357141     a    Tribune
## 2    472      AG  5.300 0.149691 27 4.992859 5.607141      b   Tribune
## 3    472      EA  5.175 0.149691 27 4.867859 5.482141      b   Tribune
## 4    579      NP  6.275 0.149691 27 5.967859 6.582141     a       Hays
## 5    579      EA  5.850 0.149691 27 5.542859 6.157141     ab      Hays
## 6    579      AG  5.750 0.149691 27 5.442859 6.057141      b      Hays
## 7    850      AG  7.900 0.149691 27 7.592859 8.207141     a  Manhattan
## 8    850      NP  6.525 0.149691 27 6.217859 6.832141      b Manhattan
## 9    850      EA  6.225 0.149691 27 5.917859 6.532141      b Manhattan

3.14.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        7.1871  3.5936     2 21.524  24.336 2.985e-06 ***
## precip         2.0640  2.0640     1 29.999  13.977 0.0007789 ***
## landuse:precip 9.1369  4.5684     2 20.759  30.938 5.926e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        1.2084 0.60418     2 20.755  4.7937  0.019430 *  
## precip         2.5870 1.29348     2 10.673 10.2627  0.003266 ** 
## landuse:precip 9.8712 2.46781     4 17.537 19.5801 2.778e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        5.60 0.203 21.5     5.18     6.02
##  EA        5.83 0.203 21.5     5.40     6.25
##  NP        6.25 0.203 21.5     5.83     6.67
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        5.92 0.203 21.5     5.50     6.35
##  EA        6.45 0.203 21.5     6.03     6.87
##  NP        6.35 0.203 21.5     5.93     6.77
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        7.92 0.203 21.5     7.50     8.35
##  EA        5.78 0.203 21.5     5.35     6.20
##  NP        6.25 0.203 21.5     5.83     6.67
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP        6.25 0.203 21.5     5.83     6.67  a    
##  EA        5.83 0.203 21.5     5.40     6.25  ab   
##  AG        5.60 0.203 21.5     5.18     6.02   b   
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  EA        6.45 0.203 21.5     6.03     6.87  a    
##  NP        6.35 0.203 21.5     5.93     6.77  ab   
##  AG        5.92 0.203 21.5     5.50     6.35   b   
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  AG        7.92 0.203 21.5     7.50     8.35  a    
##  NP        6.25 0.203 21.5     5.83     6.67   b   
##  EA        5.78 0.203 21.5     5.35     6.20   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean        SE       df lower.CL upper.CL .group
## 1      NP    472  6.250 0.2030546 21.51346 5.828338 6.671662     a 
## 2      EA    472  5.825 0.2030546 21.51346 5.403338 6.246662     ab
## 3      AG    472  5.600 0.2030546 21.51346 5.178338 6.021662      b
## 4      EA    579  6.450 0.2030546 21.51346 6.028338 6.871662     a 
## 6      NP    579  6.350 0.2030546 21.51346 5.928338 6.771662     ab
## 5      AG    579  5.925 0.2030546 21.51346 5.503338 6.346662      b
## 8      AG    850  7.925 0.2030546 21.51346 7.503338 8.346662     a 
## 9      NP    850  6.250 0.2030546 21.51346 5.828338 6.671662      b
## 7      EA    850  5.775 0.2030546 21.51346 5.353338 6.196662      b
##   precip landuse lsmean        SE       df lower.CL upper.CL .group  location
## 1    472      NP  6.250 0.2030546 21.51346 5.828338 6.671662     a    Tribune
## 2    472      EA  5.825 0.2030546 21.51346 5.403338 6.246662     ab   Tribune
## 3    472      AG  5.600 0.2030546 21.51346 5.178338 6.021662      b   Tribune
## 4    579      EA  6.450 0.2030546 21.51346 6.028338 6.871662     a       Hays
## 5    579      NP  6.350 0.2030546 21.51346 5.928338 6.771662     ab      Hays
## 6    579      AG  5.925 0.2030546 21.51346 5.503338 6.346662      b      Hays
## 7    850      AG  7.925 0.2030546 21.51346 7.503338 8.346662     a  Manhattan
## 8    850      NP  6.250 0.2030546 21.51346 5.828338 6.671662      b Manhattan
## 9    850      EA  5.775 0.2030546 21.51346 5.353338 6.196662      b Manhattan

3.14.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)    
## landuse        4.1797 2.08987     2    30   7.972 0.001673 ** 
## precip         0.4226 0.42259     1    30   1.612 0.213972    
## landuse:precip 5.2378 2.61891     2    30   9.990 0.000473 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF   DenDF F value    Pr(>F)    
## landuse        0.6350 0.31750     2 17.6191  1.7504 0.2025847    
## precip         0.3108 0.15539     2  5.3252  0.8567 0.4758018    
## landuse:precip 7.2183 1.80458     4 17.6191  9.9490 0.0002148 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        6.60 0.235 24.1     6.12     7.08
##  EA        6.60 0.235 24.1     6.12     7.08
##  NP        7.00 0.235 24.1     6.52     7.48
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        6.15 0.235 24.1     5.67     6.63
##  EA        6.90 0.235 24.1     6.42     7.38
##  NP        6.53 0.235 24.1     6.04     7.01
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL
##  AG        7.42 0.235 24.1     6.94     7.91
##  EA        5.70 0.235 24.1     5.22     6.18
##  NP        6.20 0.235 24.1     5.72     6.68
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  NP        7.00 0.235 24.1     6.52     7.48  a    
##  EA        6.60 0.235 24.1     6.12     7.08  a    
##  AG        6.60 0.235 24.1     6.12     7.08  a    
## 
## precip = 579:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  EA        6.90 0.235 24.1     6.42     7.38  a    
##  NP        6.53 0.235 24.1     6.04     7.01  ab   
##  AG        6.15 0.235 24.1     5.67     6.63   b   
## 
## precip = 850:
##  landuse lsmean    SE   df lower.CL upper.CL .group
##  AG        7.42 0.235 24.1     6.94     7.91  a    
##  NP        6.20 0.235 24.1     5.72     6.68   b   
##  EA        5.70 0.235 24.1     5.22     6.18   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip lsmean        SE       df lower.CL upper.CL .group
## 1      NP    472  7.000 0.2346744 24.07352 6.515734 7.484266      a
## 2      EA    472  6.600 0.2346744 24.07352 6.115734 7.084266      a
## 3      AG    472  6.600 0.2346744 24.07352 6.115734 7.084266      a
## 4      EA    579  6.900 0.2346744 24.07352 6.415734 7.384266     a 
## 6      NP    579  6.525 0.2346744 24.07352 6.040734 7.009266     ab
## 5      AG    579  6.150 0.2346744 24.07352 5.665734 6.634266      b
## 8      AG    850  7.425 0.2346744 24.07352 6.940734 7.909266     a 
## 9      NP    850  6.200 0.2346744 24.07352 5.715734 6.684266      b
## 7      EA    850  5.700 0.2346744 24.07352 5.215734 6.184266      b
##   precip landuse lsmean        SE       df lower.CL upper.CL .group  location
## 1    472      NP  7.000 0.2346744 24.07352 6.515734 7.484266      a   Tribune
## 2    472      EA  6.600 0.2346744 24.07352 6.115734 7.084266      a   Tribune
## 3    472      AG  6.600 0.2346744 24.07352 6.115734 7.084266      a   Tribune
## 4    579      EA  6.900 0.2346744 24.07352 6.415734 7.384266     a       Hays
## 5    579      NP  6.525 0.2346744 24.07352 6.040734 7.009266     ab      Hays
## 6    579      AG  6.150 0.2346744 24.07352 5.665734 6.634266      b      Hays
## 7    850      AG  7.425 0.2346744 24.07352 6.940734 7.909266     a  Manhattan
## 8    850      NP  6.200 0.2346744 24.07352 5.715734 6.684266      b Manhattan
## 9    850      EA  5.700 0.2346744 24.07352 5.215734 6.184266      b Manhattan

3.15 Water Moisture

3.15.1 0-5 cm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         710.43  355.22     2 19.605  13.056 0.0002487 ***
## precip         1737.68 1737.68     1 29.609  63.866 7.041e-09 ***
## landuse:precip 1125.78  562.89     2 19.605  20.688 1.476e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        1385.5  692.73     2    27  44.306 2.971e-09 ***
## precip         2044.2 1022.11     2    27  65.373 4.477e-11 ***
## landuse:precip 1373.0  343.26     4    27  21.955 3.692e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        7.99 1.98 27     3.94     12.1
##  EA       18.48 1.98 27    14.42     22.5
##  NP       23.09 1.98 27    19.04     27.2
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        9.90 1.98 27     5.85     14.0
##  EA       11.51 1.98 27     7.46     15.6
##  NP       13.40 1.98 27     9.35     17.5
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG       26.85 1.98 27    22.79     30.9
##  EA       14.15 1.98 27    10.09     18.2
##  NP       47.42 1.98 27    43.36     51.5
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       23.09 1.98 27    19.04     27.2  a    
##  EA       18.48 1.98 27    14.42     22.5  a    
##  AG        7.99 1.98 27     3.94     12.1   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       13.40 1.98 27     9.35     17.5  a    
##  EA       11.51 1.98 27     7.46     15.6  a    
##  AG        9.90 1.98 27     5.85     14.0  a    
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP       47.42 1.98 27    43.36     51.5  a    
##  AG       26.85 1.98 27    22.79     30.9   b   
##  EA       14.15 1.98 27    10.09     18.2    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip    lsmean       SE df  lower.CL upper.CL .group
## 1      NP    472 23.094419 1.977066 27 19.037815 27.15102     a 
## 2      EA    472 18.477766 1.977066 27 14.421162 22.53437     a 
## 3      AG    472  7.993733 1.977066 27  3.937129 12.05034      b
## 4      NP    579 13.403599 1.977066 27  9.346995 17.46020      a
## 5      EA    579 11.511773 1.977066 27  7.455169 15.56838      a
## 6      AG    579  9.902225 1.977066 27  5.845620 13.95883      a
## 8      NP    850 47.419133 1.977066 27 43.362528 51.47574    a  
## 7      AG    850 26.845140 1.977066 27 22.788536 30.90174     b 
## 9      EA    850 14.150962 1.977066 27 10.094358 18.20757      c
##   precip landuse    lsmean       SE df  lower.CL upper.CL .group  location
## 1    472      NP 23.094419 1.977066 27 19.037815 27.15102     a    Tribune
## 2    472      EA 18.477766 1.977066 27 14.421162 22.53437     a    Tribune
## 3    472      AG  7.993733 1.977066 27  3.937129 12.05034      b   Tribune
## 4    579      NP 13.403599 1.977066 27  9.346995 17.46020      a      Hays
## 5    579      EA 11.511773 1.977066 27  7.455169 15.56838      a      Hays
## 6    579      AG  9.902225 1.977066 27  5.845620 13.95883      a      Hays
## 7    850      NP 47.419133 1.977066 27 43.362528 51.47574    a   Manhattan
## 8    850      AG 26.845140 1.977066 27 22.788536 30.90174     b  Manhattan
## 9    850      EA 14.150962 1.977066 27 10.094358 18.20757      c Manhattan

3.15.2 5-10 cm

## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        514.37  257.19     2    30  12.184 0.0001339 ***
## precip         392.79  392.79     1    30  18.608 0.0001602 ***
## landuse:precip 662.01  331.01     2    30  15.681 2.179e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        106.94   53.47     2    27  6.9673  0.003632 ** 
## precip         677.58  338.79     2    27 44.1462 3.084e-09 ***
## landuse:precip 803.26  200.82     4    27 26.1674 6.022e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        19.8 1.39 27     17.0     22.7
##  EA        22.7 1.39 27     19.8     25.5
##  NP        20.3 1.39 27     17.5     23.2
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        16.4 1.39 27     13.6     19.3
##  EA        19.8 1.39 27     17.0     22.6
##  NP        13.4 1.39 27     10.5     16.2
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        24.5 1.39 27     21.7     27.3
##  EA        18.6 1.39 27     15.8     21.5
##  NP        38.2 1.39 27     35.4     41.1
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  EA        22.7 1.39 27     19.8     25.5  a    
##  NP        20.3 1.39 27     17.5     23.2  a    
##  AG        19.8 1.39 27     17.0     22.7  a    
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  EA        19.8 1.39 27     17.0     22.6  a    
##  AG        16.4 1.39 27     13.6     19.3  ab   
##  NP        13.4 1.39 27     10.5     16.2   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        38.2 1.39 27     35.4     41.1  a    
##  AG        24.5 1.39 27     21.7     27.3   b   
##  EA        18.6 1.39 27     15.8     21.5    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE df lower.CL upper.CL .group
## 1      EA    472 22.68240 1.385125 27 19.84036 25.52444      a
## 3      NP    472 20.32969 1.385125 27 17.48765 23.17174      a
## 2      AG    472 19.80959 1.385125 27 16.96755 22.65163      a
## 6      EA    579 19.80418 1.385125 27 16.96214 22.64623     a 
## 4      AG    579 16.44994 1.385125 27 13.60790 19.29198     ab
## 5      NP    579 13.35845 1.385125 27 10.51641 16.20050      b
## 8      NP    850 38.21451 1.385125 27 35.37247 41.05655    a  
## 7      AG    850 24.50495 1.385125 27 21.66290 27.34699     b 
## 9      EA    850 18.62572 1.385125 27 15.78368 21.46776      c
##   precip landuse   lsmean       SE df lower.CL upper.CL .group  location
## 1    472      EA 22.68240 1.385125 27 19.84036 25.52444      a   Tribune
## 2    472      NP 20.32969 1.385125 27 17.48765 23.17174      a   Tribune
## 3    472      AG 19.80959 1.385125 27 16.96755 22.65163      a   Tribune
## 4    579      EA 19.80418 1.385125 27 16.96214 22.64623     a       Hays
## 5    579      AG 16.44994 1.385125 27 13.60790 19.29198     ab      Hays
## 6    579      NP 13.35845 1.385125 27 10.51641 16.20050      b      Hays
## 7    850      NP 38.21451 1.385125 27 35.37247 41.05655    a   Manhattan
## 8    850      AG 24.50495 1.385125 27 21.66290 27.34699     b  Manhattan
## 9    850      EA 18.62572 1.385125 27 15.78368 21.46776      c Manhattan

3.15.3 10-15 cm

## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         11.20   5.598     2    27  0.9336    0.4055    
## precip         301.46 150.728     2    27 25.1375 6.837e-07 ***
## landuse:precip 730.36 182.591     4    27 30.4514 1.186e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse         11.20   5.598     2    27  0.9336    0.4055    
## precip         301.46 150.728     2    27 25.1375 6.837e-07 ***
## landuse:precip 730.36 182.591     4    27 30.4514 1.186e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        24.2 1.22 27     21.6     26.7
##  EA        24.1 1.22 27     21.6     26.6
##  NP        17.8 1.22 27     15.3     20.3
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        21.1 1.22 27     18.6     23.6
##  EA        22.6 1.22 27     20.1     25.1
##  NP        13.7 1.22 27     11.2     16.2
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL
##  AG        24.4 1.22 27     21.9     27.0
##  EA        19.7 1.22 27     17.1     22.2
##  NP        34.4 1.22 27     31.9     36.9
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95
## precip = 472:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  AG        24.2 1.22 27     21.6     26.7  a    
##  EA        24.1 1.22 27     21.6     26.6  a    
##  NP        17.8 1.22 27     15.3     20.3   b   
## 
## precip = 579:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  EA        22.6 1.22 27     20.1     25.1  a    
##  AG        21.1 1.22 27     18.6     23.6  a    
##  NP        13.7 1.22 27     11.2     16.2   b   
## 
## precip = 850:
##  landuse lsmean   SE df lower.CL upper.CL .group
##  NP        34.4 1.22 27     31.9     36.9  a    
##  AG        24.4 1.22 27     21.9     27.0   b   
##  EA        19.7 1.22 27     17.1     22.2    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   landuse precip   lsmean       SE df lower.CL upper.CL .group
## 3      AG    472 24.15175 1.224351 27 21.63959 26.66391     a 
## 2      EA    472 24.13344 1.224351 27 21.62127 26.64560     a 
## 1      NP    472 17.83541 1.224351 27 15.32325 20.34758      b
## 6      EA    579 22.58671 1.224351 27 20.07455 25.09887     a 
## 4      AG    579 21.08764 1.224351 27 18.57548 23.59981     a 
## 5      NP    579 13.68315 1.224351 27 11.17098 16.19531      b
## 8      NP    850 34.41712 1.224351 27 31.90496 36.92928    a  
## 7      AG    850 24.44505 1.224351 27 21.93288 26.95721     b 
## 9      EA    850 19.65619 1.224351 27 17.14402 22.16835      c
##   precip landuse   lsmean       SE df lower.CL upper.CL .group  location
## 1    472      AG 24.15175 1.224351 27 21.63959 26.66391     a    Tribune
## 2    472      EA 24.13344 1.224351 27 21.62127 26.64560     a    Tribune
## 3    472      NP 17.83541 1.224351 27 15.32325 20.34758      b   Tribune
## 4    579      EA 22.58671 1.224351 27 20.07455 25.09887     a       Hays
## 5    579      AG 21.08764 1.224351 27 18.57548 23.59981     a       Hays
## 6    579      NP 13.68315 1.224351 27 11.17098 16.19531      b      Hays
## 7    850      NP 34.41712 1.224351 27 31.90496 36.92928    a   Manhattan
## 8    850      AG 24.44505 1.224351 27 21.93288 26.95721     b  Manhattan
## 9    850      EA 19.65619 1.224351 27 17.14402 22.16835      c Manhattan

4 Physical

4.1 Clay

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        12.9964  6.4982     2 22.847  57.104 1.294e-09 ***
## precip          6.5823  6.5823     1 23.428  57.843 8.948e-08 ***
## landuse:precip 25.8357 12.9179     2 22.847 113.517 1.356e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        12.996   6.498     2    30  5.4125 0.0098387 ** 
## precip         90.146  90.146     1    30 75.0844 1.151e-09 ***
## landuse:precip 25.836  12.918     2    30 10.7595 0.0003001 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        149.487  74.743     2    30  12.582 0.0001076 ***
## precip          73.753  73.753     1    30  12.416 0.0013871 ** 
## landuse:precip 172.037  86.019     2    30  14.480 3.966e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "natural log"
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq   Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.017841 0.0089206     2 22.868  49.422 4.985e-09 ***
## precip         0.008716 0.0087161     1 23.368  48.289 4.026e-07 ***
## landuse:precip 0.036219 0.0181094     2 22.868 100.329 4.776e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.017841 0.008921     2    30  4.0815  0.027050 *  
## precip         0.147685 0.147685     1    30 67.5722 3.553e-09 ***
## landuse:precip 0.036219 0.018109     2    30  8.2858  0.001365 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.181244 0.090622     2    30  12.610  0.000106 ***
## precip         0.086866 0.086866     1    30  12.087  0.001571 ** 
## landuse:precip 0.209133 0.104567     2    30  14.550 3.828e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.2 20 minute MWD

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.63493 0.31746     2    30  2.7254 0.08174 .
## precip         0.00050 0.00050     1    30  0.0043 0.94830  
## landuse:precip 0.19307 0.09653     2    30  0.8287 0.44636  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value  Pr(>F)   
## landuse        1.06993 0.53496     2    30  8.5670 0.00114 **
## precip         0.01212 0.01212     1    30  0.1940 0.66274   
## landuse:precip 0.47723 0.23861     2    30  3.8212 0.03324 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        0.65137 0.32569     2    30  6.8321 0.003589 **
## precip         0.00811 0.00811     1    30  0.1702 0.682870   
## landuse:precip 0.12831 0.06416     2    30  1.3458 0.275587   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq   Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.060137 0.0300684     2    30  2.9226 0.06924 .
## precip         0.000008 0.0000083     1    30  0.0008 0.97749  
## landuse:precip 0.031352 0.0156760     2    30  1.5237 0.23429  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        0.159537 0.079768     2    30  8.2602 0.001387 **
## precip         0.005053 0.005053     1    30  0.5233 0.475051   
## landuse:precip 0.094247 0.047123     2    30  4.8797 0.014629 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.106866 0.053433     2 21.66 11.6764 0.0003627 ***
## precip         0.005291 0.005291     1 27.29  1.1563 0.2916543    
## landuse:precip 0.029561 0.014781     2 21.66  3.2299 0.0592062 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.3 20 minute 2 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        257.872 128.936     2    30  2.5984 0.09105 .
## precip          16.078  16.078     1    30  0.3240 0.57344  
## landuse:precip  67.492  33.746     2    30  0.6801 0.51422  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        355.82 177.909     2    30  7.7613 0.001921 **
## precip           0.09   0.085     1    30  0.0037 0.951753   
## landuse:precip 103.18  51.589     2    30  2.2506 0.122835   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        227.845 113.922     2    30  7.6604 0.002053 **
## precip           3.956   3.956     1    30  0.2660 0.609815   
## landuse:precip  28.670  14.335     2    30  0.9639 0.392888   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.64726 0.32363     2    30  2.1407 0.1352
## precip         0.10429 0.10429     1    30  0.6898 0.4128
## landuse:precip 0.48888 0.24444     2    30  1.6169 0.2153
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        2.1452 1.07261     2    30  8.0973 0.001542 **
## precip         0.0638 0.06380     1    30  0.4816 0.493038   
## landuse:precip 1.0203 0.51013     2    30  3.8510 0.032461 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)   
## landuse        1.28476 0.64238     2 19.249 10.0436 0.00103 **
## precip         0.01512 0.01512     1 27.831  0.2364 0.63065   
## landuse:precip 0.40826 0.20413     2 19.249  3.1916 0.06352 . 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.4 20 minute 0.25 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse          3.94    1.97     2 23.092  0.0803 0.923113   
## precip         341.74  341.74     1 24.299 13.9426 0.001013 **
## landuse:precip  62.02   31.01     2 23.092  1.2652 0.301037   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        360.27 180.134     2 27.856  1.4587 0.2498
## precip          69.43  69.427     1 30.000  0.5622 0.4592
## landuse:precip 505.52 252.758     2 27.535  2.0468 0.1483
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        221.19 110.594     2 25.052  1.1477 0.3335
## precip         241.96 241.956     1 27.240  2.5109 0.1246
## landuse:precip 103.89  51.946     2 25.052  0.5391 0.5899
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        0.05170 0.02585     2 22.977  0.4960 0.615349   
## precip         0.51807 0.51807     1 24.917  9.9398 0.004182 **
## landuse:precip 0.05467 0.02733     2 22.977  0.5244 0.598801   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.60901 0.30451     2    30  1.8762 0.1707
## precip         0.15543 0.15543     1    30  0.9576 0.3356
## landuse:precip 0.79190 0.39595     2    30  2.4396 0.1043
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        0.38261 0.19130     2 24.375  2.0485 0.15054  
## precip         0.42734 0.42734     1 28.702  4.5760 0.04106 *
## landuse:precip 0.17359 0.08680     2 24.375  0.9294 0.40832  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.5 20 minute 0.053 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)  
## landuse        367.52 183.762     2 23.449  2.8888 0.0756 .
## precip           5.75   5.748     1 24.917  0.0904 0.7662  
## landuse:precip  77.14  38.571     2 23.449  0.6064 0.5537  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        930.68  465.34     2 27.303  4.8421 0.01586 *
## precip           2.62    2.62     1 29.962  0.0273 0.86988  
## landuse:precip 845.25  422.62     2 26.965  4.3976 0.02224 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        316.069 158.034     2 23.931  1.3110 0.2882
## precip           2.318   2.318     1 25.017  0.0192 0.8908
## landuse:precip 285.469 142.735     2 23.931  1.1841 0.3233
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        0.309159 0.154580     2 24.294  2.1950 0.1330
## precip         0.014643 0.014643     1 23.593  0.2079 0.6526
## landuse:precip 0.104474 0.052237     2 24.294  0.7418 0.4867
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        0.56122 0.280610     2 27.219  3.1677 0.05797 .
## precip         0.00188 0.001883     1 29.966  0.0213 0.88507  
## landuse:precip 0.49259 0.246295     2 26.875  2.7803 0.07989 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.230347 0.115173     2    30  0.7312 0.4897
## precip         0.004269 0.004269     1    30  0.0271 0.8703
## landuse:precip 0.284333 0.142167     2    30  0.9026 0.4162

4.6 20 minute 0.02 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        2.6836  1.3418     2    30  0.5539 0.5805
## precip         6.5803  6.5803     1    30  2.7164 0.1098
## landuse:precip 1.9851  0.9926     2    30  0.4097 0.6675
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        14.0129  7.0064     2    30  1.9852 0.1550
## precip          0.2395  0.2395     1    30  0.0679 0.7963
## landuse:precip  2.6187  1.3094     2    30  0.3710 0.6932
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        19.0033  9.5016     2 19.059  1.6749 0.2137
## precip          5.9732  5.9732     1 18.225  1.0529 0.3183
## landuse:precip  3.6821  1.8411     2 19.059  0.3245 0.7268
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.106686 0.053343     2    30  1.2644 0.29702  
## precip         0.133321 0.133321     1    30  3.1602 0.08559 .
## landuse:precip 0.017007 0.008503     2    30  0.2016 0.81855  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse        0.36605 0.183023     2    30  3.4747 0.04392 *
## precip         0.04118 0.041181     1    30  0.7818 0.38361  
## landuse:precip 0.09286 0.046428     2    30  0.8815 0.42463  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.238546 0.119273     2    30  1.7265 0.1951
## precip         0.059379 0.059379     1    30  0.8595 0.3613
## landuse:precip 0.034691 0.017345     2    30  0.2511 0.7796

4.7 5 minute MWD

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        0.16741 0.08371     2 20.580  0.4895 0.6199
## precip         0.43517 0.43517     1 28.778  2.5449 0.1216
## landuse:precip 0.44009 0.22005     2 20.580  1.2869 0.2974
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        1.44609 0.72305     2 21.868  3.5546 0.04606 *
## precip         0.00573 0.00573     1 26.971  0.0282 0.86801  
## landuse:precip 0.13294 0.06647     2 21.277  0.3268 0.72480  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        0.17090 0.08545     2 18.602  0.3335 0.7206
## precip         0.60384 0.60384     1 18.081  2.3568 0.1421
## landuse:precip 0.13578 0.06789     2 18.602  0.2650 0.7701

4.8 5 minute 2 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse         25.486  12.743     2 17.893  0.1477 0.8637
## precip         144.368 144.368     1 29.716  1.6734 0.2058
## landuse:precip 138.516  69.258     2 17.893  0.8028 0.4636
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        262.953 131.477     2 20.208  1.5236 0.2419
## precip           0.614   0.614     1 28.249  0.0071 0.9334
## landuse:precip  17.133   8.566     2 19.442  0.0993 0.9060
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        204.601 102.300     2    30  0.9953 0.3815
## precip          90.142  90.142     1    30  0.8770 0.3565
## landuse:precip 184.757  92.379     2    30  0.8988 0.4177
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.08897 0.04448     2    30  0.1824 0.8341
## precip         0.08931 0.08931     1    30  0.3663 0.5496
## landuse:precip 0.65263 0.32631     2    30  1.3383 0.2775
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        0.60278 0.30139     2 17.829  1.4588 0.2590
## precip         0.18603 0.18603     1 29.187  0.9004 0.3505
## landuse:precip 0.42991 0.21495     2 16.975  1.0404 0.3748
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.49414 0.24707     2    30  0.8954 0.4191
## precip         0.41720 0.41720     1    30  1.5120 0.2284
## landuse:precip 0.41178 0.20589     2    30  0.7462 0.4828

4.9 5 minute 0.25 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse          3.94    1.97     2 23.092  0.0803 0.923113   
## precip         341.74  341.74     1 24.299 13.9426 0.001013 **
## landuse:precip  62.02   31.01     2 23.092  1.2652 0.301037   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        360.27 180.134     2 27.856  1.4587 0.2498
## precip          69.43  69.427     1 30.000  0.5622 0.4592
## landuse:precip 505.52 252.758     2 27.535  2.0468 0.1483
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        221.19 110.594     2 25.052  1.1477 0.3335
## precip         241.96 241.956     1 27.240  2.5109 0.1246
## landuse:precip 103.89  51.946     2 25.052  0.5391 0.5899
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse        0.05170 0.02585     2 22.977  0.4960 0.615349   
## precip         0.51807 0.51807     1 24.917  9.9398 0.004182 **
## landuse:precip 0.05467 0.02733     2 22.977  0.5244 0.598801   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.60901 0.30451     2    30  1.8762 0.1707
## precip         0.15543 0.15543     1    30  0.9576 0.3356
## landuse:precip 0.79190 0.39595     2    30  2.4396 0.1043
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        0.38261 0.19130     2 24.375  2.0485 0.15054  
## precip         0.42734 0.42734     1 28.702  4.5760 0.04106 *
## landuse:precip 0.17359 0.08680     2 24.375  0.9294 0.40832  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.10 5 minute 0.053 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse         45.748  22.874     2 24.548  0.3194 0.7296
## precip         114.908 114.908     1 23.550  1.6044 0.2177
## landuse:precip 188.411  94.206     2 24.548  1.3153 0.2866
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        1030.38  515.19     2 27.693  6.1854  0.006015 ** 
## precip         2384.34 2384.34     1 30.000 28.6266 8.661e-06 ***
## landuse:precip  742.17  371.09     2 27.353  4.4553  0.021140 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse         181.59   90.80     2    30  0.5894 0.560970   
## precip         1344.57 1344.57     1    30  8.7277 0.006045 **
## landuse:precip   21.70   10.85     2    30  0.0704 0.932147   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF DenDF F value Pr(>F)
## landuse        0.00423 0.002115     2    30  0.0239 0.9764
## precip         0.23393 0.233933     1    30  2.6407 0.1146
## landuse:precip 0.36298 0.181490     2    30  2.0487 0.1466
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.59609 0.29805     2    30  3.8148   0.03341 *  
## precip         2.86852 2.86852     1    30 36.7155 1.181e-06 ***
## landuse:precip 0.68207 0.34103     2    30  4.3651   0.02168 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        0.07433 0.03717     2    30  0.2268 0.798460   
## precip         1.52502 1.52502     1    30  9.3050 0.004746 **
## landuse:precip 0.00635 0.00318     2    30  0.0194 0.980815   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.11 5 minute 0.02 mm

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value  Pr(>F)  
## landuse         11.504   5.752     2    30  0.1633 0.85008  
## precip         108.763 108.763     1    30  3.0879 0.08908 .
## landuse:precip   5.499   2.749     2    30  0.0781 0.92510  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse         30.773  15.387     2 22.826  4.0730   0.03071 *  
## precip         243.082 243.082     1 25.455 64.3453 1.962e-08 ***
## landuse:precip  16.505   8.252     2 22.467  2.1844   0.13586    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse          6.022   3.011     2 21.422  0.9562    0.4002    
## precip         168.705 168.705     1 26.991 53.5781 7.145e-08 ***
## landuse:precip   6.465   3.233     2 21.422  1.0266    0.3752    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value   Pr(>F)   
## landuse        0.11786 0.05893     2    30  0.3942 0.677673   
## precip         1.33738 1.33738     1    30  8.9457 0.005514 **
## landuse:precip 0.21808 0.10904     2    30  0.7294 0.490571   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.20095 0.10047     2 21.861  1.6924    0.2073    
## precip         2.72769 2.72769     1 27.710 45.9462 2.445e-07 ***
## landuse:precip 0.16265 0.08133     2 21.215  1.3699    0.2757    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        0.15263 0.07632     2 20.772  1.9679   0.16495    
## precip         2.21443 2.21443     1 28.214 57.1017 2.978e-08 ***
## landuse:precip 0.24548 0.12274     2 20.772  3.1649   0.06312 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.12 NRCS hand method

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## landuse        2093.8 1046.87     2 18.889 21.3706 1.411e-05 ***
## precip           53.8   53.80     1 27.343  1.0983 0.3038046    
## landuse:precip 1167.9  583.95     2 18.982 11.9206 0.0004431 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        1251.63  625.82     2    29  7.6618 0.0021309 ** 
## precip         1529.48 1529.48     1    29 18.7253 0.0001635 ***
## landuse:precip  928.86  464.43     2    29  5.6860 0.0082531 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value   Pr(>F)   
## landuse         207.53  103.76     2 16.542  1.0882 0.359743   
## precip         1155.97 1155.97     1 18.062 12.1235 0.002651 **
## landuse:precip  744.85  372.43     2 16.542  3.9059 0.040800 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

4.13 Bulk Density

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value  Pr(>F)  
## landuse        0.112061 0.056031     2 23.669  3.2959 0.05461 .
## precip         0.106779 0.106779     1 29.171  6.2811 0.01803 *
## landuse:precip 0.060558 0.030279     2 23.669  1.7811 0.19030  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value Pr(>F)  
## landuse        0.037800 0.018900     2    30  0.7264 0.4920  
## precip         0.145961 0.145961     1    30  5.6099 0.0245 *
## landuse:precip 0.017941 0.008971     2    30  0.3448 0.7111  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        0.030539 0.015269     2 21.597  1.3400 0.2828
## precip         0.018907 0.018907     1 29.346  1.6592 0.2078
## landuse:precip 0.022982 0.011491     2 21.597  1.0084 0.3813

5 Ratio Graphs

5.1 Total Organic Carbon to Nitrogen

## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        2919.0 1459.50     2    27 29.2361 1.753e-07 ***
## precip         2990.0 1495.00     2    27 29.9472 1.403e-07 ***
## landuse:precip 1979.9  494.98     4    27  9.9152 4.539e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF  DenDF F value Pr(>F)
## landuse        132.975  66.488     2 11.647  1.6235 0.2388
## precip          62.711  62.711     1 29.377  1.5313 0.2257
## landuse:precip 125.403  62.701     2 10.936  1.5311 0.2593
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        165.906  82.953     2    27 13.9815 6.799e-05 ***
## precip         214.256 107.128     2    27 18.0562 1.052e-05 ***
## landuse:precip  88.463  22.116     4    27  3.7275   0.01534 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## function (...)  .Primitive("c")
## [1] "0-5 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.202409 0.101205     2    27 22.9987 1.475e-06 ***
## precip         0.071087 0.035544     2    27  8.0773  0.001781 ** 
## landuse:precip 0.108997 0.027249     4    27  6.1924  0.001138 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "5-10 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq  Mean Sq NumDF DenDF F value    Pr(>F)    
## landuse        0.102732 0.051366     2    27 13.4444 8.875e-05 ***
## precip         0.163367 0.081683     2    27 21.3797 2.721e-06 ***
## landuse:precip 0.077402 0.019350     4    27  5.0647  0.003545 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "10-15 cm"
## Type III Analysis of Variance Table with Satterthwaite's method
##                 Sum Sq  Mean Sq NumDF  DenDF F value   Pr(>F)    
## landuse        0.07008 0.035042     2 22.339  8.7260 0.001583 ** 
## precip         0.34671 0.173357     2 11.837 43.1681 3.65e-06 ***
## landuse:precip 0.02183 0.005457     4 22.339  1.3588 0.279796    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1