Investigating fold-river interactions for major rivers using a scheme of remotely sensed characteristics of river and fold geomorphology

There are frequently interactions between active folds and major rivers (mean annual water discharges > 70 m3s-1). The major river may incise across the fold, to produce a water gap across the fold, or a bevelling (or lateral planation) of the top of the fold. Alternatively, the major river may be defeated to produce a diversion of the river around the fold, with wind gaps forming across the fold in some cases, or ponding of the river behind the fold. Why a river incises or diverts is often unclear, though influential characteristics and processes have been identified. A new scheme for investigating fold-river interactions has been devised, involving a short description of the major river, climate, and structural geology, and 13 characteristics of river and fold geomorphology: 1) Channel width at location of fold axis, w, 2) Channel-belt width at location of fold axis, cbw, 3) Floodplain width at location of fold axis, fpw, 4) Channel sinuosity, Sc, 5) Braiding index, BI, 6) General river course direction, RCD, 7) Distance from fold core to location of river crossing, C-RC, 8) Distance from fold core to river basin margin, C-BM, 9) Width of geological structure at location of river crossing, Wgs, 10) Estimate of erosion resistance of surface sediments/rocks and deeper sediments/rocks in fold, ERs, ERd, 11) Channel water surface slope at location of fold axis, s, 12) Average channel migration rate, Rm, 13) Estimate of fold total uplift rate, TUR. The first 10 geomorphological characteristics should be readily determinable for nearly all major rivers using widely available satellite imagery and fine scale geological maps. The last 3 characteristics should be determinable for most major rivers where other data sources are available. This study demonstrates the methodology of this scheme, using the example of the major rivers Karun and Dez interacting with active folds in the foreland basin tectonic setting of lowland south-west Iran. For the rivers Karun and Dez (mean annual water discharges 575 m3s-1 and 230 m3s-1, respectively), it was found that geomorphological characteristics Nos. 2, 3 and 7 had statistically significant differences (p-value ≤ 0.05) between the categories of river incision across a fold and river diversion around a fold. For river incision, at the fold axis, channel-belt width was always < 2.7 km, and floodplain width was generally (80 % of cases) < 5.7 km; whereas for river diversion, at the projection of the fold axis, these two characteristics had a wide range of values. For river incision, the distance from the fold core to the location where the river channel crossed the fold axis, was generally (80 % of cases) ≤ 8.5 km; whereas for river diversion, this distance was always > 22 km. Since it is highly likely that different characteristics will be important for other major rivers interacting with other folds, it is recommended that this scheme is now used to investigate a variety of major rivers from across the globe. By comparing the same parameters for different major rivers, a better understanding of fold-river interactions should be achieved.


Introduction
Interpreting the interactions between rivers and tectonics can be challenging.Principally, this is because rivers are inherently variable and complex, influenced by a wide range of both autogenic factors that include topography, hydrology and sedimentology, and allogenic factors that include structural geology and active tectonics, plus human activities, climate and relative sea-level (or base level) changes [1][2][3][4][5][6][7][8][9].Disentangling the various internal and external factors and their influences on geomorphology can be difficult.However, for major rivers, with mean annual water discharges of 70 m 3 s −1 or more [10], interacting with active folds over horizontal spatial scales of metres to tens of kilometres (river channel dimensions to fold dimensions), the difficulties are lessened, especially at locations upstream of coastal plain-valleys [11][12][13].This is because for a single major river at such scales, climate and rates of sediment supply from the basin hinterland are likely to be similar, as climate zones typically extend over scales of hundreds of kilometres [14][15][16], and upstream of the extent of the backwater length (typically a distance of more than 150 km from the shoreline) the influences of relative sea-level changes are likely to be minimal [12,17,18].Hence, at these river reach scales, the significant allogenic factors will be limited to tectonics and human activities, with prominent human impacts being limited to the last few millennia [13,[19][20][21].
Major rivers frequently interact with active folds, particularly as transverse rivers in foreland basin systems, where folds oriented roughly parallel to the orogenic axis may form a succession of "obstacles" to river courses, particularly in the orogenic wedge and foredeep [22,23].Conceptual models of the interactions between transverse rivers and growing folds have been constructed [22,[24][25][26][27][28].Such models indicate that where rates of river aggradation exceed rates of structural uplift associated with the fold, a river will flow without impedance across the fold and may bevel off the top of the emerging fold with little or no topographic relief developing [27,28].Where a fold does develop a surface topographic expression, a river will either flow across the fold by maintaining basinward-dipping channel slopes across the fold, or it will be defeated by the growing fold.To maintain a transverse course across a fold, a river needs sufficient stream power to erode and incise into the crest and across the axis of the fold at a rate greater than the difference between the rates of structural uplift and the rates of river aggradation [22,29].Whilst the precise controls on river erosion are debated, due to factors such as bed armouring [30,31], it is likely that river erosion into bedrock and sediments will increase with stream power.If the river is defeated, then it will be diverted around the fold by channel migrations or avulsions to flow through structural low points, frequently flowing initially roughly parallel to the fold axis and thence around the nose of the fold.Alternatively, the river may be ponded in a basin upstream of the fold [22,25,27,32].
According to such conceptual models, the responses of rivers and major rivers should be fairly predictable.A river may incise across an active fold as a water gap (a river valley of a maintained river course) or it may be defeated by the fold and diverted to leave a wind gap (a dry valley of a previous river course), with the configuration of these water and wind gaps varying with a number of factors, such as the type of fold [22,27,33,34].For instance, detachment folds would be expected to have a wind gap near the centre of the fold and a water gap near the propagating fold tip, whilst fault bend folds would be expected to a have a number of wind gaps across the length of the fold, with the defeated rivers diverted parallel to the fold axis [33,34].Whilst conceptually it is clear that a major river should incise across an active fold in some cases and divert around it in other cases, in practice it is often unclear as to how and why this occurs.For instance, there is a seemingly paradoxical tendency for a number of major rivers to transect many growing anticlines in the vicinity of their greatest structural and topographic relief [35][36][37].By contrast, some rivers frequently cross a growing fold near to the laterally propagating tip or nose of the fold [27,38].Alternatively, rivers may be diverted around the fold tips of laterally propagating anticlinal fold segments until these fold segments coalesce; after which the river may divert to feed a longitudinal river or it may incise across the coalesced fold at the topographic low of the merger location [39].
These different responses are probably due to changes in the fold-river interactions with time and the variable and complex nature of river systems [8,13,40,41].There may be different reaction, relaxation and recurrence times for events [42], multiple processes may act in combination to produce a specific phenomenon [42,43], different factors may result in similar effects [41], a river system may not adjust in a progressive and systematic fashion to modifications [44], and a river system may be dominated by autogenic processes and exhibit variability independent of external factors, due to systems of non-linearity or self-organised criticality [9,44,45].Nevertheless, with such systems there may be characteristics of the river or the fold which act as thresholds which the river needs to cross for the dynamic equilibrium of river incision across an active fold to develop and be maintained [42].
The characteristics which may act as thresholds will probably include those associated with the main controlling variables for the persistence of an antecedent river across a growing fold, as shown in Table 1 [22,25,28].
Table 1.The main controlling variables for the persistence of antecedent rivers crossing growing folds (Modified from [22,25,28]).

Variable Effect
Rate of sediment aggradation and rate of structural uplift Lower rates of sediment aggradation and lower rates of structural uplift promote persistence of an antecedent river, due to less erosion of the fold hanging wall being required Erosion resistance of rocks and sediments within fold Lower erosion resistances (thick alluvial strata, poor cementation and readily erodible bedrock) mean that lower stream power are required, thus promoting persistence of an antecedent river Water discharge of river Higher water discharges and higher stream power promote persistence of an antecedent river Stream power, flow depth, channel width, channel water surface slope of river Higher stream power promote persistence of an antecedent river.Narrower channel widths and steeper channel water surface slopes promote persistence of the antecedent river, due to associated increased stream power

Sediment load
Increased sediment load decreases proportion of stream power available for bed erosion, mantling of the bed with sediment precludes erosion of bed; thus, reduced sediment load may promote persistence of an antecedent river

Width of geological structure
Widening of a geological structures causes reduced channel water surface slopes and stream power; thus, narrower geological structures promote persistence of an antecedent river

Transverse structures
Transverse structures, such as faults, provide zones of less erosion resistant rocks that cut across structures, exploited by antecedent rivers The influences of some of these controlling variables are quite intricate, particularly those associated with river hydrology and sediment load [31].For instance, a river crossing a fold will produce aggradation upstream and downstream of the fold in a dynamic equilibrium, in which sufficient foreland-dipping channel slopes for producing erosive stream power across the zone of greatest fold uplift are maintained [22,27,46,47].If upstream or downstream aggradation is insufficient, as may be the case with reduced sediment load, then the river may be defeated and diverted around the fold [22,24].If upstream aggradation is excessive, then the river may be defeated by producing slopes that promote channel migrations or avulsions to other upstream locations [27,38,48].If downstream aggradation is excessive, then the river may also be defeated by reducing channel slopes to such an extent that stream power are insufficient to maintain erosion into the fold and maintain transport away of the eroded material [47,49].Nevertheless, Table 1 still provides an adequate foundation for differentiating between river incision across a fold and river diversion around a fold.Some of the controlling variables, such as stream power, flow depth, and sediment load, involve characteristics which need to be determined by fieldwork; whereas other controlling variables, such as width of geological structure, involve characteristics which can be determined relatively easily from remote sensing imagery and fine scale geological maps.

Aim of the Study-A Scheme for Investigating Fold-River Interactions Using Remote Sensing
The aim of this study is to demonstrate a new scheme which uses a short description of the major river and 13 remotely sensed characteristics of river and fold geomorphology to investigate fold-river interactions.The short description of the major river should include river measurements (including mean annual discharge) and short descriptions of the river course, climate, and structural geology.The first 10 geomorphological characteristics should be readily determinable from widely available remote sensing imagery and fine scale geological maps.This use of remote sensing allows a large number of major rivers to be investigated relatively easily, including those in remote or inaccessible areas, without recourse to expensive fieldwork.The last three geomorphological characteristics should be determinable where additional data sources are available.This study utilises the example of the major rivers Karun and Dez interacting with folds in lowland south-west Iran to show how to apply the scheme in practice.

Selection of 13 Remotely Sensed Characteristics of River and Fold Geomorphology
Remote sensing imagery and fine-scale geological maps have the advantage of being widely available data sources which only need processing for interpretation, rather than detailed fieldwork, but have the drawback that certain parameters, such as sediment grain size, sediment load, flow velocity, and channel depth, cannot be measured accurately from them.A number of the controlling variables in Table 1 involve geomorphological characteristics which are readily determinable from remote sensing and fine scale geological and topographical maps, as are other significant geomorphological characteristics, such as channel width, that are associated with other conceptual models [22,[25][26][27][28]. Also, previous detailed studies on interactions between specific major rivers and tectonics, particularly those of Jorgensen [50] involving rivers in western U.S.A., Lavé and Avouac [51,52] involving upland rivers in Nepal, and Woodbridge [13] involving lowland rivers in south-west Iran, have identified useful characteristics determinable using remote sensing and geological maps in their investigations of such interactions.All of these data sources have been used to compile a suite of 13 useful geomorphological characteristics to be determined in investigations of fold-river interactions: (1) Channel width at location of fold axis, w (2) Channel-belt width at location of fold axis, cbw (3) Floodplain width at location of fold axis, fpw (4) Channel sinuosity, Sc (5) Braiding index, BI (6) General river course direction, RCD (7) Distance from fold core to location of river crossing, C-RC (8) Distance from fold core to river basin margin, C-BM (9) Width of geological structure at location of river crossing, Wgs (10) Estimate of erosion resistance of surface sediments/rocks and deeper sediments/rocks in fold, ERs, ERd (11) Channel water surface slope at location of fold axis, s (12) Average channel migration rate, Rm (13) Estimate of fold total uplift rate, TUR In summary, channel width at the location of the fold axis should be a useful parameter since the conceptual model of Amos and Burbank [25] and the studies of Lavé and Avouac [51,52] in Nepal indicate that channel width may act as a key characteristic of river responses, with channel narrowing to enhance incision rates, apparently taking precedence over other changes for upland rivers crossing rapidly uplifting folds [27,[52][53][54].Channel-belt width and floodplain width at the location of the fold axis should be useful parameters, since narrowing of the channel-belt and narrowing of the floodplain will increase the proportion of stream power available for vertical erosion and thus promote the maintenance of a river incising across a fold.The study of Woodbridge [13] demonstrated the importance of channel-belt width, with a narrow average channel-belt width of less than c. 2.7 km being hypothesised as a threshold needed for the rivers Karun and Dez to produce and maintain river incision across a fold in lowland south-west Iran.Channel sinuosity and braiding index should both be useful parameters, since the study of Woodbridge [13] found trends for both reduced sinuosity and braiding index for river reaches incising across a fold; though, as with the studies of Jorgensen [50] in U.S.A. and Zámolyi et al. [55] in Hungary, these trends did not always achieve statistical significance.General river course direction should be a useful parameter, as the study of Woodbridge [13] found a tendency for river incision across a fold to have a general river course direction orthogonal to the fold axis for the river reaches which crossed the fold, whereas river diversion had a general river course parallel to the fold axis upstream of the fold, followed by a change in river course bearing of about 20 • -70 • to flow around the fold.
Distance from the fold core to the location where the river crosses the fold axis should be a useful discriminative parameter since, naturally, there is very strong tendency for river incision across a fold to occur between the fold core and the fold nose, and for river diversion to occur beyond the fold nose [13,20,56].Distance from the fold core to the river basin margin, should be a useful parameter if the timing of initial fold-river interactions is important, as hypothesised by Woodbridge [13].Where a river incises across a fold due to it initially encountering the fold as a small, emerging fold, the fold core location is likely to be within the margins of the drainage basin of the river crossing the fold axis (positive measurement); whereas where a river diverts around a fold due to it initially encountering the fold as a larger, more developed fold, the fold core location is likely to be beyond the margins the drainage basin of the river crossing the fold axis, or its projection (negative measurement) [13,20].Width of geological structure should be a useful parameter, since, as shown in Table 1, the conceptual model of Burbank et al. [22] indicates that narrow geological structures promote river incision across a fold by avoiding the reduced channel slopes, stream power, and vertical erosion associated with widening geological structures.The erosion resistances of sediments and rocks in a fold can be estimated from fine scale geological maps where details of the sediments and rock types are known, and should be useful parameters since the conceptual models of Burbank et al. [22] and Bufe et al. [28] indicate that low erosion resistances promote river incision across a fold and river bevelling of the top of a fold.Also, some studies, such as that on the meandering of the River Dniester by Yeromenko and Ivanov [57], have found that variations in erosion resistances of rocks and sediments were significant in influencing river responses; though other studies, such as those on rivers and growing folds in northern Alaska reviewed by Burbank et al. [22] have found that variations in erosion resistances were not.
To determine the last three geomorphological characteristics precisely, data sources in addition to one set of remote sensing imagery and one set of fine scale geological maps are preferable.Whilst slope can be measured from a DEM, greater precision for channel water surface slope measurement will be obtained from other data sources, such as precise hydrological and topographical surveys.Channel water surface slope at the location of the fold axis should be a useful parameter, since it was found to be a key characteristic for upland rivers in studies by Lavé and Avouac [51,52] in Nepal, Yanites et al. [58] in Taiwan, and by Amos and Burbank [25] in New Zealand.Average channel migration rate over time intervals of about 20-40 years should be a useful parameter, since lateral migration rates have been found to be significant in studies of river incision across a fold [13] and river lateral planation of the top of a fold [28].Fold total uplift rate, as estimated from additional data sources, should be a useful parameter since, in a number of conceptual models, low rates of structural uplift promote the maintenance of a river incising across a fold [22,24,28].

Summary of Methods
There are three main elements to the application of the scheme to a specific major river or river system: (1) A short description of the river, including its course, and the climate and structural geology of the region through which it flows (2) Measurement of geomorphological characteristics Nos. 1 to 10 (3) Measurement of geomorphological characteristics Nos.11 to 13

Short Description of River
The short description of the river introduces the major river and the context of the fold-river interactions.It should include data on river length, drainage basin area, mean annual water discharge, seasonality of discharge, and major direct human impacts on the river, and a short description of the river course.It should also include short descriptions of the regional climate and structural geology, with some details of the tectonic setting and the types of faults and folds.The short description of the river can be supplemented by maps of the river system and structural geology.

Measurement of Geomorphological Characteristics Nos. 1 to 10
The measurement of the first 10 characteristics of river and fold geomorphology provides the main data for investigating different fold-river interactions.The only data sources needed to determine these 10 remotely sensed characteristics are: high-resolution remote sensing images, fine scale geological maps (preferably at 1:100,000 scale or finer), and maps of oil and gas fields and seismic survey sections (in cases where there are sub-surface folds).Such widespread data sources should be available for most of the major rivers of the world.
For characteristics Nos. 1 to 3 (channel width, channel-belt width, and floodplain width), the measurements are made solely at the location of the fold axis or its projection.This is because these characteristics vary continuously along the length of the river and their measurements are not dependent on how the river is sub-divided into river reaches.By contrast, characteristics Nos. 4 and 5 (channel sinuosity and braiding index) are heavily dependent on how the river is sub-divided into river reaches.Hence, for these characteristics the measurements are made for river reaches immediately upstream of the fold, across the fold axis (or its projection), and immediately downstream of the fold.This is done so that changes in these characteristics associated with the fold can be more easily differentiated from changes due to the sub-division into river reaches and other variations.Similarly, for characteristic No. 6 (general river course direction), the measurements are made for river reaches immediately upstream of the fold, across the fold axis (or its projection), and immediately downstream of the fold.For general river course direction, it is changes relative to the fold axis which are more indicative of changes associated with the fold.Hence, these measurements are also made relative to the fold axis, and there is an emphasis on changes in river course direction between river reaches immediately upstream of the fold, across the fold axis, and immediately downstream of the fold.Characteristics Nos. 7 to 10 (distances from the fold core to the river crossing and river basin margin, width of geological structure, and estimate of erosion resistance) are mainly associated with the structural geology, rocks and sediment of the fold.Hence, for these characteristics, the measurements are made relative to structures of the fold, especially the fold core, the fold axis, and the fold limbs.

Measurement of Geomorphological Characteristics Nos. 11 to 13
The measurement of the last three geomorphological characteristics provides additional data for investigating different fold-river interactions.The data sources needed for these characteristics may not be available for all major rivers worldwide, hence they may be considered as supplementary characteristics.The additional data sources could be precise hydrological or topographical surveys of the river, databases superimposing two sets of high-resolution remote sensing imagery separated by about 20-40 years, and data relating to vertical Earth surface movements, such as dating of displaced geomorphic surfaces e.g., [59], repeated precision GPS surveys, and precise levelling e.g., [60].
For characteristic No. 11 (channel water surface slope) the measurements are made solely at the location of the fold axis or its projection.This is because channel water surface slope is highly variable and the fold axis is a key location where similar conditions can be compared.Characteristic No. 12 (average channel migration rate) is heavily dependent on how the river is sub-divided into reaches.Hence, measurements for this characteristic are made for river reaches immediately upstream of the fold, across the fold axis (or its projection), and immediately downstream of the fold.Characteristic No. 13 (fold total uplift rate) is mainly associated with the structural geology of the fold.Hence, the estimates or measurements for this characteristic are made for the crest of the fold relative to the surrounding region.

Details of Methods for the 13 Geomorphological Characteristics, as Applied to the Rivers Karun and Dez
To introduce and demonstrate the use of the new scheme in practice, it has been applied to the River Karun and River Dez in the province of Khuzestan in lowland south-west Iran, as an example.As shown in Figures 1 and 2, the major rivers Karun and Dez (mean annual water discharges c. 575 m 3 s −1 and 230 m 3 s −1 , respectively) flow from the Zagros orogen in the N and NE across the Upper and Lower Khuzestan Plains into the Mesopotamian-Persian Gulf Foreland Basin to the S and SW [61].Their interactions with folds within the Upper and Lower Khuzestan Plains have been subjected to detailed investigations, as described by Woodbridge [13], Woodbridge and Frostick [56] and Woodbridge et al. [20].The data in these investigations was used to provide short descriptions of the rivers Karun and Dez, as given in Section 4.1.The data was also used to demonstrate the measurement of each of the 13 geomorphological characteristics, by using the example of the Sardarabad Anticline (SDA on Figure 2) to the north-west of Band-e Qir Figure 1, and its interactions with the River Dez (river incision across the fold) and the River Karun (Shuteyt branch) (river diversion around the fold).The Sardarabad Anticline appears to be a doubly plunging, segmented, asymmetric detachment fold which is about 58 km long × 9 km wide, and which rises to more than 70 m above the surrounding plains.The fold axis is oriented roughly ESE-WNW, curving to SE-NW at the eastern end, where it apparently merges with a roughly N-S oriented oblique lateral ramp [13,[62][63][64][65].

Measurement of Geomorphological Characteristics Nos. 1 to 10
For determining these 10 characteristics for the rivers Karun and Dez, the remote sensing images used were 30 m resolution false-colour Landsat Enhanced Thematic Mapper Plus (ETM+) images (28 July 2001 and4 August 2001) with Band 4 (near-infrared, 750-900 nm) displayed red, Band 3 (red, 630-690 nm) displayed green, Band 2 (visible green, 525-605 nm) displayed blue, and pan-sharpened with pan-chromatic Band 8 [66].The fine-scale geological maps used were mainly 1:100,000 scale geological maps, such as "Sheet 20824E Mulla Sani" of the Iranian Oil Operating Companies IOOC [67].The maps of oil and gas fields and seismic survey sections were from a variety of sources [68][69][70][71].The Landsat ETM+ images and detailed surveys of the rivers undertaken by the Dez Ab Engineering Company from 1997-2000 were used to sub-divide the main river courses of the Karun and Dez from the vicinity of Gotvand and Dezful to the Persian Gulf into a succession of straight-line river reaches.The average river reach length was 8.0 km, with an extreme range of 0.8-50.5 km.A river reach was defined as a length of river channel with a relatively homogeneous discharge and morphology [72].Significant changes in general river course direction, river planform, and river morphology were used to demarcate the end of one reach and the start of the next.This sub-division into river reaches, whilst necessarily subjective, facilitated the measurement of characteristics associated with river reaches, such as channel sinuosity and general river course direction [13,20].

Channel Width at Location of Fold Axis (or its Projection)
Symbol: w Units: m (quoted to two decimal places) Measurement location: Where river channel crosses the fold axis (or its projection) Channel width is defined as the maximum extent of the river channel water surface, as distinguished on remote sensing images (or survey records), measured orthogonal to the river thalweg.Since the channel-forming discharge is commonly taken as the bankfull discharge and channel width varies significantly with river discharge, the aim is to measure the width between the channel banks at bankfull discharge [42,73].In practice, channel width also varies with distance along the channel, local irregularities and outcrops, vegetation, human impacts, and other factors, so it is recommended that the distance between the channel banks is measured from remote sensing images of a single date, preferably at a time of relatively high flows.Whilst variations could be reduced by determining average channel width over a distance of one or two meander wavelengths [74][75][76], this is not recommended since subtle changes in channel width would be missed in the frequent cases where the zone of maximal uplift is considerably smaller than the meander wavelength of a major river.Instead, for a single-thread meandering channel pattern, the width of the channel at or very near to the fold axis should be measured, with care to avoid measuring at localised broadening or constriction of the channel.For a multi-thread braided channel pattern, the widths of all channels at the fold axis location should be measured, and the sum recorded.For anastomosing or anabranching channel patterns, the widths of all channels associated with the main branch of the river at the fold axis location should be measured, and the sum recorded [77].
For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, channel width, w = 202,19 m, at the location where the projection of the fold axis intersects with the thalweg of the main river channel, as shown in Figure 3. Channel-belt width is defined as the maximum extent of the channel-belt of the river, as distinguished on remote sensing images, measured orthogonal to the axis of the river reach.For single-thread meandering and straight channel patterns, the measurement is to the extremities of all channels, abandoned channels, meanders, levées, crevasse channels and splays, oxbows, and meander scars that are associated with the active river channel.For a multi-thread braided channel pattern, the measurement is to the extremities of all channels, bars, islands, and abandoned channels associated with the active river channel [73].For anastomosing or anabranching channel patterns, the measurement is to the extremities of the main active river channels, with any anabranches clearly separated by floodplain areas being considered as discrete channel-belts not included in the measurement [77].Where there is uncertainty, such as discriminating between extensive braided rivers and discrete channels of anastomosing rivers, the default is to use the larger channel-belt width measurement.
For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, channel-belt width, cbw = 2.051 km, as shown in Figure 3.In this figure, the channel-belt of the River Karun (Shuteyt) is highlighted in light red (that of the River Gargar is highlighted in yellow), and the channel-belt width measurement is indicated by the white and black checked straight line.The location of the measurement is the same as that for geomorphological characteristic No. 1.

Floodplain Width at Location of Fold Axis (or its Projection)
Symbol: fpw Units: km (quoted to three decimal places) Measurement location: Where river channel crosses the fold axis (or its projection)

Channel-Belt Width at Location of Fold Axis (or its Projection)
Symbol: cbw Units: km (quoted to three decimal places) Measurement location: Where river channel crosses the fold axis (or its projection) Channel-belt width is defined as the maximum extent of the channel-belt of the river, as distinguished on remote sensing images, measured orthogonal to the axis of the river reach.For single-thread meandering and straight channel patterns, the measurement is to the extremities of all channels, abandoned channels, meanders, levées, crevasse channels and splays, oxbows, and meander scars that are associated with the active river channel.For a multi-thread braided channel pattern, the measurement is to the extremities of all channels, bars, islands, and abandoned channels associated with the active river channel [73].For anastomosing or anabranching channel patterns, the measurement is to the extremities of the main active river channels, with any anabranches clearly separated by floodplain areas being considered as discrete channel-belts not included in the measurement [77].Where there is uncertainty, such as discriminating between extensive braided rivers and discrete channels of anastomosing rivers, the default is to use the larger channel-belt width measurement.
For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, channel-belt width, cbw = 2.051 km, as shown in Figure 3.In this figure, the channel-belt of the River Karun (Shuteyt) is highlighted in light red (that of the River Gargar is highlighted in yellow), and the channel-belt width measurement is indicated by the white and black checked straight line.The location of the measurement is the same as that for geomorphological characteristic No. 1.

Floodplain Width at Location of Fold Axis (or its Projection)
Symbol: fpw Units: km (quoted to three decimal places) Measurement location: Where river channel crosses the fold axis (or its projection) Floodplain width is defined as the maximum extent of the floodplain of the river, as distinguished on remote sensing images, measured orthogonal to the axis of the river valley.The floodplain width can vary from the channel-belt width to many tens of channel-belt widths [73].The margins of the floodplain are usually fairly clear due to a slight change in slope at the base of the enclosing valley walls.Interpretive difficulties with floodplain width may arise where two or more major rivers occupy a large plain, especially a large coastal plain, and, in these cases, the measurement is to the extremities of the floodplain of the streams and wetlands within the drainage basin of the major river in question [78].
For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, floodplain width, fpw = 17.603 km, as shown in Figure 3.In this figure, the floodplain width measurement is indicated by the light brown and black checked straight line.The location of the measurement is the same as that for geomorphological characteristic No. 1.

Channel Sinuosity
Symbol: Sc No units (ratio quoted to three decimal places) Measurement location: River reaches immediately upstream of fold, across the fold axis (or its projection), and immediately downstream of fold Channel sinuosity is the ratio defined by the equation Sc = L c /L v where L c is channel length (m), and L v is straight-line valley length (m) [42].The channel length is the total distance between the two ends of the river reach measured along the thalweg of the main channel.For multi-thread braided, anastomosing, and anabranching channel patterns there can be interpretive difficulties regarding the main channel thalweg, though, generally, it should be interpreted as the course of the broadest channel.The straight-line valley length is the distance between the two ends of the river reach measured in a straight line along the axis of the river reach.Measurements are made for river reaches immediately upstream of the fold, across the fold axis (or its projection), and immediately downstream of the fold, to elucidate any changes in channel sinuosity associated with the fold.
For the example of the River Dez incising across the Sardarabad Anticline, channel sinuosity, Sc = 1417 (immediately upstream of fold); 1120 (across fold axis); 1585 (immediately downstream of fold), as shown in Figure 4.

Remote Sens. 2019, 11, x FOR PEER REVIEW 12 of 39
Floodplain width is defined as the maximum extent of the floodplain of the river, as distinguished on remote sensing images, measured orthogonal to the axis of the river valley.The floodplain width can vary from the channel-belt width to many tens of channel-belt widths [73].The margins of the floodplain are usually fairly clear due to a slight change in slope at the base of the enclosing valley walls.Interpretive difficulties with floodplain width may arise where two or more major rivers occupy a large plain, especially a large coastal plain, and, in these cases, the measurement is to the extremities of the floodplain of the streams and wetlands within the drainage basin of the major river in question [78].
For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, floodplain width, fpw = 17.603 km, as shown in Figure 3.In this figure, the floodplain width measurement is indicated by the light brown and black checked straight line.The location of the measurement is the same as that for geomorphological characteristic No. 1.

Channel Sinuosity
Symbol: Sc No units (ratio quoted to three decimal places) Measurement location: River reaches immediately upstream of fold, across the fold axis (or its projection), and immediately downstream of fold Channel sinuosity is the ratio defined by the equation Sc = Lc/Lv where Lc is channel length (m), and Lv is straight-line valley length (m) [42].The channel length is the total distance between the two ends of the river reach measured along the thalweg of the main channel.For multi-thread braided, anastomosing, and anabranching channel patterns there can be interpretive difficulties regarding the main channel thalweg, though, generally, it should be interpreted as the course of the broadest channel.The straight-line valley length is the distance between the two ends of the river reach measured in a straight line along the axis of the river reach.Measurements are made for river reaches immediately upstream of the fold, across the fold axis (or its projection), and immediately downstream of the fold, to elucidate any changes in channel sinuosity associated with the fold.
For the example of the River Dez incising across the Sardarabad Anticline, channel sinuosity, Sc = 1417 (immediately upstream of fold); 1120 (across fold axis); 1585 (immediately downstream of fold), as shown in Figure 4.

Braiding Index
Symbol: BI No units (index quoted to one decimal place) Measurement location: River reaches immediately upstream of fold, across the fold axis (or its projection), and immediately downstream of fold The braiding index is a measure of the intensity of braiding, and for a river reach can be defined as the channel count index of the mean number of anabranches (or links) per river cross-section for that reach [79,80].Since the intensity of braiding varies with flow stage [81], it is recommended that measurements are undertaken from remote sensing images of a single date, preferably at a time of relatively high flows for compatibility with other measurements, such as channel width.The river reach is sub-divided into river cross-sections orthogonal to the valley axis which are approximately 1 km apart.For each river cross-section, the number of distinct anabranches is counted and the mean for the entire river reach is calculated.For single-thread meandering and straight channel patterns, the braiding index will be 1, or slightly greater than 1 where there are channel islands.For anastomosing or anabranching channel patterns, the braiding index is calculated for the main branch of the river.Measurements are made for river reaches immediately upstream of the fold, across the fold axis (or its projection), and immediately downstream of the fold, to elucidate any changes in braiding index associated with the fold.
For the example of the River Dez incising across the Sardarabad Anticline, braiding index, BI = 1.0 (immediately upstream of fold); 1.2 (across fold axis); 1.2 (immediately downstream of fold), as shown in Figure 4.In this figure, thin yellow lines indicate the sub-division of each river reach into river cross-sections orthogonal to the valley axis which are 1 km apart.

General River Course Direction
Symbol: RCD Units: degrees (quoted to the nearest 5 • , as a compass bearing in degrees relative to true north, and as a bearing in degrees relative to the fold axis) Measurement location: River reaches immediately upstream of fold, across the fold axis (or its projection), and immediately downstream of fold The general river course direction is the general overall direction towards which the river flows for the length of a river reach [13].This can be gauged "by eye" by carefully viewing the remote sensing images and drawing a straight line of that orientation on the remote sensing image -the orientation of which will be similar to the river reach axes in the vicinity -and then measuring the bearing of that line to the nearest 5 • to avoid false precision.
For the example of the River Dez incising across the Sardarabad Anticline, general river course direction, RCD = 130 • (10 • to fold axis) (immediately upstream of fold); 230 • (70 • to fold axis) (across fold axis); 135 • (15 • to fold axis) (immediately downstream of fold), as shown in Figure 5.In this figure, the general river course direction is indicated by white lines with black arrowheads.C-RC is defined as the horizontal distance from the centre of the fold core measured along the fold axis, and along the projection of the fold axis, where appropriate), to the location where the river channel thalweg crosses the fold axis or its projection [13].This is most easily measured on fine scale geological maps (typically 1:100,000 or 1:50,000 scale geological maps, depending on availability) on which the surface lithology, structural geology (including the surface extent and anticlinal axis of each fold), and river channels are accurately shown.
The river crossing location is determined simply from where the fold axis (or its projection) intersects with the thalweg of the main river channel, as indicated on the fine scale geological map or on the remote sensing image.Where the main river channel has more than one intersection with the fold axis, as may be the case with a sinuous river, the intersection that is nearest to the fold core will be considered the river crossing location.The location of the centre of the fold "core" (the centre of the main part of the fold which emerged first on the ground surface) is considerably more difficult to determine, since the detailed developmental history of a fold is usually not known.For ease of measurement, the centre of the fold core should be located on the fold axis.For sub-surface folds with little or no surface topographic expression, known principally from oil and gas field locations and seismic surveys, the centre of the fold core should be interpreted as being midway along the approximate location of the fold axis on the ground surface (with particular consideration of the dip of sub-surface structures and stratigraphy).This interpretation can be modified in cases where the sub-surface structural geology is well known.For young, emerging folds the centre of the fold core can be interpreted with more confidence and will usually be coincident with the centre of the surface topographic expression of the fold.For older, emerged folds the location of the centre of the fold core is much less certain.It can generally be interpreted to be in the vicinity of the C-RC is defined as the horizontal distance from the centre of the fold core measured along the fold axis, and along the projection of the fold axis, where appropriate), to the location where the river channel thalweg crosses the fold axis or its projection [13].This is most easily measured on fine scale geological maps (typically 1:100,000 or 1:50,000 scale geological maps, depending on availability) on which the surface lithology, structural geology (including the surface extent and anticlinal axis of each fold), and river channels are accurately shown.
The river crossing location is determined simply from where the fold axis (or its projection) intersects with the thalweg of the main river channel, as indicated on the fine scale geological map or on the remote sensing image.Where the main river channel has more than one intersection with the fold axis, as may be the case with a sinuous river, the intersection that is nearest to the fold core will be considered the river crossing location.The location of the centre of the fold "core" (the centre of the main part of the fold which emerged first on the ground surface) is considerably more difficult to determine, since the detailed developmental history of a fold is usually not known.For ease of measurement, the centre of the fold core should be located on the fold axis.For sub-surface folds with little or no surface topographic expression, known principally from oil and gas field locations and seismic surveys, the centre of the fold core should be interpreted as being midway along the approximate location of the fold axis on the ground surface (with particular consideration of the dip of sub-surface structures and stratigraphy).This interpretation can be modified in cases where the sub-surface structural geology is well known.For young, emerging folds the centre of the fold core can be interpreted with more confidence and will usually be coincident with the centre of the surface topographic expression of the fold.For older, emerged folds the location of the centre of the fold core is much less certain.It can generally be interpreted to be in the vicinity of the structurally highest part of the present-day fold, which depending on the specific fold could be near its highest topographic expression, midway along the fold axis, or near to where it merges with an older, more developed fold [13,82].
For the example of the River Dez incising across the Sardarabad Anticline, distance from fold core to location of river crossing, C-RC = 1.3 km, as shown in Figure 6.For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, distance from fold core to location of river crossing, C-RC = 32.2km, also as shown in Figure 6.In this figure, the centre of the fold core is indicated by the black and yellow circle, the C-RC measurement along the fold axis to the River Dez crossing is indicated by the solid dark green line with two black arrowheads, and the C-RC measurement along the fold axis to the River Karun crossing is indicated by the dashed dark green line with two black arrowheads.In Figures 6-8 structurally highest part of the present-day fold, which depending on the specific fold could be near its highest topographic expression, midway along the fold axis, or near to where it merges with an older, more developed fold [13,82].
For the example of the River Dez incising across the Sardarabad Anticline, distance from fold core to location of river crossing, C-RC = 1.3 km, as shown in Figure 6.For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, distance from fold core to location of river crossing, C-RC = 32.2km, also as shown in Figure 6  3.1.8.Distance From Fold Core to River Basin Margin Symbol: C-BM Units: km (quoted to one decimal place, indicating positive or negative) Measurement location: Along the fold axis, from the centre of the fold core to the nearest margin of the drainage basin of the river interacting with the fold C-BM is defined as the horizontal distance from the centre of the fold core, measured along the fold axis (and along the projection of the fold axis, where appropriate), to the nearest margin of the drainage basin of the river interacting with the fold [13].The centre of the fold core is determined from fine scale geological maps, as described in Section 3.1.7.The drainage basin margins are demarcated from remote sensing images or topographical maps, by determining which river

Width of Geological Structure at Location of River Crossing
Symbol: Wgs Units: km (quoted to one decimal place) Measurement location: Where river channel crosses the fold axis (or its projection), measured orthogonal to the fold axis (or its projection) Wgs is defined as the maximum horizontal surface extent of the geological structure at the location where the river channel thalweg crosses the fold axis (or its projection), measured orthogonal to the fold axis (or its projection) [13].For sub-surface folds with little or no surface topographic expression, known principally from oil and gas field locations and seismic surveys, this measurement is necessarily approximate.For river incision across a fold, the measurement is made orthogonal to the interpreted fold axis, between the margins of the mapped oil or gas field.
For river diversion around a fold, the measurement is made orthogonal to the projection of the interpreted fold axis, between the projected margins of the nose of the mapped oil or gas field; a measurement which is highly subjective.For emerged folds with significant surface topographic expression, this measurement is much more certain.For river incision across a fold, the measurement is made orthogonal to the fold axis, between the surface extent of the fold limbs as determined from fine scale geological maps and fine scale topography.For river diversion around a fold, the measurement is made orthogonal to the projection of the fold axis, between the projected surface extent of the fold limbs of the nose of the fold; a measurement which is moderately subjective.
For the example of the River Dez incising across the Sardarabad Anticline, width of geological structure, Wgs = 4.3 km, and for the River Karun (Shuteyt) diverting around the Sardarabad Anticline, Wgs = 4.1 km, as shown on Figure 8.In this figure, the Wgs measurement for the River Dez crossing is indicated by the thick blue-grey line, and the Wgs measurement for the River Karun crossing is indicated by the thick red line.This characteristic is defined as the resistance of sediments and rocks in a fold to river erosion, a parameter which can be difficult to quantify.It depends upon a variety of characteristics

Distance From Fold Core to River Basin Margin
Symbol: C-BM Units: km (quoted to one decimal place, indicating positive or negative) Measurement location: Along the fold axis, from the centre of the fold core to the nearest margin of the drainage basin of the river interacting with the fold C-BM is defined as the horizontal distance from the centre of the fold core, measured along the fold axis (and along the projection of the fold axis, where appropriate), to the nearest margin of the drainage basin of the river interacting with the fold [13].The centre of the fold core is determined from fine scale geological maps, as described in Section 3.1.7.The drainage basin margins are demarcated from remote sensing images or topographical maps, by determining which river channels, wadis, lakes, streams and creeks are associated with each major river and by drawing a line midway between the extents of these.The zero point for measurements is at the centre of the fold core, with +ve values where the fold core is located within the drainage basin of the river interacting with the fold, and -ve values where the fold core is located outside of the drainage basin of the river interacting with the fold.
For the example of the River Dez incising across the Sardarabad Anticline, distance from fold core to river basin margin, C-BM = +3.8km, as shown in Figure 7.For the example of the River Karun (Shuteyt) diverting around the Sardarabad Anticline, distance from fold core to river basin margin, C-BM = −25.7 km, also as shown in Figure 7.In this figure, the centre of the fold core is indicated by the black and yellow circle, drainage basin margins are indicated by dashed blue lines, the C-BM measurement along the fold axis to the nearest River Dez basin margin is indicated by the solid dark purple line with one black arrowhead, and the C-BM measurement along the fold axis to the nearest River Karun basin margin is indicated by the dashed dark purple line with one black arrowhead.

Width of Geological Structure at Location of River Crossing
Symbol: Wgs Units: km (quoted to one decimal place) Measurement location: Where river channel crosses the fold axis (or its projection), measured orthogonal to the fold axis (or its projection) Wgs is defined as the maximum horizontal surface extent of the geological structure at the location where the river channel thalweg crosses the fold axis (or its projection), measured orthogonal to the fold axis (or its projection) [13].For sub-surface folds with little or no surface topographic expression, known principally from oil and gas field locations and seismic surveys, this measurement is necessarily approximate.For river incision across a fold, the measurement is made orthogonal to the interpreted fold axis, between the margins of the mapped oil or gas field.For river diversion around a fold, the measurement is made orthogonal to the projection of the interpreted fold axis, between the projected margins of the nose of the mapped oil or gas field; a measurement which is highly subjective.For emerged folds with significant surface topographic expression, this measurement is much more certain.For river incision across a fold, the measurement is made orthogonal to the fold axis, between the surface extent of the fold limbs as determined from fine scale geological maps and fine scale topography.For river diversion around a fold, the measurement is made orthogonal to the projection of the fold axis, between the projected surface extent of the fold limbs of the nose of the fold; a measurement which is moderately subjective.
For the example of the River Dez incising across the Sardarabad Anticline, width of geological structure, Wgs = 4.3 km, and for the River Karun (Shuteyt) diverting around the Sardarabad Anticline, Wgs = 4.1 km, as shown on Figure 8.In this figure, the Wgs measurement for the River Dez crossing is indicated by the thick blue-grey line, and the Wgs measurement for the River Karun crossing is indicated by the thick red line.

Estimate of Erosion Resistance of Surface Sediments/Rocks and Deeper Sediments/Rocks in Fold
Symbols: ERs (surface); ERd (deeper) No units (estimate quoted on a relative scale from 1 to 8) Measurement location: Where river channel crosses the fold axis (or its projection) This characteristic is defined as the resistance of sediments and rocks in a fold to river erosion, a parameter which can be difficult to quantify.It depends upon a variety of characteristics including structural geology, rock type, sediment type, strength of intact rock (especially rock compressive strength, rock tensile strength, and rock mass strength), resistance to weathering, jointing and fracturing (especially width, spacing, orientation, continuity, and infilling of joints), degree of movement of water through the rock mass, porosity, grain size, type and degree of cementation; as well as characteristics of the river, such as discharge, nature and frequency of floods, river sediment supply, suspended sediment concentration, and river bed roughness.Many of these characteristics are difficult to measure and their relative importance in determining the general erosion resistance of a fold is not fully known [22,31,[83][84][85][86][87][88].
Hence, for each case an estimate is made that is quoted as an integer on a scale, accompanied by a short description of the lithology or sedimentology, where known.The estimate of the erosion resistance of sediments and rocks is according to this scale: Very low (Unlithified floodplain sediments-predominantly sands) 2.
Low/Moderate (Mainly unlithified floodplain sediments-predominantly sands and silts; some quite poorly consolidated bedrock-such as Agha Jari Formation bedrock (quite poorly consolidated sandstones)-and other similar rocks-such as mudstones, evaporites and poorly consolidated limestones) 4.
Moderate (Mainly quite poorly consolidated bedrock-such as Agha Jari Formation bedrock (quite poorly consolidated sandstones)-and other similar rocks-such as mudstones, evaporites and poorly consolidated limestones; some unlithified floodplain sediments) 5.
Moderate/High (Mainly well consolidated bedrock-such as Bakhtyari Formation bedrock (very well consolidated conglomerates)-and other similar rocks-such as well consolidated limestones, marbles, sandstones and schists; some unlithified floodplain sediments and rocks of relatively low erosion resistance) 6.
High (Mainly well consolidated bedrock-such as Bakhtyari Formation bedrock (very well consolidated conglomerates)-and other similar rocks-such as well consolidated limestones, marbles, sandstones and schists) 7.
Very high (Very erosion resistant bedrock-basalts, gabbros, metasandstones and other very erosion resistant igneous and metamorphic rocks) 8.
Extremely high (Extremely erosion resistant bedrock-quartzite, cherts, granites, andesites, gneisses and other extremely erosion resistant igneous and metamorphic rocks) The position on this scale can be determined by careful interpretation of fine scale geological maps and remote sensing images, plus fieldwork and work on the properties of rocks and sediments, where available.The surface erosion resistance of the fold, ERs, is that of the surface lithology and sedimentology of the fold; especially that in the general vicinity of a river channel at the upstream location where the river first encounters the limb of the fold.The deeper erosion resistance of the fold, ERd, is that of the deeper lithology and sedimentology of the fold; especially that exposed in the general vicinity of an incising river channel at the location of the fold axis.With emerging folds, ERd may be unknown in some cases where the sub-surface geology is only poorly known.
For the example of the River Dez incising across the Sardarabad Anticline, ERs = 4, ERd = 5, and for the River Karun (Shuteyt) diverting around the Sardarabad Anticline, ERs = 4, ERd = 5, as shown on Figure 8.For the location of the River Dez crossing, ERs = Moderate (surface of unlithified floodplain sediments, with outcrops of Bakhtyari Formation bedrock (well consolidated conglomerates) and Agha Jari Formation bedrock (quite poorly cemented sandstones) at SW, W and E edges of floodplain), and ERd = Moderate/High (Bakhtyari Formation bedrock (well consolidated conglomerates) overlying Agha Jari Formation bedrock (quite poorly consolidated sandstones)).For the location of the River Karun (Shuteyt) crossing, ERs = Moderate (surface of unlithified floodplain sediments, with outcrops of Bakhtyari Formation bedrock (well cemented conglomerates) at SW and W edges of floodplain), and ERd = Moderate/High (assuming Bakhtyari Formation bedrock overlying Agha Jari Formation bedrock).

Measurement of Geomorphological Characteristics Nos. 11 to 13
For determining these three characteristics for the rivers Karun and Dez, the precise hydrological and topographical surveys used were those of the Dez Ab Engineering Company from 1997-2000, supplemented by geomorphological fieldwork [13].This data facilitated the measurement of channel water surface slopes.The superimposed database used included false-colour Landsat ETM+ satellite images (dated 2001), fine-scale geological maps, and CORONA satellite images (dated 1966 and 1968) which had been geo-referenced, orthorectified and enhanced in a unified database using ArcGIS ® software [85,89].This database facilitated the measurement of channel migrations with time, and facilitated easier measurement of characteristics associated with the fold core and fold axis.The data relating to vertical Earth surface movements were of various types, including radiocarbon dating of marine terrace sediments and Optically Stimulated Luminescence (OSL) dating of river terrace sediments, and enabled estimates of fold total uplift rate to be made [13,56].

Channel Water Surface Slope at Location of Fold Axis (or its Projection)
Symbol: s Units: m m −1 (quoted in standard form to nearest 1 × 10 −7 ) Measurement location: Where river channel crosses the fold axis (or its projection), measured for the river reach crossing the fold axis (or its projection) Channel water surface slope is determined by the equation, s = H c /L c where H c is change in channel water surface elevation (m) and L c is change in channel length (m) [90].Vertical accuracy should be at least of the order decimetres or better, ideally of the order of centimetres, or it will not be possible to discriminate the fine changes in slope associated with Earth surface movements, especially in lowland areas with very gentle slopes.
It is very difficult to determine sufficiently accurate channel water surface slopes from satellite remote sensing and fine scale geological maps, mainly due to their relatively poor vertical accuracy.Whilst Digital Elevation Models constructed from Shuttle Radar Topography Mission-30 m (SRTM-30m) data and Advanced Land Observing Satellite World 3D-30 m (AW3D30) data can be useful, they only have a vertical accuracy (Root Mean Square Errors) of about 5.7 m and 8.3 m, respectively [91].Greater accuracy may be obtainable from fine-scale topographical maps (especially maps of 1:25,000 scale and finer, with contour intervals of 5 m or less and frequent spot heights) [76], to determine river bank (and thus channel water surface) elevations.However, in general, more accurate additional data sources, such as precise hydrological or topographical surveys of the river, or hydrological and geomorphological fieldwork, will be needed to sufficient vertical accuracy.Even with precise surveys, there will be a variety of factors influencing water surface elevation measurements, such as local vegetation and obstructions, human modifications, levées, pools and riffles, eddies, and daily variations in discharge.These factors frequently induce appreciable errors, particularly in lowland areas with very gentle slopes.Thus, to reduce the influence of these errors, the channel water surface slope should be measured for the entirety of the river reach crossing the fold axis, or its projection.With all channel pattern types, the water surface of the main channel thalweg is used for the measurement.
For the example of the River Dez incising across the Sardarabad Anticline, channel water surface slope, s = 2.999 × 10 −4 mm −1 , and for the River Karun (Shuteyt) diverting around the Sardarabad Anticline, s = 3.5 × 10 −6 mm −1 .Precise hydrological and topographical surveys undertaken by the Dez Ab Engineering Company from 1997-2000 were used for these measurements.

Average Channel Migration Rate
Symbol: Rm Units: m yr −1 (quoted to three decimal places) Measurement location: River reaches immediately upstream of fold, across the fold axis (or its projection), and immediately downstream of fold Rm, the average channel migration rate over a specified period, can be defined by the equation Rm = (A/L c )/yr where A is total area of "migration polygons" drawn as shape files in a river reach between corresponding points of a river bank between remote sensing images of different dates (m 2 ); L c is channel length of reach (m), and yr is number of years between the remote sensing images [92].
To determine average channel migration rate, it is necessary to have access to high-resolution remote sensing images separated by a time interval of c. 20-40 years and Geographic Information System (GIS) software (such as ArcGIS ® ) to orthorectify and superimpose the two sets of remote sensing images.A time interval of c. 20-40 years should be long enough for significant channel migration to have taken place, though not so long that a channel may have migrated back to its original location.Where possible, one set of images should be high resolution aerial photographs or satellite images from the 1960 s or earlier, so that the time interval includes periods prior to major dam building and other major human impacts.
One of the river banks-the left bank when facing downstream -is manually digitised for each image set, the "migration polygons" created by their intersections are highlighted and saved as shape files, and the total area of these "migration polygons" for the river reach is calculated.For single-thread meandering and straight channel patterns, the left bank of the channel is used for the measurement.For multi-thread braided channel patterns, the left bank of the outermost braid channel is used for the measurement.For anastomosing or anabranching channel patterns, the left bank of the main active river channels is used, with any anabranches clearly separated by floodplain areas being considered as discrete channel-belts not included in the measurement.Using this value of A for the total area of the "migration" polygons, the channel length of the reach, L c , and the mean time interval in years, the average channel migration rate, Rm, can be calculated.
For the example of the River Dez incising across the Sardarabad Anticline, average channel migration rate, Rm = 11,129 m yr −1 (immediately upstream of fold); 1578 m yr −1 (across fold axis); 4502 m yr −1 (immediately downstream of fold), as shown in Figure 9

Estimate of Fold Total Uplift Rate
Symbol: TUR No units (estimate quoted on a relative scale from 0 to 8, roughly equivalent to ranges of rates of uplift in mm yr −1 ) Measurement location: At, or near, the fold crest The fold total uplift rate is defined as the rate at which a fold is rising above the surrounding region; that is, the single fold uplift rate less the sum of the regional subsidence rate and the sediment aggradation rate [29].Generally, it is estimated or measured at, or near to, the fold crest because, in most cases, that is the part of the fold undergoing the greatest uplift relative to the surrounding region [93].
Fold total uplift rate cannot be determined solely from remote sensing images, remote sensing data, topographical maps, and geological maps.Other data sources are needed, which may be precision topographic survey (recurrent surveys over several decades to determine vertical surface movements) e.g., [90,94], or precision GPS survey (recurrent measurements from GPS stations over several years to determine horizontal and vertical surface movements) e.g., [60,95,96].Alternatively, the data sources may be the measurement and dating of uplifted geomorphic markers, especially marine terraces and river terraces e.g., [51,56,59,97,98], the measurement and dating of archaeological structures, especially disused ancient canals e.g., [13,99], and the measurement and dating of structural geology, especially the development and erosion of fold growth strata e.g., [32,100].Where such data is available for a fold, either by direct measurement or by careful interpretation, the estimated fold total uplift rate can be quoted as an integer on this relative scale: 1.
Net subsidence (less than 0 mm yr −1 , the fold uplift rate is less than the sum of the regional subsidence rate and the sediment aggradation rate) 2.
Extremely high (more than 12.0 mm yr −1 ) For the example of the River Dez incising across the Sardarabad Anticline, TUR = 3, and for the River Karun (Shuteyt) diverting around the Sardarabad Anticline, TUR = 3.The TUR for the Sardarabad Anticline was estimated to be Low/Moderate (about 0.2-0.5 mm yr −1 ) because OSL dating of river terrace sediments indicated uplift of the back-limb of the Sardarabad Anticline at a rate of 0.23-0.29 mm yr −1 [13,20].

Results for the Rivers Karun and Dez
The results from applying the scheme to the River Karun and River Dez interacting with folds and emerging folds in lowland south-west Iran are given as a short description of the two rivers in Section 4.1, followed by tables of the results for the 13 geomorphological characteristics in Section 4.2.It is recommended that a similar format is used when applying the scheme to other major rivers in different parts of the world.The findings of Analysis of Variance (ANOVA) between river incision across a fold and river diversion around a fold, applied to the 13 geomorphological characteristics for the rivers Karun and Dez, are given in Section 4.3.

Short Description of River
River Karun (Iran) Length: 890 km Drainage basin area: 45,230 km 2 Mean annual water discharge: 575 m 3 s −1 (at Ahvaz in the Khuzestan Plains) Seasonality of discharge: Maximum in April (c.850 m 3 s −1 , or more than 2,000 m 3 s −1 prior to major dam construction); Minimum in October (c.280 m 3 s −1 ) (at Ahvaz in the Khuzestan Plains).
Short description of river course: Source in central/eastern Zagros on slopes of Zardeh Kuh, elevation c. 4200 m-very winding, roughly west course through Zagros Mountains, often in accordance with general NW-SE structural grain and folding -generally west course across Zagros foothills-generally south course from Gotvand onwards across Upper Khuzestan Plains, with bifurcation at Shushtar into River Shuteyt (larger branch) to the west and River Gargar to the east, which re-unite at Band-e Qir in vicinity of confluence with River Dez -generally south-west course from Ahvaz across Lower Khuzestan Plains-joins Tigris-Euphrates-Karun delta at Khorramshahr and fans out in south-east direction into Persian Gulf [13,35,54,102].

Geomorphological Characteristics
The geomorphology of the study area is complex and it was affected by tectonic processes over time.The results for the 13 geomorphological characteristics when applying the scheme to the River Karun and Dez interacting with folds in lowland south-west Iran are given in Table 2.This has been sub-divided into: Table 2 for the Turkalaki, Shushtar and Qal'eh Surkheh Anticlines and the River Karun; Table 3 for the Sardarabad, Qal'eh Surkheh and Kupal Anticlines and the River Karun and River Gargar; Table 4 for the Dezful Uplift and Sardarabad and Shahur Anticlines and the River Dez; Table 5 for the Ramin Oilfield, Ahvaz and Ab-e Teymur Oilfield Anticlines and the River Karun; and Table 6 for the Dorquain Oilfield Anticline and the River Karun.

Statistical Analysis of Geomorphological Characteristics
In tabulated form, the results for the geomorphological characteristics can be readily subjected to statistical analyses for investigating fold-river interactions.For the river Karun and Dez results given in Tables 2-6, Analysis of Variance (ANOVA) was applied between the categories of river incision across a fold and river diversion around a fold for each of the 13 geomorphological characteristics [13].The ANOVA findings are summarised in Table 7, in which F = Obtained F value (mean sums of squares due to between-group differences/mean sums of squares due to within-group differences), F crit = Critical F value needed to reject the null hypothesis, p-value = Level of significance of F value.In Table 7, bold text and yellow shading is used to highlight statistical significance, that is where p-value ≤ 0.05 (equivalent to a 5% significance level or a 95% confidence level or better) [74,121,122].

Discussion
The new scheme for investigating fold-river interactions was successfully applied to the major rivers Karun and Dez in lowland south-west Iran.It was found to be relatively easy to use in practice, with the notable exception of determining geomorphological characteristics Nos.11, 12 and 13.The acquisition and interpretation of precise hydrological and topographical surveys of the rivers, the creation of a superimposed database of satellite images and fine scale geological maps, subdivision of the Karun and Dez into river reaches, creation of river "migration polygons", and acquisition of data relating to fold uplift rates, was time-consuming.Hence, especially since data for these last 3 geomorphological characteristics may not be available for all major rivers, they can be considered as supplementary characteristics.

Significant Geomorphological Characteristics in Fold-River Interactions for the Rivers Karun and Dez in Lowland South-West Iran
As is frequently the case in geomorphology, the measurements of the 13 geomorphological characteristics will vary according to where and how the observer take the measurements.Thus it is necessary to follow the directions given with the descriptions of each of the 13 characteristics, so that the measurements are standardised.This is especially the case with the subdivision into river reaches and the determination of the location of the fold core, where a greater degree of subjectivity is involved.Hence, it may be useful to include error estimates with the measurements of the geomorphological characteristics.
Considering this, and the natural variability and complexity of major rivers, it might be expected that none of the 13 geomorphological characteristics would show statistically significant differences between the categories of river incision across a fold and river diversion around a fold.Nevertheless, the Analysis of Variance (ANOVA) findings in Table 7 for the rivers Karun and Dez interacting with folds in lowland south-west Iran, show three geomorphological characteristics with statistically significant differences at the 95% confidence level (p-value ≤ 0.05): channel-belt width at location of fold axis, floodplain width at location of fold axis, and distance from fold core to location of river crossing (geomorphological characteristics Nos. 2, 3 and 7).
Both channel-belt width and floodplain width at the location of the fold axis are significantly narrower for river incision across a fold compared with river diversion around a fold.In cases of river diversion, channel-belt width and floodplain width at the projection of the fold axis may have a wide range of values.By contrast, in cases of river incision, channel-belt width is always (100% of cases) less than 2.7 km, and floodplain width is generally (80% of cases) less than 5.7 km, at the location of the fold axis.A narrow channel-belt and a narrow floodplain at the location of the fold axis are indicative of a reduction in the lateral migration of the river at the fold axis to increase vertical incision of the river to keep pace with fold uplift.The general scenario is one of broader channel-belts and floodplains immediately upstream and downstream of the fold due to increased aggradation to maintain channel slopes across the fold, and narrow channel-belts and floodplains across the fold due to increased erosion and incision to keep pace with fold uplift [32,46,47].A narrow channel-belt is present in all cases of river incision, probably because a channel-belt is a relatively small feature that typically develops over time intervals of several decades or more [123].Indeed, an average channel-belt width of 2.7 km or less may be a threshold for the rivers Karun and Dez in the Khuzestan Plains that needs to be maintained if a major river is incise across a fold in the long-term [13].By contrast, a narrow floodplain is not present in all cases or river incision, probably because a floodplain is a significantly larger feature that typically develops over time intervals of centuries [63].Thus, in extensive, relatively flat areas, such as the Lower Khuzestan Plains, the streams and wetlands of the floodplain of a major river may extend far beyond the surface expression of a small fold and thus be unaffected by it.
It is not unexpected that the distance from the fold core to the location of the river crossing should discriminate between the two categories of fold-river interactions, since river incision occurs between the fold core and the fold nose, and river diversion occurs beyond the fold nose.Indeed, in all but one case of river incision the distance is less than 16 km (the exception of 43.6 km for the River Gargar incising across the Kupal Anticline is associated with pronounced human influences on the development of the River Gargar), whereas in all cases of river diversion the distance is greater than 22 kilometres.Interestingly, there is a strong tendency for a river to incise across the fold at locations near to that of the fold core (8.5 km or less in 80% of cases).Since the folds in the Khuzestan Plains are relatively young folds, this suggests that river incision across a fold at, or near to, the fold core is initiated at a very early stage in fold development, probably when the fold is initially emerging on the ground surface [13].
These findings can help to explain the seemingly paradoxical tendency of rivers to transect both young and old anticlines at or near to locations of their greatest structural and topographic relief.It can be considered that a fold initially emerges on the ground surface as a fold core, which in plan form may be an "oval", a "sausage", or another similar form, depending on the type of fold [33,82,[122][123][124][125].Where a major river initially encounters the fold as an emerging fold core, then the river may flow across the uplifting fold for sufficient time (at least several decades [123]) for the development of a narrow channel-belt; thus producing an incising river course across the fold in the vicinity of the fold core.As the fold grows vertically and laterally, depending on the size and nature of the river, the incising river course may be maintained and become "fixed" to produce a water gap in the fold in the vicinity of the subsequent structural culmination, or the river may be subsequently defeated to produce a wind gap and a diverted river course [124].By contrast, where a major river initially encounters a fold as a larger, emerged fold, the river may not flow across the uplifting fold for sufficient time for a narrow channel-belt to develop, due to repeated channel migration in response to lateral fold growth; thus producing a river course diverting around the fold nose [13,20].
Also, whilst other geomorphological characteristics are not discriminatory at the 95% confidence level, they do show some trends which support this model.River incision across a fold frequently has a general river course direction across the fold orthogonal to the fold axis, and river reaches across the fold axis have low channel sinuosities (generally < 1.4), steep channel water surface slopes (generally > 1 × 10 −4 mm −1 ), and low average channel migration rates (generally < 2 m yr −1 ).These trends might be related to the river initially encountering the fold as an emerging fold core, with reductions in lateral migration and increases in specific stream power as the river incises vertically in response to fold uplift.By contrast, river diversion around a fold frequently has a general river course direction upstream of the fold parallel to the fold axis and a course change of about 20 • -70 • to flow around the fold, and river reaches across the projection of the fold axis that have quite widely ranging channel sinuosities, gentle channel water surface slopes (generally < 1 × 10 −4 mm −1 ), and quite widely ranging average channel migration rates.These trends might be related to the river initially encountering the fold later in its development as a relatively large "obstacle", with a diverted river course, with frequent lateral migration in which there is only limited time for any increases in specific stream power to develop [13,20].

Significant Geomorphological Characteristics in Fold-River Interactions for Other Major Rivers
These observed changes apply for the major rivers Karun and Dez in lowland south-west Iran.To investigate whether similar or different changes apply with other major rivers and other folds, the scheme should now be applied to a variety of major rivers across the globe.For other fold-river interactions, it is highly likely that different changes will be found, and that other characteristics of river and fold geomorphology may discriminate between river incision across a fold and river diversion around a fold.
For instance, for the rivers Karun and Dez interacting with active folds in lowland south-west Iran, it was found that channel-belt width was a key discriminative characteristic, whereas channel width and channel water surface slope were not.By contrast, an investigation of two side-by-side upland rivers crossing rapidly uplifting folds (rates of uplift exceeding 10 mm yr −1 ) in the Himalayan foreland of central Nepal, found that both of the rivers exhibited a significant reduction in channel width across the zone of rock uplift [51,53,126].The smaller Bakeya River became steeper across the zone of rapid uplift, whereas the larger Bagmati River showed no significant profile steepening across the same zone [51,52].The research indicated that channel width acted as a key characteristic of river responses, and that if structural uplift should become sufficiently great, the channel width would reduce to less than a certain threshold width value to maintain an incising river course across a zone of uplift.Channel narrowing to enhance incision rates appeared to take precedence over other changes, such as channel steepening and reduced river profile concavity [27,52,53]; a scenario which has also been found with upland rivers elsewhere in the world.In central Taiwan, in response to increasing rates of differential uplift, upland rivers in studies were found to have progressively narrower channel widths until a channel width:depth ratio of about 10 was reached, after which they also steepened [58].In southern New Zealand, surveys of small upland channels indicated that 1 m-2 m of uplift resulted in a five-to ten-fold narrowing of river channels [25].Such findings enabled Amos and Burbank (2007) [25] to produce a conceptual model for a given river discharge, in which decreased channel width produced sufficient increased erosion to keep pace with uplift for small folds; whereas decreased channel width to a minimum value followed by subsequent channel steepening was needed to keep pace with uplift for larger folds [27].Hence, it is highly likely that there are significant differences between fold-river interactions in upland and lowland river catchments, with the geomorphological characteristics of channel width and channel water surface slope probably being more significant with upland rivers.This should be investigated by extending the database for the scheme to a variety of upland rivers.
Also, it has been hypothesised that the seemingly paradoxical tendency for the rivers Karun and Dez in lowland south-west Iran to transect anticlines near to locations of their greatest structural and topographic relief, is primarily due to the nature and timing of the initial fold-river interactions [13,20].However, there are other mechanisms that may account for this, which apply after the initial stages of fold development.It may arise by the drainage network being superimposed from above via a structurally conformable more easily eroded horizon [35,36].It may arise in areas where the crust is deforming plastically in response to regional compression, as a consequence of focussed rock uplift in response to significant differences between net erosion along major rivers and the surrounding regions [127], or in response to significant unloading of the crust by river erosion that amplifies the background deformation to produce a doubly plunging anticline with a river valley at its centre [128,129].Alternatively, with continued crustal shortening and thickening, it may arise with amplification of a regional slope that produces higher erosion rates in transverse catchments than in longitudinal catchments, and which creates a new organisation of the drainage system following the regional slope [130].It is likely that there will be notable differences in the relative significance of the geomorphological characteristics with each of these mechanisms, which should be investigated by extending the database for the scheme to a wide variety of rivers across the globe.

Conclusions
This study has introduced and demonstrated a new scheme using remote sensing for investigating fold-river interactions for major rivers.This scheme involved a short description of the river, climate, and structural geology, and 13 geomorphological characteristics.
The scheme was successfully applied to the major rivers Karun and Dez in lowland south-west Iran, using widely available satellite imagery and fine scale geological maps.It was relatively easy to use in practice, though geomorphological characteristics Nos.11, 12 and 13 involved additional data sources and additional processing, which was more difficult and time-consuming.Since the data needed for these last three geomorphological characteristics may not be available for all major rivers, they can be considered as supplementary characteristics.
For the major rivers Karun and Dez (mean annual water discharges 575 m 3 s −1 and 230 m 3 s −1 , respectively) interacting with folds in lowland south-west Iran, it was found that geomorphological characteristics Nos. 2, 3 and 7 (channel-belt width, floodplain width, and distance from fold core to location of river crossing) had statistically significant differences (p-value ≤ 0.05) between the categories of river incision across a fold and river diversion around a fold.These findings suggest that the nature and timing of initial fold-river interactions is important in determining whether a river incises across a fold or diverts around it, and that the formation and maintenance of a narrow channel-belt and a narrow floodplain are necessary for a major river to incise across a fold, with this incision frequently being in the vicinity of the fold core and subsequent structural culmination.
The scenario in the foreland basin tectonic setting of lowland south-west Iran involves major rivers (of which the Karun and Dez are the largest) interacting with relatively young, emerging, thrust-related folds, with gradual earth surface movements predominating due to lubricated décollements on evaporite layers.In this scenario, the new scheme was found to be useful and identified channel-belt width, floodplain width, and distance from fold core to river crossing as important characteristics in the interactions between the major rivers and the folds.The scheme should now be applied to a wide variety of major rivers across the globe, to determine its usefulness in other scenarios and improve our knowledge of fold-river interactions.By comparing the same parameters for different major rivers, a better understanding of fold-river interactions should be achieved.

Figure 1 .
Figure 1.The River Karun, River Dez, and other main rivers of Khuzestan province and its environs (Modified from Heyvaert et al., 2013) [63].Centred on 31°33'N 49°02'E.HM Huwayzah marshes SM Shadegan marshes International border Border of Khuzestan province.The Sardarabad Anticline is a 58 km long fold that is oriented roughly ESE-WNW and located to the north-west of the settlement of Band-e Qir

Figure 1 .
Figure 1.The River Karun, River Dez, and other main rivers of Khuzestan province and its environs (Modified from Heyvaert et al., 2013) [63].Centred on 31 • 33'N 49 • 02'E.HM Huwayzah marshes SM Shadegan marshes International border Border of Khuzestan province.The Sardarabad Anticline is a 58 km long fold that is oriented roughly ESE-WNW and located to the north-west of the settlement of Band-e Qir.

Figure 3 .
Figure 3.The measurement of w, cbw and fpw (False-colour Landsat image (2001) of the River Karun (Shuteyt branch) diverting around the Sardarabad Anticline, centred on c. 31°52′N 48°53′E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.3.1.2.Channel-Belt Width at Location of Fold axis (or its Projection) Symbol: cbw Units: km (quoted to three decimal places) Measurement location: Where river channel crosses the fold axis (or its projection)

Figure 3 .
Figure 3.The measurement of w, cbw and fpw (False-colour Landsat image (2001) of the River Karun (Shuteyt branch) diverting around the Sardarabad Anticline, centred on c. 31 • 52 N 48 • 53 E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

Figure 4 .
Figure 4.The measurement of Sc and BI (False-colour Landsat image (2001) of the River Dez incising across the Sardarabad Anticline, centred on c. 31 • 57 N 48 • 36 E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

Figure 5 .
Figure 5.The measurement of RCD (False-colour Landsat image (2001) of the River Dez incising across the Sardarabad Anticline, centred on c. 31°57′N 48°37′E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

3. 1 . 7 .
Distance From Fold Core to Location of River Crossing Symbol: C-RC Units: km (quoted to one decimal place) Measurement location: Along the fold axis, from the centre of the fold core to where the river channel crosses the fold axis (or its projection)

Figure 5 .
Figure 5.The measurement of RCD (False-colour Landsat image (2001) of the River Dez incising across the Sardarabad Anticline, centred on c. 31 • 57 N 48 • 37 E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.3.1.7.Distance From Fold Core to Location of River Crossing Symbol: C-RC Units: km (quoted to one decimal place) Measurement location: Along the fold axis, from the centre of the fold core to where the river channel crosses the fold axis (or its projection)C-RC is defined as the horizontal distance from the centre of the fold core measured along the fold axis, and along the projection of the fold axis, where appropriate), to the location where the river channel thalweg crosses the fold axis or its projection[13].This is most easily measured on fine scale geological maps (typically 1:100,000 or 1:50,000 scale geological maps, depending on availability) on which the surface lithology, structural geology (including the surface extent and anticlinal axis of each fold), and river channels are accurately shown.The river crossing location is determined simply from where the fold axis (or its projection) intersects with the thalweg of the main river channel, as indicated on the fine scale geological map or on the remote sensing image.Where the main river channel has more than one intersection with the fold axis, as may be the case with a sinuous river, the intersection that is nearest to the fold core will be considered the river crossing location.The location of the centre of the fold "core" (the centre of the main part of the fold which emerged first on the ground surface) is considerably more difficult to determine, since the detailed developmental history of a fold is usually not known.For ease of measurement, the centre of the fold core should be located on the fold axis.For sub-surface folds with little or no surface topographic expression, known principally from oil and gas field locations and seismic surveys, the centre of the fold core should be interpreted as being midway along the approximate location of the fold axis on the ground surface (with particular consideration of the dip of sub-surface structures and stratigraphy).This interpretation can be modified in cases where the sub-surface structural geology is well known.For young, emerging folds the centre of the fold core can be interpreted with more confidence and will usually be coincident with the centre of the surface topographic expression of the fold.For older, emerged folds the location of the centre of the fold core : White = Quaternary Alluvium and Recent Deposits (c. 1 Ma-Present; generally unconsolidated alluvial sands, muds, gravels, and marls).Yellow (Bk) = Bakhtyari Formation (Middle Pliocene to Pleistocene, c. 3 Ma-1 Ma; well-consolidated conglomerates, sandstones, and mudstones).Dark orange (Aj) = Agha Jari Formation (Middle Miocene to Middle Pliocene, c. 10 Ma-3 Ma; sandstones, marls, and mudstones).Light orange (Lbm) = Lahbari Member of Agha Jari Formation (Early to Middle Pliocene, c. 5.5 Ma-3 Ma; mudstones, marls, and sandstones) [69].Remote Sens. 2019, 11, x FOR PEER REVIEW 15 of 39 . In this figure, the centre of the fold core is indicated by the black and yellow circle, the C-RC measurement along the fold axis to the River Dez crossing is indicated by the solid dark green line with two black arrowheads, and the C-RC measurement along the fold axis to the River Karun crossing is indicated by the dashed dark green line with two black arrowheads.In Figures 6 to 8: White = Quaternary Alluvium and Recent Deposits (c. 1 Ma-Present; generally unconsolidated alluvial sands, muds, gravels, and marls).Yellow (Bk) = Bakhtyari Formation (Middle Pliocene to Pleistocene, c. 3 Ma-1 Ma; well-consolidated conglomerates, sandstones, and mudstones).Dark orange (Aj) = Agha Jari Formation (Middle Miocene to Middle Pliocene, c. 10 Ma-3 Ma; sandstones, marls, and mudstones).Light orange (Lbm) = Lahbari Member of Agha Jari Formation (Early to Middle Pliocene, c. 5.5 Ma-3 Ma; mudstones, marls, and sandstones) [69].

Figure 6 .
Figure 6.The measurement of C-RC (1:100,000 fine scale geological map (IOOC, 1969 [67]) of the rivers Dez and Karun (Shuteyt branch) interacting with the Sardarabad Anticline, centred on c. 31°54′N 48°42′E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

Figure 6 .
Figure 6.The measurement of C-RC (1:100,000 fine scale geological map (IOOC, 1969 [67]) of the rivers Dez and Karun (Shuteyt branch) interacting with the Sardarabad Anticline, centred on c. 31 • 54 N 48 • 42 E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

Figure 7 .
Figure 7.The measurement of C-BM (1:100,000 fine scale geological map (IOOC, 1969 [67]) of the rivers Dez and Karun (Shuteyt branch) interacting with the Sardarabad Anticline, centred on c. 31°53′N 48°43′E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

Figure 7 .
Figure 7.The measurement of C-BM (1:100,000 fine scale geological map (IOOC, 1969 [67]) of the rivers Dez and Karun (Shuteyt branch) interacting with the Sardarabad Anticline, centred on c. 31 • 53 N 48 • 43 E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

Figure 8 .
Figure 8.The measurement of Wgs, ERs and ERd (1:100,000 fine scale geological map (IOOC, 1969 [67]) of the rivers Dez and Karun (Shuteyt branch) interacting with the Sardarabad Anticline, centred on c. 31°53′N 48°43′E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.3.1.10.Estimate of Erosion Resistance of Surface Sediments/Rocks and Deeper Sediments/Rocks in Fold Symbols: ERs (surface); ERd (deeper) No units (estimate quoted on a relative scale from 1 to 8) Measurement location: Where river channel crosses the fold axis (or its projection)

Figure 8 .
Figure 8.The measurement of Wgs, ERs and ERd (1:100,000 fine scale geological map (IOOC, 1969 [67]) of the rivers Dez and Karun (Shuteyt branch) interacting with the Sardarabad Anticline, centred on c. 31 • 53 N 48 • 43 E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.
. In this figure, false-colour Landsat ETM+ images (28 July 2001 and 4 August 2001) make up the background; thin red lines indicate the location of the river channel banks on CORONA satellite images (23 September 1966 and 5 February 1968); yellow "migration polygons" indicate left bank channel migration inwards (or to the right) over a mean time interval of 34.2 years; green "migration polygons" indicate left bank channel migration outwards (or to the left) over 34.2 years.migration rate, Rm = 11,129 m yr (immediately upstream of fold); 1578 m yr (across fold axis); 4502 m yr −1 (immediately downstream of fold), as shown in Figure 9.In this figure, false-colour Landsat ETM+ images (28 July 2001 and 4 August 2001) make up the background; thin red lines indicate the location of the river channel banks on CORONA satellite images (23 September 1966 and 5 February 1968); yellow "migration polygons" indicate left bank channel migration inwards (or to the right) over a mean time interval of 34.2 years; green "migration polygons" indicate left bank channel migration outwards (or to the left) over 34.2 years.

Figure 9 .
Figure 9.The measurement of Rm (False-colour Landsat image (2001) of the River Dez incising across the Sardarabad Anticline, centred on c. 31 • 57 N 48 • 36 E).The axis of the anticline is shown as a red line with cross-bar, and straight-line river reaches are shown as thin green lines, with roughly orthogonal thin green lines demarcating successive reaches.

Table 2 .
Results for 13 geomorphological characteristics in fold-river interactions for the Turkalaki, Shushtar and Qal eh Surkheh Anticlines.

Table 3 .
Results for 13 geomorphological characteristics in fold-river interactions for the Sardarabad, Qal eh Surkheh and Kupal Anticlines.

Table 4 .
Results for 13 geomorphological characteristics in fold-river interactions for the Dezful Uplift and the Sardarabad and Shahur Anticlines.

Table 5 .
Results for 13 geomorphological characteristics in fold-river interactions for the Ramin Oilfield Anticline, Ahvaz Anticline and Ab-e Teymur Oilfield Anticline.

Table 6 .
Results for 13 geomorphological characteristics in fold-river interactions for the Dorquain Oilfield Anticline.

Table 7 .
Analysis of Variance (ANOVA) between river incision across a fold and river diversion around a fold, applied to 13 geomorphological characteristics in fold-river interactions for the rivers Karun and Dez in lowland south-west