SciELO - Scientific Electronic Library Online

 
vol.45 número4FTIR-ATR spectra as fingerprints of forest floor compositionEstrategia de cierre con Tecnosoles para la recuperación ambiental de la Mina de São Domingos índice de autoresíndice de assuntosPesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Revista de Ciências Agrárias

versão impressa ISSN 0871-018Xversão On-line ISSN 2183-041X

Rev. de Ciências Agrárias vol.45 no.4 Lisboa dez. 2022  Epub 01-Dez-2022

https://doi.org/10.19084/rca.28627 

Artigo

How does climate change affect soil salinity? The case study of a Mediterranean pasture with saline groundwater

O impacto das alterações climáticas na salinidade do solo: caso de estudo de uma pastagem mediterrânica com água subterrânea salina

A. Paz1 

N. Castanheira1 

P. Matos Soares2 

M. Farzamian1 

M.C. Paz3 

M.C. Gonçalves1 

1 Instituto Nacional de Investigação Agrária e Veterinária, I.P. (INIAV), Quinta do Marquês, 2780-157, Oeiras, Portugal

2 Instituto Dom Luiz (IDL), Faculdade de Ciências da Universidade de Lisboa, Campo Grande, Ed. C1, Piso 1, 1749-016, Lisbon, Portugal

3 ESTBarreiro and CIQuiBio, Polytechnic Institute of Setúbal, Rua Américo da Silva Marinho, 2839-001, Lavradio, Portugal


Abstract

The case study is located in an agricultural area of alluvial origin, approximately 10 km northeast of Lisbon (Portugal). The area is a peninsula surrounded by River Tejo, its estuary, and by River Sorraia. The peninsula lies 1 to 2 m above sea level, with saline groundwater arising from the estuarine tides. The soil is typically clayey and homogeneous. The future climate projections for the Mediterranean region show a decrease in precipitation, increase in mean temperature, and increase in sea water level. These conditions aggravate the processes leading to salt-affected soils (increased evapotranspiration, decrease in salt leaching during the rainy season, and increase of the groundwater salinity). This work evaluates the influence of climate change on soil salinity and water content in a rainfed spontaneous pasture, using Hydrus-1D. Three climate scenarios, comprising series of 30 years obtained with a regional climate model, were used: reference scenario (referent to 1971-2000), RCP4.5, and RCP 8.5 (both for 2011-2040). The results show an increase in the soil salinity in the RCP scenarios compared to the reference and a decrease in the mean yearly water content of the topsoil. Both effects limit the productivity of the pasture and even compromise the plants’ growth.

Keywords: soil salinity; soil water; soil modelling; climate change.

Resumo

O caso de estudo situa-se na Lezíria de Vila Franca de Xira, uma zona agrícola de origem aluvial, a cerca de 10 km a nordeste de Lisboa. A Lezíria é uma península rodeada pelo rio Tejo, pelo seu estuário e pelo rio Sorraia, entre 1 a 2 m acima do nível do mar, com água subterrânea salina resultante da influência marinha no estuário. As projeções climáticas para a região indicam uma diminuição da precipitação, aumento da temperatura média e aumento do nível da água do mar. Estas condições agravam os processos de salinização do solo (aumento da evapotranspiração, diminuição da lixiviação dos sais pela chuva e aumento da salinidade das águas subterrâneas). Este trabalho pretende analisar a influência das alterações climáticas na salinidade do solo e teor de água numa pastagem espontânea de sequeiro, utilizando Hydrus-1D. Foram simulados três cenários climáticos, compreendendo séries de 30 anos obtidas a partir de um modelo climático regional: cenário de referência (referente a 1971-2000), RCP4.5 e RCP 8.5 (ambos para 2011-2040). Os resultados mostram um aumento da salinidade do solo nos cenários RCP em relação à referência e uma diminuição no teor de água médio anual na camada superficial do solo. Ambos os efeitos limitam a produtividade da pastagem, podendo comprometer totalmente o crescimento das plantas.

Palavras-chave: salinidade do solo; água do solo; modelação; alterações climáticas

INTRODUCTION

The peninsula of Lezíria de Vila Franca de Xira is an agricultural area of alluvial origin, approximately 10 km northeast of Lisbon (Portugal). It is surrounded by River Tejo, its estuary, and by River Sorraia. The peninsula lies 1 to 2 m above sea level, with saline groundwater arising from the estuarine tides, which salinity is close to that of the sea, with an electrical conductivity (EC) of 40-50 dSm-1 (Rodrigues et al., 2019). Climate projections for the Mediterranean region show a decrease in precipitation, increase in mean temperature, and increase in sea water level. These conditions potentiate the accumulation of salts in the soil profile, which may compromise its agricultural use.

Hydrus-1D is a modelling tool which has been found to accurately model the soil water and salts dynamics (Kanzari et al., 2018; Ramos et al., 2011). In this study, Hydrus-1D was used to evaluate the influence of climate change projections on soil salinity and water content, in a rainfed spontaneous pasture in the south of Lezíria.

MATERIAL AND METHODS

Study site

The site is a rainfed spontaneous pasture dominated by Poaceae. The soil is classified as a Solonchak, according to the Harmonized World Soil Database (Fischer et al., 2012). The soil texture is relatively uniform along the soil profile and is classified as silty-clay. The site is located 2 m above the sea and at about 2.5 km from the estuary of the River Tejo. Primary salinization occurs due to the influence of the marine water from the estuary in groundwater.

Calibration of the Hydrus-1D model

Hydrus-1D is a public domain modelling tool for analysis of one dimensional water flow and solute transport in variably saturated porous media (Simunek et al., 2018).

Five soil layers to a depth of 1.5 m were considered (average depths of 0.15, 0.45, 0.75, 1.05, and 1.35 m). For the model calibration, soil samples were collected from the five soil layers, between May 2017 and October 2018, in a total of 20 samples. The soil volumetric water content (θ) was determined from the soil gravimetric water content (θ w ) and the dry bulk density. The electrical conductivity of the soil saturated paste (EC e ) was determined with a conductivity meter in the extract collected with suction filters from the saturated soil paste. The electrical conductivity of the soil solution, EC sw , was calculated according to the following equation,

where SP, is the gravimetric water content of the saturated soil paste (Rhoades et al., 1989). In Hydrus-1D the van Genuchten-Mualem model was selected to calculate the water transport, using the saturated and residual water content (θ s and θ r ), the saturated hydraulic conductivity (K s ), the inverse of the air-entry value (α), the pore-size distribution index (n), and the pore connectivity (l). θ s and θ r were obtained from the particle size distribution using pedotransfer functions developed by Ramos et al. (2014), which were found to be well-suited for the soils in this region. The equilibrium model was used for determination of the solute transport. The inverse modelling tool was used to calibrate K s , α and n, and the solute transport parameter, dispersivity (D). The final soil parameters used in Hydrus-1D are presented in Table 1.

Table 1 Soil parameters used in Hydrus-1D 

Avg. depth θ r θ s α n K s l D
0.15 m 0.09 0.57 0.03 1.53 17 0.5 15
0.45 m 0.09 0.57 0.02 1.4 65 0.5 14
0.75 m 0.08 0.56 0.015 1.3 10 0.5 14
1.05 m 0.08 0.55 0.01 1.1 65 0.5 25
1.35 m 0.08 0.55 0.01 1.1 25 0.5 25

The crop evapotranspiration (ET c ) was calculated according to the single crop coefficient approach proposed by Allen et al. (1998), in which ET c , is calculated by multiplying the reference crop evapotranspiration (ET o ), by a crop coefficient, K c . ET o was calculated using the Penman-Monteith equation with data from a nearby meteorological station. The K c was set to 0.4 for the initial development and 0.85 in the developed phase, with the cycle length adapted to site conditions. As required by Hydrus-1D, ET c was fractioned in the components of soil evaporation and potential crop transpiration, which were estimated as a function of the Leaf Area Index (LAI) and the corresponding Soil Cover Factor, according to the approach proposed by Ritchie (1972). The LAI of the developed pasture was 3.5. In Hydrus-1D, T was integrated over the root domain using the Feddes model for root water uptake with the parameters for pasture proposed by Wesseling et al. (1991) and the solute stress was considered using the multiplicative threshold model.

The calibration of the model with the field data resulted in a RMSE of calibration of 0.064 for a θ range between 0.2 and 0.55 and a RMSE of 6.5 for a EC sw range between 5 and 70 dSm-1.

Simulation of the climate scenarios

Hydrus-1D was also used to simulate θ and EC sw for three scenarios comprised by 30-years periods: reference, RCP 4.5, and RCP 8.5. The reference scenario corresponds to the model data for the historical period of 1971-2000. RCP 4.5 and RCP 8.5 represent different levels of greenhouse gas concentration in the atmosphere during the period 2041-70, corresponding to an intermediate and a severe scenario, respectively. The climate data was obtained from the CORDEX experiment, with the Regional Climate Model, RACMO22E (Meijgaard et al., 2008), forced by the Global Climate Model, MOHC-HadGEM2-ES. This combination was found to have better results in climate description for this region of Portugal (Soares et al., 2017).

ET c was calculated as described for the calibration process, but ET 0 was now calculated using the climate data from the projection model, using the Hargreaves-Samani equation (Hargreaves & Samani, 1985):

where C is the radiation adjustment coefficient [°C-0.5], Ra is extraterrestrial solar radiation [MJ], λ is the latent heat of vaporization equal to 2.45 MJkg-1, tavg is the average daily temperature, and t min and t max are respectively, the minimum and maximum daily temperatures [°C]. For ‘interior’ locations, where land mass dominates and air masses are not strongly influenced by a large water body (Allen et al., 1998) proposed C=0.16. t min , t max , and the precipitation, were obtained from the climate projections. The depth and EC of ground water was considered content during the simulations (1.2 m), which is an approximation, as its depth varies during the year and is also likely to rise as a consequence of the projected increase in sea water level.

RESULTS

Because EC e allows the analysis of the soil salinity independently of its actual θ, it was calculated using EC sw and θ simulated with Hydrus-1D in eq. (1), where SP is θ s . Figures 1 and 2 show that the median of EC e increases in both RCP scenarios, when compared to the reference scenario. Furthermore, there is a very large increment in the maximum EC e in RCP scenarios, reaching values only tolerable to few plants. EC e shows a similar increase for the deeper soil layers. Figure 3 show the yearly average θ at 15 cm. In the topsoil a decrease in average θ in the RCP scenarios compared to the reference scenario is evident. The decrease in θ does not occur at the layers below, which is likely due to the effect of capillary rise of the groundwater to these layers.

Figure 1 Daily EC e at the average soil depth of 15 cm. 

Figure 2 Daily EC e at the average soil depth of 45 cm. 

Figure 3 Yearly average θ at the average soil depth of 15 cm. 

CONCLUSIONS

The simulations show that soil salinity increases (median and maximum of daily EC e ) for all soil layers, in the future climate scenarios comparatively to the reference scenario. Also, the average yearly θ of the topsoil decreases in the future scenarios, compared to the reference. Both changes in salinity and water content may limit the productivity of the pasture and even compromise the ability of plants to grow. This evidences the need for adaptation strategies to avoid degradation of these soils due to salinity.

Acknowledgements

This work was carried out with support from the projects Saltfree (ARIMNET 2/000572015) and Soil4Ever (PTDC/ASP-SOL/28796/2017), both financed by FCT.

References

Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. & FAO (1998) - Crop evapotranspiration: Guidelines for computing crop water requirements. FAO. [ Links ]

Fischer, G.; Nachtergaele, F.O.; Prieler, S.; Teixeira, E.; Toth, G.; van Velthuizen, H.; Verelst, L. & Wiberg, D. (2012) - Global Agro-ecological Zones (GAEZ v3.0)- Model Documentation. FAO. [ Links ]

Hargreaves, G.H. & Samani, Z.A. (1985) - Reference Crop Evapotranspiration from Temperature. Applied Engineering in Agriculture, vol. 1, n. 2, p. 96-99. https://doi.org/10.13031/2013.26773Links ]

Kanzari, S.; Ben Nouna, B.; Ben Mariem, S. & Rezig, M. (2018) - Hydrus-1D model calibration and validation in various field conditions for simulating water flow and salts transport in a semi-arid region of Tunisia. Sustainable Environment Research, vol. 28, n. 6, p. 350-356. https://doi.org/10.1016/j.serj.2018.10.001Links ]

Meijgaard, E.; Bosveld, F. & Lenderink, G. (2008) - The KNMI regional atmospheric climate model RACMO version 2.1. Undefined. [ Links ]

Ramos, T.; Gonçalves, M.; Martins, J. & Pereira, L. (2014) - Comparison of different pedotransfer function for estimating soil hydraulic properties in Portugal. [ Links ]

Ramos, T.; Šimůnek, J.; Gonçalves, M.C.; Martins, J.C.; Prazeres, A.; Castanheira, N.L. & Pereira, L.S. (2011) - Field evaluation of a multicomponent solute transport model in soils irrigated with saline waters. Journal of Hydrology, vol. 407, n. 1-4, p. 129-144. https://doi.org/10.1016/j.jhydrol.2011.07.016Links ]

Rhoades, J.D.; Manteghi, N.A.; Shouse, P.J. & Alves, W.J. (1989) - Estimating Soil Salinity from Saturated Soil-Paste Electrical Conductivity. Soil Science Society of America Journal, vol. 53, n. 2, p. 428-433. https://doi.org/10.2136/sssaj1989.03615995005300020019xLinks ]

Ritchie, J.T. (1972) - Model for predicting evaporation from a row crop with incomplete cover. Water Resources Research, vol. 8, n. 5, p. 1204-1213. https://doi.org/10.1029/WR008i005p01204Links ]

Rodrigues, M.; Fortunato, A.B. & Freire, P. (2019) - Saltwater Intrusion in the Upper Tagus Estuary during Droughts. Geosciences, vol. 9, n. 9, art. 400. https://doi.org/10.3390/geosciences9090400Links ]

Simunek, J.; Sejna, M. & van Genuchten, M.Th. (2018) - Hydrus 1D (4.17.0140) [Computer software]. PC-Progress. [ Links ]

Soares, P.M.M.; Cardoso, R.M.; Lima, D.C.A. & Miranda, P.M.A. (2017) - Future precipitation in Portugal: High-resolution projections using WRF model and EURO-CORDEX multi-model ensembles. Climate Dynamics, vol. 49, p. 2503-2530. https://doi.org/10.1007/s00382-016-3455-2Links ]

Wesseling, J.G.; Elbers, J.A.; Kabat, P. & Van den Broek, B.J. (1991) - SWATRE: Instructions for input. Internal Note, Winand Staring Centre, Wageningen, the Netherlands. [ Links ]

(*E-mail: ana.paz@iniav.pt)

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License