Next Article in Journal
Distribution of Subsurface Nitrogen and Phosphorus from Different Irrigation Methods in a Maize Field
Previous Article in Journal
Impacts of Climate Change on Groundwater in the Al-Badan Sub-Catchment, Palestine: Analyzing Historical Data and Future Scenarios
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrated Approach to the Regional Estimation of Soil Moisture

by
Luis Pastor Sánchez-Fernández
1,*,
Diego Alberto Flores-Carrillo
1 and
Luis Alejandro Sánchez-Pérez
2
1
Centro de Investigación en Computación, Instituto Politécnico Nacional, Mexico City CP 07738, Mexico
2
Department of Electrical and Computer Engineering, University of Michigan, Dearborn, MI 48128, USA
*
Author to whom correspondence should be addressed.
Hydrology 2024, 11(10), 170; https://doi.org/10.3390/hydrology11100170
Submission received: 19 September 2024 / Revised: 8 October 2024 / Accepted: 9 October 2024 / Published: 11 October 2024

Abstract

:
Automatic or smart irrigation systems benefit irrigation water management. However, measurement sensor networks in automatic irrigation systems are complex, and maintenance is essential. Regional soil moisture estimation avoids the multiple measurements necessary when deploying an irrigation system. In this sense, a fuzzy estimation approach based on decision-making (FEADM) has been used to obtain soil moisture point estimates. However, FEADM requires intelligent weather adjustment based on spatial features (IWeCASF) to perform regional soil moisture estimation. The IWeCASF-FEADM integrated approach for regional soil moisture estimation is developed in this work. IWeCASF provides the inputs for FEADM. FEADM is performed R times; R is the number of checkpoints at which a point estimate is obtained. In this way, regional estimation is achieved when the set of R soil moisture point estimates is completed. Additionally, IWeCASF-FEADM considers the irrigation water records, which are not included in either method individually. This method can detect when the soil moisture is deficient in a region, allowing actions to prevent water stress. This regional estimation reduces an irrigation system’s operational and maintenance complexity. This integrated approach has been tested over several years by comparing the results of regional soil moisture estimation with measurements obtained at many points in the study region.

1. Introduction

Currently, enhancing irrigation water consumption is vital in appropriately exploiting agricultural resources. As a result, improved automatic irrigation systems [1,2,3,4], or intelligent systems for water management, have been developed [5,6,7,8].
A set of checkpoints is commonly deployed throughout an irrigation region when an automatic irrigation system is used [9,10]. At each checkpoint, soil moisture measurements determine the irrigation water requirements. Therefore, a measurement sensor network is needed [11]. The soil moisture determined within the irrigation region based on measurements can provide a very accurate representation; nevertheless, the inherent complexity of its implementation and maintenance is prominent [12,13].
An intelligent model for the estimation of soil moisture can reduce the complexity involved in measuring soil moisture and, consequently, the costs of automatic irrigation systems. Hydrological models can estimate soil moisture based on land surface models and the weather conditions [14,15,16]. These models can accurately estimate soil moisture; however, their development complexity is high due to the difficulties in acquiring the input data and the assumptions performed during model development. Moreover, soil moisture can be obtained with soil water balance calculations, which are easy to implement while being very accurate, because these calculations require regular recalibration [17,18,19].
Soil moisture remote sensing can be performed at any time without being affected by cloud cover, vegetation or night conditions; it is achieved using satellite microwave sensors [20,21]. Nevertheless, the scale (~25 to 50 km) and the complexity in implementing these microwave sensors render soil moisture remote sensing an alternative that is only suitable for huge farmland systems. Earth observation technologies are essential in monitoring sustainable application development at different scales; these include soil moisture estimations [22], soil moisture mapping services [23], and soil moisture estimation, combining satellite radar measurements and land surface modeling [24], among others.
Artificial neural networks (ANNs) [25,26,27] and support vector machines (SVMs) [28,29] are used for soil moisture estimation, and they have different advantages. These models enable accurate estimation but present restrictions for implementation in small farms. An SVM performs better than an ANN when an ensemble Kalman filter (EnKF) is utilized [30,31]. However, both approaches require many tests to obtain an accurate soil moisture estimate [32,33].
Fuzzy logic is a powerful technique for analyzing hydrologic elements and decision making in water resources. Many issues connected to hydrology frequently deal with imprecision and vagueness, and they can be very well handled by fuzzy logic-based models [34]. Other works are related to optimizing irrigation systems; the authors propose using several water microturbines to generate sustainable energy without affecting crop irrigation efficiency [35]. In this sense, sustainable and energy-efficient agricultural practices require integrating different systems that utilize water adequately, including small irrigation regions or greenhouses [36].
A model was developed to estimate the soil point moisture based on decision-making using a fuzzy inference system (FEADM) [32]. FEADM obtains point estimates of the soil moisture from the measurements of certain variables, such as the air temperature, rain, solar radiation, wind speed, and evapotranspiration, using a wireless environmental monitoring station that uses radio frequency signals to transmit data (details in [32]).
FEADM models the relationships among the soil moisture, weather conditions, and checkpoint features where the estimation is performed. FEADM is suitable for regional soil moisture estimation when obtaining point soil moisture estimates at several regional checkpoints. Nevertheless, FEADM requires the weather conditions at every checkpoint, but measuring the weather conditions at several checkpoints within a region is impractical.
The weather conditions remain similar within a specific region [37,38], although its landscape features can lead to inconsistencies. Many analyses and experiments support the significance of such discrepancies in the regional weather conditions. In this sense, an intelligent weather condition adjustment based on spatial features (IWeCASF) [39] was developed. This model permits the weather conditions to be obtained anywhere within a region, since IWeCASF considers the weather condition inconsistencies mentioned above. IWeCASF requires the measured weather conditions at a single checkpoint, and, from these measurements, IWeCASF determines the weather conditions at any other checkpoint by performing intelligent adjustment. This intelligent adjustment considers the relationships between the weather conditions and the checkpoints’ spatial features.
Moreover, both methods consider the weather conditions and soil and crop features to be essential to determine the soil moisture content.
Nevertheless, as with conventional irrigation systems, the integrated approach aims to determine the soil moisture content at several checkpoints scattered across an irrigation region. Thus, it eliminates the need for the full implementation of the conventional irrigation system and mitigates its maintenance complexity [12,13] by estimating the regional soil moisture instead of measuring it.
Measuring the weather conditions at several checkpoints located in a region can be more complex than measuring the soil moisture as in conventional irrigation systems. For example, in a region, as depicted in Figure 1, measuring the weather conditions at each R = 12 checkpoints, P r = 1 , 2 , , R , is more complex than simply measuring the soil moisture at each checkpoint, P r , due to the requirements of the measuring devices. However, the weather conditions remain similar within a region [37,38]. For example, given an irrigation region, as shown in Figure 1, the air temperature can be highly similar at checkpoints P 0 , P 4 , and P 8 . The same can be true for other weather conditions, such as solar radiation, which can even be similar at checkpoints P0, P 4 , P 8 , and P 11 , which are located relatively far apart.
Thus, the weather conditions anywhere within a specific region can be determined by measuring only one checkpoint (the primary checkpoint). Moreover, despite the similarity of the weather conditions, there can be minor variations due to the spatial features of the region. For example, in the region in Figure 1, the air temperature is almost the same at most of the checkpoints, as described in the previous example; however, there could be a noticeable air temperature difference between checkpoint P 4 , which is a grassland area, and the checkpoints P 2 or P 10 , which are mainly tree-covered areas. The tree-covered areas at checkpoints P 2 or P 10 partially block solar radiation at the soil level, which can result in a cooler air temperature. Therefore, considering that the spatial features influence the weather conditions, these variations within a specific region can be determined if a spatial analysis is performed. This is the aim of the “intelligent weather condition adjustment based on spatial features” (IWeCASF) model [39].
In brief, the IWeCASF-FEADM integrated approach utilizes IWeCASF to determine the weather conditions at any checkpoint where a point soil moisture estimation is performed with FEADM. Figure 1 demonstrates an application of the integrated approach for regional soil moisture estimation. Firstly, the weather conditions C i = 1 , 2 , , I 0 are measured at the primary checkpoint ( P 0 ), and then the measured weather conditions C i = 1 , 2 , , I 0 and the rest of the inputs required by IWeCASF are processed to determine the adjusted weather conditions C i = 1 , 2 , , I r at each checkpoint P r . The adjusted weather conditions C i = 1 , 2 , , I r for checkpoint P r and the irrigation water record are introduced to FEADM to perform the point soil moisture estimation. In this way, the set of point soil moisture estimates at each checkpoint P r = 1 , 2 , , R is defined as the regional soil moisture estimation.
The locations of checkpoints P r = 1 , 2 , , R are determined by the number of zones identified in conventional irrigation systems, but, if a more detailed regional soil moisture estimate is required, more checkpoints must be defined. This approach is suitable for use in different regions; it is only necessary to modify the input parameters.
The theoretical framework presents mathematical formulations and symbols supporting the analyses. However, they do not limit other readers since the work includes many descriptions, detailed figures, diagrams, and tables.
This work hypothesizes that linking FEADM with ICASF allows regional soil moisture estimation (in a specific region). Likewise, the IWeCASF-FEADM approach for regional soil moisture estimation can reduce the implementation and maintenance complexity of conventional automatic irrigation systems and their inherent costs.
In this study, the integrated approach to regional soil moisture estimation (called the integrated approach) is developed. This approach utilizes the joint IWeCASF-FEADM to perform a regional analysis (set of point estimates) of soil moisture; additionally, the irrigation water supply is introduced to FEADM to complement the rainfall for the processing of a broader set of water supply records.
Based on mathematical formulations, the integrated approach determines a regional estimation of the soil moisture anywhere within a region through image geoprocessing algorithms and fuzzy inference systems. The intelligent weather condition adjustment through spatial features (IWeCASF) [39] and the soil moisture fuzzy point estimation model based on decision-making (FEADM) [32] exploit the advantages of evaluating soil moisture with FEADM but eliminate its limitations.

2. Materials and Methods

Figure 2 depicts an overview of the integrated IWeCASF-FEADM approach for the regional estimation of soil moisture. The integrated method developed herein accomplishes regional soil moisture estimation ( R e ) from a set of soil moisture point estimates ( S m r = 1 , 2 , , R   ). The integrated approach requires as inputs the satellite imagery of the estimation region, as well as the data of the soil and crop zones ( M t = 1 , 2 , , T ) identified within this region. Moreover, the locations of the checkpoints ( P r = 0 , 1 , 2 , , R ), the weather conditions measured ( C i = 1 , 2 , , I 0 ) at one of the checkpoints (primary checkpoint P 0 ) and the water irrigation record ( Δ r = 1 , 2 , , R ) are needed. The integrated approach links the IWeCASF model [39] and the fuzzy estimation approach based on decision-making (FEADM) [32]. IWeCASF determines the weather conditions ( C i = 1 , 2 , , I r ) at every checkpoint ( P r ), which are utilized by FEADM to obtain each soil moisture point estimate ( S m r ). The set of soil moisture point estimates ( R e = S m 1 , S m 2 , ,   S m R ) of all checkpoints ( P r ) represents the regional estimation of soil moisture ( R e ).

2.1. Inputs

The integrated approach conforms with the IWeCASF-FEADM linkage; therefore, certain inputs are required by these two models. As defined in [39], the information necessary to develop IWeCASF includes the satellite imagery, soil and crop zones, checkpoint locations, and measured weather conditions; meanwhile, FEADM [32] only needs the weather conditions as inputs.
The satellite imagery corresponds to images (with a scale of 1:10,000) obtained with the USGS National Map Viewer from Landsat in March 2017. QGIS [40] is used for processing. The images are divided into S sectors, resulting in pairwise sectors s x , y , where x = 1 , 2 , , X , y = 1 , 2 , , Y , X is the image width, and Y is the image height.
The soil and crop zones constitute a geographic data set organized in two databases. The first has the region’s geo-information divided into zones concerning the soil texture, field capacity and available water content [39,41,42]. The second database includes the region’s geo-information divided into zones according to the crop type and its development stage [39,43]. Both databases are utilized to define the feature matrices M t = 1 , 2 , , T , where T is the total number of features identified; in this case, T = 3 ( M t = 1 soil type, M t = 2 crop type, and M t = 3 development stage). The size of the feature matrices is X x Y , defined by the S sectors established in the satellite imagery. Additionally, elements m x , y t of feature matrices M t contain the number of pixels at sector s x , y that belong to a feature [39].
Moreover, the checkpoint location s x r , y r is another input; it denotes the sector s x , y where checkpoint P r is located, and it is used to identify the particular features at this location. Additionally, the input weather conditions C i = 1 , 2 , , I 0 are measured at the primary checkpoint P 0 with an integrated sensor suite [32], as defined in Table 1.
The input irrigation water is a matrix Δ = δ x r , y r X Y obtained from a geographic database; it denotes the amount of irrigation water supplied in each sector s x , y from the irrigation region. This input, as a complement to the rain data C 2 r , is used in the integrated approach to determine the actual water supply in a sector s x r , y r . The external subscripts X and Y are as follows: X is the image width and Y is the image height.

2.2. IWeCASF-FEADM Integrated Approach

The IWeCASF-FEADM integrated approach is composed of two stages. The first stage is IWeCASF [39]. Firstly, IWeCASF’s results are combined with the irrigation water record; afterwards, they are supplied as input to the second stage. The second stage is FEADM [32], which determines the soil moisture point estimates at each checkpoint. The linkage of IWeCASF and FEADM gives a regional soil moisture estimation as a point estimate set. Section 2.2.1 and Section 2.2.2 present the details of the integration of IWeCASF and FEADM, respectively.

2.2.1. IWeCASF

IWeCASF [39] determines the weather conditions within a region, providing FEADM [32] with the inputs required for a point estimate. This goal is achieved by performing three tasks: landscape feature extraction, selecting specific landscape features at each checkpoint, and adjusting the weather conditions at each checkpoint. IWeCASF’s inputs are the satellite imagery of the estimation region, the soil and crop zones, the checkpoint locations, and the weather conditions measured at the primary checkpoint. As mentioned above, the results derived from IWeCASF (the weather conditions at each checkpoint) are utilized as inputs by FEADM. The five stages (A to F) are described as follows:
A. Landscape feature extraction: This is the stage of IWeCASF in which an image is processed to define the landscape feature matrices F l = 1 , 2 , , L of the irrigation region’s satellite imagery. The landscape feature matrices F l = 1 , 2 , , L are necessary to obtain a steady adjustment factor. L is the total number of landscape features; in this case, L = 5 (grassland, tree cover area, buildings, elevation, and spatial configuration). First, the image is enhanced through color space conversion from red, green, and blue (RGB) to the International Commission of Lighting Luminosity band A/band B color space (CIELAB) [44]. Afterwards, a decorrelation process is implemented to simplify the image color segmentation. Then, the landscape feature matrices F l are defined, where each matrix represents a landscape feature identified from the image. Color segmentation assigns each pixel of the image to one prototype cluster K l = 1 , 2 , , L , which represents a landscape feature, using the objective function G defined in Equation (1), as described in [44].
G = l = 1 L a d ϵ K l a d b l 2
When color segmentation is completed, the image is separated pixel by pixel into L = 5 landscape features, which consist of the grassland ( l = 1), tree cover area ( l = 2 ), buildings ( l = 3 ), elevation ( l = 4 ), and spatial configuration ( l = 5). However, to define the landscape feature matrices, the image is divided into sectors s ( x , y ) , where x = 1 , 2 , ,   X and y = 1 , 2 , , Y . X is the width and Y is the height of the satellite image divided into S sectors. Thus, each element f x , y l = ( n u m b e r   o f   p i x e l s ) of a landscape feature matrix F l = f x , y l X Y denotes the number of pixels at sector s ( x , y ) corresponding to the landscape feature l .
B. Checkpoint landscape feature selection: In this stage, the input locations s x r , y r of checkpoints P r = 1 , 2 , , R and the feature matrices M t = 1 , 2 , , T = m x , y t X x Y are utilized. Additionally, the landscape feature matrices F l = 1 , 2 , , L = f x , y l X Y extracted from the image are required. The checkpoint location s x r , y r selects from its respective matrices the elements f x = x r ,   y = y r l and m x = x r ,   y = y r t that represent the landscape feature l or the feature t at sector s x = x r ,   y = y r . As a result of this selection, the particular feature vector Φ r F l , M t , P r is obtained, as shown in Equation (2). The elements of this vector Φ r F l , M t , P r are defined as the number of pixels at sector s x = x r , y = y r . Then, the elements of the particular feature vector Φ r F l , M t , P r are represented as a percentage of the particular feature vector Φ r F l , M t , P r , as described in Equation (3).
Φ r F l , M t , P r = [ f x r , y r 1 , f x r , y r 2 , , f x r , y r L , m x r , y r 1 , m x r , y r 2 , , m x r , y r T ]
Φ r F l , M t , P r = [ ϕ e = 1 r , ϕ e = 2 r , , ϕ e = L r , ϕ e = L + 1 r , ϕ e = L + 2 r , , ϕ e = E r ]
C. Fuzzy adjustment: This is performed to determine the weather conditions C i = 1 , 2 , , I 0 at checkpoint P r after the particular feature vector Φ r F l , M t , P r is identified. Fuzzy adjustment begins with the fuzzification of the weather conditions measured, C i = 1 , 2 , , I 0 , and of the landscape feature vector Φ r F l , M t , P r . A sigma function, an L-shaped function, and a triangle function are used for fuzzifying, as defined in [39]. The fuzzy weather conditions C ~ i 0 ( v i ) , where v i is the current value of weather condition C i 0 , are of the form defined in Equations (4)–(8).
C ~ 1 0 ( v i ) = c ~ 1 , 1 0 , c ~ 1 , 2 0 , , c ~ 1 , 5 0         Temperature
C ~ 2 0 ( v i ) = c ~ 2 , 1 0 , c ~ 2 , 2 0 , , c ~ 2 , 5 0         Rain
C ~ 3 0 ( v i ) = c ~ 3 , 1 0 , c ~ 3 , 2 0 , c ~ 3 , 5 0         Solar   radiation
C ~ 4 0 v i = c ~ 4 , 1 0 , c ~ 4 , 2 0 , , c ~ 4 , 5 0         Wind   speed
C ~ 5 0 ( v i ) = c ~ 5 , 1 0 , c ~ 5 , 2 0 , , c ~ 5 , 5 0         Evapotranspiration
Moreover, to fuzzify the landscape feature vectors Φ r F l , M t , P r , it is necessary to compare the landscape features Φ r F l , M t , P r at the checkpoint where the adjustment is performed and the landscape features at the primary checkpoint Φ 0 F l , M t , P 0 . This comparison Φ 0 , r is executed using Equation (9).
Φ 0 , r = Φ 0 F l , M t , P 0 Φ r F l , M t , P r = [ ϕ e = 1 0 ϕ e = 1 r , ϕ e = 2 0 ϕ e = 2 r , , ϕ e = E 0 ϕ e = E r ]
Then, the elements of comparison Φ 0 , r are fuzzified using membership functions, i.e., the form triangle, L-shaped or sigma-shaped function, as described in [39]. In this way, the fuzzy particular feature matrix Φ ~ r is obtained, as shown in Equation (10).
Φ ~ r = ϕ ~ 1 , 1 r ϕ ~ 1 , 2 r ϕ ~ 1 , 3 r ϕ ~ 2 , 1 r ϕ ~ 2 , 2 r ϕ ~ e , 3 r ϕ ~ E , 1 r ϕ ~ E , 1 r ϕ ~ E , 3 r
D. Landscape adjustment: This stage models the influence of the particular checkpoint features Φ r F l , M t , P r over the weather conditions measured C i = 1 , 2 , , I 0 . Landscape adjustment uses a fuzzy inference system (FIS), which receives as inputs the fuzzy particular feature matrix Φ ~ r . This FIS utilizes a compound of n = 1 , 2 , , N  IF THEN rules, where N = 245 when modeling the adjustment factor according to the actual particular feature differences between the checkpoint P r and the primary checkpoint P 0 . The landscape adjustment FIS presents as output a steady adjustment factor element α r , i |   i = 1 , 2 , , I , where I = 5 is the number of weather conditions. Each output is separated into λ i , o   |   o = 1 , 2 , , O output sets, with O = 5 according to [39]. The result of landscape adjustment is a weather condition steady adjustment factor α r , i for each weather condition, which is gathered as Λ r = α r , i 1 , I .
E. Variable adjustment: Before developing the variable adjustment, it is necessary to determine the certainty of the weather condition replication η i = 1 , 2 , , I , which improves the modeling of inconsistent weather conditions such as rain ( C 2 0 ). This certainty η i = 1 , 2 , , I is obtained using a distribution function, as in [39]. Afterwards, the certainty of the weather condition replication η i = 1 , 2 , , I and the weather conditions measured C i = 1 , 2 , , I 0 are processed through the Q = 445  IF THEN rules of a second FIS (variable adjustment). Its result is associated with the output membership functions’ variables β i   |   i = 1 , 2 , , I to model the influence of the weather conditions over the rest. As a result of defuzzification, a variable adjustment B factor is obtained.
F. Final adjustment: This gives the result of the weather conditions C i = 1 , 2 , , I r at checkpoint P r . Both the factors steady adjustment Λ r and variable adjustment B are used to determine the suggested adjustment of each weather condition C i r at checkpoint P r . According to [39], the adjusted weather conditions C i r are obtained as shown in Equation (11).
C i r = ( C i 0 ) ( 1 + α r , i + β i )     |   r = 1 , 2 , , R ; i = 1 , 2 , , I

2.2.2. Addition of Irrigation Water

IWeCASF’s final result is the weather conditions C i = 1 , 2 , , I r at any checkpoint within a region. In the first stage of the integrated IWeCASF-FEADM approach developed herein, these weather conditions C i = 1 , 2 , , I r , are used as the inputs of FEADM (second stage). A particular case is the weather condition C i = 2 r (rain), which is complemented before being introduced to FEADM. Thus, the actual weather conditions C ˇ i r supplied to FEADM are defined in Equation (12), in which the irrigation water record δ x r , y r of checkpoint P r is added to weather condition C i = 2 r (rain) at checkpoint P r . Meanwhile, the rest of the weather conditions remain the same.
C ˇ i r = C i r + δ x r , y r i = 2 C i r i 2

2.2.3. FEADM

The integrated IWeCASF-FEADM approach continues with the second FEADM stage [32]. FEADM is utilized to estimate the soil moisture content (point estimation of soil moisture S m r ) in a delimited area (checkpoint P r ) as from the weather conditions determined with IWeCASF.
Furthermore, unlike conventional decision-making methods, FEADM not only selects the best decision alternative (qualitative result) but also estimates the actual value of an alternative when available (quantitative result). The aim is to select, from a set of soil moisture levels (defined in FEADM as alternatives A k = 1 , 2 , , K ), the best-assessed alternative according to the actual weather conditions C i = 1 , 2 , , I (defined in FEADM as decision criteria C i = 1 , 2 , , I ) and then estimate its actual value. Therefore, FEADM is executed at each checkpoint P r , in which the actual weather conditions C ˇ i = 1 , 2 , , I r derived from IWeCASF and modified with the addition of irrigation water are required as inputs. FEADM’s result is a point soil moisture estimate S m r at checkpoint P r . The three stages (A to C) are described as follows:
A. FEADM begins with a fuzzy analytical hierarchy process [32,45,46,47], which compares the relevance of each decision criterion (adjusted weather conditions derived from IWeCASF) over the rest. This comparison is performed for each decision criterion. The method developed by [45] and adapted by [46] is used in FEADM. As a result, a weight set W = w 1 , w 2 , ,   w I is obtained.
Firstly, a judgment matrix P ~ is defined as in Equation (13). The elements p ~ x , y are triangle fuzzy numbers (TFN) that represent the relative importance values, which are derived from criteria pairwise comparison C a x , C a y , where x = 1 , 2 , , I and y = 1 , 2 , , I .
P ~ = C 1     C 2         C I     C 1 C 2 C I p ~ 1,1 p ~ 1 , 2 p ~ 1 , I p ~ 2 , 1 p ~ 2 , 2 p ~ 2 , I p ~ I , 1 p ~ I , 2 p ~ i , I
Then, the relative importance values p ~ x , y are utilized in Equation (14) to obtain the synthetic extent value S ~ i = 1 , 2 , , I of criteria C i = 1 , 2 , , I [45].
S ~ i = 1 , 2 , , I = n = 1 I p ~ i , n i = 1 I n = 1 I p ~ i , n 1 i = 1 , 2 , , I
Afterwards, the total integral value with the index of optimism I T S ~ i = 1 , 2 , , I [32,46,47,48] is used in Equation (15) to calculate the weight set W = w 1 , w 2 , , w I .
w i = I T S ~ i n = 1 I I T S ~ n i = 1 , 2 , , I
B. The expert evaluation starts with the fuzzification of the criteria C i = 1 , 2 , , I . J = 5 linguistic labels are associated with membership functions μ C i , j z i   |   j = 1 , 2 , , J , where z i Z i is the adjusted weather condition value at the checkpoint under analysis, and Z i is its universe of discourse. The membership functions μ C i , j z i   |   j = 1 , 2 , , J are used to fuzzify each criterion C i . In this way, each fuzzified criterion is defined as C ~ i = μ C i , 1 ( z i ) , , μ C i , 5 ( z i ) .
Afterwards, the fuzzified criteria vectors C ~ i = 1 , 2 , , I are gathered in the fuzzy criteria matrix C ~ , as shown in Equation (16). Each element c i , j of matrix C ~ represents a membership function μ C i , j z i . The rows of the fuzzy criteria matrix C ~ are a combination of the membership values of criteria C ~ i = 1 , 2 , , I .
C ~ = c 1 , 1 c 1 , J c I , 1 c I , J
After obtaining the fuzzy criteria matrix C ~ , weighting is developed to assign greater relevance to those criteria that are more important. Subsequently, the weight set W = w 1 , w 2 , , w I obtained in FAHP is used to weight elements c i , j of fuzzy criteria matrix C ~ , as defined in Equation (17). All weighted elements c i , j w are combined in the fuzzy weighted criteria matrix C ~ w shown in Equation (18).
c i , j w = c i , j w i
C ~ w = c 1,1 w c 1 , J w c I , 1 w c I , J w
The next step of the expert evaluation is the selective assessment, which is responsible for modeling the relationships among the fuzzy weighted criteria c i , j w . However, the selective assessment determines the criteria c i , j w , which are relevant for the assessment of each alternative A k [32].
Furthermore, the selective assessment defines a set of evaluations E k = e 1 k , e 2 k , , e N k k , which contains the result e n = 1 , 2 , , N k k defined in Equation (19). The number of evaluations N k = 1 , 2 , , K is resolved by the decision makers according to the fuzzy weighted criteria c i , j w selected for the assessment of each alternative A k = 1 , 2 , , K .
e n = 1 , 2 , , N k = f k n , C ~ w
where the selective assessment function f k n , C ~ w   |   n = 1 , 2 , , N k , defined in Equation (20), assesses the fuzzy weighted criteria c i , j w   ϵ   C ~ w selected for the evaluations e n = 1 , 2 , , N k k of alternative A k = 1 , 2 , , K .
f k n , C ~ w = min D n k = d i , j k | i = 1 , 2 , , I ; j = 1 , 2 , , J
in which the elements d i , j k of the relationship modeling matrix D n k determine whether a fuzzy weighted criterion c i , j w is used ( d i , j k = c i , j w ) or not used ( d i , j k = 1 ) in the assessment n of alternative A k . Moreover, the relationship modeling matrix, defined by the decision makers, depends on the relationships among the weather conditions, the spatial features of the checkpoint, the crop features, and the soil moisture content.
After assessing the whole set E k , the best assessment for each alternative A k is obtained with Equation (21).
a k = m a x E k
Once the best assessments a k are obtained for every alternative A k , they are combined in a matrix a 1 a K T . The best choice  is achieved through Equation (22).
 = max a 1 a K
C. Fuzzy inference: As FEADM does not only select a qualitative value (soil moisture level) but can also estimate a qualitative value (soil moisture value), a fuzzy inference process is included in the second stage.
An output membership function μ A k = 1 , 2 , , K is defined for each alternative A k , as in [32]. Moreover, the best assessment a k determines the maximum value for each output membership function μ A k = 1 , 2 , , K ; in this way, the actual output membership function μ A k is obtained, as shown in Equation (23).
μ A k = μ A k ; μ A k < a k a k ; μ A k a k
Then, the function μ S m e is obtained using Equation (24). This function is utilized in the centroid method, shown in Equation (25). As a result, the soil moisture estimates S m r at checkpoint P r are achieved.
μ S m e = m a x μ A 1 , μ A 2 , , μ A 5
S m r = l = 0 L μ S m e z s m l z l l = 0 L μ S m e z s m l

2.3. Regional Estimation (Point Estimates of Soil Moisture)

The final result of FEADM, the second stage of the integrated approach, is a point estimate of the soil moisture. Additionally, in the integrated IWeCASF-FEADM approach, the set of checkpoints is scattered within a region, as in conventional systems. Thus, the soil moisture at these checkpoints must be determined. As a consequence, FEADM, which is responsible for obtaining the soil moisture, is executed at every checkpoint P r = 1 , 2 , , R . Consequently, FEADM is implemented a total of R times, giving, as a result, soil moisture point estimates S m r = 1 , 2 , , R   . The soil moisture point estimates S m r = 1 , 2 , , R   are combined in a set, which is named the regional soil moisture estimation R e = S m 1 , S m 2 , ,   S m R . If a more accurate regional soil moisture estimation is required, more checkpoints P r must be defined to reduce the area without a point estimate.

3. Experimental Case Study

As presented above, the regional soil moisture estimation R e = S m 1 , S m 2 , ,   S m R is a set of soil moisture point estimates S m r = 1 , 2 , , R   , obtained based on two stages. The first stage is IWeCASF, which is necessary to supply the weather conditions C i = 1 , 2 , , I r to FEADM, and the second stage is the integrated approach, IWeCASF-FEADM.

3.1. Inputs

As described in Section 2.1 and [39], the satellite imagery of the estimation region is divided into S = 2150 sectors s ( x , y ) . The total number of sectors, S = 2150 , is determined from the experimental results so that the extension of each sector s ( x , y ) is 40 x 40   m . The estimation region is part of a basin with a length of almost 2 km. The different areas (grassy, tree-covered, and building areas) within the estimation region make it suitable for the exploration of diverse spatial features.
The satellite image in Figure 3 depicts the region of interest, the locations s ( x r , y r ) of some checkpoints P r and the feature matrix M 2 (crop type). As described in [39] and Section 3.1, the feature matrices M t are given in the input soil and crop zones. The feature matrix M 2 is composed of elements ( m x , y 1 ) X x Y . In this case, the feature matrix M 2 denotes two different crop types. The crop type influences the soil moisture estimation approach utilized in the integrated IWeCASF-FEADM approach.
Additionally, with the integrated IWeCASF-FEADM approach, there can be as many checkpoints P r as required, because this approach does not require the physical deployment of measurement equipment to obtain the soil moisture. Therefore, the more checkpoints there are, the more accurate the soil moisture estimation is, because fewer sectors (x,) remain without a soil moisture estimate. In this experiment, each sector s x , y is considered a checkpoint. The irrigation water matrix δ x r , y r is utilized to complement the rain record and corresponds to the irrigation water supplied at each sector s ( x , y ) within the estimation region.
Moreover, IWeCASF [39], the first stage of the integrated approach, requires the weather conditions measured C i = 1 , 2 , , I 0 at a single checkpoint (primary checkpoint P 0 ), which can also be identified in Figure 3. The weather condition C i = 1 , 2 , , I 0 and values v i = 1 , 2 , , I used in this experiment are presented in Table 2 and are in line with [39].

3.2. Applying the IWeCASF-FEADM Integrated Approach

3.2.1. IWeCASF

IWeCASF processes the input satellite imagery to determine the landscape feature matrices F l . Then, the feature matrices M t , which are given as input, as well as the landscape feature matrices F l , are used to determine the particular features Φ r F l , M t , P r at a checkpoint P r . Next, two adjustment factors for the weather conditions are obtained from two fuzzy inference systems. The first is the steady adjustment α r , i determined by processing the influence of particular features Φ r F l , M t , P r over the weather conditions C i = 1 , 2 , , I r . The second is the variable adjustment B , which is the result of modeling the interactions among weather conditions C i = 1 , 2 , , I r . Later, the steady adjustment α r , i and the variable adjustment B are used to establish the weather conditions C i = 1 , 2 , , I r at checkpoint P r , which is called the final adjustment. The two stages (A and B) are described as follows:
A. Landscape features: As described in Section 2.2.1 A, B, the landscape features are extracted from the satellite imagery; then, the landscape features of the checkpoint under review are selected for further analysis.
Firstly, a color space conversion is performed. The image of the estimation region is taken from the RGB color space, as shown in Figure A1a, to the CIELAB color space. Afterwards, the converted image is decorrelated; in this way, the prototype clusters K l = 1 , 2 , , L are more distinguishable. The prototype clusters K l = 1 , 2 , , L are required to perform image segmentation using the objective function G , as defined in Equation (1). Therefore, a more discernible prototype cluster K l = 1 , 2 , , L encourages better image segmentation. Figure A1b presents the decorrelated image of the estimation zone. As described in [39], in Figure A1b, tree-covered areas are shown in blue, buildings are shown in pink, and grasslands are shown in green, while they appear yellow when the grass is dry. Moreover, in Figure A1b, the prototype clusters K l = 1 , 2 , 3 , which are used to define the landscape feature matrices F l = 1 , 2 , 3 , are shown. The color segmentation utilizes a merge algorithm to assign each pixel of the image to a prototype cluster k l .
After performing the color segmentation of the estimation zone image, the original image is separated into L = 5 landscape feature matrices F l = 1 , 2 , , L , whose elements f x , y l contain the number of pixels of sector s ( x , y ) that belong to landscape feature l . In Figure 4, the landscape feature matrix F l = 2 (tree-covered area) is depicted. Furthermore, a magnified image of sector s ( 22 , 26 ) is presented to illustrate the element f 22 , 26 2 = 48 , i.e., there are 48 of a total of 156 pixels that belong to landscape feature l = 2 at sector s ( 22 , 26 ) .
After determining the landscape feature matrices F l = 1 , 2 , , L , the next step is selecting the particular features at a checkpoint. As mentioned previously, every sector within the estimation region can be considered a checkpoint; thus, the checkpoint P r under consideration in the experiment described in this section is the sector s ( 22 , 26 ) . In this way, using Equations (2) and (3), the particular feature vector Φ r F t , M t , P r at checkpoint P r = s ( 22 , 26 ) is as presented in Equations (26) and (27), and they are depicted in Figure 5.
In Equation (26), the particular feature vector Φ r F l , M t , P r is expressed as a number of pixels; for example, there are 21 pixels that belong to landscape feature l = 1 (grassland), and 69 pixels belong to landscape feature l = 3 (buildings). Furthermore, the particular feature vector Φ r F l , M t , P r is expressed as a percentage in Equation (27); for example, 65.38% of the pixels at the checkpoint P r = s ( 22 , 26 ) correspond to l = 1 (grassland), and 35.1 % of the pixels belong to l = 3 (buildings).
Φ r F l , M t , P r = [ 21 , 105 , 69 , ( 196 / A ) , 574 , ( 196 / I I ) , ( 196 / 1 ) , ( 196 / 1 ) ]
Φ r F l , M t , P r = [ 65.38 , 53.57 , 35.1 , ( 100 / A ) , 35.2 , ( 100 / I I ) , ( 100 / 1 ) , ( 100 / 1 ) ]
B. Weather condition adjustment: As explained previously, there are two adjustment factors for the weather conditions. The first is derived from steady adjustment, which analyzes the landscape features of the estimation region and their influence over the weather conditions. The second is obtained from variable adjustment, which considers the interactions between the weather conditions.
The steady and variable adjustments are developed using a fuzzy inference system (FIS). For this reason, a fuzzification stage is performed. According to Equations (4)–(8), the weather conditions C i = 1 , 2 , , I 0 are measured at the primary checkpoint P 0 .
Five membership functions ( J = 5 ) are defined to fuzzify each weather condition, as in [39]; the parameters used to fuzzify the weather conditions C i = 1 , 2 , , I 0 and their actual values are also described in [39]. The resulting fuzzy weather condition vectors C ~ i 0 ( v i ) are shown in Equations (28)–(32).
C ~ 1 0 ( 19.5 ) = 0 , 0 , 0.9 , 0 , 0                   Temperature
C ~ 2 0 ( 5.8 ) = 0 , 0 , 0.0909 , 0.2667 , 0     Rain
C ~ 3 0 ( 265 ) = 0 , 0 , 0.25 , 0 , 0                 Solar   radiation
C ~ 4 0 13 = 0 , 0 , 0.6 , 0 , 0                     Wind   speed
C ~ 5 0 ( 3.28 ) = 0 , 0 , 0.98 , 0 , 0               Evapotranspiration
Moreover, a comparison Φ 0 , r between the particular features Φ r F l , M t , P r at checkpoint P r and the particular features Φ 0 F l , M t , P 0 at the primary checkpoint P 0 is performed using Equation (9); this is shown in Equations (33) and (34) and finally in Equation (35). The particular feature comparison Φ 0 , r is intended to assess the weather condition disparities regarding a variation in the landscape features.
Φ 0 F l , M t , P 0 = [ 100 , 0 , 0 , ( 100 / A ) , 22.76 , ( 100 / I I ) , ( 100 / 1 ) , ( 100 / 1 ) ]
Φ 0 , r = [ 100 10.2 , 0 53.2,0 35.1 , ( 100 100 ) / A , 22.76 35.6 , ( 100 100 ) / I , ( 100 100 ) / 1 , ( 100 100 ) / 1 ]
Φ 0 , r = [ 89.8 , 53.2 , 35.1 , ( 0 / A ) , 12.84 , ( 0 / I ) , ( 0 / 1 ) , ( 0 / 1 ) ]
Afterwards, the fuzzy particular feature matrix Φ ~ r , shown in Equation (36), is derived by fuzzifying the comparison Φ 0 , r presented in Equation (35) according to Equation (10). The resulting fuzzy particular feature matrix Φ ~ r presented in Equation (36) shows the great similarity between primary checkpoint P 0 and checkpoint P r ; as a result, a barely significant weather condition adjustment is expected at checkpoint P r .
Φ ~ r = 0 0 1 0.9388 0 0 0.2812 0.125 0 0 1 0 0 0.68 0 0 1 0 0 1 0 0 1 0
As defined in [39], the steady adjustment FIS utilizes IF THEN rules as the form of the rule n = 23 , which is described as follows:
  • IF the grassland is high ( ϕ ~ 1 , 2 r ), the tree cover area is equal ( ϕ ~ 2 , 2 r ), the building area is less ( ϕ ~ 3 , 1 r ) , the elevation is equal ( ϕ ~ 4 , 2 4 ) , and the spatial configuration is less ( ϕ ~ 5 , 2 r ), THEN the adjustment factor for the temperature ( α 4 , 1 ) is null ( λ 1 , 4 ), the adjustment factor for rain ( α 4 , 2 ) is null ( λ 2 , 3 ), the adjustment factor for solar radiation ( α 4 , 3 ) is barely high ( λ 3.4 ), the adjustment factor for the wind speed ( α 4 , 4 ) is null ( λ 4 , 3 ), and the adjustment factor for evapotranspiration ( α 4 , 5 ) is barely high ( λ 5 , 4 ).
After assessing the N = 245 rules as described in [39] and Section 3.2.1 B, the result is the steady adjustment factor α r , i , presented as the vector Λ r = α r , i 1 , I , shown in Equation (37). For example, the adjustment factor for solar radiation α r , 3 shows that, at checkpoint P r , the solar radiation adjustment factor is approximately + 4.99 % of its value.
Λ r = [ 0.0288 , 0.003 , 0.0499 , 0.225 , 0.288 ]
Moreover, the certainty of the weather condition replication η i = 1 , 2 , , I improves the modeling of inconsistent weather conditions, as described in Section 3.2.1 B. According to [39], the certainty η i = 1 , 2 , , I is calculated with a distribution function. Then, the variable adjustment is developed to determine the second adjustment factor. In the same way as for the steady adjustment FIS, a set of IF THEN rules is utilized. However, in the variable adjustment, Q = 445 rules are assessed; for example, the following rule corresponds to q = 412 .
  • IF the temperature is high c ~ 1 , 4 0  and rain is low ( c ~ 2 , 2 0 )  and solar radiation is higher c ~ 3 , 5 0  and evapotranspiration is higher, THEN outputs β 1 , β 5 , β 2   β 3 , β 4 are equal ( ω 1 , 3 , ω 2 , 3 , ω 3 , 3 , ω 4 , 3 , ω 5 , 3 ).
The variable adjustment factor B is unique for the whole estimation region. In Equation (38), the adjustment factor B resulting from the variable adjustment of the weather conditions C i = 1 , 2 , , I 0 is presented. For example, the variable adjustment factor for rain β 2 is 0.0167 , which means that there is an adjustment of 1.67 % in the rain C 2 0 .
B = [ 0 , 0.0167 , 0.0498 , 0.0498 , 0.0214 ]
Finally, IWeCASF determines the weather conditions at the checkpoint P r by using Equation (11), where the steady adjustments α r , i 1 , I and variable B , which are obtained in Equations (37) and (38), are utilized. The resulting adjusted weather conditions C i = 1 , 2 , , I r are shown in Table 3. The value of rain is C 2 0 = 5.8 , and its adjustment factors are α r , 1 = 0.003 and β 2 = 0.0167 , according to Equations (37) and (38). As a consequence, the measured value C 2 0 = 5.8 is decreased in a combined factor of 0.3 % 1.67 % = 1.97 % of the measured value. Thus, the adjusted value of rain is C 2 r = 5.68 .

3.2.2. Irrigation Water Addition

Afterwards, the weather conditions C i = 2 r are applied, as shown in Equation (12), with the irrigation water record δ x r , y r at checkpoint P r . In this experiment, the irrigation water record δ x r , y r at checkpoint P r is δ x r , y r = 4.12 . Therefore, according to Equation (12), the actual weather conditions C ˇ i = 1 , 3 , 4 , 5 r are the same as the weather conditions given in Table 3, as shown in Equations (39) and (41)–(43); meanwhile, the actual weather condition C ˇ i = 2 r is calculated using Equation (40). They are introduced to FEADM.
C ˇ i = 1 r = C 1 r = 18.93             Temperature
C ˇ i = 2 r = C 2 r + δ x r , y r = 5.8 + 4.12 = 9.92       Rain
C ˇ i = 3 r = C 3 r = 238.57           Solar   radiation
C ˇ i = 4 r = C 4 r = 9.42               Wind   speed
C ˇ i = 5 r = C 5 r = 2.62               Evapotranspiration

3.2.3. FEADM

The second stage of the integrated approach is the “fuzzy estimation approach based on decision making” (FEADM) [32]. Unlike IWeCASF [39], which is performed only one time for each regional estimation R e , FEADM is executed as many times as there are checkpoints P r . FEADM is implemented to obtain a soil moisture point estimate S m r at each checkpoint P r . Consequently, FEADM receives as inputs the actual weather conditions C ˇ i = 1 , 2 , , I r at checkpoint P r , obtained from IWeCASF in the previous stage of the integrated approach. Then, FEADM determines the soil moisture level that best fits the current weather conditions. Afterwards, a quantitative value of the chosen soil moisture level is obtained. This quantitative value is the soil moisture point estimate S m r at checkpoint P r . The two stages (A and B) are described as follows:
A. FEADM decision-making stage: This begins with the fuzzy analytical hierarchy process (FAHP), as described in [32] and Section 3.2.3 A. In this stage, the criteria C i are compared to obtain the judgment matrix P ~ defined in Equation (44). Each element p ~ x , y of judgement matrix P ~ denotes the relevance of criteria C x regarding criteria C y . The assignment of the relative importance values is performed as in [32].
P ~ =   C 1 C 2 C 3 C 4 C 5 C 1 C 2 C 3 C 4 C 5 1 , 1 , 1 1 , 1 , 1 1 , 2 , 3 1 , 2 , 3 1 , 2 , 3 1 , 1 , 1 1 , 1 , 1 2 , 4 , 6 4 , 6 , 8 4 , 6 , 8 1 3 , 1 2 , 1 1 3 , 1 2 , 1 1 , 1 , 1 4 , 6 , 8 1 , 2 , 3 1 3 , 1 2 , 1 1 8 , 1 6 , 1 4 1 8 , 1 6 , 1 4 1 , 1 , 1 1 , 1 , 1 1 3 , 1 2 , 1 1 8 , 1 6 , 1 4 1 3 , 1 2 , 1 1 , 1 , 1 1 , 1 , 1
Afterwards, the extent of synthetic values S ~ i = 1 , 2 , , I is calculated with Equation (13), as in [32]. Furthermore, the total integral value with the index of optimism I T S ~ i = 1 , 2 , , I and the extent of synthetic values S ~ i = 1 , 2 , , I are used in Equation (14) to obtain the weight set W , which is shown in Equation (45). According to Equation (45), the most relevant criterion is rain ( C 2 ), as was expected.
W = 0.1834 , 0.3805 , 0.2555 , 0.0791 , 0.1015
The FEADM decision-making stage continues with the expert evaluation. Firstly, the assessment criteria C i = 1 , 2 , , I are fuzzified, as performed in [32]. As shown in Table 3, the actual weather conditions C ˇ i = 1 , 2 , , I r are utilized as assessment criteria C i = 1 , 2 , , I . Thus, in this experiment, the values C ˇ i = 1 , 2 , , I r shown in Equations (39)–(43) are the values z i , which, in addition to the membership functions μ C i , j z i defined in [32], are used to obtain the fuzzy criterion C ~ i = μ C i , 1 ( z i ) , , μ C i , 5 ( z i )   . After fuzzifying the five criteria ( I = 5 ), the fuzzy criteria matrix C ~ can be obtained using Equation (15). In this experiment, the resulting fuzzy criteria matrix C ~ is as shown in Equation (46).
C ~ = 0 0.0333 0.36 0 0 0 0 0 0 1 0 0.6667 0 0 0 0 0 0.8 0 0 0 0.3 0 0 0
Next, the weight set W from Equation (45) and the fuzzy criteria matrix C ~ from Equation (46) are used in Equation (47) to calculate the elements of the fuzzy weighed criteria matrix C ~ w , shown in Equation (18). The C ~ w resulting from this experiment is given in Equation (47).
C ~ w = 0 0.0061 0.066 0 0 0 0 0 0 0.385 0 0.1703 0 0 0 0 0 0.0633 0 0 0 0.0305 0 0 0
The fuzzy weighed criteria matrix C ~ w is required in the selective assessment [32]. The number of assessments N k = 1 , 2 , , K used in this experiment is the same as defined in [32]. To accomplish each set of evaluations E k = e 1 k , e 2 k , , e N k k , according to Equations (19) and (20), the selective assessment function f k n , C ~ w evaluates the fuzzy weighted criteria c i , j w defined in the relationship modeling matrix D n k , which only contains those c i , j w relevant for evaluation e n k . The matrices D n = 1 , 2 , , 25 5 shown in Equations (48)–(50) correspond to alternative very high ( A 5 ) instances; meanwhile, the matrices depicted in Equations (51)–(53) are three of the matrices used to assess alternative very low ( A 1 ) instances. All matrices D n k presented in this experiment correspond to evaluations e n k performed at checkpoint P r = s ( 22 , 26 ) , as mentioned before. The matrices D n k are determined by the particular features Φ r F l , M t , P r of the checkpoint P r , where the evaluation e n k is conducted. Thus, given a checkpoint P r = x and a checkpoint P r = y with particular features Φ r = x F l , M t , P r = x = Φ r = y F l , M t , P r = y , the matrices D n k used to assess checkpoint P r = x can also be used to assess checkpoint P r = y ; i.e., if the spatial features of one or more checkpoints are the same, the assessment can be conducted using the same matrices D n k .
D 2 5 = 1 0.0061 1 1 1 1 1 1 1 0.385 1 1 1 1 1 1 1 1 1 1 1 0.0305 1 1 1
D 5 5 = 1 1 0.066 1 1 1 1 1 1 0.385 1 1 1 1 1 1 1 1 1 1 1 0.0305 1 1 1
D 17 5 = 1 1 0.066 1 1 1 1 1 0 1 0 1 1 1 1 1 1 1 1 1 1 0.0305 1 1 1
D 3 1 = 1 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 1 0.0633 1 1 1 1 1 1 0
D 32 1 = 1 1 1 0 1 0 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 0 1
D 37 1 = 1 1 1 1 1 0 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1
After assessing the N k matrices D n = 1 , 2 , , N k k of each alternative A k , the set of evaluations E k = e 1 k ,   e 2 k , , e N k k is performed. The resulting sets E k are of the form described in Equations (54) and (55), which correspond to sets E 1 and E 5 of very low ( A 1 ) and very high ( A 5 ) alternatives.
E 1 = e n 1 = 0 | n = 1 , 2 , , 37
E 5 = . 0385 , e n 5 = 0 | n = 2 , 3 , , 25
Afterwards, the best assessment a k is obtained for each evaluation set E k = e 1 k ,   e 2 k , , e N k k using Equation (22). The best assessment a k of each alternative in this experiment is shown in Table 4.
A 1 A 2 A 3 A 4 A 5 The FEADM decision-making stage ends with the selection of the best-assessed alternative A ^ , according to Equation (22) and Table 4. In this experiment, the best-assessed alternative A ^ is high ( A 4 ), as depicted in Equation (56).
A ^ = max 0 , 0 , 0 , 0.41 , 0 T
Thus far, FEADM has determined that the soil moisture level high ( A 4 ) is the best alternative according to the actual weather conditions ( C ˇ i = 1 , 2 , , I r ) at checkpoint P r . Next, the crisp value of soil moisture S m r is calculated in the FEADM fuzzy inference stage.
B. Fuzzy inference: In the fuzzy inference proposed in [32], a membership function μ A k = 1 , 2 , , K z s m assigned to every alternative A k = 1 , 2 , , K is modified using the best assessment a k of each alternative A k , as shown in Equation (23); i.e., the actual output membership function μ A k z s m of alternative A k is delimited by the best assessment a k of the same alternative. For the soil moisture estimate S m r , the range Z s m = [ 0 ,   10 ] is defined considering the measurement range of the soil moisture sensor utilized. Finally, the soil moisture estimate S m e is obtained using Equation (24) for fuzzy aggregation and Equation (25) for defuzzification. The soil moisture estimate S m r at checkpoint P r is shown in Equation (57).
S m r = 7.14

3.3. Point Estimates of Soil Moisture

As depicted in Figure 2, when regional soil moisture estimation R e = S m 1 , S m 2 , ,   S m R is conducted, IWeCASF, the first stage of the integrated approach, is completed only once, whereas FEADM, the second stage of the integrated approach, is performed as many times as required for soil moisture point estimate S m r . Furthermore, the previous analysis developed in Section 3.2.1, Section 3.2.2 and Section 3.2.3 has to be performed for every checkpoint P r . As a result, a soil moisture point estimate is obtained at each sector s ( x , y ) .
Figure 6 illustrates the regional estimation R e = S m 1 , S m 2 , ,   S m R of the area of interest. The irrigation water δ x r , y r increases the soil water content. This is why irrigation water causes higher soil moisture levels in specific sectors ( 9 S m < 10 ) . Moreover, there are areas in which the soil moisture estimation S m is unavailable ( N / A ) due to the existence of buildings.

4. Results and Discussion

The integrated approach based on linking FEADM with ICASF allowed regional soil moisture estimation, reducing the implementation and maintenance complexity of conventional automatic irrigation systems and their inherent costs. Although the experimental case study was for the region in Figure 1, this integrated approach can be applied in any estimation region with the required input data (satellite imagery, crop and soil data, weather conditions measured at a single checkpoint, and irrigation water record).
The method to obtain the regional soil moisture estimate combines two models published by our research team. In addition, it considers the irrigation water records, which were not previously used. The integrated approach is functional and has been tested for several years. The IWeCASF model [39] provides the required inputs for FEADM [32]. In the previous section, IWeCASF’s results were complemented by the irrigation water records of the region of interest. However, if the irrigation water record is not considered and FEADM is applied, the results would be as depicted in Figure 7.
In this case, the soil moisture S m r is higher in mostly tree-covered checkpoints, such as the checkpoint at the sector P r = s ( 13 ,   22 ) with S m r = 6.32 , compared to mainly grassland checkpoints, such as P r = s ( 28 ,   28 ) , where the soil moisture content is S m r = 5.09. This is because tree-covered areas preserve more soil moisture when performing this test under the actual weather conditions. This fact can be relevant if the crop type is a fruit that is ripened on trees. Furthermore, it is worth noting that the integrated approach can be performed without an irrigation water record. Therefore, if required, this integrated approach is suitable for the determination of soil moisture using only rain in sowing regions. Thus, the integrated approach’s application is comprehensive.
The regional estimation R e = S m 1 , S m 2 , ,   S m R is presented as a set of soil moisture point estimates S m r . In Figure 6 and Figure 7, each sector s ( x , y ) is a checkpoint P r so that FEADM is executed at every sector s ( x , y ) of the region of interest. In this way, the regional estimation R e = S m 1 , S m 2 , ,   S m R is composed of R = 1366   point estimates S m r = 1 , 2 , , 1366 .
A total of 921 measured weather conditions C i = 1 , 2 , , 5 0 and the same number of irrigation water records Δ r were utilized to develop the IWeCASF-FEADM integrated approach. A total of 645 records were used to improve the performance of the approach, whereas 276 records were assigned to validate it. These 921 records were collected over almost three years to test the soil moisture estimation approach, considering a wide range of irrigation and weather conditions.
The estimation errors of four checkpoints from the estimation region are presented in Figure 8. The graphs correspond to checkpoints P r = ( 17 ,   8 ) , P r = 22 ,   26 , P r = 14 ,   35 , and P r = 33 ,   47 . The estimations at checkpoint P r = 17 ,   8 are the most accurate; the normalized error E R M S for this checkpoint is E R M S = 0.0361 . Meanwhile, the normalized error for checkpoint P r = 22 ,   26 is E R M S = 0.0479 ; for checkpoint P r = 14 ,   35 , the error is E R M S = 0.0519 , and, for checkpoint P r = 33 ,   47 , the error is E R M S = 0.0544 .
According to these results, and recalling the location of the primary checkpoint P 0 = 16 ,   6 , the performance of the integrated approach depends on the distance between the checkpoint P r and the primary checkpoint P 0 ; i.e., the point estimate of soil moisture is more accurate when obtained closer to the checkpoint where the weather conditions C i = 1 , 2 , , I 0 are measured.
Nevertheless, the results obtained at checkpoints that are distant from the primary checkpoint present a normalized error that does not imply a limitation in obtaining the soil moisture to determine the irrigation water supply.

Soil Moisture Behavior under Landscape Variability

The soil moisture is highly dependent on the landscape features and the weather conditions. An additional advantage of the method developed in this work (the integrated approach) is that it allows simulations to observe and verify these dependencies. In this sense, we simulate and present in Figure 9 the results when the landscape features F t vary within a determined section of the estimation region. The influence of a change in a landscape feature F t over the weather conditions C i = 1 , 2 , , I r and soil moisture S m r is outstanding. A determined area (I, II, and III) of the region of interest is depicted. Each colored rectangle is a sector s ( x , y ) .
In Figure 9a, changes in the landscape features are simulated: image I presents an area with the original landscape; image II shows the modified original landscape so that the grassland prevails; and, in image III, after the modifications, the tree-covered areas prevail. Moreover, Figure 9b–f present the temperature ( C 1 r ) , rain ( C 2 r ) , solar radiation ( C 3 r ) , wind speed ( C 4 r ) , and evapotranspiration ( C 5 r ) for each category: the original landscape, grassland, and tree-covered areas. For these categories, the weather conditions are compared with the weather conditions at the primary checkpoint P 0 . Finally, Figure 9g presents the soil moisture regional estimation ( R e ) for the simulated areas.
The simulation highlights how changes in the landscape features impact the weather conditions and soil moisture. Section 2.2.1 D and Section 3.2.1 B presented the landscape adjustment and weather conditions adjustment, respectively, explaining that there are two adjustment factors for the weather conditions. The first is derived from steady adjustment, which analyzes the landscape features of the estimation region and its influence over the weather conditions. The second is obtained from variable adjustment, which considers the interactions among the weather conditions. Using a fuzzy inference system, the adjustment factor is the result of the comparison of the characteristics of the sector under analysis and the primary checkpoint P 0 (where the weather conditions are measured). Therefore, the linguistic variables refer to the values of the weather conditions ( C i r ) analyzed in the sector s ( x , y ) and whether they are lower, equal to or higher than the values of the weather conditions at the primary checkpoint P 0 .
For example, in Figure 9b, for the temperature ( C 1 r ) , column II (mainly grassland) has more sectors with equal adjustments than the original landscape (column I). On the other hand, column III (tree-covered areas) has more sectors with lower adjustments than the original landscape (column I).
Similarly, in Figure 9d, for solar radiation ( C 3 r ) , column II (mainly grassland) has more sectors with equal adjustments than the original landscape (column I); likewise, there are sectors with higher adjustments. Column III (tree-covered areas) has more sectors with lower adjustments than the original landscape (column I). Reasonable results can be observed in Figure 9c,e,f. For example, in Figure 9f, column II has more sectors with evapotranspiration ( C 5 r ) with equal and higher adjustments than the original landscape (column I), since this area is mainly grassland. We note again that all lower, equal or higher adjustments are with respect to the primary checkpoint P 0 (where the weather conditions are measured).
Finally, Figure 9g compares the soil moisture regional estimation ( R e ) concerning the primary checkpoint P 0 . Column III (tree-covered areas) has more sectors with higher humidity.
The results of soil moisture were compared indirectly with interpolation methods; however, these interpolation methods show limitations when used in a regional estimation as proposed in this study because spatial interpolation models use a network of measurement points; rainfall is not remarkably homogenous in a region, thus basing its interpolation value only on spatial measurements or statistical records can be misleading.
This new regional soil moisture estimation method uses an intelligent adjustment of weather conditions based on spatial features. The regional estimation of soil moisture depends on the accuracy of the intelligent adjustment of weather conditions based on spatial features [39] and the fuzzy point estimation based on decision-making, validated in [32]. In [39], Tables 10–14 compared weather conditions adjusted, measured, and interpolated for temperature, rain, solar radiation, wind speed, and evapotranspiration, respectively. All cases highlight the advantage of using an intelligent adjustment of weather conditions based on spatial features; therefore, the regional estimation of soil moisture will depend on the accuracy of those results and will be best.

5. Conclusions

A new method, named the integrated approach, to obtain a regional soil moisture estimation is developed in this study. This method has been tested over several years by comparing the results of soil moisture regional estimation with measurements at many points in the study region. The integrated approach achieves a regional soil moisture estimation as a set of soil moisture point estimates. The method can detect when the soil moisture is deficient in a region, allowing actions to prevent water stress. A regional estimate can reduce an irrigation system’s operational and maintenance complexity.
This new method exploits the advantages of FEADM, such as obtaining a quantitative estimate of the soil moisture and avoiding the complexity involved in measuring it. Likewise, by combining IWeCASF and FEADM, this integrated approach addresses the input acquisition for FEADM; consequently, FEADM can obtain soil moisture point estimates at several checkpoints scattered within an estimation region.
FEADM and IWeCASF are two computer methods already presented in the literature, but they had not been integrated as part of a single method, as proposed in this work.
Furthermore, the integrated approach allows the more accessible display of the relationships among the weather conditions, landscape features and soil moisture content. Some weather conditions, such as rain C i = 2 r , exert a greater influence over the soil moisture content; for example, when there is rainfall, the soil moisture content is high. In addition, the landscape features, such as tree-covered areas ( F t = 2 ) , can also modify the soil moisture content because solar radiation does not fully reach the soil.
This new method based on linking FEADM with IWeCASF allowed regional soil moisture estimation, reducing the implementation and maintenance complexity of conventional automatic irrigation systems and their inherent costs. Although the experimental case study was for a specific region, this integrated approach can be applied in any estimation region with the required input data (satellite imagery, crop and soil data, weather conditions measured at a single checkpoint, and irrigation water record).
IWeCASF justifies the advantage of using an intelligent adjustment of weather conditions based on spatial features. Therefore, linking FEADM-IWeCASF for the regional soil moisture estimation will be better than interpolation methods. All of the above verify compliance with the hypotheses proposed for this work.
However, as a limitation, we point out that the experimental case study was for the region of Figure 1. Other experimental case studies in different regions are already being performed, and the analysis of their results is considered future work.

Author Contributions

Conceptualization, methodology, investigation, supervision, validation, funding acquisition, project administration, and review and editing, L.P.S.-F. Conceptualization, investigation, methodology, data acquisition, formal analysis, validation, software, visualization, and writing—original draft preparation, D.A.F.-C. Conceptualization, methodology, investigation, supervision, validation, and review and editing, L.A.S.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Instituto Politécnico Nacional of Mexico and the National Council of Humanities, Science and Technology (CONAHCYT) of Mexico. SIP 20241747 and SIP 20231285.

Data Availability Statement

The software used is MATLAB 2022, and it is available under license. The same applies to QML (Qt Modeling Language) and Limbo. Online data and codes are available: https://github.com/Lsanchez2024/An-integrated-approach-for-regional-estimations-of-soil-moisture- (accessed on 18 September 2024).

Acknowledgments

We thank the Instituto Politécnico Nacional, Mexico and the National Council of Humanities, Science and Technology (CONAHCYT) of Mexico for their financial support.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Color Space Image Conversion

A color space conversion, as shown in Figure A1, highlights the landscape features of an estimation region. The colors that represent each landscape feature F l are more distinguishable. As a result, better color segmentation is performed.
Figure A1. Estimation region. (a) The original image is in RGB color space. (b) Image in CIELAB color space under decorrelation used to obtain landscape matrices F l = 1 , 2 , 3 .
Figure A1. Estimation region. (a) The original image is in RGB color space. (b) Image in CIELAB color space under decorrelation used to obtain landscape matrices F l = 1 , 2 , 3 .
Hydrology 11 00170 g0a1

References

  1. Goumopoulos, C.; O’Flynn, B.; Kameas, A. Automated zone-specific irrigation with wireless sensor/actuator network and adaptable decision support. Comput. Electron. Agric. 2014, 105, 20–33. [Google Scholar] [CrossRef]
  2. Adeyemi, O.; Grove, I.; Peets, S.; Norton, T. Advanced Monitoring and Management Systems for Improving Sustainability in Precision Irrigation. Sustainability 2017, 9, 353. [Google Scholar] [CrossRef]
  3. Lloret, J.; Sendra, S.; Garcia, L.; Jimenez, J.M. A Wireless Sensor Network Deployment for Soil Moisture Monitoring in Precision Agriculture. Sensors 2021, 21, 7243. [Google Scholar] [CrossRef]
  4. El Mezouari, A.; El Fazziki, A.; Sadgal, M. Smart Irrigation System. IFAC-PapersOnLine 2022, 55, 3298–3303. [Google Scholar] [CrossRef]
  5. Phasinam, K.; Kassanuk, T.; Shinde, P.P.; Thakar, C.M.; Sharma, D.K.; Mohiddin, M.K.; Rahmani, A.W. Application of IoT and Cloud Computing in Automation of Agriculture Irrigation. J. Food Qual. 2022, 2022, 8285969. [Google Scholar] [CrossRef]
  6. Sánchez-Sutil, F.; Cano-Ortega, A. Smart Control and Energy Efficiency in Irrigation Systems Using LoRaWAN. Sensors 2021, 21, 7041. [Google Scholar] [CrossRef] [PubMed]
  7. Abioye, E.A.; Abidin, M.S.Z.; Mahmud, M.S.A.; Buyamin, S.; Ishak, M.H.I.; Rahman, M.K.I.A.; Otuoze, A.O.; Onotu, P.; Ramli, M.S.A. A review on monitoring and advanced control strategies for precision irrigation. Comput. Electron. Agric. 2020, 173, 105441. [Google Scholar] [CrossRef]
  8. Raouhi, E.M.; Zouizza, M.; Lachgar, M.; Zouani, Y.; Hrimech, H.; Kartit, A. AIDSII: An AI-based digital system for intelligent irrigation. Softw. Impacts 2023, 17, 100574. [Google Scholar] [CrossRef]
  9. Moradkhani, H. Hydrologic remote sensing and land surface data assimilation. Sensors 2008, 8, 2986–3004. [Google Scholar] [CrossRef]
  10. Crow, W.T.; Chen, F.; Reichle, R.H.; Liu, Q. L band microwave remote sensing and land data assimilation improve the representation of prestorm soil moisture conditions for hydrologic forecasting. Geophys. Res. Lett. 2017, 44, 5495–5503. [Google Scholar] [CrossRef]
  11. Phillips, A.J.; Newlands, N.K.; Liang, S.H.L.; Ellert, B.H. Integrated sensing of soil moisture at the field-scale: Measuring, modeling and sharing for improved agricultural decision support. Comput. Electron. Agric. 2014, 107, 73–88. [Google Scholar] [CrossRef]
  12. Romero, R.; Muriel, J.L.; García, I.; Muñoz de la Peña, D. Research on automatic irrigation control: State of the art and recent results. Agric. Water Manag. 2012, 114, 59–66. [Google Scholar] [CrossRef]
  13. Millán, S.; Casadesús, J.; Campillo, C.; Moñino, M.J.; Prieto, M.H. Using Soil Moisture Sensors for Automated Irrigation Scheduling in a Plum Crop. Water 2019, 11, 2061. [Google Scholar] [CrossRef]
  14. Aubert, D.; Loumagne, C.; Oudin, L. Sequential assimilation of soil moisture and streamflow data in a conceptual rainfall–runoff model. J. Hydrol. 2003, 280, 145–161. [Google Scholar] [CrossRef]
  15. Perrin, C.; Michel, C.; Andréassian, V. Does a large number of parameters enhance model performance? Comparative assessment of common catchment model structures on 429 catchments. J. Hydrol. 2001, 242, 275–301. [Google Scholar] [CrossRef]
  16. Agyeman, B.T.; Bo, S.; Sahoo, S.R.; Yin, X.; Liu, J.; Shah, S.L. Soil moisture map construction by sequential data assimilation using an extended Kalman filter. J. Hydrol. 2021, 598, 126425. [Google Scholar] [CrossRef]
  17. Jones, H.G. Irrigation scheduling: Advantages and pitfalls of plant-based methods. J. Exp. Bot. 2004, 55, 2427–2436. [Google Scholar] [CrossRef]
  18. Pan, F.; Nieswiadomy, M.; Qian, S. Application of a soil moisture diagnostic equation for estimating root-zone soil moisture in arid and semi-arid regions. J. Hydrol. 2015, 524, 296–310. [Google Scholar] [CrossRef]
  19. Gu, Z.; Qi, Z.; Burghate, R.; Yuan, S.; Jiao, X.; Xu, J. Irrigation Scheduling Approaches and Applications: A Review. J. Irrig. Drain. Eng. 2020, 146, 04020007. [Google Scholar] [CrossRef]
  20. Kumar, S.V.; Dirmeyer, P.A.; Peters-Lidard, C.D.; Bindlish, R.; Bolten, J. Information theoretic evaluation of satellite soil moisture retrievals. Remote Sens. Environ. 2018, 204, 392–400. [Google Scholar] [CrossRef]
  21. Beck, H.E.; Pan, M.; Miralles, D.G.; Reichle, R.H.; Dorigo, W.A.; Hahn, S.; Sheffield, J.; Karthikeyan, L.; Balsamo, G.; Parinussa, R.M.; et al. Evaluation of 18 satellite- And model-based soil moisture products using in situ measurements from 826 sensors. Hydrol. Earth Syst. Sci. 2021, 25, 17–40. [Google Scholar] [CrossRef]
  22. Jarray, N.; Ben Abbes, A.; Rhif, M.; Dhaou, H.; Ouessar, M.; Farah, I.R. SMETool: A web-based tool for soil moisture estimation based on Eo-Learn framework and Machine Learning methods. Environ. Model. Softw. 2022, 157, 105505. [Google Scholar] [CrossRef]
  23. Pulvirenti, L.; Squicciarino, G.; Cenci, L.; Boni, G.; Pierdicca, N.; Chini, M.; Versace, C.; Campanella, P. A surface soil moisture mapping service at national (Italian) scale based on Sentinel-1 data. Environ. Model. Softw. 2018, 102, 13–28. [Google Scholar] [CrossRef]
  24. Tischler, M.; Garcia, M.; Peters-Lidard, C.; Moran, M.S.; Miller, S.; Thoma, D.; Kumar, S.; Geiger, J. A GIS framework for surface-layer soil moisture estimation combining satellite radar measurements and land surface modeling with soil physical property estimation. Environ. Model. Softw. 2007, 22, 891–898. [Google Scholar] [CrossRef]
  25. Souissi, R.; Al Bitar, A.; Zribi, M. Accuracy and Transferability of Artificial Neural Networks in Predicting in Situ Root-Zone Soil Moisture for Various Regions across the Globe. Water 2020, 12, 3109. [Google Scholar] [CrossRef]
  26. Ghasemloo, N.; Matkan, A.A.; Alimohammadi, A.; Aghighi, H.; Mirbagheri, B. Estimating the Agricultural Farm Soil Moisture Using Spectral Indices of Landsat 8, and Sentinel-1, and Artificial Neural Networks. J. Geovisualization Spat. Anal. 2022, 6, 1–12. [Google Scholar] [CrossRef]
  27. Han, H.; Choi, C.; Kim, J.; Morrison, R.R.; Jung, J.; Kim, H.S.; Han, H.; Choi, C.; Kim, J.; Morrison, R.R.; et al. Multiple-Depth Soil Moisture Estimates Using Artificial Neural Network and Long Short-Term Memory Models. Water 2021, 13, 2584. [Google Scholar] [CrossRef]
  28. Yu, Z.; Liu, D.; Lü, H.; Fu, X.; Xiang, L.; Zhu, Y. A multi-layer soil moisture data assimilation using support vector machines and ensemble particle filter. J. Hydrol. 2012, 475, 53–64. [Google Scholar] [CrossRef]
  29. He, B.; Jia, B.; Zhao, Y.; Wang, X.; Wei, M.; Dietzel, R. Estimate soil moisture of maize by combining support vector machine and chaotic whale optimization algorithm. Agric. Water Manag. 2022, 267, 107618. [Google Scholar] [CrossRef]
  30. Liu, D.; Mishra, A.K.; Yu, Z. Evaluating uncertainties in multi-layer soil moisture estimation with support vector machines and ensemble Kalman filtering. J. Hydrol. 2016, 538, 243–255. [Google Scholar] [CrossRef]
  31. Zhu, Q.; Wang, Y.; Luo, Y. Improvement of multi-layer soil moisture prediction using support vector machines and ensemble Kalman filter coupled with remote sensing soil moisture datasets over an agriculture dominant basin in China. Hydrol. Process. 2021, 35, e14154. [Google Scholar] [CrossRef]
  32. Flores-Carrillo, D.A.; Sánchez-Fernández, L.P.; Sánchez-Pérez, L.A.; Carbajal-Hernández, J.J. Soil moisture Fuzzy Estimation Approach based on Decision-Making. Environ. Model. Softw. 2017, 91, 223–240. [Google Scholar] [CrossRef]
  33. Wang, S.; Li, R.; Wu, Y.; Wang, W. Estimation of surface soil moisture by combining a structural equation model and an artificial neural network (SEM-ANN). Sci. Total Environ. 2023, 876, 162558. [Google Scholar] [CrossRef] [PubMed]
  34. Kambalimath, S.; Deka, P.C. A basic review of fuzzy logic applications in hydrology and water resources. Appl. Water Sci. 2020, 10, 1–14. [Google Scholar] [CrossRef]
  35. Rodríguez-Pérez, A.M.; Rodríguez-Gonzalez, C.A.; López, R.; Hernández-Torres, J.A.; Caparrós-Mancera, J.J. Water Microturbines for Sustainable Applications: Optimization Analysis and Experimental Validation. Water Resour. Manag. 2024, 38, 1011–1025. [Google Scholar] [CrossRef]
  36. Rodríguez-Pérez, Á.M.; García-Chica, A.; Caparros-Mancera, J.J.; Rodríguez, C.A. Turbine-Based Generation in Greenhouse Irrigation Systems. Hydrology 2024, 11, 149. [Google Scholar] [CrossRef]
  37. Elshorbagy, A.; Parasuraman, K. On the relevance of using artificial neural networks for estimating soil moisture content. J. Hydrol. 2008, 362, 1–18. [Google Scholar] [CrossRef]
  38. Ondieki, J.; Laneve, G.; Marsella, M.; Mito, C.; Ondieki, J.; Laneve, G.; Marsella, M.; Mito, C. Enhancing Surface Soil Moisture Estimation through Integration of Artificial Neural Networks Machine Learning and Fusion of Meteorological, Sentinel-1A and Sentinel-2A Satellite Data. Adv. Remote Sens. 2023, 12, 99–122. [Google Scholar] [CrossRef]
  39. Sánchez-Fernández, L.P.; Flores-Carrillo, D.A.; Sánchez-Pérez, L.A. Computer Model for an Intelligent Adjustment of Weather Conditions Based on Spatial Features for Soil Moisture Estimation. Mathematics 2024, 12, 152. [Google Scholar] [CrossRef]
  40. Welcome to the QGIS Project! Available online: https://qgis.org/en/site/ (accessed on 22 December 2023).
  41. Sridhar, V.; Hubbard, K.G.; You, J.; Hunt, E.D. Development of the Soil Moisture Index to Quantify Agricultural Drought and Its “User Friendliness” in Severity-Area-Duration Assessment. J. Hydrometeorol. 2008, 9, 660–676. [Google Scholar] [CrossRef]
  42. Saha, A.; Patil, M.; Goyal, V.C.; Rathore, D.S. Assessment and Impact of Soil Moisture Index in Agricultural Drought Estimation Using Remote Sensing and GIS Techniques. Proceedings 2018, 7, 2. [Google Scholar] [CrossRef]
  43. WCA 2020|World Programme for the Census of Agriculture|Food and Agriculture Organization of the United Nations. Available online: https://www.fao.org/world-census-agriculture/wcarounds/wca2020/en/ (accessed on 28 November 2023).
  44. CIE Publications—Premium Source for Knowledge on Light and Lighting|CIE. Available online: https://cie.co.at/publications (accessed on 11 December 2023).
  45. Chang, D.-Y. Applications of the extent analysis method on fuzzy AHP. Eur. J. Oper. Res. 1996, 95, 649–655. [Google Scholar] [CrossRef]
  46. Şen, C.G.; Çınar, G. Evaluation and pre-allocation of operators with multiple skills: A combined fuzzy AHP and max–min approach. Expert Syst. Appl. 2010, 37, 2043–2053. [Google Scholar] [CrossRef]
  47. Wang, Y.M.; Luo, Y.; Hua, Z. On the extent analysis method for fuzzy AHP and its applications. Eur. J. Oper. Res. 2008, 186, 735–747. [Google Scholar] [CrossRef]
  48. Liou, T.-S.; Wang, M.-J.J. Ranking fuzzy numbers with integral value. Fuzzy Sets Syst. 1992, 50, 247–255. [Google Scholar] [CrossRef]
Figure 1. Example of application of integrated soil moisture estimation approach for 12 checkpoints. In the image, * FEADM [32]; ** IWeCASF [39]; and *** addition of irrigation water record δ x r , y r of checkpoint P r .
Figure 1. Example of application of integrated soil moisture estimation approach for 12 checkpoints. In the image, * FEADM [32]; ** IWeCASF [39]; and *** addition of irrigation water record δ x r , y r of checkpoint P r .
Hydrology 11 00170 g001
Figure 2. Overview of integrated approach for regional estimation of soil moisture. In the image, ** IWeCASF [39] and * FEADM [32].
Figure 2. Overview of integrated approach for regional estimation of soil moisture. In the image, ** IWeCASF [39] and * FEADM [32].
Hydrology 11 00170 g002
Figure 3. Feature matrix M t = 2 (crop type).
Figure 3. Feature matrix M t = 2 (crop type).
Hydrology 11 00170 g003
Figure 4. Landscape feature matrix F 2 (tree-covered area).
Figure 4. Landscape feature matrix F 2 (tree-covered area).
Hydrology 11 00170 g004
Figure 5. Landscape feature matrices F l = 1 , 2 , , 5 .
Figure 5. Landscape feature matrices F l = 1 , 2 , , 5 .
Hydrology 11 00170 g005
Figure 6. Soil moisture regional estimation = S m 1 , S m 2 , , S m R , considering irrigation water record δ x r , y r , which increases soil water content.
Figure 6. Soil moisture regional estimation = S m 1 , S m 2 , , S m R , considering irrigation water record δ x r , y r , which increases soil water content.
Hydrology 11 00170 g006
Figure 7. Soil moisture regional estimation R e = S m 1 , S m 2 , ,   S m R without considering irrigation water record δ x r , y r .
Figure 7. Soil moisture regional estimation R e = S m 1 , S m 2 , ,   S m R without considering irrigation water record δ x r , y r .
Hydrology 11 00170 g007
Figure 8. Soil moisture estimation error ( E R M S ) for four checkpoints in the estimation region.
Figure 8. Soil moisture estimation error ( E R M S ) for four checkpoints in the estimation region.
Hydrology 11 00170 g008
Figure 9. Landscape features F t , weather conditions C i = 1 , 2 , , I r , and soil moisture S m r variability. (a) Changes in landscape features. (bf) Weather conditions of temperature ( C 1 r ) , rain ( C 2 r ) , solar radiation C 3 r ) , wind speed C 4 r ) , and evapotranspiration ( C 5 r ) for each category: original landscape, grassland, and tree-covered areas. For these categories, the weather conditions are compared with those at the primary checkpoint. (g) Soil moisture regional estimation (Re) for simulated areas.
Figure 9. Landscape features F t , weather conditions C i = 1 , 2 , , I r , and soil moisture S m r variability. (a) Changes in landscape features. (bf) Weather conditions of temperature ( C 1 r ) , rain ( C 2 r ) , solar radiation C 3 r ) , wind speed C 4 r ) , and evapotranspiration ( C 5 r ) for each category: original landscape, grassland, and tree-covered areas. For these categories, the weather conditions are compared with those at the primary checkpoint. (g) Soil moisture regional estimation (Re) for simulated areas.
Hydrology 11 00170 g009aHydrology 11 00170 g009bHydrology 11 00170 g009c
Table 1. Weather conditions C i = 1 , 2 , , I r , I = 5 measured at primary checkpoint r = 0 .
Table 1. Weather conditions C i = 1 , 2 , , I r , I = 5 measured at primary checkpoint r = 0 .
C i = 1 , 2 , , I 0 Variable
C 1 0 Temperature ( T )
C 2 0 Rain ( R )
C 3 0 Solar radiation ( S r )
C 4 0 Wind speed ( W s )
C 5 0 Evapotranspiration ( E t )
Table 2. Weather conditions (measured) at the primary checkpoint P 0 .
Table 2. Weather conditions (measured) at the primary checkpoint P 0 .
C i = 1 , 2 , , I 0 Weather ConditionValue ( v i )
C 1 0 Temperature ( ° C ) 19
C 2 0 Rain ( m m ) 5.8
C 3 0 Solar radiation ( W m 2 ) 265
C 4 0 Wind speed ( k m h 1 ) 13
C 5 0 Evapotranspiration ( m m ) 3.28
Table 3. Weather condition comparison at checkpoint P 4 .
Table 3. Weather condition comparison at checkpoint P 4 .
Weather ConditionAdjusted
C a i 4
Measured
C m i 4
i = 1 Temperature   ( ° C ) 18.9319.5
i = 2 Rain (mm)5.685.8
i = 3 Solar radiation (W m−2) 238.57265
i = 4 Wind speed (Km h−1)9.4213
i = 5 Evapotranspiration (mm)2.623.8
Table 4. Best assessment a k .
Table 4. Best assessment a k .
Soil Moisture Level Alternative   A k = 1 , 2 , , K Best   Assessment   a k = 1 , 2 , , K
Very Low A 1 0
Low A 2 0
Medium A 3 0
High A 4 0.41
Very High A 5 0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sánchez-Fernández, L.P.; Flores-Carrillo, D.A.; Sánchez-Pérez, L.A. An Integrated Approach to the Regional Estimation of Soil Moisture. Hydrology 2024, 11, 170. https://doi.org/10.3390/hydrology11100170

AMA Style

Sánchez-Fernández LP, Flores-Carrillo DA, Sánchez-Pérez LA. An Integrated Approach to the Regional Estimation of Soil Moisture. Hydrology. 2024; 11(10):170. https://doi.org/10.3390/hydrology11100170

Chicago/Turabian Style

Sánchez-Fernández, Luis Pastor, Diego Alberto Flores-Carrillo, and Luis Alejandro Sánchez-Pérez. 2024. "An Integrated Approach to the Regional Estimation of Soil Moisture" Hydrology 11, no. 10: 170. https://doi.org/10.3390/hydrology11100170

APA Style

Sánchez-Fernández, L. P., Flores-Carrillo, D. A., & Sánchez-Pérez, L. A. (2024). An Integrated Approach to the Regional Estimation of Soil Moisture. Hydrology, 11(10), 170. https://doi.org/10.3390/hydrology11100170

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