Load data

Read HWSD NetCDF files into a data frame

df <- rbeni::nc_to_df("/alphadata01/bstocker/data/soil/hwsd/hwsd_wieder/data/T_SAND.nc4","T_SAND") %>%
  rename( fsand = myvar ) %>%
  left_join(
    rbeni::nc_to_df("/alphadata01/bstocker/data/soil/hwsd/hwsd_wieder/data/T_CLAY.nc4","T_CLAY") %>%
      rename( fclay = myvar ),
    by = c("lon", "lat")
  ) %>%
  left_join(
    rbeni::nc_to_df("/alphadata01/bstocker/data/soil/hwsd/hwsd_wieder/data/T_OC.nc4","T_OC") %>%
      rename( forg = myvar ),
    by = c("lon", "lat")
  ) %>%
  left_join(
    rbeni::nc_to_df("/alphadata01/bstocker/data/soil/hwsd/hwsd_wieder/data/T_GRAVEL.nc4","T_GRAVEL") %>%
      rename( fgravel = myvar ),
    by = c("lon", "lat")
  )

For testing

df_sub <- df %>% 
  tidyr::drop_na() %>%
  sample_n(10000)

Calculate soil parameters

Calculate soil parameters following Saxton & Rawls 2006

df_sr <- df_sub %>% calc_soilparams(method = "saxtonrawls")
## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

## Warning in log(pwp): NaNs produced

Calculate soil parameters following Ballard et al. 2008

df_bl <- df_sub %>% calc_soilparams(method = "ballard")
## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

## Warning in log(silt): NaNs produced

Field capacity

out <- df_sr %>% 
  select(fc_sr = fc) %>% 
  bind_cols(
    select(df_bl, fc_bl = fc)
  ) %>% 
  analyse_modobs2(mod = "fc_bl", obs = "fc_sr", type = "heat")
## Loading required package: LSD
## Loading required package: ggthemes
## Loading required package: RColorBrewer
out$gg +
  labs(x = "FC Ballard 2008 (m3 m-3)", y = "FC Saxton & Rawls 2006 (m3 m-3)",
       title = "FC")

Permanent wilting point

out <- df_sr %>% 
  select(pwp_sr = pwp) %>% 
  bind_cols(
    select(df_bl, pwp_bl = pwp)
  ) %>% 
  analyse_modobs2(mod = "pwp_bl", obs = "pwp_sr", type = "heat")
out$gg +
  labs(x = "PWP Ballard 2008 (m3 m-3)", y = "PWP Saxton & Rawls 2006 (m3 m-3)",
       title = "PWP")

Water holding capacity

out <- df_sr %>% 
  select(whc_sr = whc) %>% 
  bind_cols(
    select(df_bl, whc_bl = whc)
  ) %>% 
  analyse_modobs2(mod = "whc_bl", obs = "whc_sr", type = "heat")
out$gg +
  labs(x = "WHC Ballard 2008 (m3 m-3)", y = "WHC Saxton & Rawls 2006 (m3 m-3)",
       title = "WHC")