Comparing Two Saline-Gypseous Wetland Soils in NE Spain

: Small (<1 km 2 ) saline wetlands scattered across the landscape often go unnoticed or are threatened by urbanization or other interventions, despite their role as biodiversity shelters. This study is needed to show methods for monitoring this speciﬁc kind of wetland, and to guide the selection of analytical techniques. We provide data and comparisons for salient soil traits of two quasi-pristine gypsiferous and saline wetlands named Farrachuela (FA) and Agust í n (AG). The soil characteristics presented in this article are a more sensitive indicator of their ecological status than some of the most used indicators, such as birds and plants. We found signiﬁcant differences between the two saladas in percent water saturation, equivalent calcium carbonate, gypsum content, and soil salinity expressed as electrical conductivity both of 1:5 soil-to-water ratio and of saturation extracts. The differences were also signiﬁcant in the concentrations of Mg 2+ , Na + , and Cl − , while they were non-signiﬁcant for Ca 2+ , HCO 32 − , and SO 42 − . The mean contents of the six ions were lower in FA than in AG. Both pH and sodium adsorption ratios were signiﬁcantly different between the two wetlands. The data are mainly examined and plotted by displaying their non-parametric statistics, a synoptic approach that will allow us to monitor the evolution of the wetlands against both traditional agricultural pressures and emerging green energy infrastructures. Last but not least, we discuss the shortcomings of some standard laboratory methods when applied to gypsum-rich soils.


Introduction
European inland saline habitats have been quasi-eliminated [1], with many of the surviving ones reduced in size or badly degraded.This is also the case in the arid Central Ebro Basin, Spain [2].The scarcity of regulations or incentives for protecting soil diversity occurs in several European countries [3] and is more pronounced for inland saline wetlands in general.
This scarcity also happens for the lands whose composition, color, vegetation or other characteristics are due to an abundance of gypsum (CaSO 4 •2H 2 O).These whitish gypsum soil landscapes, known by local people in NE Spain as chesas, occur mainly in warm countries around the world [4], but also in cold regions with gyprock outcrops, e.g., [5,6].Chesas have also been traditionally disdained due to their low agricultural production and geotechnical characteristics.Irrigation of these lands, needed in dry areas for agricultural production, is technically unfeasible by inundation due to the solubility of gypsum, and only the pressurized application of water has alleviated the drawbacks.The disdain is counteracted by the agricultural pressure and by the increasing attention of scientists to the gypsophilous plants protected by environmental rules [7].Thus, saline wetlands located in gypseous soils pose a double challenge to life, and are therefore of twofold interest to science.Moreover, the study and description of these saline wetlands face methodological issues in their lab analyses-often overlooked in the literature-derived from the abundance of gypsum [4,8].
The old judgment on both the wetlands and the gypseous soils based on their little or no agricultural value is evolving towards legal protection due to the biodiversity they Land 2023, 12,1990 2 of 17 harbor and their role in ecological balances.In addition, insights into biodiversity and ecosystem services raise the need to move beyond single-lake thinking [9] and to include the small wetlands [10], a line of reasoning that also applies to gypseous soils.Together with the coexistence with agriculture, including irrigation and intensive animal husbandry, there now arise threats from the "green energy" sources with several kinds of impacts, such as those recently reported in [11,12].
The Ebro Basin, located in northeastern Spain (Figure 1), has a dry climate and a high presence of gypsum rocks 1.9 10 6 hectares of gypsiferous lithofacies [13] together with many landscape evidences of salinity, such as saline ravines associated with springs or oozing sources, and especially many endorheic saline wetlands and saline depressions.Some of these wetlands had artisanal salterns until one century ago.Nowadays, these athalassohaline wetlands -locally named saladas-are hot spots for nature conservation due to the rare or singular habitats they host.Scientific interest in these wetlands can be traced back to the 19th century.The data set [14] lists many works describing the physical conditions, geology, soils, vegetation, and other features of the saladas of Sástago-Bujaraloz in the Central Ebro Basin, an endorheic complex included in the RAMSAR list.Salada Agustín is included in this complex.By contrast, salada Farrachuela located in the chesas of the Barbastro Gypsum Formation, about 70 km NE of salada Agustín (Figure 1), is isolated and went unnoticed by scientists until recently.This kind of wetland is found in areas with arid or Mediterranean climates, but is also comparable to wetlands developed in gypsum karst depressions in cold climates [6].
Land 2023, 11, x FOR PEER REVIEW 2 of 17 The old judgment on both the wetlands and the gypseous soils based on their little or no agricultural value is evolving towards legal protection due to the biodiversity they harbor and their role in ecological balances.In addition, insights into biodiversity and ecosystem services raise the need to move beyond single-lake thinking [9] and to include the small wetlands [10], a line of reasoning that also applies to gypseous soils.Together with the coexistence with agriculture, including irrigation and intensive animal husbandry, there now arise threats from the "green energy" sources with several kinds of impacts, such as those recently reported in [11,12].
The Ebro Basin, located in northeastern Spain (Figure 1), has a dry climate and a high presence of gypsum rocks 1.9 10 6 hectares of gypsiferous lithofacies [13] together with many landscape evidences of salinity, such as saline ravines associated with springs or oozing sources, and especially many endorheic saline wetlands and saline depressions.Some of these wetlands had artisanal salterns until one century ago.Nowadays, these athalassohaline wetlands -locally named saladas-are hot spots for nature conservation due to the rare or singular habitats they host.Scientific interest in these wetlands can be traced back to the 19th century.The data set [14] lists many works describing the physical conditions, geology, soils, vegetation, and other features of the saladas of Sástago-Bujaraloz in the Central Ebro Basin, an endorheic complex included in the RAMSAR list.Salada Agustín is included in this complex.By contrast, salada Farrachuela located in the chesas of the Barbastro Gypsum Formation, about 70 km NE of salada Agustín (Figure 1), is isolated and went unnoticed by scientists until recently.This kind of wetland is found in areas with arid or Mediterranean climates, but is also comparable to wetlands developed in gypsum karst depressions in cold climates [6].Protection of the European saladas should be implemented throughout the Natura 2000 network of Special Conservation Areas (SCA).Farrachuela is located within SCA "ES2410074 Yesos de Barbastro", but the salada went unnoticed in the data form.Agustín, located within SCA "ES2430082 Monegros", has been protected by excluding it and its surrounding lands from irrigation.
Several investigations, e.g., [15], ratify the traditional knowledge that salinity together with intermittent flooding are the main conditioners of life in the saladas.On the other hand, many published reports and investigations on saline wetlands overlook the shortcomings of some routine concepts and methods unsuited for analyzing gypsum-rich materials, like the soils or sediments of the saladas studied here.
This paper aims to: (i) compare certain compositional features, especially the contents of the major ions relevant for life, at two gypseous and hypersaline wetlands located in the Ebro Basin, NE Spain (Figure 1); and (ii) discuss the adequacy of some common basic operations in the analytical determinations of the gypseous soils.Protection of the European saladas should be implemented throughout the Natura 2000 network of Special Conservation Areas (SCA).Farrachuela is located within SCA "ES2410074 Yesos de Barbastro", but the salada went unnoticed in the data form.Agustín, located within SCA "ES2430082 Monegros", has been protected by excluding it and its surrounding lands from irrigation.
Several investigations, e.g., [15], ratify the traditional knowledge that salinity together with intermittent flooding are the main conditioners of life in the saladas.On the other hand, many published reports and investigations on saline wetlands overlook the shortcomings of some routine concepts and methods unsuited for analyzing gypsum-rich materials, like the soils or sediments of the saladas studied here.
This paper aims to: (i) compare certain compositional features, especially the contents of the major ions relevant for life, at two gypseous and hypersaline wetlands located in the Ebro Basin, NE Spain (Figure 1); and (ii) discuss the adequacy of some common basic operations in the analytical determinations of the gypseous soils.
Beyond the scientific interest of the above objectives, we emphasize their usefulness for defining the living conditions of the protected organisms harbored by these wetlands, and monitoring for undesirable alterations.

The Context of the Wetlands Studied
We compare the isolated wetland Farrachuela (FA) to Agustín (AG) which belongs to the Sástago-Bujaraloz complex.The two hypersaline wetlands lie about 70 km apart (Figure 1), and both of them occupy Quaternary depressions formed by dissolution of rocks rich in gypsum and more soluble salts, with the remaining residuum in their bottoms.A conspicuous network of surficial cracks appearing after each rain episode demonstrates the inundation, while the fresh plant residues visible to the naked eye evidence the recent biological activity in the soil.Contrariwise, the microcrystalline gypsum, i.e., the flour-like gypsum described elsewhere [4], is not observed within the two saladas studied, but does occur in the areas surrounding them.This fact mirrors the contrasting soil hydric regime of the almost permanently dry soils in the Central Ebro Basin against the persistently moist or wet soils in the enclaved saladas.
The two saladas occupy flat-bottomed valleys, which allows for preferential water circulation and accumulation.FA stands isolated in the landscape at the outcrops of the Barbastro Gypsum Formation from the Eocene-Oligocene age [16], whose dipping strata facilitate deep percolation of rainfall.AG is part of the largest group of saline wetlands in the Ebro Basin, with a total of 149 after the inventory of [17], a number that is in stark contrast with the absence of wetlands around FA.The near-horizontal Miocene strata of gypsum and limestone interbedded with lutites of the Bujaraloz Formation in Monegros [18] allow the occurrence of karstic depressions aligned along flat-bottomed valleys and low-lying areas hosting wetlands.Other differentiating factors are the presence of shallow saline aquifers (aquitards) in Monegros [19], responsible for groundwater discharge into the intermittently flooded depressions, and the wind regime that influences the excavation and shaping of the playa lakes in this area.
Winter cereals were cultivated in the watersheds of both saladas (Figures A1 and A2 in the Appendix A) at the time of the samplings, i.e., July 2013 for Farrachuela; and February, March, April, July, and August of the years 1979-1980 and 1999-2000 for Agustín.The plants collected (Table 1) are well known as salinity and waterlogging tolerant or as gypsophilous, and they illustrate the ecology of these wetlands.The floristic composition has remained stable up to the present.The aerial photographs from 1956 of both saladas attest that irrelevant anthropic interventions occurred prior to the years of the soil samplings presented in this study.We deem AG representative of those saladas of Sástago-Bujaraloz that remained in an acceptable state of conservation, even threatened by flooding with fresh water from forthcoming irrigation.The land around FA, cleared more than two hundred years ago, bears barley and sparse olive trees, with no irrigation planned.

Climate
Climate data were obtained from the nearest automatic weather stations, in Tamarite, located 11.5 km south of Farrachuela, and Valfarta located 8.5 km NNW of Agustín.Both stations belong to the Spanish SIAR network.For the period 2004-2021, the mean annual rainfall in the Farrachuela area was 350 mm, the mean temperature was 14 • C, and the mean evapotranspiration rate was 1041 mm.The same mean temperature and precipitation of 359 mm were obtained in Agustín area for the same period, with a higher mean evapotranspiration rate (1258 mm).
Using the terminology of the Soil Survey Staff [22] and the data from the two weather observatories, the soil moisture regime at both saladas is aquic and the most likely soil temperature regimes are mesic for Farrachuela and thermic for Agustín.

Sampling
We collected soil samples by hand augering down to a depth of 2 m in Farrachuela, and until reaching an impenetrable layer in Agustín.All the samples were examined visually, including with a magnifying glass, and by touch both in the field and in the laboratory.Then, the samples were air-dried and passed through a 2 mm sieving mill to obtain the fine earth fraction, i.e., <2 mm Ø, used for further analyses.Table 2 shows details of the augerings, and the number of analyses.

Analyses
The official methods of the Agricultural Ministry of Spain [23] provided the methodology for the chemical determinations.Calcium carbonate (CaCO 3 ) equivalent (CCE) was measured by gasometry, and gypsum was titrated as per [24].
All the determinations of electrical conductivity are reported as dS m −1 at 25 • C. We measured the electrical conductivity (EC1:5) of the aqueous extracts at an earth-to-water ratio weight:weight of 1:5.Moreover, we prepared water-saturated pastes of the fine earth [25], recording the water saturation percentage (SP, %) and measuring the electrical conductivity of the extract (ECe), and its pH (pHe).The major ions in the saturation extracts were determined using the classical methods specified in [20].
We report the ionic contents (ionic-C, mmolc L −1 ) as the sum of cations + anions divided by two in order to counterbalance the deviations to electroneutrality, frequent when analyzing highly concentrated solutions.The sodium adsorption ratio SAR (mmolc L −1 ) 0.5 is calculated from the concentrations in mmolc L −1 of Na + , Ca 2+ , and Mg 2+ by the expression SAR = Na + /[(Ca 2+ + Mg 2+ )/2] 0.5 .

Statistical Procedures
We display and analyze most of our data through graphical methods, representing the measures of exploratory data analysis by means of the intuitive boxplots proposed by [26].Our boxplots, drawn following [27], include the 95% confidence intervals for the medians estimated as per [28].Figure A3 in Appendix A sketches the meanings of all the parts of the boxplot diagrams presented in this paper.
We calculate the regression lines by the Ordinary Least Squares (OLS) method with p = 0.05, and we check their equations with a non-parametric simple regression using Theil's method [29,30] in order to cope with the habitual non-Gaussian distributions of our data, the frequent outliers, and the limited effectiveness of some data batches.

Saturation Percentages
The mean SP in the 16 samples from FA with this determination is 59.0%, far from the 39.3% in the 133 samples from AG.The other statistics and the distributions of SP also are very distinct between the two saladas (Figure 2).The use of SP to compare soil characteristics is sounder than the classical particle size distribution (PSD) determinations, which, in some materials, produce artifactual results, as argued in Section 3.9.

Statistical Procedures
We display and analyze most of our data through graphical methods, representing the measures of exploratory data analysis by means of the intuitive boxplots proposed by [26].Our boxplots, drawn following [27], include the 95% confidence intervals for the medians estimated as per [28].Figure A3 in Appendix A sketches the meanings of all the parts of the boxplot diagrams presented in this paper.
We calculate the regression lines by the Ordinary Least Squares (OLS) method with p = 0.05, and we check their equations with a non-parametric simple regression using Theil's method [29,30] in order to cope with the habitual non-Gaussian distributions of our data, the frequent outliers, and the limited effectiveness of some data batches.

Saturation Percentages
The mean SP in the 16 samples from FA with this determination is 59.0%, far from the 39.3% in the 133 samples from AG.The other statistics and the distributions of SP also are very distinct between the two saladas (Figure 2).The use of SP to compare soil characteristics is sounder than the classical particle size distribution (PSD) determinations, which, in some materials, produce artifactual results, as argued in Section 3.9.The saturated paste extraction approximates the field capacity of the soil much better than other, more diluted, soil-to-water ratios, while making water extraction feasible with simple equipment [31,32].Thus, the extract of the saturated paste is broadly used [22,23,25] to assess the salt tolerance of plants by measuring the electrical conductivity of the extract (ECe).Also, eventually, its ionic composition is determined to appraise the effects of individual ions or their co-occurrence, as is the case of the structural stability of the soil.Many studies of soil salinity report SP because of its interest in understanding the hydric behavior of the soil.The simplicity of the preparation of the saturated paste, based on the skill of the operator, and the saturation point determined by feel and experience, together with the use of few-if any-sophisticated laboratory tools, probably contribute to SP being overlooked in some studies.
However, even if the composition of saturated extracts was not to be analyzed, the SP is worth recording because SP can be a surrogate for textural composition, or "a The saturated paste extraction approximates the field capacity of the soil much better than other, more diluted, soil-to-water ratios, while making water extraction feasible with simple equipment [31,32].Thus, the extract of the saturated paste is broadly used [22,23,25] to assess the salt tolerance of plants by measuring the electrical conductivity of the extract (ECe).Also, eventually, its ionic composition is determined to appraise the effects of individual ions or their co-occurrence, as is the case of the structural stability of the soil.Many studies of soil salinity report SP because of its interest in understanding the hydric behavior of the soil.The simplicity of the preparation of the saturated paste, based on the skill of the operator, and the saturation point determined by feel and experience, together with the use of few-if any-sophisticated laboratory tools, probably contribute to SP being overlooked in some studies.
However, even if the composition of saturated extracts was not to be analyzed, the SP is worth recording because SP can be a surrogate for textural composition, or "a quantitative expression of soil texture" as noted by [33], and represents key functional soil properties, like the available water capacity [34].This is an old issue in soil science when "dealing with single-valued expressions instead of trying to interpret the complex series of numbers represented by the mechanical analysis" [35], as in the class-moisture equivalent diagram of [36].When dealing with soil salinity, SP has been used in recent attempts [37] at estimating ECe from extracts at fixed dilution ratios.Also, SP is a relevant soil characteristic for precision agriculture, as noted in [38], with SP performing better than organic matter, CaCO 3 , ECe, and pH for yield prediction.Furthermore, SP can help to check consistency over time in sample preparation [32], which we also did successfully to compare two operators [39].
The surrogate role of SP for texture is especially useful in gypsum-rich soils where the standard methods for PSD determination are unsuitable [4,40] and many substitute methods are unfeasible in routine labs.In [8], hand textural class estimates were used to check the proposed method for PSD determination.After the above considerations, we did not determine PSD, and we claim SP as a strong indicator of hydric behavior and a reliable surrogate for textural composition [40-42].

Major Mineral Components
The boxplots of calcium carbonate equivalent (CCE) and gypsum for each salada (Figure 3) show their distribution along the samples, and the prominence of both components.The mean percent rounded to the nearest integer for FA and AG are 24% and 10% for CCE while for gypsum they are 30% and 62%, respectively.Then, the remaining material, i.e., silica and silicates-including the mineralogical clays-plus organic components is <50% on average.These compositions agree with the contrasting average SP of 59% in FA, and 39% in AG (Figure 2), since gypsum and calcium carbonate have a much lower water-holding capacity than the other minority components, like clay and organic matter.
quantitative expression of soil texture" as noted by [33], and represents key functional soil properties, like the available water capacity [34].This is an old issue in soil science when "dealing with single-valued expressions instead of trying to interpret the complex series of numbers represented by the mechanical analysis" [35], as in the class-moisture equivalent diagram of [36].When dealing with soil salinity, SP has been used in recent attempts [37] at estimating ECe from extracts at fixed dilution ratios.Also, SP is a relevant soil characteristic for precision agriculture, as noted in [38], with SP performing better than organic matter, CaCO3, ECe, and pH for yield prediction.Furthermore, SP can help to check consistency over time in sample preparation [32], which we also did successfully to compare two operators [39].
The surrogate role of SP for texture is especially useful in gypsum-rich soils where the standard methods for PSD determination are unsuitable [4,40] and many substitute methods are unfeasible in routine labs.In [8], hand textural class estimates were used to check the proposed method for PSD determination.After the above considerations, we did not determine PSD, and we claim SP as a strong indicator of hydric behavior and a reliable surrogate for textural composition [40-42].

Major Mineral Components
The boxplots of calcium carbonate equivalent (CCE) and gypsum for each salada (Figure 3) show their distribution along the samples, and the prominence of both components.The mean percent rounded to the nearest integer for FA and AG are 24% and 10% for CCE while for gypsum they are 30% and 62%, respectively.Then, the remaining material, i.e., silica and silicates-including the mineralogical clays-plus organic components is <50% on average.These compositions agree with the contrasting average SP of 59% in FA, and 39% in AG (Figure 2), since gypsum and calcium carbonate have a much lower water-holding capacity than the other minority components, like clay and organic matter.The abundance of gypsum in the soils and geological materials governs the chemistry of the soil solution and the exchange complex due to their saturation in Ca 2+ .The contents of gypsum (Figure 3) in the samples analyzed-with a minimum of 12.2% in FA and 9.5% in AG-largely guarantee Ca 2+ saturation, preventing clay dispersion.The abundance of gypsum in the soils and geological materials governs the chemistry of the soil solution and the exchange complex due to their saturation in Ca 2+ .The contents of gypsum (Figure 3) in the samples analyzed-with a minimum of 12.2% in FA and 9.5% in AG-largely guarantee Ca 2+ saturation, preventing clay dispersion.

Soil Salinity
In the present study, a common visual symptom of soil salinity is the white bright efflorescences mostly produced by evapoconcentration during dry periods at the bottoms and margins of both FA and AG.These efflorescences are often composed mainly of crystals of gypsum and more soluble salts, as evidenced in the field by naked-eye and hand-lens observations, as well as by taste.The halophilous vegetation also attests the soil salinity.
The visual symptoms do not quantify the salt stock in the soils.Such stock is a key indicator of the ecological status of the soil, which can be tracked by recording the changes in the salt content through time.To establish this stock, the soil salt must be extracted at fixed-ratio dilution, higher than the saturated paste [25,31] as is evident for hypersaline soils [43].To appraise the salt stock of the soil, we use the easy, unsophisticated, and reproducible determinations of EC1:5 [25].Such attributes make EC1:5 well-suited for future comparative assessments.
The mean of EC1:5 in 32 extracts of FA was 7.18 dS m −1 , versus 8.78 dS m −1 in the 72 extracts from AG for which EC1:5 was analyzed (Figure 4).We consider these conductivities a proxy for the total content of highly soluble salts that can allow tracking the evolution of the soil's salt stock.The 1.68 dS m −1 of the difference between the means of EC1:5 seems moderate.However, AG is a more salt-stressing environment than FA, as shown by the vegetation (Table 1) and by the mean ECe of 37.04 dS m −1 for FA against 68.58 dS m −1 for AG.The difference is more pronounced than in EC1:5, and the same is true for the other statistics graphed in Figure 4.The longer range of the distributions in AG seems sound if considering the 20-fold surface area of AG versus FA (Table 2), which allows for more variable conditions in AG.

Soil Salinity
In the present study, a common visual symptom of soil salinity is the white bright efflorescences mostly produced by evapoconcentration during dry periods at the bottoms and margins of both FA and AG.These efflorescences are often composed mainly of crystals of gypsum and more soluble salts, as evidenced in the field by naked-eye and hand-lens observations, as well as by taste.The halophilous vegetation also attests the soil salinity.
The visual symptoms do not quantify the salt stock in the soils.Such stock is a key indicator of the ecological status of the soil, which can be tracked by recording the changes in the salt content through time.To establish this stock, the soil salt must be extracted at fixed-ratio dilution, higher than the saturated paste [25,31] as is evident for hypersaline soils [43].To appraise the salt stock of the soil, we use the easy, unsophisticated, and reproducible determinations of EC1:5 [25].Such attributes make EC1:5 well-suited for future comparative assessments.
The mean of EC1:5 in 32 extracts of FA was 7.18 dS m −1 , versus 8.78 dS m −1 in the 72 extracts from AG for which EC1:5 was analyzed (Figure 4).We consider these conductivities a proxy for the total content of highly soluble salts that can allow tracking the evolution of the soil's salt stock.The 1.68 dS m −1 of the difference between the means of EC1:5 seems moderate.However, AG is a more salt-stressing environment than FA, as shown by the vegetation (Table 1) and by the mean ECe of 37.04 dS m −1 for FA against 68.58 dS m −1 for AG.The difference is more pronounced than in EC1:5, and the same is true for the other statistics graphed in Figure 4.The longer range of the distributions in AG seems sound if considering the 20-fold surface area of AG versus FA (Table 2), which allows for more variable conditions in AG.The much higher values of EC1:5 and ECe in AG than in FA (Figure 4) agree with the species of halophytes recorded at both saladas (Table 1).The absence of Arthrocnemum macrostachyum in FA is relevant because this plant lives in the more often inundated areas of AG where ECe is >80 dS m −1 , a soil salinity not reached in FA.The difference is also pronounced if only the upper soil layer is considered, with about 37 dS m −1 in FA versus 85 dS m −1 in AG.Accordingly, A. macrostachyum, the most salinity perennial tolerant plant, lives in AG and not in FA.
The scatterplot of ECe on EC1:5 (Figure 5) does not show the inflection associated with gypsum [43]   The much higher values of EC1:5 and ECe in AG than in FA (Figure 4) agree with the species of halophytes recorded at both saladas (Table 1).The absence of Arthrocnemum macrostachyum in FA is relevant because this plant lives in the more often inundated areas of AG where ECe is >80 dS m −1 , a soil salinity not reached in FA.The difference is also pronounced if only the upper soil layer is considered, with about 37 dS m −1 in FA versus 85 dS m −1 in AG.Accordingly, A. macrostachyum, the most salinity perennial tolerant plant, lives in AG and not in FA.
The scatterplot of ECe on EC1:5 (Figure 5) does not show the inflection associated with gypsum [43] because all samples have EC > 2.25 dS m −1 in both EC1:5 and ECe.The wide scattering of the values-seen in the plot-hampers the estimation of ECe from the plain EC1:5.If desired, estimates of ECe can be made from EC1:5 by the method in [44].In this article, we apply the easy transformation of EC1:5 by SP as per [45], i.e., by calculating regressions of the shape ECe = a + b × EC1:5 × (500/SP) q with empirical powers of q.The best correlation coefficients (Table A1, in Appendix A) were 0.863 for FA, and 0.928 for AG.
scattering of the values-seen in the plot-hampers the estimation of ECe from the plain EC1:5.If desired, estimates of ECe can be made from EC1:5 by the method in [44].In this article, we apply the easy transformation of EC1:5 by SP as per [45], i.e., by calculating regressions of the shape ECe = a + b × EC1:5 × (500/SP) q with empirical powers of q.The best correlation coefficients (Table A1, in Appendix A) were 0.863 for FA, and 0.928 for AG.For the present article, the equations allowed us to check the coherence between ECe and EC1:5.The future studies of these wetlands will have to decide whether or not to estimate ECe from EC1:5 with this kind of equation.The reason is that ECe is the standard indicator of saline stress on plants, used to express a plant's tolerance to soil salinity, a key feature for crop feasibility [32].The equations in Table A1 allowed us to check the consistency between ECe and EC1:5 and will help to decide in future studies of these wetlands whether or not to estimate ECe from EC1:5 with these kinds of equations.

Soil Classification
The inspection of the auger samples and their analyses for ECe and gypsum content enabled us to identify the diagnostic horizons defined by the Soil Survey Staff [22] without the need to open a trench, a much more soil disturbing procedure.According to the data of the ECe and gypsum content, and following step-by-step the statements in [22], the soils of FA and AG have the Salic and Gypsic horizon in conjunction.Thus, the soils of both wetlands are Gypsic Aquisalids.

Ionic Contents of the Saturation Extracts
Bicarbonate, titrated in all the extracts, was always below the limits of detection.Potassium contents are irrelevant.The pairs of boxplot diagrams in Figure 6 compare the content of each of the ions Ca 2+ , Mg 2+ , Na + , HCO3 − , SO4 2− , and Cl − in the saturation extracts of FA to AG.The range of all the distributions is broader for AG, which also happens for the interquartile ranges.
Of note are: (i) the greater surface area of AG and the greater number of samples than in FA, and (ii) the greater surface area of AG allowing for more diverse hydric conditions including the redistribution of water by the wind.These conditions are mirrored by the occurrence of the highly halophilic A. macrostachyum in AG.Both saladas are similar with regard to the predominance of SO4 2− and Mg 2+ , but the concentrations of Cl − and Na + are significantly lower in Farrachuela, hence its lower salinity.
Figure 6 shows that the median concentrations are significantly lower in FA than in AG for the six ions in question, but the difference is insignificant for sulfate.The mean For the present article, the equations allowed us to check the coherence between ECe and EC1:5.The future studies of these wetlands will have to decide whether or not to estimate ECe from EC1:5 with this kind of equation.The reason is that ECe is the standard indicator of saline stress on plants, used to express a plant's tolerance to soil salinity, a key feature for crop feasibility [32].The equations in Table A1 allowed us to check the consistency between ECe and EC1:5 and will help to decide in future studies of these wetlands whether or not to estimate ECe from EC1:5 with these kinds of equations.

Soil Classification
The inspection of the auger samples and their analyses for ECe and gypsum content enabled us to identify the diagnostic horizons defined by the Soil Survey Staff [22] without the need to open a trench, a much more soil disturbing procedure.According to the data of the ECe and gypsum content, and following step-by-step the statements in [22], the soils of FA and AG have the Salic and Gypsic horizon in conjunction.Thus, the soils of both wetlands are Gypsic Aquisalids.

Ionic Contents of the Saturation Extracts
Bicarbonate, titrated in all the extracts, was always below the limits of detection.Potassium contents are irrelevant.The pairs of boxplot diagrams in Figure 6 compare the content of each of the ions Ca 2+ , Mg 2+ , Na + , HCO 3 − , SO 4 2− , and Cl − in the saturation extracts of FA to AG.The range of all the distributions is broader for AG, which also happens for the interquartile ranges.
Of note are: (i) the greater surface area of AG and the greater number of samples than in FA, and (ii) the greater surface area of AG allowing for more diverse hydric conditions including the redistribution of water by the wind.These conditions are mirrored by the occurrence of the highly halophilic A. macrostachyum in AG.Both saladas are similar with regard to the predominance of SO 4 2− and Mg 2+ , but the concentrations of Cl − and Na + are significantly lower in Farrachuela, hence its lower salinity.
Figure 6 shows that the median concentrations are significantly lower in FA than in AG for the six ions in question, but the difference is insignificant for sulfate.The mean concentrations are also lower in FA than in AG for the ions analyzed.This happens when computing all the analyses from AG extracts (Figure 6), and also if computing only the extracts with the six ions analyzed (Table 3).This Table also shows numerically the differences between the mean contents of each individual ion in the two saladas.The greatest difference occurs in Cl − , followed by Na + and SO 4 2− , and then by Mg 2+ , with 594, 410, 317, and 224 mmolc L −1 , respectively.We deem irrelevant the differences in Ca 2+ and Land 2023, 12,1990 9 of 17 HCO 3 − .The differences between the two wetlands for Mg 2+ , Na + , SO 4 2− , and Cl − agree with their distinct physical characteristics outlined in the second paragraph of Section 2.1.Of note is the detection of magnesium by taste in the efflorescences of AG, agreeing with the boxplot in Figure 6.The abundance of magnesium should be due to evapoconcentration from the parental rocks of gypsum that contain this element.
concentrations are also lower in FA than in AG for the ions analyzed.This happens when computing all the analyses from AG extracts (Figure 6), and also if computing only the extracts with the six ions analyzed (Table 3).This Table also shows numerically the differences between the mean contents of each individual ion in the two saladas.The greatest difference occurs in Cl − , followed by Na + and SO4 2− , and then by Mg 2+ , with 594, 410, 317, and 224 mmolc L −1 , respectively.We deem irrelevant the differences in Ca 2+ and HCO3 − .The differences between the two wetlands for Mg 2+ , Na + , SO4 2− , and Cl − agree with their distinct physical characteristics outlined in the second paragraph of Section 2.1.Of note is the detection of magnesium by taste in the efflorescences of AG, agreeing with the boxplot in Figure 6.The abundance of magnesium should be due to evapoconcentration from the parental rocks of gypsum that contain this element.Figure 6.Pairs of boxplots for the six main ions in all the analyzed saturation extracts from Farrachuela (FA) and Agustín (AG).The red squares mark the means, for the meaning of all parts, see Figure A3.Table 3. Mean contents (mmolc L −1 ) of six major ions for: (a) the 16 saturation extracts from FA, (b) all the analytical determinations from AG as shown in Figure 6, (c) differences between these means, (d) means of the 76 extracts from AG with the six ions titrated, and (e) differences between the means of all determinations from AG and the means of FA.

Ca 2+
Mg 2+ Na + Cl − SO4 2− HCO3 − mmolc L −1 (a) Farrachuela Table 3. Mean contents (mmolc L −1 ) of six major ions for: (a) the 16 saturation extracts from FA, (b) all the analytical determinations from AG as shown in Figure 6, (c) differences between these means, (d) means of the 76 extracts from AG with the six ions titrated, and (e) differences between the means of all determinations from AG and the means of FA.Table 3 also shows the consistency of the results obtained by computing all the available titrations or only the ones from the samples with the six ions analyzed.This feature is interesting for long-term tracking, when the homogeneity of the number of ions analyzed along successive campaigns cannot be guaranteed.

SAR and pH
Due to their salinity and frequent inundation, the saladas are either bald or seasonally populated by halophytes, and are unusable for the cash crops feasible in the Ebro Basin.Notwithstanding, we present (Table 4) the statistics of SAR and pH of the saturation extracts in both wetlands, showing their significant difference.The values of SAR in Farrachuela are well below the classical threshold of SAR = 13 for clay dispersion in soils.By contrast, the values of SAR in Agustín are largely in excess of 13, denoting the unfeasibility of cash crops.The pH values are significantly different between the two wetlands (Table 4); however, pH by itself should not be a problem for most plants, as it ranges from 7.28 to 8.44 in the two saladas.These basic pH values are allowable for most plants.For the infiltration behavior, SAR must be considered jointly with pH and salinity [45].For SAR, it does not make sense to apply the commonly used threshold SAR < 13 (mmol c L −1 ) 0.5 to judge the sustainability of irrigation, since for most soil samples the ECe is >>20 dS m −1 (Table 4).Then, the SAR/EC ratios do not affect the clay dispersion and no associated imperviousness would happen [45,46].Moreover, the contents of mineralogical clay are low, as the sum of the means of gypsum and calcium carbonate equivalent are 54% for Farrachuela, and 72% for Agustín.The mineralogical clays are thus in the minority, reducing the value of SAR to predict the hydrological behavior of these gypsum-rich soils often saturated in water.More relevant to the hydraulic behavior is the horizontal microstructure in FA and AG, similar to the layering described in other saladas in Monegros.

Relationship between Ionic Concentration and ECe
Figure 7A shows the concordance between both saladas regarding the relation of ionic contents (ionic-C, mmolc L −1 ) in the extracts, as defined in Section 2.4., versus ECe.The tusk-shaped distribution in this figure illustrates the increase in the ratio of ions to ECe for the high ionic concentrations, due to the occurrence of neutral and other ionic pairs [47].The dilution of the more concentrated extracts, required for the analytical methods, also contributes to the scattering of values in the high concentrations.
For the two saladas, the distributions approach a straight line after transforming both variables using decimal logarithms (Figure 7B), as proposed in [48].This reinforces the consistency of the data.Table A2 in the Appendix A shows the regression equations obtained from these transformations.

Relationship between the Concentration of Individual Ions and ECe
The saturation extracts of AG showed a good correlation coefficient (R) of the ions between them as well as between the ECe and the ions (Table 5).These coefficients in AG are higher than their counterparts in FA, except the correlation of ECe with SO 4 2− , which is slightly higher in FA than in AG.

Relationship between the Concentration of Individual Ions and ECe
The saturation extracts of AG showed a good correlation coefficient (R) of the ions between them as well as between the ECe and the ions (Table 5).These coefficients in AG are higher than their counterparts in FA, except the correlation of ECe with SO4 2− , which is slightly higher in FA than in AG.For the saturation extracts of FA, the correlation of ECe was significant (p < 0.001) with Mg 2+ , Na + , and SO4 2− , but not with Cl − .By contrast, in AG the correlation ECe/Cl − was the highest, as shown in Table 5.
The coefficients of determination (R 2 ) attained by the OLS regression of the average ionic content versus ECe (Table 5) are the same as when using Theil calculations, while the standard errors using the OLS regression are smaller than with Theil regressions.The  For the saturation extracts of FA, the correlation of ECe was significant (p < 0.001) with Mg 2+ , Na + , and SO 4 2− , but not with Cl − .By contrast, in AG the correlation ECe/Cl − was the highest, as shown in Table 5.
The coefficients of determination (R 2 ) attained by the OLS regression of the average ionic content versus ECe (Table 5) are the same as when using Theil calculations, while the standard errors using the OLS regression are smaller than with Theil regressions.The OLS method is thus recommended for estimating the ionic content from the ECe in the saturation extracts of both saladas within the ranges indicated in Table 6.These equations could help to streamline the environmental tracking of the saladas.Table 6.Regression equations of the shape "ionic-C = a + b × ECe" with the coefficient of determination (R 2 , %) and standard error (S) calculated for 16 and 76 soil samples from Farrachuela and Agustín, respectively.The methods of calculation are ordinary least squares (OLS), and Theil using the median of the interceptors.

The Specificity in the Analytics of Gypseous Materials
In fresh soil samples of the gypseous horizons of the soils studied, the abundant gypsum crystals are bound together by fine material, the strength of which increases as the moisture content decreases.This is a well-known characteristic of many soils [49].The coarse gypsum crystals are easily separated, intact, using the fingers when the sample is wet, but when air dried, these crystals break easily -mostly along cleavage plains-before becoming detached from the mass.The breaking happens even if a dry soil sample is gently crushed between a wooden board and a wooden rolling pin, and, of course, in mechanical shakers and sieving mills, until they eventually pass through the sieve.Milling may result in the entire sample passing through the 2mm sieve, making the time at which the operation is interrupted arbitrary.For this reason, the content of coarse gypsum crystals and other gypsum fragments must be appraised before drying and sieving in the lab.
In gypsum-rich soils, the low hardness and cohesion of the gypsum crystals lessens the value and significance of the particle size separates obtained after milling.Moreover, the abundance of gypsum and its solubility (≈2.4 g L −1 ) immediately cause the sedimentation of clays when being dispersed for particle size separations [4,8].The judgments about other methods for PSD determinations based on X-ray fluorescence, e.g., [50] or its combination with FT-NIR spectroscopy, e.g., [51], or with Vis-NIR and pXRF [52] are unsound if these methods are calibrated against either pipette or wet sieving and hydrometer, which do not work in gypseous soils.The method proposed by [8], based on sonication in a 7:3 ethanol to water solution-a procedure available in few labs-allows dispersion of fine earth, but the significance of the coarse and fine earth remains arguable in the soils studied.In the common circumstances of gypsum-rich soils, the lab determination of PSD should not be taken as representing the behavior of the undisturbed soil.These shortcomings led to coining two terms in lieu of texture: "Coarse Gypsum Material" and "Fine Gypsum Material" [49].
The issues discussed in the above two paragraphs indicate the non-sense of using some routine analytical methods with gypseous sediments and soils.Then, [4] proposed the elimination of PSD for assessing the hydric behavior of gypseous soils.Specifically, we support considering saturation percentage as a sound and useful surrogate for textural composition.Accordingly, we did not separate the clay fraction, whose determination would also be impractical, cumbersome and probably irrelevant for future monitoring.

Conclusions
The article provides a snapshot of the main limiters to life in two inland saline marshes.We hope it provides a baseline for protecting of these shelters of biodiversity.The intuitive graphical procedures used in this article allow for a quick and detailed statistical comparison of the soil characteristics that are decisive for life in the saline wetlands studied.The features compared are: (i) contents of gypsum and calcium carbonate; (ii) soil salinity as EC1:5; (iii) saturation percentage; and (iv) pH, electrical conductivity and ionic contents in the saturation extracts.The data will constitute a baseline -the only one at present-for tracking possible compositional changes in the two saladas studied.The calculation methods and graphical presentation are easily transposable to other saline wetlands, and may allow for multiple comparisons in the future.
The characteristics studied herein ought to be contrasted with other saline wetlands in order to delimit the environmental conditions for the life of specific plants, microbes, or other organisms.Hopefully, the correlations presented in this article will support the management and conservation measures to be implemented by environmental authorities in charge of the wetlands.
In situ recording of variable soil parameters is advisable.For some of them, the technology is available and is continuously improving.However, the implementation of recording seems unlikely in the near future due to the limited capacity of the conservation agencies responsible.
The misuse of analytical techniques that are wrong or unsound for gypsum-rich soils, like the ones discussed in the present article, is a key point to be: (i) avoided when conducting future analyses, and (ii) checked when scrutinizing heritage data.

Figure 1 .
Figure 1.The shaded area is the Ebro Basin, where the two saline wetlands compared, Farrachuela (FA) and Agustín (AG), are located on gypsiferous substrates.

Figure 1 .
Figure 1.The shaded area is the Ebro Basin, where the two saline wetlands compared, Farrachuela (FA) and Agustín (AG), are located on gypsiferous substrates.

Figure 2 .
Figure 2. Distribution and main statistics of SP in Farrachuela and Agustín.The red squares mark the means; for the meaning of all parts, see Figure A3.

Figure 2 .
Figure 2. Distribution and main statistics of SP in Farrachuela and Agustín.The red squares mark the means; for the meaning of all parts, see Figure A3.

Figure 3 .
Figure 3. Boxplots of CCE and gypsum in the analyzed samples from Farrachuela and Agustín, with the number (n) of samples computed.The red squares mark the means; for the meaning of all parts, see Figure A3.

Figure 3 .
Figure 3. Boxplots of CCE and gypsum in the analyzed samples from Farrachuela and Agustín, with the number (n) of samples computed.The red squares mark the means; for the meaning of all parts, see Figure A3.

Figure 4 .
Figure 4. Boxplots of EC1:5 and ECe for the two saladas compared.The red squares mark the means; for the meaning of all parts, see Figure A3.
because all samples have EC > 2.25 dS m −1 in both EC1:5 and ECe.The wide

Figure 4 .
Figure 4. Boxplots of EC1:5 and ECe for the two saladas compared.The red squares mark the means; for the meaning of all parts, see Figure A3.

7 *Figure 6 .
Figure 6.Pairs of boxplots for the six main ions in all the analyzed saturation extracts from Farrachuela (FA) and Agustín (AG).The red squares mark the means, for the meaning of all parts, see Figure A3.

Figure 7 .
Figure 7. (A) Scatterplot of the ionic contents over ECe in the saturation extracts; (B) Log-log scatter plot of the same parameters.

Figure 7 .
Figure 7. (A) Scatterplot of the ionic contents over ECe in the saturation extracts; (B) Log-log scatter plot of the same parameters.

Figure A2 .
Figure A2.Salada Agustín, 12 October 1979 at 4p.m.In the foreground, a plowed field, followed by vegetation of Suaeda vera.Behind it the area covered by Arthrocnemum macrostachyum, often flooded, around the bottom occupied by water.
Figure A3.Sketch of a boxplot.

Figure A2 . 17 Figure A2 .
Figure A2.Salada Agustín, 12 October 1979 at 4p.m.In the foreground, a plowed field, followed by vegetation of Suaeda vera.Behind it the area covered by Arthrocnemum macrostachyum, often flooded, around the bottom occupied by water.

Table 2 .
Location and size of the two wetlands and data of their sampling.

Table 4 .
Statistics of SAR and pH in the saturation extracts of samples from Farrachuela (FA) and Agustín (AG).

Table 5 .
Correlation coefficients between ECe and the main ions, and between them.The p-value is < 0.001, except if non-significance (n.s.) is indicated.

Table 5 .
Correlation coefficients between ECe and the main ions, and between them.The p-value is < 0.001, except if non-significance (n.s.) is indicated.

Table A1 .
Pearson's correlation coefficients (R) obtained with OLS regression equations of the shape ECe = a + b × EC1:5 × (500/SP) q for several arbitrary values of q.

Table A1 .
Pearson's correlation coefficients (R) obtained with OLS regression equations of the shape ECe = a + b × EC1:5 × (500/SP) q for several arbitrary values of q.
Figure A3.Sketch of a boxplot.

Table A1 .
Pearson's correlation coefficients (R) obtained with OLS regression equations of the shape ECe = a + b × EC1:5 × (500/SP) q for several arbitrary values of q.

Table A2 .
Regression equations of the shape log ionic-C = a + b × log ECe calculated by OLS and Theil methods for 16 and 76 soil samples from Farrachuela and Agustín, respectively.