Drivers of Long-Term Land-Use Pressure in the Merguellil Wadi, Tunisia, Using DPSIR Approach and Remote Sensing

: Increasing land use pressure is a primary force for degradation of agricultural areas. The drivers for these pressures are initiated by a series of interconnected processes. This study presents a novel methodology to analyze drivers of changing land use pressure and the effects on society and landscape. The focus was on characterizing these drivers and relate them to land use statistics obtained from geospatial data from the important semiarid Merguellil Wadi between 1976 and 2016. Cause-and-effect relationships between different drivers of land use change were analyzed using the DPSIR approach. Results show that during the 40-year period cultivated land increased and wetland areas decreased substantially. Drivers for change were pressure from economic development, cultivation practices, and hydro-agricultural techniques. This leads to stress on water and soil resulting in soil erosion, poverty increase, and rural exodus. We show that hydro-agricultural techniques adapted to the semiarid climate, allocation of land property rights, resource allocation, and improved marketing of agricultural products can help rural residents to diversify their economy, and thus better preserve the fragile semiarid landscape. Results of this study can be used to ensure sustainable management of water and soil resources in areas with similar climate and socio-economic conditions.


Introduction
Land use has significant environmental and socio-economic impacts. Information on land use changes is among other things necessary for questions regarding planning and protection of environment for a sustainable management of natural resources [1]. The dynamics of land use are to a great extent influenced by human and natural processes, and results from complex interactions between social, economic, and biophysical developments [2,3]. These processes operate at various temporal and spatial scales and are significantly affected by agricultural development, technological advances, and population pressure [1,4,5].
The study of land use dynamics, as a function of several physical and socio-economic factors or drivers of land use change, has attracted much attention from many researchers during the latest decades [6]. Human-induced land use changes reflect various policies of land use development and management. Although, these changes may be economically beneficial for humans, the stability of the natural environment may be jeopardized (e.g., [7]). Understanding of underlying drivers for land use changes becomes important to mitigate

Study Area
The Merguellil Wadi Basin is located in central Tunisia. The Merguellil Wadi is the second most important non-perennial river in central Tunisia. It covers an area of 1183 km 2 . During the Roman period the area was intensely farmed and today irrigation is used to produce mainly olive and apricot crops. The El Haouareb dam retains all water from the Merguellil Wadi (Figure 1) [24]. The shape of the watershed is almost triangular with a downstream base (Figure 1). The bed of the watercourse is classified as very sensitive to erosion [25]. Various hydraulic and soil structures have been built in the upstream of Merguellil Wadi to protect against floods and supply agricultural areas with irrigation. These structures have significantly altered surface runoff feeding the wadi. To prevent El Haouareb dam from silting, the upstream area is equipped with more than 200 km 2 of soil benches, lakes, and dams (38 lakes and five earth dams) to collect more than one million cubic meters of water. The altitude varies between 200 and 1200 m with an average of about 500 m (Figure 1). The basin is located in a semiarid climate, characterized by a very high spatial-temporal variability of precipitation [26]. Annual rainfall varies between 220 mm downstream and 420 mm in the upstream. Occurring rainfall is often intense, especially in spring and autumn, which cause flooding of the wadis [27]. Intermittent flow in the Merguellil Wadi supplies the El Haouareb reservoir with water. The reservoir, in turn, provides water to the downstream Kairouan plain and recharges the groundwater table here. The groundwater is used for irrigation and drinking water supply.

Study Area
The Merguellil Wadi Basin is located in central Tunisia. The Merguellil Wadi is the second most important non-perennial river in central Tunisia. It covers an area of 1183 km 2 . During the Roman period the area was intensely farmed and today irrigation is used to produce mainly olive and apricot crops. The El Haouareb dam retains all water from the Merguellil Wadi ( Figure 1) [24]. The shape of the watershed is almost triangular with a downstream base (Figure 1). The bed of the watercourse is classified as very sensitive to erosion [25]. Various hydraulic and soil structures have been built in the upstream of Merguellil Wadi to protect against floods and supply agricultural areas with irrigation. These structures have significantly altered surface runoff feeding the wadi. To prevent El Haouareb dam from silting, the upstream area is equipped with more than 200 km² of soil benches, lakes, and dams (38 lakes and five earth dams) to collect more than one million cubic meters of water. The altitude varies between 200 and 1200 m with an average of about 500 m ( Figure 1). The basin is located in a semiarid climate, characterized by a very high spatial-temporal variability of precipitation [26]. Annual rainfall varies between 220 mm downstream and 420 mm in the upstream. Occurring rainfall is often intense, especially in spring and autumn, which cause flooding of the wadis [27]. Intermittent flow in the Merguellil Wadi supplies the El Haouareb reservoir with water. The reservoir, in turn, provides water to the downstream Kairouan plain and recharges the groundwater table here. The groundwater is used for irrigation and drinking water supply. The basin annual average temperature is 19.2 °C with a minimum of 10.7 °C in January and a maximum of 28.6 °C in August and the potential evapotranspiration is about 1600 mm/year [28]. The area has a strong rural character with about 85% rural and 15% urban residents and these live mainly from livestock farming and rainfed and irrigated agriculture [21]. The population density varies between 50 and 100 inhabitants per km² with an average annual population growth of about 1.5%. About half of the watershed surface is occupied by annual crops (wheat and barley) and arboriculture (olive, apricot, The basin annual average temperature is 19.2 • C with a minimum of 10.7 • C in January and a maximum of 28.6 • C in August and the potential evapotranspiration is about 1600 mm/year [28]. The area has a strong rural character with about 85% rural and 15% urban residents and these live mainly from livestock farming and rainfed and irrigated agriculture [21]. The population density varies between 50 and 100 inhabitants per km 2 with an average annual population growth of about 1.5%. About half of the watershed surface is occupied by annual crops (wheat and barley) and arboriculture (olive, apricot, and almond trees). Since about 2000, the central watershed has experienced an unprecedented increase in fruit growing. Irrigated agriculture is the economic engine of the area. The area's agricultural products are oriented towards the regional and national market [23]. Selective changes in land use, coupled with certain poor agricultural practices (over-tillage, ploughing against the slope, poor control of irrigation and spate) have resulted in increased erosion. Farmers have to a great extent moved from multi-cropping to intensive monoculture [23]. This has deteriorated the land further and, thus, the degradation of soil fertility continues. As a result, the loss of agricultural land has had large-scale impacts on food security of the area, health, and quality of the soil, which in the long-term will hamper its productivity.

Data Used
We used Landsat 1 MSS (Multi-Spectral Sensor) acquired on 21 July 1976, Landsat 5 TM (Thematic Mapper) taken on 19 July 1996, and Landsat 8 OLI-TIRS (Operational Land Imager-Thermal Infrared Sensor) taken on 21 July 2016 with a spatial resolution of 30 m. July coincides with the start of the dry season and the end of the rainy season and thus, images were clear of clouds and haze. This also made it possible to avoid the effect of annual vegetation and to better understand the behavior and distribution of the plant cover. Landsat images were downloaded free of charge from the USGS website (https://earthexplorer.usgs.gov/, access date: 3 March 2020) to produce land use maps and identify the different classes of land use. A digital elevation model (DEM) was derived from the SRTM Shuttle Radar Topography Mission data (SRTM), from 1 arc-second for global coverage (30 m) and downloaded via the USGS website (https://earthexplorer.usgs.gov/, access date: 30 November 2019) ( Figure 1).

Land Use-Land Cover Classification
The land use/land cover classification for the three years (1976, 1996, and 2016) was based on pixels supervised with the Support Vector Machine (SVM) whose algorithm is integrated with the OTB tool (Orfeo ToolBox) of QGIS (Quantum GIS) ( Figure 2). The adopted land use classes in this study are based on the classification used by the Food and Agriculture Organization of the United Nations (FAO) [29]. The land use/land cover classes were grouped into seven main categories depending on soil conditions and nature of land cover in the study area ( Table 1). The classification was performed for each acquisition date using reference data and the normalized vegetation index. Table 1. Different classes adopted for the diachronic study of land use degradation [29].

Land Use/Land Cover Class Description
Cultivated area Herbaceous crops; woody crops; mixed herbaceous and woody crops.

Forest
Tree plants with a cover of 10% or more. Other types of plants (shrubs and/or grasses) may be present, even at a higher density than trees.

Arboriculture
Woody plants (trees and/or shrubs) may be present assuming that their cover is less than 10%.

Rangeland
Natural herbaceous plants (grasslands, steppes) with a cover of 10% or more, independent of different human and/or animal activities, such as grazing.
Bare land Natural abiotic surfaces (bare soil, sand, rocks, etc.) where natural vegetation is absent or almost absent (cover less than 2%).

Accuracy Assessment
Classification results were evaluated using the QGIS software to verify the supervised classification with field reality using randomly selected reference sampling points [1]. Benchmark data for each class were in the form of classification samples for training and validation. These samples were selected via Google Earth and the land use/land cover map of the Merguellil Wadi derived from an inventory (1995-2010) published in 2010 by the Directorate General of Forests [30]. The total number of benchmarks collected was 20,549 pixels, of which 15,807 were for classification and 4742 for accuracy rating. The total number of pixels represents approximately 0.77% of the pixels for the entire study area. Precision measures, such as total accuracy, kappa coefficient, user (U), and producer (P) accuracy were calculated as percentages. Thus, an error matrix of the land use classification was generated.

Household Surveys
Household surveys were conducted to understand the changing land use situation and gain insight of effects on farming using the DPSIR approach. A total of 116 household surveys were conducted during the period from November to December 2017 [31]. The population targeted by the survey was all individuals (farmers, herders, traders, etc.). The survey included socio-economic questions on the farming life and farmers' use of the natural resources. The aim was to analyze the functioning of the farm and to identify the factors driving the land use change. Interviews were semi-structured together with participating observations [31]. The surveys were designed to allow expression of farmers´

Accuracy Assessment
Classification results were evaluated using the QGIS software to verify the supervised classification with field reality using randomly selected reference sampling points [1]. Benchmark data for each class were in the form of classification samples for training and validation. These samples were selected via Google Earth and the land use/land cover map of the Merguellil Wadi derived from an inventory (1995-2010) published in 2010 by the Directorate General of Forests [30]. The total number of benchmarks collected was 20,549 pixels, of which 15,807 were for classification and 4742 for accuracy rating. The total number of pixels represents approximately 0.77% of the pixels for the entire study area. Precision measures, such as total accuracy, kappa coefficient, user (U), and producer (P) accuracy were calculated as percentages. Thus, an error matrix of the land use classification was generated.

Household Surveys
Household surveys were conducted to understand the changing land use situation and gain insight of effects on farming using the DPSIR approach. A total of 116 household surveys were conducted during the period from November to December 2017 [31]. The population targeted by the survey was all individuals (farmers, herders, traders, etc.). The survey included socio-economic questions on the farming life and farmers' use of the natural resources. The aim was to analyze the functioning of the farm and to identify the factors driving the land use change. Interviews were semi-structured together with participating observations [31]. The surveys were designed to allow expression of farmers views regarding long-term changes. However, people s opinions on change are often related to changes in the short-term and there was no real way to collect data on long-term changes in farmers' opinions.
The El Alaa region with about 32,000 inhabitants was selected as a representative area located at the foothills of the upstream Merguellil Wadi Basin at an altitude of 463 m ( Figure 1). The area was selected after a comprehensive pre-study of the socio-economic conditions of this area [31]. The region of El Alaa provides a link between the forest environment and other landscape units forming the Merguellil Basin. In addition, the socioeconomic problems experienced by the villagers in the management of natural resources in connection with exploitation methods are representative for the larger basin area. The El Alaa region constitutes a delegation (mutamadiyah) and a second level administrative division of Tunisia between governorates and sectors (imadats). There are 24 governorates in Tunisia that are divided into 264 delegations. The El Alaa inhabitants are mostly rural, living mainly from agriculture. Most other inhabitants have left the delegation to work elsewhere in Tunisia or Libya. According to the local authorities, 80% of the population are unemployed. Among the unemployed people, 1700 have higher education degrees and have generally been unemployed for nearly 10 years. More than 85% of the rural area are agro-pastoral [21]. Despite difficult conditions and the extension of the arboriculture, the cereal farming continues to be an important resource [31]. Most of the water is pumped from the downstream groundwater table through a high density of wells [23]. The basin has more than 5000 wells while there were only about a hundred in the 1960s [21]. In the upstream basin, farmers practice dry farming and arboriculture of olive trees (71% of plantations), almond trees (12%), and cereals dominate.

Indicator-Based Approach
Land use dynamics can impact land conditions through pressures from various driving forces that cause changes. The satellite images were used to determine the state of land use changes in the area for the different time periods. These only explain "how" these changes occurred and the "why" requires careful study of various variables or indicators of this change. The DPSIR approach establishes the relationship between "how" and "why". It presents a cause-and-effect framework of driving force-pressure-state-impact-response. Each parameter is impacted by the precedent that ultimately affects political decisions at the highest level. The adaptation of DPSIR in the assessment of land use dynamics is summarized in Figure 3. views regarding long-term changes. However, people´s opinions on change are often related to changes in the short-term and there was no real way to collect data on long-term changes in farmers' opinions. The El Alaa region with about 32,000 inhabitants was selected as a representative area located at the foothills of the upstream Merguellil Wadi Basin at an altitude of 463 m ( Figure 1). The area was selected after a comprehensive pre-study of the socio-economic conditions of this area [31]. The region of El Alaa provides a link between the forest environment and other landscape units forming the Merguellil Basin. In addition, the socioeconomic problems experienced by the villagers in the management of natural resources in connection with exploitation methods are representative for the larger basin area. The El Alaa region constitutes a delegation (mutamadiyah) and a second level administrative division of Tunisia between governorates and sectors (imadats). There are 24 governorates in Tunisia that are divided into 264 delegations. The El Alaa inhabitants are mostly rural, living mainly from agriculture. Most other inhabitants have left the delegation to work elsewhere in Tunisia or Libya. According to the local authorities, 80% of the population are unemployed. Among the unemployed people, 1700 have higher education degrees and have generally been unemployed for nearly 10 years. More than 85% of the rural area are agro-pastoral [21]. Despite difficult conditions and the extension of the arboriculture, the cereal farming continues to be an important resource [31]. Most of the water is pumped from the downstream groundwater table through a high density of wells [23]. The basin has more than 5000 wells while there were only about a hundred in the 1960s [21]. In the upstream basin, farmers practice dry farming and arboriculture of olive trees (71% of plantations), almond trees (12%), and cereals dominate.

Indicator-Based Approach
Land use dynamics can impact land conditions through pressures from various driving forces that cause changes. The satellite images were used to determine the state of land use changes in the area for the different time periods. These only explain "how" these changes occurred and the "why" requires careful study of various variables or indicators of this change. The DPSIR approach establishes the relationship between "how" and "why". It presents a cause-and-effect framework of driving force-pressure-state-impactresponse. Each parameter is impacted by the precedent that ultimately affects political decisions at the highest level. The adaptation of DPSIR in the assessment of land use dynamics is summarized in Figure 3.  The DPSIR indicators were assessed using the probability of selection and based on the 5-point Likert scale (level 1-very low; level 5-very high). The surveys were constructed and processed according to the probability of occurrence, mean, standard deviation, and the consensus for the 5 points of the Likert scale. The X i denotes the degree of satisfaction of the arguments according to category i of the Likert scale and P i the probability of occurrence of X i . Thus, the average of the overall agreement for the 5-point Likert scale is defined by: The mean values can be classified according to five levels; very low (1.00-1.49); low (1.50-2.49); medium (2.50-3.49); high (3.50-4.49); very high (4.50-5.00) [32,33].
The standard deviation for the 5-point Likert scale is defined by: The consensus is defined as the agreement of an argument between individuals in a sample of the group [33,34]. It is defined according to the 5-point Likert scale by: The Pearson product moment correlation aims to determine the relationships between land use/land cover dynamics and the DPSIR parameters. This tests the linear association between two quantitative variables X and Y [35]: where R is the correlation coefficient; X corresponds to the land use/land cover type; Y is the dependent variable, which corresponds to the probability of occurrence for the DPSIR parameters as perceived by households, and N is the number of indicators for each DPSIR component. The coefficient of determination (R 2 ) is calculated to test the correlation strength [36].  The DPSIR indicators were assessed using the probability of selection and based on the 5-point Likert scale (level 1-very low; level 5-very high). The surveys were constructed and processed according to the probability of occurrence, mean, standard deviation, and the consensus for the 5 points of the Likert scale. The Xi denotes the degree of satisfaction of the arguments according to category i of the Likert scale and Pi the probability of occurrence of Xi. Thus, the average of the overall agreement for the 5-point Likert scale is defined by:
The standard deviation for the 5-point Likert scale is defined by: The consensus is defined as the agreement of an argument between individuals in a sample of the group [33,34]. It is defined according to the 5-point Likert scale by: The Pearson product moment correlation aims to determine the relationships between land use/land cover dynamics and the DPSIR parameters. This tests the linear association between two quantitative variables X and Y [35]: where R is the correlation coefficient; X corresponds to the land use/land cover type; Y is the dependent variable, which corresponds to the probability of occurrence for the DPSIR parameters as perceived by households, and N is the number of indicators for each DPSIR component. The coefficient of determination (R²) is calculated to test the correlation strength [36].    Table 2 shows the confusion matrix for the three investigated years. The confusion matrix is calculated by comparing land covers derived from the Landsat images against ground truth land use data. Each column of the confusion matrix represents a ground truth class that corresponds to the image's labeling of the ground truth pixels [37]. Evaluation of the supervised classification is illustrated in Table 3. The overall precision was 79.3% for 1976, 79.5% for 1996, and 81.1% for 2016. The rate of the overall accuracy precision is related to the resolution of the satellite images used. The kappa coefficient for 1976, 1996, and 2016 was 0.73-0.75 (Table 3). The classification of cultivated area yielded high producer's and user's accuracies. However, confusions were mainly observed between cultivated area and rangeland (Table 3), indicating the difficulty of optical data in separating land covers of similar spectral signatures. The resulting land use classification of Merguellil Wadi Basin for 1976, 1996, and 2016 is shown in Figure 5 and Table 4. Therefore, the results can be said to have met the requirement for precision and it is possible to use them for an in-depth analysis of changes during the 40-year study.       Table 5 shows that the basin area has gone through substantial changes in land use during the 40-year study period. Especially, wetland areas decreased and cultivated areas increased. During the last twenty years, arboriculture, and forestry suffered a significant decrease (Table 5).

Adaptation of the DSPIR Approach in the Merguellil Wadi Basin
The DPSIR methodology was used to link land cover maps to the indicator-based approach. Based on the responses of farming households in the El Alaa region, regarding their perception of the different factors concerning land use change, a DSPIR model was developed ( Figure 6). The figure lists all responses from farm households in the questionnaire survey, but also includes the factors mentioned in the group discussions to encompass the entire Merguellil Wadi Basin.   Table 5 shows that the basin area has gone through substantial changes in land use during the 40-year study period. Especially, wetland areas decreased and cultivated areas increased. During the last twenty years, arboriculture, and forestry suffered a significant decrease (Table 5).

Adaptation of the DSPIR Approach in the Merguellil Wadi Basin
The DPSIR methodology was used to link land cover maps to the indicator-based approach. Based on the responses of farming households in the El Alaa region, regarding their perception of the different factors concerning land use change, a DSPIR model was developed ( Figure 6). The figure lists all responses from farm households in the questionnaire survey, but also includes the factors mentioned in the group discussions to encompass the entire Merguellil Wadi Basin.   Table 6 shows driving forces of land use/land cover changes as perceived through the 116 household surveys. The climate change appeared to be one of the main factors responsible for land use change (CnS = 0.68). The cultural practices implemented by farmers clearly show that drought presents a major climatic risk. Changes in the seasonal distribution of rainfall threaten the sustainability of existing production systems [31]. The plots of the El Ala region are cultivated using rain-fed systems (97.7% of the agricultural area). The most important crop distribution (51%) is a mixture between cereals and arboriculture of which olive and almond trees are the main crops due to their water stress resistance [31]. The population growth received a large consensus (CnS = 0.6) ( Table 6). Indeed, the number of inhabitants is about 29,900 people. The average size of a farming family in the study area is about 6.6 people, with a maximum of 12 people. The structure of land organization represents a main problem. The creation of agricultural plots and the division of land represent high mean driving force of 4.65 and 3.03, respectively. About half (48%) of farmers does not have ownership of the cultivated land. This adds to the fragmentation of land that is an obvious and important issue. In fact, the average number of plots per farmer is 6 with a minimum of 1 and a maximum of 20. The fragmentation of agricultural land (CnS = 0.34) and the creation of agricultural plots (CnS = 0.62) present two important driving forces. The average number of plots is 4.75 with a maximum of 22 plots. The total agricultural area is of the order of 1891.5 ha of which 1383.5 ha are for agriculture, with an average of 9.6 ha, a minimum of 1.5 ha, and a maximum of up to 200 ha per plot.

Pressures Due to Land Use/Land Cover Change
The pressures exerted due to the changing land use as perceived by agricultural households are mainly overexploitation of groundwater (CnS = 0.54) by increase in the number of illegally drilled wells, demand for more agricultural land, theft of livestock, degradation of rangelands, agroforestry, overexploitation of forest resources, and increased demand for forest products (CnS = 0.58) ( Table 7).

States Perceived Due to Land Use/Land Cover Change
In most of the region, the condition currently observed due to change in land use by households is change in forest cover (CnS = 0.76). Indeed, the survey region is characterized by traces of Aleppo pine and rosemary with diversified plantations (red juniper, etc.). In addition, water and soil resources are degraded and threatened by water erosion (CnS = 0.70). During field visits, we noticed the presence of expanding gullies and rock outcrops (Table 8).

Impacts of Land Use/Land Cover Change
The environmental and economic impacts reported by farmers are the increase in migration from rural to urban areas, work in the field of masonry, tourism and trade, and craft activities (margoum tapestry, sewing), declining land productivity, and increasing consumption of natural resources (Table 9). Table 9. Impacts of land use/land cover changes perceived by households.

Responses of Land Use/Land Cover Change
The response of farmers to changes in land use is development of diversification [31]. As a result, the population is turning towards diversification such as arboriculture (olive and almond trees), cereals, breeding (poultry, sheep, cattle, and goats) and the promotion of cactus crops (Table 10).  Table 11 illustrates the strength of the relationships between each land use land cover type and the probability of occurrence of DPSIR parameters as perceived by households. The determination coefficient is on average equal to 0.76. This shows that there is a strong correlation between the land cover class and the different components of the DPSIR approach.

Land Use Change
The land use change between 1970 and 2016 showed a continuous increase in cultivated areas. The share of crops occupied 6.3, 18.4, and 19.3% in 1976, 1996, and 2016, respectively (Table 3). Thus, over a period of 40 years between 1976 and 2016, the crop cover almost tripled, i.e., from 6.3 to 19.3% (Table 3). This is due to the conversion of denuded areas into rangeland and arboriculture. The construction of the El Houarreb dam in 1989 has helped to develop irrigation. These results are consistent with the description of the state of the soil cover in the Forest Investment Program in 2016 for the second national forest and pastoral inventory (IFPN) [38]. The inventory documented a gradual extension of crops, particularly for the practice of cereals and arboriculture. The increase of crops serves to improve selfsufficiency for certain products (vegetables, fruits, and cereals, including barley) and the provision of agricultural by-products to feed the livestock. Barley remains an important crop in the production system integrating sheep and cattle farming in the study region. In addition, the increasing cropping area is linked to the development of irrigation with an increase in drip irrigation between 1999 and 2005 [39]. Bare soil and rangeland have also experienced an important change (Table 5). Overgrazing and overexploitation of rangelands are explained by the increase in number sheep and goats and illegal logging of shrubs and trees, commonly practiced by the poorest agro-pastoral population [38]. A major land use change is the decreasing wetlands. These include springs, hill reservoirs, and tributaries of Merguellil Wadi (Table 5). This is explained by the variability of the rainfall regime and the increase in temperature during the study period [26] as well as soil erosion that causes siltation of the hill reservoirs and the El Houarreb dam. [40].

Drivers Linked to Land Use Change
Land use change is the result of a combined association between demographic, socioeconomic, biophysical, and institutional processes. Population growth has led to higher demand for natural resources, resulting in changes in land cover over time [41]. Indeed, according to the general population and housing census in 2014, the total population in the governorate of Kairouan is about 187,000 people with a growth rate of 1.8% [42]. According to forecasts by the National Statistics Institute, this population will more than double until 2030 [43]. In addition, the rugged mountainous topography of the Merguellil Basin makes agricultural production difficult resulting in large pressure on existing farmland.
Another important driver of land use change is the fragmentation of agricultural land. Indeed, the plots are small areas that usually do not exceed 10 ha [44]. This has encouraged farmers to increase their agricultural production, to use the most available space and, therefore, to remove ecological zones, namely hedges, and banks of the river. In addition, agricultural land use has been developed without regard to the capacity of the land. For example, livestock and animal production are not well integrated with cereal production.
Another driver of land use degradation is the insufficient/inefficient use of animal manure, inappropriate agricultural techniques, including among other things the repeated cycles of barley, wheat, and fallow in crop rotations. Inappropriate crops and mechanization (disc ploughs are not good for the soil structure). In addition, one of the main factors responsible for land use change is the impact of climate change. This factor is crucial for all of Tunisia. Indeed, recent trends show that Tunisia is becoming warmer and drier with intense extreme events (droughts and heavy rains). In fact, global warming will further disrupt rainfall patterns, as well as increase temperature, which will have even more disastrous consequences for water erosion and land management [28,45].

Pressures Linked to Land Use Dynamics
The drivers of land use change are exerting pressure on the sustainability of water and soil. In the Merguellil Basin, geomorphology has conditioned creation of the douars. These constitute groups of dwellings linked according to family relationships based on a common ancestry in the paternal line. The douars are grouped around water sources. Agricultural activities are characterized by diversification and extension of irrigated land; namely arboriculture and a frequent practice of market gardening especially in the center of the basin. This is due to the overexploitation of groundwater by drilling and pumping [26]. In addition, unregulated livestock grazing on public pastures and woodlands eliminates any possibility of sustainable management. There are at present no control mechanism for the number and type of animals and their grazing time on these pastures [21,23].

Status Linked to Land Use Dynamics
The degradation of the agricultural landscape is a visible and tangible phenomenon that has a direct impact on the livelihoods of the douars in the study area [23]. Indeed, soil loss presents the main type of land degradation since the most fertile upper part of the soil is lost, leading to exposure of the less fertile and more erosive part of the soil. For shallow soils, water erosion can strip the soil down to bedrock. Therefore, erosion leads to a total loss of the productive potential of agricultural land. The structure of soils can be severely damaged, and they usually become more acidic and much less fertile or non-productive [46].

Impacts Linked to the Dynamics of Land Use
Several studies have documented the impact of land use dynamics on the biodiversity (e.g., [47][48][49]). In the study area, the farmers are extremely poor due to the nature of dry farming and the use of land often less than an area of 5 ha. Income from agriculture is often insufficient to support the family. In addition, apart from the periods of sowing and harvesting, small farms do not generate labor. This, favors rural exodus and mountainous areas, traditionally the most popular for farming, olive picking, and forestry activities, have become less attractive as compared to plains and towns. In addition, the improved urban infrastructure has increased the differences between towns and the countryside and accelerated the rural exodus [23,50].

Responses Linked to Land Use Dynamics
Responses are measures to be taken by the government to mitigate the negative impacts of land use change. Tunisia has implemented different forms of adaptation strategies in view of environmental change and climate variability [51]. This strategy was developed following a period of severe drought between 1999 and 2001 that became the national strategy for the development and sustainable management of forests and rangelands (2015-2024). As well, the authorities have implemented the integrated rural agricultural development project (PADRI), the mountainous areas development policy since 1987, and the MERGUSIE project [52]. Water and soil conservation developments in the Merguellil Basin are divided into two categories, land use planning and hydrographic network development. The objective is to help increase water and soil reserves and to protect dams and hill lakes located downstream against siltation. They are essentially in the center of the basin, covering a total area of 285.2 km 2 [52].
The development of watercourses in the study area consists of stabilizing structures for riverbanks and longitudinal sections, lakes, and hill-side dams. In total, there are 44 dams and hill-side lakes in the Merguellil Basin [52]. More than half of these lakes (24 hillside lakes) are entirely for agricultural use, 9 lakes are used to protect dams against siltation, 2 to supply groundwater, and 13 for various uses.
The response of the authorities is further to encourage farmers to transform and sell their products at national and regional markets. Cereals and olive oil are collected and sold by specialized offices and organizations, namely the Cereals Office from 1992 and the National Oil Office from 1993.
However, despite continuous State interventions and significant investment, farmers still do not consider the problems to be solved. Decision-making must not only meet the needs of the people but as well constitute an appropriate land use policy to end harmful land use changes. The State must propose an appropriate management strategy based on international land resource management conventions.

Conclusions
Analyses of land use changes over a 40-year period at the watershed scale, using geospatial techniques and the DPSIR methodology, have helped to improve the understanding of drivers, pressures, conditions, impacts, and responses to adverse changes in land cover. The spatial distribution of land cover for the investigated period indicated increasing cultivated land (19.6%) and decreasing wetland area (−17.5%). The DPSIR model provided qualitative means regarding the land use changes and included an explanatory platform to understand the complexity of these changes. The novelty of this study lies in the combination of quantitative techniques for the detection of land use changes, such as remote sensing and ground truth observations, and qualitative studies, namely diagnostics, surveys, and group interviews. The purpose was to improve the understanding and link motivation of farmers to the observed change in land cover/use.
Despite State interventions in land use planning and management, the level of satisfaction of farmers is low and their social involvement is very low or even absent. A State policy in the 90s that was used as development strategy for coastal areas supplying resources and labor from "interior territories" has reinforced the marginalization and vulnerability of these rural inland territories that initiated the movement leading to the Tunisian revolution in 2011. Today, in the still brittle post-revolutionary context, the Tunisian State must thus, face the socio-economic challenges (rural exodus, unemployment), environmental degradation (water and soil resources), and climatic change for a sustainable development in these agricultural landscapes.
The results of the study have important implications for similar water-scarce and semiarid areas. They show that allocation of land property rights, resource allocation, and improved marketing of agricultural products are important instruments to help rural residents and avoid rural exodus. Results of this study can be used to ensure sustainable management of water and soil resources in areas with similar climate and socio-economic conditions.