Next Article in Journal
Influences on the Seismic Response of a Gravity Dam with Different Foundation and Reservoir Modeling Assumptions
Previous Article in Journal
Proof-of-Concept of a Quasi-2D Water-Quality Modelling Approach to Simulate Transverse Mixing in Rivers
Previous Article in Special Issue
Assessment of the Environmental Impact of Acid Mine Drainage on Surface Water, Stream Sediments, and Macrophytes Using a Battery of Chemical and Ecotoxicological Indicators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying Risk Indices to Assess and Manage Soil Salinization and Sodification in Crop Fields within a Mediterranean Hydro-Agricultural Area

1
Escola Superior Agrária, Instituto Politécnico de Beja, R. Pedro Soares S/N, 7800-295 Beja, Portugal
2
GeoBioTec, NOVA School of Science and Technology, Campus da Caparica, 2829-516 Caparica, Portugal
3
Earth Remote Sensing Laboratory (EaRSLab), Universidade de Évora, 7000-671 Évora, Portugal
4
Instituto de Ciências da Terra (ICT), Universidade de Évora, 7000-671 Évora, Portugal
5
Departamento de Física, Escola de Ciências e Tecnologia, Universidade de Évora, 7000-671 Évora, Portugal
6
Centro de Química, Universidade de Trás-os-Montes e Alto Douro, Quinta de Prados, 5000-801 Vila Real, Portugal
7
Centro Operativo e de Tecnologia de Regadio, Quinta da Saúde, Apartado 354, 7800-999 Beja, Portugal
*
Author to whom correspondence should be addressed.
Water 2021, 13(21), 3070; https://doi.org/10.3390/w13213070
Submission received: 29 September 2021 / Revised: 27 October 2021 / Accepted: 28 October 2021 / Published: 2 November 2021

Abstract

:
Irrigation-induced or secondary salinization can occur when salts are introduced by irrigation water and accumulate within the root zone due to insufficient leaching. Mediterranean regions are especially susceptible, given the predominant climate and the expanding of irrigation areas. In this study, two indices to assess the risk of salinization (RSA) and sodification (RSO), previously applied at a regional scale, were used in a hydro-agricultural area (AHA) in Southern Portugal, in ten crop fields. Information on climate, irrigation water quality, soil characteristics, and land use was obtained from large databases and from local data. The results revealed the feasibility of using the RSA and RSO indices both on large and smaller scales, seeing as most of the area in the monitored crop fields presented the same risk classes (62% in RSA and 78% in RSO). Deviations were due to the reduction in scores for drainage and, in the case of RSO, the assigned irrigation method based on the land occupation class. Considering that different spatial scales of risk assessment are associated with different objectives and management options, a risk management framework was outlined following a multi-scale perspective for mitigation actions in salt-sensitive areas, ranging from territorial planning to the adoption of on-farm practices that can contribute to the sustainability of irrigated agriculture.

1. Introduction

Soil quality declines when affected by salinization processes which reduce the soil capacity to provide ecosystem services, such as physical support for plants, nutrient cycling, and the regulation of the water supply or habitat for microorganisms [1,2]. Soil salinization problems are related with excessive concentrations of either soluble salts or exchangeable sodium (Na) [3]. An excess of soluble salts in the soil, or soil salinity, results in increased osmotic stress, affecting water uptake and causing ion toxicities or an imbalanced accumulation of specific ions in plants; the dominance of Na in the soil exchange complex, or soil sodicity, may lead to the dispersion of soil colloids that affects the movement of water and air, with soils becoming prone to the formation of surface crusts, thereby hindering the emergence of seedlings, causing water stagnation, reducing infiltration, and creating anoxic conditions [4].
Salt-affected soils are widely distributed in the world, typically in areas where potential evapotranspiration is higher than precipitation or in areas with high water tables [5]. Irrigation-induced or secondary salinization can occur when salts are introduced by irrigation water and accumulate within the root zone due to insufficient leaching, a process that is intensified by factors such as poor-quality irrigation water, low hydraulic conductivity of soil layers, or high evaporative conditions [6]. It is estimated that more than 20% of the world’s irrigated area is affected by salinization caused by irrigation [7]. In Mediterranean regions, 25% of irrigated cropland is affected by moderate to high salinization processes [8]. An increased temperature and changes in regional precipitation patterns resulting from climate change, along with water quality degradation, could contribute to the increase in salt-affect areas [9,10,11,12,13]. Additionally, the expansion of irrigation areas is characterized by rapid changes in land use and the redistribution of water resources, therefore reshaping the spatial heterogeneity of soil salinity [14,15].
In irrigated agroecosystems, salt-affected soils vary spatially and temporally, depending, among other factors, on climate, soil physicochemical properties, topography, irrigation water quality, and irrigation systems [6]. Thereby, any risk assessment methodology must be based on robust indicators to assess whether a given soil could be affected by secondary salinization. Risk assessment methods can be performed at different levels of spatial and temporal scales and associated levels of detail, in each case using different types of approaches, namely, a qualitative, quantitative, or model approach [16,17]. Tomaz et al. [18] proposed two indices to assess the risk of soil salinization (RSA) and soil sodification (RSO) at a regional scale, based on easily available data. The RSA index was obtained from a set of four indicators, namely, the aridity index, slope, soil drainage, and irrigation water quality. In the case of RSO, the type of irrigation method was added as a fifth indicator. The aim of the RSA and RSO indices was to provide a straightforward tool to help identify susceptible areas, to select appropriate agronomical strategies and management options, and to prevent soil degradation in irrigated lands. It is useful, from the standpoint of the differentiated and sustainable management of soil and water resources, to understand whether these indices have the potential to be used with greater spatial resolution at the scale of farms and agricultural plots.
Crop type, water quality, and soil properties determine to a large extent the management practices required to optimize production in salt-sensitive areas [19]. In order to achieve a long-term “salt” sustainability of an irrigated field or an irrigation plan, the monitoring and assessment of water quality for irrigation should be complemented with the recommendation of appropriate agricultural water management options for preventing those impacts [10]. Furthermore, control measures of soil salinity/sodicity must be implemented to sustain irrigated agriculture and to prevent the pollution of associated water resources. Such measures must be chosen with recognition of the natural processes operative in irrigated, geohydrologic systems, not just those on-farm, and with an understanding of how they affect the quality of soil and water resources in addition to crop production limitations [20]. The strategies that can be used to prevent, reduce, or adapt to soil salinity or sodicity vary as a function of the level of risk assessed but its implementation must also take into account the spatial scale of assessment. Thereby, the risk management needs to be adapted to different scales, from the cultivation field level to the regional/landscape level, using integrated approaches [21,22].
Taking into account the above considerations, with this study we intended to: (i) compare the salinization and sodification risk assessment based on easily available data at the hydro-agricultural area (AHA) scale, and, based on more detailed data obtained from in situ monitoring, at the crop field scale; (ii) evaluate the feasibility of using the RSA and RSO indices at the crop field scale; (iii) outline a multi-scale risk management framework for salt-sensitive areas in irrigated agroecosystems.

2. Materials and Methods

2.1. Study Area

2.1.1. Hydro-Agricultural Area

The AHA of Brinches-Enxoé was considered for the irrigation plan scale assessment. This AHA, located on the left bank of the Guadiana Basin (South Portugal), is one of the 21 areas of the Alqueva irrigation plan, part of the Multipurpose Development of Alqueva (EFMA—Empreendimento de Fins Múltiplos de Alqueva), centered in the Alqueva reservoir and presently covering 120,000 ha. The AHA of Brinches-Enxoé has both pressurized and gravity conveyance networks, with origins in the Laje reservoir and the Montinhos reservoir, respectively. The equipped area is 5061 ha with 3976 ha under exploration in 2018. The climate in the region is mainly temperate with a hot and dry summer (Mediterranean) with an annual precipitation and average mean monthly temperature of, respectively, 558 mm and 16.9 °C (long-term means for the 1981–2010 period, [23]). During 2018, data from an automatic meteorological station (37.97° N; 7.55° W) of the Irrigation Operating and Technology Center (Centro Operativo e de Tecnologia do Regadio—COTR) located in the AHA, near the monitored crop fields, showed that the annual precipitation was 603 mm, and the mean temperature was 16.7 °C. The year was characterized by a very rainy spring, particularly in the months of March (200 mm) and April (70 mm). Predominant soils are Vertisols, Cambissols, and Luvisols (Figure 1).

2.1.2. Irrigated Crop Fields

For the crop field level risk assessment, ten agricultural fields, accounting for an area of 125 ha, within the Brinches-Enxoé area were evaluated. The ten farm plots were distributed over 88 ha of field, grass, and forage crops; 22 ha of intensive olive groves; and 15 ha of vineyards (Table 1). These fields were assessed under the project Fitofarmgest that aimed at the study of the sustainable management of pesticides in olive groves, vineyards, and field crops in the EFMA’s area. Soils in the farm plots were conventionally tilled, as commonly performed by the farmers, except for maize, which was planted as a second crop in rotation shortly after the harvest of an Autumn-Winter cereal.
Solid fertilizers used were mainly binary or ternary formulations of Nitrogen, Phosphorus, and Potassium (kg N ha−1; kg P2O5 ha−1; kg K2O ha−1) applied at sowing or at the beginning of the vegetative cycle using localized placement and broadcasting (Table 2). Water-soluble and liquid fertilizers, applied over the crops’ cycle through irrigation water, were mostly nitrogen fertilizers but also included other fertilizers containing Sulfur (kg SO2 ha−1), Iron (kg Fe chelates ha−1), or Calcium (kg CaO ha−1). Foliar applications of Boron (kg B ha−1) were employed in some of the crops.

2.2. Risk Assessment

The Salinization Risk (RSA) index and the Sodification Risk (RSO) index, proposed by Tomaz et al. [18], were used. The RSA and RSO indices incorporate sets of 4 and 5 indicators, respectively, and their value is obtained by their score sum, as presented in Equations (1) and (2):
R S A = i = 1 4 S i = S A I + S S l + S S D + S E C w ,
R S O = i = 1 5 S i = S A I + S S l + S S D + S I M + S E C w S A R ,
where Si is the rating score assigned to each of the indicators’ class, as given in Table 3, namely: SAI—Aridity Index class score (Si = 1, 2, …, or 5); SSl—Slope (%) class score (Si = 1, 3, or 5); SSD—Soil drainage class score (Si = 1, 3, or 5); SECw—Water salinity (electrical conductivity; dS m−1) class score (Si = 1, 3, or 5); SIM—Irrigation method score (Si = 1 or 5); SECw˄SAR—Water sodicity (electrical conductivity and Sodium adsorption ratio) class score (Si = 1, 3, or 5). The sodium adsorption ratio was given by Equation (3):
S A R = [ N a + ] [ C a 2 + ] + [ M g 2 + ] 2 ,
where [Na+], [Ca2+], and [Mg2+] are the concentrations of Sodium, Calcium, and Magnesium in the irrigation water, in meq L−1.
From the RSA values, five risk classes were assigned, namely: Very Low, RSA < 7; Low, 7 ≤ RSA ≤ 10; Moderate, 11 ≤ RSA ≤ 14; High, 15 ≤ RSA ≤ 18; Very High, RSA > 18. In the case of RSO values, the classes considered were: Very Low, RSO < 8; Low, 8 ≤ RSO ≤ 12; Moderate, 13 ≤ RSO ≤ 17; High, 18 ≤ RSO ≤ 22; Very High, RSO > 22.

2.3. Data Collection, Sampling, and Physical-Chemical Analyses

The RSA and RSO indices were applied at two different spatial scales of assessment, respectively, the hydro-agricultural area level and crop field level. The assessment was performed during 2018. The collection of data and map building for some indicators varied from large databases and maps to in situ measurements, as presented in Table 4.
Soil sampling was carried out at the beginning (March–April) and at the end (October–November) of the crops’ cycle and consisted of collecting a composed sample of every approximate 5 ha area, obtained from a mixture of sub-samples collected at randomly selected points following a zigzag trajectory, totaling 175 georeferenced points and 27 composed samples [29]. Soil texture was obtained by ≤2 mm-sieving and a particle size analysis, following the ISO 11277:2020 method by sieving and sedimentation [30].
Table 4. Scales of evaluation, data sources, and estimation methods for the application of the RSA and RSO indices.
Table 4. Scales of evaluation, data sources, and estimation methods for the application of the RSA and RSO indices.
IndicatorData Sources and Estimation Methods
Hydro-Agricultural AreaCrop Fields
AIRatio P/ETP [27]; P and ETP obtained from the ERA5 dataset of Monthly-mean averages values. ERA5 is the fifth-generation reanalysis product provided by the European Centre for Medium-Range Weather Forecasts (ECMWF; [31])
SlPortuguese soil and land use capacity chart [32]
SDEstimated from textural classes of the surface horizons/layers of the predominant soil units in the Portuguese soil and land use capacity chart [32]Estimated from textural classes determined from particle size distribution of soil samples
ECwEvaluation of the reservoir water quality for irrigation according to the FAO guidelines [28]Evaluation of the hydrants’ water quality for irrigation according to the FAO guidelines [28]
IMA probable irrigation method is assigned based on the chart of Land Use and Occupation (COS2018) in Continental Portugal, made available from the Portuguese Directorate General for Territory (DGT—Direcção-Geral do Território) [33]Actual data of crop occupation made available from the Alqueva Development and Infrastructures Company (EDIA—Empresa de Desenvolvimento e Infraestruturas de Alqueva) [34]
SAR ˄ ECWEvaluation of the reservoir water quality for irrigation according to the FAO guidelines [28]Evaluation of the hydrants’ water quality for irrigation according to the FAO guidelines [28]
AI—Aridity Index; Sl—Slope; SD—Soil Drainage class; ECW—Irrigation water Electrical Conductivity; IM—Irrigation Method; SAR ˄ ECW—Irrigation water Sodium Adsorption Ratio and Electrical Conductivity.
The study covered four water sampling campaigns, performed in April (Apr), July (Jul), October (Oct), and November (Nov). The sampling sites at the reservoirs included two in the Laje reservoir—Laje (L) and Laje Saída (LS), and at the Laje Estação Elevatória (LE)—and one in the Montinhos reservoir (M). For the local assessment, water samples were collected from eight irrigation hydrants, namely, H5, H6, H6, H21, H22, H23, H33 (water from Laje), and H16 (water conveyed from Montinhos). At each sampling location in the reservoirs, 2 L of surface water were collected, at a depth of 50 cm, using a Van Dorn bottle, and then stored in polyethylene bottles. At the hydrants, 2 L water samples were collected in polyethylene bottles. The water samples were transported to the laboratory, in a cooler at 4 °C, and stored, following the requisites for water conservation for each parameter [35]. The concentrations of Ca2+, Mg2+, and Na+ were determined by the ionic chromatography methodology [35]. Water ECW in the reservoirs was measured in situ at the same dates at a 50 cm depth using a multiparametric probe, YSI 6820 MPS probe. In the case of the hydrants’ water samples, ECW was determined in the laboratory by conductometry. Values of SAR were calculated after Equation (3).
For the AHA scale assessment, a spatial combination of maps of the different classes of indicators was made to obtain the final risk maps. The risk assessment of the crop fields was incorporated in these maps. For this purpose, the 125 ha area of the ten monitored plots was that of the polygon delimited by the sampling point distribution. The areas were obtained using the QGIS (version 3.12.3) built-in function. The comparison analysis between the results obtained for the area of study applying the general methodology at the regional scale and those obtained from the Fitofarmgest sampling points was conducted considering that the points were homogeneously distributed in the polygon with a known area. In this way, the area that maintained (or changed) the risk class in each case was obtained considering that it was directly proportional to the number of points that maintained (or changed) the classification, as given in Equation (4):
A m / c = A p × S P m / c S P
where Am/c is the area that maintained (m) or changed (c) the classification; Ap is the area of the polygon; SP is the number of sampling points in the polygon; and SPm/c is the number of sampling points that maintained (m) or changed (c) the classification.

3. Results and Discussion

3.1. Indicators

The aridity index ranged between 0.28 and 0.29; therefore, the entire area of the AHA, including the monitored plots, presented a semi-arid climate (0.20–0.50) according to the classification of UNEP (1993; [27]). An extensive part of the area (76.7%) presented slopes between 3% and 8%, included in class II, according to the risk methodology calculation (Figure 2). Steeper terrain (class I) only occurred in 12.7% of the AHA, and flat terrain (class III) occurred on 10.6% of the total area.
The drainage class evaluation at the farm plots, obtained from the texture determined by the soil particle size analysis, showed that about 69% of the monitored area had different classes of drainage when compared with the drainage class estimated for the AHA (Figure 3). Mainly, the variation occurred from poorly drained soils (grouped in Class III) in the AHA scale assessment, to moderately drained (grouped in Class II) in the crop field scale assessment, therefore indicating a prevalence of medium textures determined in the soil samples. Such differences point to a spatial variability in the soil texture comprised by each soil unit area that is not translated by the Portuguese Soil Chart (scale 1:25,000), which is the source for the drainage estimation at a regional scale.
Regarding the land use and irrigation method, there was a change in the classification in 24% of the monitored area when using actual data (Figure 4). At the AHA scale, when assigning a probable irrigation method from the land occupation in COS2018, the class of temporary crops was considered as potentially irrigated with sprinkler systems. This assumption implies that fields actually occupied by permanent crops, irrigated by drip irrigation systems were wrongly assumed to be under sprinkler irrigation.
The temporal variation observed for ECW and SAR in the reservoirs (Laje and Montinhos) is presented in Figure 5. In the case of ECW, the highest average values were measured in April in both reservoirs. As for SAR, the highest values were registered at Laje Estação Elevatória (LE) in April and July. These higher values in Spring and early Summer could be due to the 2017 and early 2018 drought, with low precipitation rates that led to low water storage and, possibly, reduced salt dissolution.
In the case of the irrigation hydrants, H23 (water supplier to farm plots S1—sunflower and G1—grapevine) presented the higher values of ECW throughout the year, while H33, which supplied irrigation water to Alfalfa (A), showed the highest SAR value, in October (Figure 6). Regardless of these slight variations, and according to the methodology of salinization and sodification risk prediction, the irrigation water quality indicators fall in the same classes for both the regional and local scales of assessment, namely, for RSA estimation, the ECW values showed no degree of restriction of use and for RSO estimation, the combined ECW and SAR values indicated a slight-to-moderate degree of restriction of use in all water samples.

3.2. Maps of Salinization and Sodification Risk

The risk map of soil salinization in the AHA of Brinches-Enxoé shows that there was a predominance of Moderate risk areas (55%), followed by Low (43%), and Very Low risk (2%). There were no areas within the higher risk classes (Figure 7). However, when evaluating at the crop field scale, the Low-risk areas prevail (71%), and the remaining areas of the plots are of Moderate risk of salinization (29%). The moderate risk of salinization zones was mainly associated with fine textured soils and flat to low slope terrain (slopes ≤ 2%). In general, 38.2% of the total crop fields’ area changed their risk class of salinization, with the differences resulting mainly from the reduction in scores for soil drainage, which, as explained in Section 3.1, were lower due to an actual prevalence of moderately drained terrain in the monitored fields.
The risk map of soil sodification can be observed in Figure 8. Considering the first tier evaluation, in the AHA of Brinches-Enxoé, the predominant risk class is Moderate (68.0%), followed by High (25.3%), and Low risk (6.6%) classes. In the second tier, Moderate risk areas represent 84.0%, and the remaining area is of High risk of sodification. Areas with high RSO resulted mostly from under sprinkler irrigated crops in the low slope terrain, with fine textured soils. The differences for the two scales of assessment occurred in 21.0% of the fields area and are mostly deviations in a “positive” direction, i.e., from higher to lower risk. Again, the explanation for the differences relies on the reduction in scores for drainage but also, in the case of this index, on the assigned probable irrigation method based on the land occupation classes in COS2018, as referred to in Section 3.1.
Overall, given that most of the area in the monitored crop fields presents the same risk classes (62% in RSA and 78% in RSO), it is feasible to apply both the RSA and RSO for regional/landscape and local/farm scale assessment. For this purpose, a tolerable increase in error can be accepted when the spatial resolution becomes lower and assuming that the identification of soil salinization and sodification risk “Hot Spots” is the primary aim of the assessment [36].

3.3. Hints for Risk Management of Salinization and Sodification in Irrigated Agroecosystems

According to the International Organization for Standardization (ISO) guideline for risk management (ISO 31000:2018), the risk management process involves the systematic application of policies, procedures, and practices to the activities of communicating and consulting, establishing the context and assessing, treating, monitoring, reviewing, recording, and reporting risk [37]. Risk assessment is the overall process of risk identification, risk analysis, and risk evaluation, and it is followed by the risk treatment stage, where options for addressing risk, through its mitigation, should be selected and implemented. This common approach to managing any type of risk may be applied to the risk management of salinization and sodification in irrigated agroecosystems (Figure 9). In this case, risk assessment methods may allow for the identification of sensitive areas where preventive measures or adaptative agronomic options are needed [38]. Different spatial scales are associated with different objectives and measures whose sequence and integration are necessary as a function of the risk assessed. Therefore, for the prevention and mitigation of soil salinization and sodification related to anthropogenic activities, such as irrigation, a multi-scale and holistic perspective should be considered.
The following scales are proposed:
  • Regional scale: management and planning of water and soil resources—long term irrigation water and soil quality assessment regarding its uses (involves the strengthening of actions and links between institutions in the agricultural and environmental sectors responsible for the implementation of resource-planning policies).
  • Irrigation plan/perimeter scale: land use planning—selection and classification of land according to the vulnerability of the areas to salinization and sodification processes; selection and classification of land according to irrigation feasibility; development of guidelines for land use.
  • Farm/crop scale: crop management options—maintain favorable levels of salts in the root zone through the use of leaching fractions; adapted irrigation and field drainage strategies; selection of tolerant crops or varieties; other agronomic practices for prevention and mitigation.
Methods of salinity prevention and mitigation also vary as a function of the level of risk assessed. For instance, in high-risk areas, it may be necessary to monitor the groundwater levels, the amount of salt in the soils, and/or the irrigation water quality, while at moderate risk areas, the adoption of cropping practices that minimize water logging can be recommended [39]. From a farm/crop field management viewpoint, differentiation can also be established for each risk class, from the appliance of salt amendments and leaching fractions in higher risk areas, to the adoption of improved irrigation practices, such as methods with higher irrigation efficiency and uniformity, or increased irrigation frequency, in the case of moderate to low risk zones [21].
Under a regional/landscape perspective, soil salinization mitigation in irrigation areas should involve long-term irrigation water quality monitoring and assessment programs, including sampling campaigns covering different periods in the wet and dry seasons, particularly under drought conditions [10,21]. For smaller areas, as is the case of farms and agricultural fields, higher detail can be applied for the control of soil and water salinity and sodicity. The recommendation of preventive irrigation practices, such as maintenance leaching applications and properly designed drainage systems, should be considered [20,40,41]. In salt-sensitive areas, irrigation should aim to achieve a sustainable salt balance in the soil profile. For this, frequent irrigation events and schedules providing crops’ leaching requirements and adapted to the crop salt tolerance and selected irrigation method are usually encouraged [4,42]. The judicious use of multi-quality waters’ blending is also indicated by Rhoades et al. [20], Qadir and Oster [43], and Phogat et al. [44]. Since crop fertilization can be a source of salts in the soil, the appropriate timing and fractioning of fertilizers’ application should be considered [45]. When salinity is associated with shallow and saline groundwater, provisions for surface/sub-surface drainage systems further help in aeration, lateral water flows for salt removal, and trafficability [20,42].

4. Conclusions

This study used two risk indices to predict salinization and sodification in a hydro-agricultural area in Southern Portugal and ten farm plots within that area, aiming at the comparison and adaptability of the risk indices to different spatial scales of evaluation. Low and moderate salinization risk areas prevailed in the study area. Regarding sodification, moderate risk was the predominant class. The results revealed the feasibility of using the RSA and RSO indices both on large and smaller scales, as most of the area in the monitored crop fields presented the same risk classes (62% in RSA and 78% in RSO). Deviations were due to the reduction in scores for drainage and, in the case of RSO, the assigned irrigation method based on the land occupation class. Nevertheless, these deviations were mostly in a “positive” direction, from higher to lower risk classes of salinization and sodification risk.
The study outlined a risk management framework applied to soil salinization and sodification in irrigated areas adapting the ISO guidelines. The framework considered that different spatial scales are associated with different objectives and measures whose sequence and integration are necessary. Therefore, a multi-scale perspective for mitigation actions in salt-sensitive areas was delineated, ranging from territorial planning to the adoption of on-farm practices that can contribute to the sustainability of irrigated agriculture.

Author Contributions

Conceptualization, A.T., M.J.C. and P.P.; methodology, A.T., M.J.C. and P.P.; validation, A.T., J.D., I.G., M.M.P., M.F., L.B., M.P. and P.P.; investigation, A.T., M.J.C., J.C., J.D., A.C., I.M., C.M., I.G., M.M.P., M.F., L.B., M.P. and P.P.; writing—original draft preparation, A.T.; writing—review and editing, A.T., M.J.C., J.C., I.G., M.M.P., M.F., L.B., M.P. and P.P.; visualization, M.J.C.; supervision, P.P.; project administration, P.P.; funding acquisition, A.T., M.J.C., J.C., M.P. and P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research is co-funded by the European Union through the European Regional Development Fund, included in the COMPETE 2020 (Operational Program Competitiveness and Internationalization) through the ICT project (UIDB/04683/2020), with the reference POCI-01-0145-FEDER-007690, through Geobiotec (UIDB/04035/2020) and Centro de Química (UIBD/00616/2020)funded by FCT—Fundação para a Ciência e a Tecnologia, Portugal, and through the FitoFarmGest Operational Group (PDR2020-101-030926).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors thank the General Directorate of Agriculture and Rural Development of Portugal (DGADR—Direção Geral de Agricultura e Desenvolvimento Rural) and the Alqueva Development and Infrastructures Company (EDIA—Empresa de Desenvolvimento e Infraestruturas do Alqueva) for providing, respectively, the Portuguese soil and land use capacity maps, and the Alqueva crop occupation map in digital form.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Nuñez, M.; Tóth, T. Ecosystem Services and Salinisation. In Final Report of the EIP-AGRI Focus Group: Soil Salinisation; European Comission: Brussels, Belgium, 2020; p. 10. [Google Scholar]
  2. Adhikari, K.; Hartemink, A.E. Linking Soils to Ecosystem Services—A Global Review. Geoderma 2016, 262, 101–111. [Google Scholar] [CrossRef]
  3. Richards, L.A. Diagnosis and Improvement of Saline and Alkali Soils; Agriculture Handbook; USDA: Washington, DC, USA, 1954; Volume 60.
  4. Minhas, P.S.; Ramos, T.B.; Ben-Gal, A.; Pereira, L.S. Coping with Salinity in Irrigated Agriculture: Crop Evapotranspiration and Water Management Issues. Agric. Water Manag. 2020, 227, 105832. [Google Scholar] [CrossRef]
  5. Weil, R.R.; Brady, N.C. The Nature and Properties of Soils, 15th ed.; Pearson: Columbus, OH, USA, 2016; ISBN 978-0-13-325448-8. [Google Scholar]
  6. Rengasamy, P. Soil Processes Affecting Crop Production in Salt-Affected Soils. Funct. Plant. Biol. 2010, 37, 613. [Google Scholar] [CrossRef]
  7. Qadir, M.; Quillérou, E.; Nangia, V.; Murtaza, G.; Singh, M.; Thomas, R.J.; Drechsel, P.; Noble, A.D. Economics of Salt-Induced Land Degradation and Restoration. Nat. Resour. Forum. 2014, 38, 282–295. [Google Scholar] [CrossRef]
  8. Mateo-Sagasta, J.; Burke, J. Agriculture and Water Quality Interactions: A Global Overview; SOLAW Background Thematic Report-TR08; Food and Agriculture Organization of the United Nations: Rome, Italy, 2010. [Google Scholar]
  9. Stolte, J.; Tesfai, M.; Øygarden, L.; Kværnø, S.; Keizer, J.; Verheijen, F.; Panagos, P.; Ballabio, C.; Hessel, R. Soil Threats in Europe; JCR Technical Report; European Comission: Brussels, Belgium, 2016. [Google Scholar]
  10. Tomaz, A.; Palma, P.; Fialho, S.; Lima, A.; Alvarenga, P.; Potes, M.; Salgado, R. Spatial and Temporal Dynamics of Irrigation Water Quality under Drought Conditions in a Large Reservoir in Southern Portugal. Environ. Monit. Assess. 2020, 192, 93. [Google Scholar] [CrossRef] [PubMed]
  11. Tomaz, A.; Palma, P.; Alvarenga, P.; Gonçalves, M.C. Chapter 13-Soil salinity risk in a climate change scenario and its effect on crop yield. In Climate Change and Soil Interactions; Prasad, M.N.V., Pietrzykowski, M., Eds.; Elsevier: Amsterdam, The Netherlands, 2020; pp. 351–396. ISBN 978-0-12-818032-7. [Google Scholar]
  12. Yeo, A. Predicting the Interaction between the Effects of Salinity and Climate Change on Crop Plants. Sci. Hortic. 1998, 78, 159–174. [Google Scholar] [CrossRef]
  13. Zalidis, G.; Stamatiadis, S.; Takavakoglou, V.; Eskridge, K.; Misopolinos, N. Impacts of Agricultural Practices on Soil and Water Quality in the Mediterranean Region and Proposed Assessment Methodology. Agric. Ecosyst. Environ. 2002, 88, 137–146. [Google Scholar] [CrossRef]
  14. Singh, A. Soil Salinization and Waterlogging: A Threat to Environment and Agricultural Sustainability. Ecol. Indic. 2015, 57, 128–130. [Google Scholar] [CrossRef]
  15. Muyibul, Z.; Jianxin, X.; Muhtar, P.; Qingdong, S.; Run, Z. Spatiotemporal Changes of Land Use/Cover from 1995 to 2015 in an Oasis in the Middle Reaches of the Keriya River, Southern Tarim Basin, Northwest China. Catena 2018, 171, 416–425. [Google Scholar] [CrossRef]
  16. Bloem, E.; van der Zee, S.E.A.T.M.; Tóth, T.; Hagyo, A. Risk Assessment Methods of Salinity; Report EUR 24097 EN; Publications Office of the European Union: Luxembourg, 2008; ISBN 978-92-79-14291-8. [Google Scholar]
  17. Shahid, S. Developments in Soil Salinity Assessment, Modeling, Mapping, and Monitoring from Regional to Submicroscopic Scales; Shahid, S., Abdelfattah, M., Taha, F., Eds.; Springer: Dordrecht, The Netherlands, 2013; pp. 3–43. ISBN 978-94-007-5683-0. [Google Scholar]
  18. Tomaz, A.; Palma, P.; Fialho, S.; Lima, A.; Alvarenga, P.; Potes, M.; Costa, M.; Salgado, R. Risk Assessment of Irrigation-Related Soil Salinization and Sodification in Mediterranean Areas. Water 2020, 12, 3569. [Google Scholar] [CrossRef]
  19. Chartzoulakis, K.S. The Use of Saline Water for Irrigation of Olives: Effects on Growth, Phisiology, Yield and Oil Quality. Acta Hortic. 2011, 888, 97–108. [Google Scholar] [CrossRef]
  20. Rhoades, J.D.; Kandiah, A.; Mashali, A.M. The Use of Saline Waters for Crop. Production; FAO Irrigation and Drainage Paper; Food and Agriculture Organization of the United Nations: Rome, Italy, 1992; ISBN 978-92-5-103237-4. [Google Scholar]
  21. Hillel, D. Salinity Management for Sustainable Irrigation: Integrating Science, Environment, and Economics; The World Bank: Washington, DC, USA, 2000; ISBN 978-0-8213-4773-7. [Google Scholar]
  22. Ren, D.; Wei, B.; Xu, X.; Engel, B.; Li, G.; Huang, Q.; Xiong, Y.; Huang, G. Analyzing Spatiotemporal Characteristics of Soil Salinity in Arid Irrigated Agro-Ecosystems Using Integrated Approaches. Geoderma 2019, 356, 113935. [Google Scholar] [CrossRef]
  23. IPMA. Climate Normals-1981-2010-Beja. Available online: https://www.ipma.pt/en/oclima/normais.clima/1981-2010/002/ (accessed on 6 July 2019).
  24. IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015; FAO: Rome, Italy, 2014. [Google Scholar]
  25. Madeira, M.; Monteiro, F. A Base de Referência Para Os Solos Do Mundo e a Classificação Dos Solos de Portugal. Rev. De Cienc. Agrárias 2004, 27, 13–23. [Google Scholar]
  26. Lopes, R.; Madeira, M.; Arsénio, P. Mapping of Land Units and Land Capability Classification in Portugal. The Case of Lourinhã Municipality. Finisterra 2017, 52, 103–127. [Google Scholar] [CrossRef] [Green Version]
  27. UNEP. World Atlas of Desertification; The United Nations Environmental Programme: London, UK, 1993. [Google Scholar]
  28. Ayers, R.S.; Westcot, D.W. Water Quality for Agriculture; FAO Irrigation and Drainage Paper; Food and Agriculture Organization of the United Nations: Rome, Italy, 1985; ISBN 978-92-5-102263-4. [Google Scholar]
  29. Varennes, A.D. Produtividade Dos Solos e Ambiente; Escolar Editora: Lisboa, Portugal, 2003. [Google Scholar]
  30. International Organization for Standardization. Soil Quality—Determination of Particle Size Distribution in Mineral. Soil Material—Method by Sieving and Sedimentation (ISO Standard No. 11277:2020); International Organization for Standardization: London, UK, 2020; p. 16. [Google Scholar]
  31. Hersbach, H.; Dee, D. ERA5 Reanalysis Is in Production. In ECMWF Newsletter; European Centre for Medium-Range Weather Forecasts: Reading, UK, 2016; 147p. [Google Scholar]
  32. SROA/CNROA/IEADR Cartas de Solos e Capacidade de Uso (1:25000); Direção Geral de Agricultura e Desenvolvimento Rural: Lisbon, Portugal, 1960.
  33. DGT Carta de Uso e Ocupação Do Solo (COS2018); Direção Geral do Território: Lisbon, Portugal, 2018.
  34. EDIA Crop Occupation 2018; Empresa de Desenvolvimento e Infraestruturas de Alqueva: Beja, Portugal, 2018.
  35. APHA. Standard Methods for the Examination of Water and Wastewater, 20th ed.; American Public Health Association, American Water Works Association and Water Environmental Federation: Washington, DC, USA, 1998. [Google Scholar]
  36. Eckelmann, W.; Baritz, R.; Bialousz, S.; Bielek, P.; Carre, F.; Houskova, B.; Jones, R.J.A.; Kibblewhite, M.; Kozak, J.; le Bas, C.; et al. Common Criteria for Risk Area Identification According to Soil Threats; Office for Official Publications of the European Communities: Luxembourg, 2006. [Google Scholar]
  37. International Organization for Standardization. Risk Management—Guidelines (ISO Standard No. 31000:2018); International Organization for Standardization: London, UK, 2018; p. 16. [Google Scholar]
  38. Thompson, R.; Amezketa, E.; Nuñez, M.; Falsone, G.; Toth, T.; Paz, A.; Gonçalves, M.; Castanheira, N.; Gould, I.; Tamvakidis, S. EIP-AGRI Focus Group Soil Salinisation. Measuring, Mapping and Monitoring of Soil Salinity; The European Innovation Partnership: Brussels, Belgium, 2020. [Google Scholar]
  39. Department of Environment and Resource Management. Salinity Management Handbook, 2nd ed.; Department of Environment and Resource Management: Queensland, Australia, 2011; ISBN 0-7242-7412-X. [Google Scholar]
  40. Corwin, D.L.; Rhoades, J.D.; Šimůnek, J. Leaching Requirement for Soil Salinity Control: Steady-State versus Transient Models. Agric. Water Manag. 2007, 90, 165–180. [Google Scholar] [CrossRef]
  41. Wallender, W.W.; Tanji, K.K. Agricultural Salinity Assessment and Management, 2nd ed.; American Society of Civil Engineers: Reston, VA, USA, 2011; ISBN 978-0-7844-1169-8. [Google Scholar]
  42. Hanson, B.R.; Grattan, S.R.; Fulton, A. Agricultural Salinity and Drainage; Water Management Series publication; University of California Irrigation Program: Davis, CA, USA, 2006. [Google Scholar]
  43. Qadir, M.; Oster, J.D. Crop and Irrigation Management Strategies for Saline-Sodic Soils and Waters Aimed at Environmentally Sustainable Agriculture. Sci. Total Environ. 2004, 323, 1–19. [Google Scholar] [CrossRef]
  44. Phogat, V.; Mallants, D.; Cox, J.W.; Šimůnek, J.; Oliver, D.P.; Awad, J. Management of Soil Salinity Associated with Irrigation of Protected Crops. Agric. Water Manag. 2020, 227, 105845. [Google Scholar] [CrossRef]
  45. Machado, R.; Serralheiro, R. Soil Salinity: Effect on Vegetable Crop Growth. Management Practices to Prevent and Mitigate Soil Salinization. Horticulturae 2017, 3, 30. [Google Scholar] [CrossRef]
Figure 1. Reference Soil Groups (RSG) of the AHA of Brinches-Enxoé (based on [24,25,26]) and position of Laje and Montinhos reservoirs in the area (line pattern filled in the map). The bottom-left inset shows the location of the AHA under study.
Figure 1. Reference Soil Groups (RSG) of the AHA of Brinches-Enxoé (based on [24,25,26]) and position of Laje and Montinhos reservoirs in the area (line pattern filled in the map). The bottom-left inset shows the location of the AHA under study.
Water 13 03070 g001
Figure 2. Slope classes in the AHA of Brinches-Enxoé. I: >9%; II: 3–8%; III: 0–2%. The black open circles represent the monitored crop fields.
Figure 2. Slope classes in the AHA of Brinches-Enxoé. I: >9%; II: 3–8%; III: 0–2%. The black open circles represent the monitored crop fields.
Water 13 03070 g002
Figure 3. Soil drainage classes in the AHA of Brinches-Enxoé (a) and in the monitored crop fields (b), displaying the similarities and differences of the two scales of evaluation. I: Good; II: Moderate; III: Restricted.
Figure 3. Soil drainage classes in the AHA of Brinches-Enxoé (a) and in the monitored crop fields (b), displaying the similarities and differences of the two scales of evaluation. I: Good; II: Moderate; III: Restricted.
Water 13 03070 g003
Figure 4. Irrigation methods in the AHA of Brinches-Enxoé (a) and in the monitored crop fields (b), displaying the similarities and differences of the two scales of evaluation.
Figure 4. Irrigation methods in the AHA of Brinches-Enxoé (a) and in the monitored crop fields (b), displaying the similarities and differences of the two scales of evaluation.
Water 13 03070 g004
Figure 5. Spatial and temporal variation of ECW (a) and SAR (b), during 2018, at three sites in the Laje reservoir and one site at the Montinhos reservoir. Bars represent the mean, and whiskers represent the standard deviation for three repetitions. Dashed purple lines correspond to the limit values from FAO water quality standards as presented in Table 3. L—Laje; LS—Laje Saída; LE—Laje Estação Elevatória; M—Montinhos. Apr—April; Jul—July; Oct—October; Nov—November.
Figure 5. Spatial and temporal variation of ECW (a) and SAR (b), during 2018, at three sites in the Laje reservoir and one site at the Montinhos reservoir. Bars represent the mean, and whiskers represent the standard deviation for three repetitions. Dashed purple lines correspond to the limit values from FAO water quality standards as presented in Table 3. L—Laje; LS—Laje Saída; LE—Laje Estação Elevatória; M—Montinhos. Apr—April; Jul—July; Oct—October; Nov—November.
Water 13 03070 g005
Figure 6. Spatial and temporal variation of ECW (a) and SAR (b), during 2018, at eight irrigation hydrants. Bars represent the mean, and whiskers represent the standard deviation for three repetitions Dashed purple lines correspond to the limit values from FAO water quality standards as presented in Table 3. Apr—April; Jul—July; Oct—October; Nov—November.
Figure 6. Spatial and temporal variation of ECW (a) and SAR (b), during 2018, at eight irrigation hydrants. Bars represent the mean, and whiskers represent the standard deviation for three repetitions Dashed purple lines correspond to the limit values from FAO water quality standards as presented in Table 3. Apr—April; Jul—July; Oct—October; Nov—November.
Water 13 03070 g006
Figure 7. Salinization risk in the AHA of Brinches-Enxoé (a) and in the monitored crop fields (b) based on the Salinization Risk index (RSA). Similarities and differences of the two scales of evaluation are displayed. Very Low risk: RSA < 7; Low risk: 7 ≤ RSA ≤ 10; Moderate risk: 11 ≤ RSA ≤ 14; High risk: 15 ≤ RSA ≤ 18; Very High risk: RSA > 18.
Figure 7. Salinization risk in the AHA of Brinches-Enxoé (a) and in the monitored crop fields (b) based on the Salinization Risk index (RSA). Similarities and differences of the two scales of evaluation are displayed. Very Low risk: RSA < 7; Low risk: 7 ≤ RSA ≤ 10; Moderate risk: 11 ≤ RSA ≤ 14; High risk: 15 ≤ RSA ≤ 18; Very High risk: RSA > 18.
Water 13 03070 g007
Figure 8. Sodification risk in the AHA of Brinches-Enxoé (a) and in the Fitofarmgest crop fields (b) based on the Sodification Risk index (RSO). Similarities and differences of the two scales of evaluation are displayed. Very Low risk: RSO < 8; Low risk: 8 ≤ RSO ≤ 12; Moderate risk: 13 ≤ RSO ≤ 17; High risk: 18 ≤ RSO ≤ 22; Very High risk: RSO > 18.
Figure 8. Sodification risk in the AHA of Brinches-Enxoé (a) and in the Fitofarmgest crop fields (b) based on the Sodification Risk index (RSO). Similarities and differences of the two scales of evaluation are displayed. Very Low risk: RSO < 8; Low risk: 8 ≤ RSO ≤ 12; Moderate risk: 13 ≤ RSO ≤ 17; High risk: 18 ≤ RSO ≤ 22; Very High risk: RSO > 18.
Water 13 03070 g008
Figure 9. Outline of a salinization and sodification risk management process based on ISO 31000:2018 [37] with risk treatment stage adapted to different spatial levels of assessment.
Figure 9. Outline of a salinization and sodification risk management process based on ISO 31000:2018 [37] with risk treatment stage adapted to different spatial levels of assessment.
Water 13 03070 g009
Table 1. General characterization of the crop fields monitored in the study.
Table 1. General characterization of the crop fields monitored in the study.
Crop 1Monitored Area (ha)Predominant Soil Types 2Irrigation MethodIrrigation Hydrant
Alfalfa (Medicago sativa L.) (A)27.6Calcaric Cambisols; RegosolsSprinklerH33
Permanent Pasture (grasses (70%), legumes (18%) and others (12%)) (P)21.7Calcaric Cambisols; Chromic Luvisols; RegosolsSprinklerH22
Grapevine (Vitis vinifera L. Cv. ‘Antão Vaz) (G1)4.5Chromic VertisolsDripH23
Sunflower (Helianthus annus L.) (S1)13.1Calcaric Cambisols; Chromic VertisolsSprinklerH23
Olive (Olea europaea L. Cv. ‘Cobrançosa’) (O1)8.0Chromic Vertisols; Pellic VertisolsDripH21
Sunflower (Helianthus annus L.) (S2)15.0Pellic Vertisols; Calcaric VertisolsSprinklerH7
Grapevine (Vitis vinifera L. Cv. ‘Aragonez) (G2)9.0Calcaric Cambisols; Chromic VertisolsDripH5
Maize (Zea mays L.) (M)10.3Calcaric Cambisols; Chromic VertisolsSprinklerH6
Grapevine (Vitis vinifera L. Cv. ‘Antão Vaz) (G3)1.5Chromic Luvisols; Calcaric CambisolsDripH16
Olive (Olea europaea L. Cv. ‘Cordovil’) (O2)14.0Calcaric Cambisols; Chromic Vertisols; Vertic LuvisolsDripH16
1 Crop and attributed symbol. 2 RSG—Reference Soil Group and principal qualifiers, when applicable (according to the World reference base for soil resources classification [24]).
Table 2. Main crop management data in 2018.
Table 2. Main crop management data in 2018.
CropSowing/Beginning of Cycle (dd/mm)Seasonal Irrigation Water (m3 ha−1)First Irrigation (dd/mm)Last Irrigation (dd/mm)N (kg N ha−1)P (kg P2O5 ha−1)K (kg K2O ha−1)Other Fertilizers (kg ha−1)Harvest (dd/mm)
A26/04400001/0530/091010--- 1
P- 2100015/0515/1043--6 CaO-
G115/03396028/0630/10505056-31/08
S118/04251719/0401/0812734-16 SO2; 0.2 B27/08
O107/02 3330001/0430/098384920.3 B06/12
S227/04460628/0426/08109401216 SO218/09
G202/04260004/0610/09961021798 Fe; 53 CaO24/09
M18/07480018/0704/1020214421627 SO217/01 4
G315/03300015/0515/10662151-23/08
O215/01 3170015/0523/10692435-21/11
1 Pluriannual forage (2 to 3 cuts per year). 2 Permanent pasture (>10 years) with continuous grazing. 3 Date of the first fertilization 4 Delayed harvest due to the occurrence of a long rainy period following the physiological maturity stage.
Table 3. Indicators, classification method, and rating scores considered in the RSA and RSO indices (adapted from [18]).
Table 3. Indicators, classification method, and rating scores considered in the RSA and RSO indices (adapted from [18]).
IndicatorRating Score
12345
AI—Aridity Index 1≥0.65 (Humid)0.50–0.65 (Dry sub-humid)0.20–0.50 (Semi-arid)0.05–0.20 (Arid)<0.05 (Hyper-arid)
Sl—Slope (%)>9-3–8-0–2
SD—Soil Drainage class 2Good-Moderate -Restricted
ECw (dS m−1) —water Electrical Conductivity<0.7 (No restriction)-0.7–3.0 (Slight to moderate restriction)->3.0 (Severe restriction)
IM—Irrigation MethodDrip---Surface or Sprinkler
SAR ˄ ECW (dS m−1) 3—water Sodium Adsorption Ratio and Electrical ConductivitySAR 0–3 and ECW > 0.7 (No restriction)-SAR 0–3 and 0.2 < ECW < 0.7 (Slight to moderate restriction)-SAR 0–3 and ECW < 0.2 (Severe restriction)
1 After the United Nations Environmental Programme (UNEP) classification [27]. 2 Good: Sand and Loamy Sand; Moderate: Sandy Loam, Sandy Clay Loam, Loam, Clay Loam, Silty Clay Loam, Silty Loam, and Silt; Restricted: Sandy Clay, Silty Clay, and Clay; a restricted drainage class was assigned whenever soil use capacity was severely to very severely limited due to excess water [18]. 3 For different SAR values, different ECW intervals apply [28].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tomaz, A.; Costa, M.J.; Coutinho, J.; Dôres, J.; Catarino, A.; Martins, I.; Mourinha, C.; Guerreiro, I.; Pereira, M.M.; Fabião, M.; et al. Applying Risk Indices to Assess and Manage Soil Salinization and Sodification in Crop Fields within a Mediterranean Hydro-Agricultural Area. Water 2021, 13, 3070. https://doi.org/10.3390/w13213070

AMA Style

Tomaz A, Costa MJ, Coutinho J, Dôres J, Catarino A, Martins I, Mourinha C, Guerreiro I, Pereira MM, Fabião M, et al. Applying Risk Indices to Assess and Manage Soil Salinization and Sodification in Crop Fields within a Mediterranean Hydro-Agricultural Area. Water. 2021; 13(21):3070. https://doi.org/10.3390/w13213070

Chicago/Turabian Style

Tomaz, Alexandra, Maria João Costa, João Coutinho, José Dôres, Adriana Catarino, Inês Martins, Clarisse Mourinha, Isabel Guerreiro, Maria Margarida Pereira, Marta Fabião, and et al. 2021. "Applying Risk Indices to Assess and Manage Soil Salinization and Sodification in Crop Fields within a Mediterranean Hydro-Agricultural Area" Water 13, no. 21: 3070. https://doi.org/10.3390/w13213070

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop