Input files and model results for a schematized morphodynamic model of channel response of the lower Rhine River to climate change and human intervention. 

These model files are related to the manuscript Climate change vs human intervention: Centennial Channel Response to Climate Change and Human Intervention in an Engineered River.  

Input files (for modeling tool SOBEK-RE) include:
- cross-section data, with main channel and floodplain width at each cross section (cross-sections.txt)
- initial mean bed elevation (bedlevel_ini.xlsx)
- initial bed surface and substrate composition per fraction (gs_ini.txt)
- boundary conditions for the reference case and (climate) scenarios: upstream water discharge, upstream sediment flux, downstream water level (#scenario#_QKoln.txt, #scenario#_BaseLevel.txt, #scenario#_SedFlux.txt)

Model results (ModelRuns_2050-2100.xlsx) include for the reference case and climate scenarios: 
- Main channel bed elevation at 2050 and 2100 
- Geometric mean grain size at 2050 and 2100 

In addition, the field data collected by Rijkswaterstaat on bed surface grain size in the Dutch Rhine in 2020 (grainsize_field_DutchRhine.xlsx) is also included. These data has been used as information to derive the model input grain size fraction content.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
.txt file specs are explained below
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% cross-sections.txt %%%%%

description "#cross_section_name" MainChannelWidth(m) FloodplainWidth(m) SedimentTransportWidth(m)
{
          elevation_LowestLevel (m)     flowWidthLowestLevel (m)         totalWidthLowestLevel (m)
          .				.				 .	
	  .				.				 .
			
          	        
         elevation_HighestLevel (m)     flowWidthHighestLevel (m)         totalWidthHighestLevel (m)
}

.
.
.

crosssection "#profile_name" "#associated_cross_section_name_as_tagged_in_description" "branch" chainage_relative_to_branch(m) vertical_translation(m)
.
.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% gs_ini.txt %%%%%

$GSINIT BRANCH num_branch AT chainage_relative_to_branch(m)
LAYNUM = 1 fraction_content_finest_fraction ..  .. .. ..  fraction_content_coarsest_fraction -------- %this is the layer closest to the surface, each layer is 0.5 m thick
.
.
.
LAYNUM = 20 fraction_content_finest_fraction ..  .. .. ..  fraction_content_coarsest_fraction -------- %this is the deapest layer


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% #scenario#_QKoln.txt %%%%%

"YYYY/MM/DD; hh:mm:ss" discharge(m3/s)



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% #scenario#_BaseLevel.txt %%%%%

"YYYY/MM/DD; hh:mm:ss" water_level(m)



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% #scenario#_sedflux.txt %%%%%

%QsF#number_fraction#  ----%number fraction from fine to coarse
X Q1  ... Qn ---> discharge values (m3/s)
Y Qs1 ... Qsn---> sediment flux values associated to discharge values (m3/s)


