diff --git a/code/05_ExtractEnvironment.Rmd b/code/05_ExtractEnvironment.Rmd
index b834e31efac5b54f73f88fa15d7e21e1a5d05eb9..45db8744eb9bb00f5d8d0c0e1e3fdaeabfa04dc3 100644
--- a/code/05_ExtractEnvironment.Rmd
+++ b/code/05_ExtractEnvironment.Rmd
@@ -274,12 +274,21 @@ knitr::kable(tmp.sum,
                   full_width = F, position = "center")
 ```
 
+## 4 Elevation variability
+Reimport output (calculated in `04_buildHeader`)
+```{r, message=F}
+elevation.out <- read_csv("../_derived/elevatr/elevation.out.csv")
+```
+
 
-## 4 Create output and export
+## 5 Create output and export
 
 ```{r}
 soilclim <- header %>% 
   dplyr::select(PlotObservationID) %>% 
+  left_join(elevation.out %>% 
+              dplyr::select(PlotObservationID, Elevation_median, Elevation_q2.5, Elevation_q97.5, Elevation_DEM.res=DEM.res), 
+            by="PlotObservationID") %>% 
   left_join(header.shp@data %>% 
               dplyr::select(PlotObservationID) %>% 
               bind_cols(chelsa.out) %>%
@@ -288,7 +297,6 @@ soilclim <- header %>%
               bind_cols(isric.sd.out) %>% 
               distinct(), 
             by="PlotObservationID")
-
 ```
 
 ```{r, echo=F}