Plots and analyses in Feeney et al. (2022): https://doi.org/10.1016/j.crm.2022.100439

Christopher Feeney 1st of June 2022 
Any queries, please email chrfee@ceh.ac.uk:

The R scripts in this folder contain all the necessary code to complete analysis and visualisation of model results post-use of CAESAR-Lisflood.
In alphabetical order, these scripts do the following:

R script:		Description of script's function:
---------		---------------------------------

Calibration_DoDs	Maps the simulated river channel changes between 1976 & 2018 for the Low and High Erosion models (Fig. 4 in the main text of the paper)
			The data for this can be found in the 'Calibrated Runs' folder which should be included with the R scripts.

Calibration_TimeSeries	Plots some of the figures on the annual frequencies of geomorphologically effective flows (as determined by JBA, 2009) of 75 & 105 m3/s.
			These plots form parts B & C of Fig. S2 in the Supplementary Materials of the Feeney et al. (2022) paper linked at the top.
			The data for this can be found in the 'Flows' folder which should be included with the R scripts.

DoD_histograms		Plots histograms and maps comparing the vertical erosion & deposition amounts during the 2006-2018 part of the model calibration period.
			Data are derived from DEMs of Difference (DoDs), with model outputs compared with data derived from EA LiDAR data (resampled to 15m grids).
			These plots are shown in Fig. S4 in the Supplementary Materials of the Feeney et al. (2022) paper linked at the top.
			The data for this can be found in the 'DoD' folder which should be included with the R scripts.

Flow_Stats		Computes a series of annual flow series statistics (annual maximum, median and flows exceeded 5 & 10% of the time) for the 33 years of each
			flow magnitude scenario (2018-2050). This creates Fig. 3a in the main text of the paper (3b is created using the next script below). Also, a
			analyses and plots (not included in the paper) are performed, including the probability of a geomorphologically effective flow threshold (as
			defined by JBA, 2009) being exceed in a given year within each of the 6 flow magnitude scenarios.
			The data for this can be found in the 'Flows' folder which should be included with the R scripts.

Flow_TimeSeries		Plots the annual frequency of geomorphologically effective flows under each flow magnitude scenario (Fig. 3b in the main text of the paper).
			The data for this can be found in the 'Flows' folder which should be included with the R scripts.

Scenario_erodep_pdfs	Similar to the function of 'DoD_histograms' only this covers the scenario simulations for the Low & High Erosion models under increasing
			flow magnitudes. Produces the plots in Fig. 6 in the main text of the paper.
			The data for this can be found in the 'Scenarios' subfolder of the 'DoD' folder which should be included with the R scripts.

Scenario_TimeSeries	Plots time-series of erosion & deposition occurring at each of the 40 transmission towers present in the Mersey Valley. Two main figures are
			produced: Time-series for the Low Erosion model (Fig. 7 in the main text of the paper) and for the High Erosion model (Fig. 8). Each figure is
			split into 40 panels (one for each tower), with 6 lines each representing the effects of a flow magnitude scenario. Extra figures (not included
			in the paper) focus on the 4 most-threatened towers in more detail.
			The data for this can be found in the 'Towers' folder which should be included with the R scripts.