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
checkpoints,
, is more complex than simply measuring the soil moisture at each checkpoint,
, 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
,
, and
. The same can be true for other weather conditions, such as solar radiation, which can even be similar at checkpoints
P0,
,
, and
, 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
, which is a grassland area, and the checkpoints
or
, which are mainly tree-covered areas. The tree-covered areas at checkpoints
or
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
are measured at the primary checkpoint (
), and then the measured weather conditions
and the rest of the inputs required by IWeCASF are processed to determine the adjusted weather conditions
at each checkpoint
. The adjusted weather conditions
for checkpoint
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
is defined as the
regional soil moisture estimation.
The locations of checkpoints 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 (
) from a set of soil moisture
point estimates (
). The integrated approach requires as inputs the satellite imagery of the estimation region, as well as the data of the soil and crop zones (
) identified within this region. Moreover, the locations of the checkpoints (
), the weather conditions measured (
) at one of the checkpoints (primary checkpoint
) and the water irrigation record (
) 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 (
) at every checkpoint (
), which are utilized by FEADM to obtain each soil moisture point estimate (
). The set of soil moisture point estimates (
) of all checkpoints (
) represents the
regional estimation of soil moisture (
).
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
sectors, resulting in pairwise sectors
, where
,
is the image width, and
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
, where
is the total number of features identified; in this case,
(
soil type,
crop type, and
development stage). The size of the feature matrices is
, defined by the
sectors established in the satellite imagery. Additionally, elements
of feature matrices
contain the number of pixels at sector
that belong to a feature [
39].
Moreover, the checkpoint location
is another input; it denotes the sector
where checkpoint
is located, and it is used to identify the particular features at this location. Additionally, the input weather conditions
are measured at the primary checkpoint
with an integrated sensor suite [
32], as defined in
Table 1.
The input irrigation water is a matrix obtained from a geographic database; it denotes the amount of irrigation water supplied in each sector from the irrigation region. This input, as a complement to the rain data , is used in the integrated approach to determine the actual water supply in a sector . The external subscripts X and Y are as follows: is the image width and 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
of the irrigation region’s satellite imagery. The landscape feature matrices
are necessary to obtain a steady adjustment factor.
is the total number of landscape features; in this case,
(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
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
, which represents a landscape feature, using the objective function
defined in Equation (1), as described in [
44].
When color segmentation is completed, the image is separated pixel by pixel into landscape features, which consist of the grassland (1), tree cover area (), buildings (), elevation (), and spatial configuration (5). However, to define the landscape feature matrices, the image is divided into sectors , where and . is the width and is the height of the satellite image divided into sectors. Thus, each element of a landscape feature matrix denotes the number of pixels at sector corresponding to the landscape feature .
B. Checkpoint landscape feature selection: In this stage, the input locations
of checkpoints
and the feature matrices
are utilized. Additionally, the landscape feature matrices
extracted from the image are required. The checkpoint location
selects from its respective matrices the elements
and
that represent the landscape feature
or the feature
at sector
. As a result of this selection, the particular feature vector
is obtained, as shown in Equation (2). The elements of this vector
are defined as the number of pixels at sector
. Then, the elements of the particular feature vector
are represented as a percentage of the particular feature vector
, as described in Equation (3).
C. Fuzzy adjustment: This is performed to determine the weather conditions
at checkpoint
after the particular feature vector
is identified. Fuzzy adjustment begins with the fuzzification of the weather conditions measured,
, and of the landscape feature vector
. A sigma function, an L-shaped function, and a triangle function are used for fuzzifying, as defined in [
39]. The fuzzy weather conditions
, where
is the current value of weather condition
, are of the form defined in Equations (4)–(8).
Moreover, to fuzzify the landscape feature vectors
, it is necessary to compare the landscape features
at the checkpoint where the adjustment is performed and the landscape features at the primary checkpoint
. This comparison
is executed using Equation (9).
Then, the elements of comparison
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
is obtained, as shown in Equation (10).
D. Landscape adjustment: This stage models the influence of the particular checkpoint features
over the weather conditions measured
. Landscape adjustment uses a fuzzy inference system (FIS), which receives as inputs the fuzzy particular feature matrix
. This FIS utilizes a compound of
IF THEN rules, where
when modeling the adjustment factor according to the actual particular feature differences between the checkpoint
and the primary checkpoint
. The landscape adjustment FIS presents as output a steady adjustment factor element
, where
= 5 is the number of weather conditions. Each output is separated into
output sets, with
according to [
39]. The result of landscape adjustment is a weather condition steady adjustment factor
for each weather condition, which is gathered as
.
E. Variable adjustment: Before developing the variable adjustment, it is necessary to determine the certainty of the weather condition replication
, which improves the modeling of inconsistent weather conditions such as rain (
). This certainty
is obtained using a distribution function, as in [
39]. Afterwards, the certainty of the weather condition replication
and the weather conditions measured
are processed through the
IF THEN rules of a second FIS (variable adjustment). Its result is associated with the output membership functions’ variables
to model the influence of the weather conditions over the rest. As a result of defuzzification, a variable adjustment
factor is obtained.
F. Final adjustment: This gives the result of the weather conditions
at checkpoint
. Both the factors steady adjustment
and variable adjustment
are used to determine the suggested adjustment of each weather condition
at checkpoint
. According to [
39], the adjusted weather conditions
are obtained as shown in Equation (11).
2.2.2. Addition of Irrigation Water
IWeCASF’s final result is the weather conditions
at any checkpoint within a region. In the first stage of the integrated IWeCASF-FEADM approach developed herein, these weather conditions
, are used as the inputs of FEADM (second stage). A particular case is the weather condition
(rain), which is complemented before being introduced to FEADM. Thus, the actual weather conditions
supplied to FEADM are defined in Equation (12), in which the irrigation water record
of checkpoint
is added to weather condition
(rain) at checkpoint
. Meanwhile, the rest of the weather conditions remain the same.
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
) in a delimited area (checkpoint
) 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 ), the best-assessed alternative according to the actual weather conditions (defined in FEADM as decision criteria ) and then estimate its actual value. Therefore, FEADM is executed at each checkpoint , in which the actual weather conditions derived from IWeCASF and modified with the addition of irrigation water are required as inputs. FEADM’s result is a point soil moisture estimate at checkpoint . 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
is obtained.
Firstly, a judgment matrix
is defined as in Equation (13). The elements
are triangle fuzzy numbers (TFN) that represent the relative importance values, which are derived from criteria pairwise comparison
, where
and
.
Then, the relative importance values
are utilized in Equation (14) to obtain the synthetic extent value
of criteria
[
45].
Afterwards, the total integral value with the index of optimism
[
32,
46,
47,
48] is used in Equation (15) to calculate the weight set
.
B. The expert evaluation starts with the fuzzification of the criteria . linguistic labels are associated with membership functions , where is the adjusted weather condition value at the checkpoint under analysis, and is its universe of discourse. The membership functions are used to fuzzify each criterion . In this way, each fuzzified criterion is defined as .
Afterwards, the fuzzified criteria vectors
are gathered in the fuzzy criteria matrix
, as shown in Equation (16). Each element
of matrix
represents a membership function
. The rows of the fuzzy criteria matrix
are a combination of the membership values of criteria
.
After obtaining the fuzzy criteria matrix
, weighting is developed to assign greater relevance to those criteria that are more important. Subsequently, the weight set
obtained in FAHP is used to weight elements
of fuzzy criteria matrix
, as defined in Equation (17). All weighted elements
are combined in the fuzzy weighted criteria matrix
shown in Equation (18).
The next step of the expert evaluation is the selective assessment, which is responsible for modeling the relationships among the fuzzy weighted criteria
. However, the selective assessment determines the criteria
, which are relevant for the assessment of each alternative
[
32].
Furthermore, the selective assessment defines a set of evaluations
, which contains the result
defined in Equation (19). The number of evaluations
is resolved by the decision makers according to the fuzzy weighted criteria
selected for the assessment of each alternative
where the selective assessment function
, defined in Equation (20), assesses the fuzzy weighted criteria
selected for the evaluations
of alternative
.
in which the elements
of the relationship modeling matrix
determine whether a fuzzy weighted criterion
is used (
) or not used (
) in the assessment
of alternative
. 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
, the best assessment for each alternative
is obtained with Equation (21).
Once the best assessments
are obtained for every alternative
, they are combined in a matrix
. The best choice
is achieved through Equation (22).
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
is defined for each alternative
, as in [
32]. Moreover, the best assessment
determines the maximum value for each output membership function
; in this way, the actual output membership function
is obtained, as shown in Equation (23).
Then, the function
is obtained using Equation (24). This function is utilized in the centroid method, shown in Equation (25). As a result, the soil moisture estimates
at checkpoint
are achieved.
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 . Consequently, FEADM is implemented a total of times, giving, as a result, soil moisture point estimates . The soil moisture point estimates are combined in a set, which is named the regional soil moisture estimation . If a more accurate regional soil moisture estimation is required, more checkpoints must be defined to reduce the area without a point estimate.
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 is higher in mostly tree-covered checkpoints, such as the checkpoint at the sector with , compared to mainly grassland checkpoints, such as , where the soil moisture content is 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
is presented as a set of soil moisture
point estimates
. In
Figure 6 and
Figure 7, each sector
is a checkpoint
so that FEADM is executed at every sector
of the region of interest. In this way, the regional estimation
is composed of
point estimates
.
A total of 921 measured weather conditions and the same number of irrigation water records 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
,
,
, and
. The estimations at checkpoint
are the most accurate; the normalized error
for this checkpoint is
. Meanwhile, the normalized error for checkpoint
is
; for checkpoint
, the error is
and, for checkpoint
, the error is
.
According to these results, and recalling the location of the primary checkpoint , the performance of the integrated approach depends on the distance between the checkpoint and the primary checkpoint ; i.e., the point estimate of soil moisture is more accurate when obtained closer to the checkpoint where the weather conditions 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
vary within a determined section of the estimation region. The influence of a change in a landscape feature
over the weather conditions
and soil moisture
is outstanding. A determined area (I, II, and III) of the region of interest is depicted. Each colored rectangle is a sector
.
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
(, rain
, solar radiation
, wind speed
, and evapotranspiration
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
. Finally,
Figure 9g presents the soil moisture
regional estimation
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
(where the weather conditions are measured). Therefore, the linguistic variables refer to the values of the weather conditions
analyzed in the sector
and whether they are lower, equal to or higher than the values of the weather conditions at the primary checkpoint
.
For example, in
Figure 9b, for the temperature
(, 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
, 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
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
(where the weather conditions are measured).
Finally,
Figure 9g compares the soil moisture
regional estimation
concerning the primary checkpoint
. 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 , 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 , 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.