Forestry roads in the Purapel fluvial catchment and related changes in sediment connectivity

This dataset contains georeferenced data of forestry roads and sediment connectivity in the Purapel catchment, which drains the Chilean Coastal Range. The forestry road network consists of all the dirt and gravel roads mapped in QGIS by observing open satellite images and vectorial data available during January 2021. The observed data are maps that were listed in the QGIS OpenLayers plugin (https://github.com/sourcepole/qgis-openlayers-plugin), such as Google Satellite (Map data ©2015 Google) and OpenStreetMap ¹, the road network of the Chilean Congress National Library (https://www.bcn.cl/siit/mapas_vectoriales) and compositions of Sentinel 2 images (European Space Agency, courtesy of the U.S. Geological Survey) of the post-2017 fire period. Sediment Connectivity maps were calculated on a 5 m resolution LiDAR DTM using the Connectivity Index². The maps were derived from the stand-alone, free and open-source executable SedInConnect 2.3³ using the Weighting factor of ² and two different targets, which are available as tif files: ICs.tif contains <em>IC_s</em>, the Connectivity Index to the stream network. ICrs.tif contains <em>ICr_s</em>, the Connectivity Index to the road and the stream network. Here, the Road Connectivity, <em>RC </em>(dimensionless) is defined as the difference between both previous maps, with the aim to describe the change in sediment connectivity due to forestry road network: <em>RC = IC_rs - IC_s</em> It is available as RC.tif file. The area of<em> high RC </em>was defined using the percentile 95 (3.12). File RC95.tif is a mask of <em>RC </em><em>≥</em><em> 3.12</em>. The contributing area <em>CA </em>(m²) was calculated using the multiple flow D-infinity approach ⁴ using TauDEM (https://hydrology.usu.edu/taudem/taudem5/downloads.html). The file CA_RC95.tif contains the contributing area (m²) of the surfaces with highest changes in sediment connectivity due to the road network. That is: <em>CA_RC95 = </em>{<em>CA </em>|<em> RC </em><em>≥</em><em> 3.12</em>} The landscape distribution of those surfaces, in terms of proximity to the hilltops and valleys, is described by the density plot of the raster file CA_RC95.tif in R: <pre><code>library("raster") library("ggplot2") CA_RC95&lt;-raster("CA_RC95.tif") CA_RC95&lt;-CA_RC95*0.0025 df = as.data.frame(CA_RC95) df = na.omit(df) ggplot(df,aes(CA_RC95)) + geom_histogram(aes(y=..count..*25),binwidth = 50)+ geom_density(aes(y=50 * ..count..*25), col="blue",size=2, adjust=10000)+ xlab("Contributing Area [ha] \n Hilltop Valley") + ylab("Area [m2]")+ theme(axis.text.x = element_text(face="bold", size=30), plot.title = element_text(color="black", size=40, face="bold",hjust=0.5), axis.title.x=element_text(color="blue", size=40, face="bold"), axis.text.y = element_text(face="bold", size=30), axis.title.y=element_text(color="blue", size=40, face="bold"))+ scale_y_continuous(trans = 'log10')+ ggtitle("Upstream area of surfaces with \n High Road Connectivity (RC &gt; 3.12)") </code></pre> Bibliography 1. OpenStreetMap contributors. Planet dump retrieved from https://planet.osm.org. https://www.openstreetmap.org/ (2017). 2. Cavalli, M., Trevisani, S., Comiti, F. &amp; Marchi, L. Geomorphometric assessment of spatial sediment connectivity in small Alpine catchments. <em>Geomorphology</em> <strong>188</strong>, 31–41 (2013). 3. Crema, S. &amp; Cavalli, M. SedInConnect: a stand-alone, free and open source tool for the assessment of sediment connectivity. <em>Computers and Geosciences</em> <strong>111</strong>, 39–45 (2018). 4. Tarboton, D. G. A new method for the determination of flow directions and upslope areas in grid digital elevation models. <em>Water Resources Research</em> <strong>33</strong>, 309–319 (1997).