In this paper, we have developed three algorithms, namely hybrid weighted
particle swarm optimization (wPSO) with the gravitational search algorithm (GSA),
known as wPSOGSA; GSA; and PSO in MATLAB to interpret one-dimensional
magnetotelluric (MT) data for some corrupted and non-corrupted synthetic
data, as well as two examples of MT field data over different geological
terrains: (i) geothermally rich area, island of Milos, Greece, and (ii) southern Scotland due to the occurrence of a significantly high electrical
conductivity anomaly under crust and upper mantle, extending from the Midland
Valley across the Southern Uplands into northern England. Even though the
fact that many models provide a good fit in a large predefined search space,
specific models do not fit well. As a result, we used a Bayesian statistical
technique to construct and assess the posterior probability density function
(PDF) rather than picking the global model based on the lowest misfit error.
The study proceeds using a 68.27 % confidence interval for selecting a region
where the PDF is more prevalent to estimate the mean model which is more
accurate and close to the true model. For illustration, correlation matrices
show a significant relationship among layer parameters. The findings
indicate that wPSOGSA is less sensitive to model parameters and produces
more stable and reliable results with the least uncertainty in the
model, compatible with existing borehole samples. Furthermore, the
present methods resolve two additional geologically significant layers, one
highly conductive (less than 1.0

The magnetotelluric (MT) method is a natural source electromagnetic method that explores various natural resources, namely hydrocarbon, minerals, geothermal prospects, groundwater, and metalliferous ores, etc. (Nabighian and Asten, 2002; Simpson and Bahr, 2005). Due to its instability, non-unique solution, and algorithm sensitivity, MT data interpretation is thought-provoking. Many researchers have attempted and developed various inversion algorithms to interpret and improve the model accuracy, convergence speed, and stability and reduce the uncertainty of the solutions (Kirkpatrick et al., 1983; Constable et al., 1987; Rodi and Mackie, 2001; Li et al., 2018; Zhang et al., 2019; Khishe and Mosavi, 2020). There are mainly two categories of the inversion algorithm: first, the local optimization methods, namely the conjugate gradient, Levenberg–Marquardt/ridge regression, the Gauss-Newton algorithm, the steepest descent, and Occam's inversion, which require a good initial guess (Shaw and Srivastava, 2007; Wen et al., 2019; Roy and Kumar, 2021). Another category is global optimization techniques (i.e., ant colony optimization, genetic algorithm, particle swarm optimization, gravitational search algorithm, and simulated annealing, etc.), which do not require an initial guess. Many researchers have carried out numerous metaheuristic optimization algorithms to invert MT data (Dosso and Oldenburg, 1991; Pérez-Flores and Schultz, 2002; Miecznik et al., 2003; Sen and Stoffa, 2013). These algorithms are inspired by natural phenomena and have various geophysical applications, including particle swarm optimization (Kennedy and Eberhart, 1995; Essa et al., 2023), the genetic algorithm (Whitley, 1994), Bat algorithm (Yang, 2010a; Essa and Diab, 2023), differential evolution (Storn and Price, 1997), biogeographically based optimization (Simon, 2008), the firefly algorithm (Yang, 2010b), grey wolf optimizer (Mirjalili et al., 2014), ant colony optimization (Colorni et al., 1991), the gravitational search algorithm (Rashedi et al., 2009), and the novel barnacles mating optimizer algorithm (Ai et al., 2022).

However, unique characteristics, namely exploration and exploitation, persist in any global optimization algorithm. For example, the particle swarm optimization (PSO) algorithm has a very high potential for exploitation, which implies that the algorithm performs well in a local search but is inferior in exploration (S̨enel et al., 2019). This suggests that the algorithm has a limited capacity to estimate the best model in an extensive search range. Because of low exploration characteristics, it gets trapped at the local minima (Mirjalili and Hashim, 2010). So, integrating the two algorithms with opposite characteristics is the best way to balance exploration and exploitation characteristics to achieve better solutions than the results obtained from an individual algorithm.

Here, we utilized weighted hybrid PSOGSA (wPSOGSA), a new global optimization method that takes into account the algorithm based on natural behavior seen in birds, fish, and insects, known as PSO, and gravity-based Newton's law (with high exploration capability), known as the gravity search algorithm (GSA). Researchers interested in artificial intelligence and developing effective optimization algorithms for comparative analysis of different metaheuristic algorithms (Pace et al., 2022) have been drawn to notable characteristics in such social behavior. The wPSOGSA, PSO, and GSA are used to estimate the resistivity distribution of 1D multi-layered Earth models using synthetic (noise-free and noisy) data for three- and four-layer cases, taken from Shaw and Srivastava (2007) and Xiong et al. (2018), respectively, and field MT sounding data for six- and four-layer cases, taken from Hutton et al. (1989) and Jones and Hutton (1979), respectively.

Furthermore, numerous (here 10 000) models that fit well are optimized to obtain the mean model, which is preceded by calculating the posterior probability density function (PDF) based on Bayesian concepts using all accepted models to find the optimal mean solution with the least uncertainty, as well as a correlation matrix to determine the relationships among the layer parameters. Thus, the research reveals that the wPSOGSA method may be utilized to provide a more accurate and reliable model with superior stability, a quick rate of convergence, and the least amount of model uncertainty.

Different MT datasets are utilized to evaluate the proposed wPSOGSA's effectiveness, sensitivity, stability, and robustness in outlining the genuine subsurface structure. These datasets are noise-free and Gaussian noise synthetic data produced for several geological formations, and two MT field data have been optimized for analysis.

To demonstrate and evaluate the robustness of the present algorithms, we
have generated a synthetic MT apparent resistivity and apparent phase dataset
without noise and with noise levels (10 % and 20 % noise) considering
a three-layer typical continental crustal model with a total thickness of
33 000 m (i.e., 33.0 km) and a resistivity of the upper crust of 30 000.0

For the second example of the synthetic data, a typical four-layer HK-type of Earth model taken from Xiong et al. (2018) is generated by forward modeling equations for the demonstration of the wPSOGSA, PSO, and GSA, and their performance is compared with improved differential evolution (IDE) results obtained by Xiong et al. (2018).

We utilized the first example of field data taken from Hutton et al. (1989) on the island of Milos, Greece. Milos is a part of the South Aegean active volcanic arc, an example of an emergent volcanic edifice (Stewart and McPhie, 2006) formed by monogenetic effusive and explosive magmatism pulses. Milos is the world's biggest exporter of bentonite, and it also has a diverse variety of metalliferous and non-metalliferous mineral reserves. It is a conserved on-land laboratory for studying shallow underwater hydrothermal ore-forming processes. The accompanying shallow subsurface hydrothermal venting fields have had significantly less attention. Dawes (1986) used magnetotelluric data to assess the resistivity structure of the geothermal area on Milos' west side. With around 3.0 km spacing, 37 MT probes in the bandwidth of 100–0.01 Hz and 12 investigations in the bandwidth of 0.01–0.0001 Hz were installed, along with various profiles that were perpendicular to the Zephyria graben in the W–E direction, as well as along the graben in the S–N direction (Hutton et al., 1989).

Another field example of MT data from Jones
and Hutton (1979) was picked to illustrate our technique from Newcastleton
(55.196

The ability to formulate an effective inversion method requires a thorough
understanding of the forward modeling technique for the issue of interest.
Factors like frequency range, actual resistivity, and layer thickness are
used to create datasets of synthetic MT apparent resistivity,

If displacement currents are not taken into account, Eq. (1) becomes

The techniques that we have used for the joint modeling of metaheuristic global optimization, namely PSO, GSA, and wPSOGSA in Step 1 and posterior Bayesian probability density function technique in Step 2, obtain the global model by utilizing the synthetic data generated by using forward modeling and field MT apparent resistivity and phase curves, depicted in the schematic diagram (Fig. 1).

Schematic diagram demonstrating the essential processes considered for the joint modeling of metaheuristic global optimization (Step 1) and the posterior PDF technique (Step 2) to obtain the global model by utilizing synthetic and field MT data.

In the present study, we have implemented a new innovative global
optimization technique known as wPSOGSA, in which swarm particles and mass
particles provide the best particle, i.e., the best model. The best model is
chosen based on the fitness of the particles, and the cost function or
objective function is used to estimate this fitness. Thus the magnetotelluric
(MT) inverse problem can be formulated through the forward modeling
operator,

The particle swarm optimization (PSO) technique is a widespread evolutionary
optimization approach for determining the optimal global solution to a
nonlinear inverse problem (Kennedy and Eberhart, 1995). This
technique is analogous to the particle's natural behavior in search of food
with the help of collaborative support from the population represented by
geophysical resistivity solutions/models (known as particles) in a swarming
group. The best model/position obtained among the particles so far is stored
for each iteration, which helps in search of the global best solution,
defined by the fitness of each particle estimated using Eq. (8). The
velocity and location of the

The gravitational search algorithm (GSA) is a metaheuristic algorithm based
on Newton's gravitational law (Rashedi et al., 2009),
which states that mass particles attract each other with a gravitational
force that is directly proportional to the product of their masses and
inversely proportional to the square of the distance between them. It
signifies that massive particles (here, particle represents the resistivity
layer model/solution) attract the neighboring lighter particles. Similar
to PSO, the gravitational search optimizer works with a population of
particles known as mass particles in the universe. Thus the best
model, solution, or particle is achieved among the mass particles. The best model
is defined by each particle's capability (i.e., the fitness), calculated
using Eq. (8). The initialization of their position in the search spaces is
given by

During execution time, the gravitational acting force on agent

Here, the gravitational constant

The following equations are used to update the particle's velocity and
location:

The weighted hybrid of PSO and GSA, known as wPSOGSA, integrates two essential characteristics, exploration (i.e., the ability of an algorithm to search the whole range of a given parameter) and exploitation (i.e., the ability to converge the solution nearest to the best solution) of the global optimization algorithm that increases its efficiency and converges the objective function to achieve global minima. The velocity and location of the particles updated in the wPSOGSA are illustrated in the schematic diagram (see Fig. 2).

The wPSOGSA combines the characteristic of social thinking of PSO and the
searching capability of GSA; thus, the particle's velocity is defined as

Flow chart of the weighted hybrid particle swarm optimization and gravity search algorithm, known as wPSOGSA (after Mirjalili and Hashim, 2010).

In a Bayesian framework, the probability distribution of the model
parameters (known as posterior probability distribution) is computed using
given observed data and models obtained from inversion. The posterior for a
model is calculated using Bayes' theorem and previous model space
information. Individual model parameter ranges are incorporated into the prior
knowledge. The two fundamental stages in the Bayesian statistics method are
the representation of previous knowledge as a probability density function
and calculation of the likelihood functional derived from data misfit
(Tarantola and Valette, 1982). Specific
characteristics, such as the best fitting model, mean model, and correlation
matrix, may be determined from posterior distribution of models. According to
Bayes' theorem,

After the application of the PDF, the study further proceeds by choosing a
confidence interval (CI) of 68.27 % that is based on the empirical rule,
known as the 68–95–99.7 rule (Ross, 2009). The model parameters below
68.27 % CI are discarded, and the remaining parameters are used for
determining the mean model and uncertainty. Thus, the mean model (

The effectiveness, sensitivity, stability, and robustness of the proposed wPSOGSA in identifying the authentic subsurface structure are evaluated using various MT datasets. These datasets consist of synthetic data with no noise and Gaussian noise, which simulate different geological formations. Additionally, two MT field datasets have been optimized for analysis.

The first example of synthetic MT data was executed for 10 000 runs, keeping the same lower and upper bounds as given in Table 1 and iterations to 1000. Figure 3 shows (a) the observed apparent resistivity with the computed data; (b) the observed apparent phase with the computed data; (c) the 1D inverted model by wPSOGSA (red color), GSA (green color), and PSO (blue color) with a true model (black color); and (d) the relation between the misfit and iterations for the noise-free synthetic data.

The inverted MT response by PSO (blue color), GSA (green color),
and hybrid wPSOGSA (red color) with a true model (black color) over
three-layer synthetic data as shown in

The misfit curve as shown in Fig. 3d is gradually decreasing with
increasing iterations and becomes constant where the algorithm converges.
The PSO, GSA, and wPSOGSA converge at iterations 492, 35, and 316, with
associated errors of

The 20 % noisy synthetic MT data were executed for 10 000 runs, keeping the same lower bound, upper bound, and iterations. The well-fitted inverted MT response (see Fig. 4) was as follows: (a) the corrupted synthetic and calculated apparent resistivity data, (b) the corrupted synthetic and calculated apparent phase data, (c) the inverted 1D depth model, and (d) convergence response in terms of misfit error versus iterations. We analyzed Fig. 4d and found that the PSO, GSA, and wPSOGSA using noisy synthetic data converge at iterations 236, 7, and 73, with associated errors of 0.0394, 0.0408, and 0.0393, respectively.

The inverted MT response by PSO (blue color), GSA (green color),
and hybrid wPSOGSA (red color) with a true model (black color) over
three-layer synthetic data with 20 % random noise as shown in

The posterior Bayesian probability density function (PDF) with 68.27 % CI for wPSOGSA for three-layered synthetic data.

Two methods are used to estimate a mean solution and uncertainty: one method is the mean solution for all accepted best-fitted solutions acquired from 10 000 runs for all three global optimization techniques; another method is the model derived from all approved solutions using the posterior Bayesian PDF within 1 standard deviation. To get the global best solutions in our study, we incorporated the posterior PDF based on the Bayesian approach to enhance the efficacy of the inverted model and minimize the uncertainty in the model. The process for obtaining the mean solution proceeds by selecting an initial threshold error, which is essential because the smaller the threshold value, the more significant the number of models with less uncertainty in the model parameters (Sharma, 2012). Thus, a more considerable threshold gives a lower number of models with enormous uncertainty in the model parameter (Sen and Stoffa, 1996; Sharma, 2012). The study proceeds by calculating the PDF for each parameter value using Eq. (22). In order to select values of each parameter that having a higher posterior PDF, a 68.27 % CI is used. The mean model obtained from selected model parameters is near to the actual model.

Figure 5 shows the output of the posterior Bayesian PDF, which selects model
parameters with less error. The straight lines (dashed lines) present the
actual value of the respective layer parameters. The first-layer thickness,
second-layer thickness, and first-layer resistivity have higher
uncertainties, i.e., 61.25 m, 51.47 m, and 210.61

True model, search range, and inverted layer parameters by wPSOGSA, GSA, and PSO using different amounts of noisy (0 %, 10 %, and 20 %) synthetic MT apparent resistivity and apparent phase datasets for a three-layer Earth model.

Table 1 shows the inverted layer parameters using wPSOGSA, GSA, and PSO for noise-free and noisy synthetic MT data based on the posterior Bayesian PDF, as well as the actual model and the search range. In addition, layered properties of synthetic data corrupted with 10 % and 20 % random noise are compared and statistically analyzed. Our findings, as shown in Table 1, were compared to those obtained using the genetic algorithm (GA), ridge regression (RR), and PSO by Shaw and Srivastava (2007), which consistently outperform GA and RR, which is closer to the genuine model.

Histogram of selected models for misfit error below a defined threshold error of wPSOGSA.

The mean values of the accepted model parameters are

As well as being based on the histograms (see Fig. 6), the posterior PDF and
uncertainty of the inverted layer parameters resistivity (

Cross-plots of thickness and resistivity of first layer for the
three-layered synthetic resistivity model using

The accepted models, which have a posterior PDF value within 68.27 % CI, are used to calculate the correlation matrix. This correlation matrix gives the relationship among model parameters. Thus, the lower correlation value gives weak relation among the parameters and vice versa. The correlation matrix of PSO, GSA, and wPSOGSA was examined on one set of synthetic data, as shown in Figs. 8, 9, and 10, demonstrating the sensitivity among inverted model parameters. The value of correlation matrix, 1.0, indicates that the two parameters are strongly correlated.

Figure 8 shows that first-layer resistivity is correlated highly positively
with a first-layer thickness (0.97) and second-layer thickness (0.98), while
the second-layer resistivity (

Figure 9 indicates that first-layer resistivity is highly associated with
second-layer thickness (1.00) and weakly with second-layer resistivity
(

Correlation matrix calculated from the PSO inverted model using three-layer noise-free synthetic MT apparent resistivity and apparent phase data.

Correlation matrix calculated from the GSA inverted model using three-layer noise-free synthetic MT apparent resistivity and apparent phase data.

Figure 10 shows the correlation matrix of wPSOGSA. The analyses reveal that
the first-layer resistivity is strongly negative with the second-layer
resistivity, substantially negative (

Correlation matrix calculated from the wPSOGSA inverted model using three-layer noise-free synthetic MT apparent resistivity and apparent phase data.

We used two different search ranges for stability evaluation of the proposed wPSOGSA and executed the algorithms over three layers of synthetic MT data, one of which is expanded, and the other is contracted by 10 % of the initial search range. We infer from three layers of synthetic data that results fluctuate by approximately 3 % from the true value when the search range is changed. This variation is about 10 % on average for synthetic data corrupted with 30 % random noise, as shown in Table 2.

Stability analysis of a hybrid algorithm for three layers of synthetic data.

For the second example of the synthetic data, a typical four-layer HK-type of Earth model to analyze the performance of the present algorithm with improved differential evolution (IDE) results obtained by Xiong et al. (2018). Analysis over noisy synthetic data is done by corrupting synthetic data with 10 % and 20 % Gaussian random noise to mimic the real field data because different types of noises influence apparent resistivity data. Following that, all three optimization methods are run using the noisy synthetic data. As the misfit error increases with the noise in the data, the Bayesian PDF of 68.27 % CI is calculated with respect to the threshold misfit error of 0.01, and thus the mean model is calculated.

Enormous uncertainty is shown in the inverted results; hence, we calculated the mean model for 68.27 % CI using the posterior Bayesian PDF to reduce the uncertainty and produce the global best solution. The optimized results obtained from the posterior PDF and the true model are shown in Table 3. Figure 11 illustrates that the inverted responses for PSO, GSA, and wPSOGSA are well-fitting as follows (a) observed and calculated apparent resistivity data, (b) observed and calculated apparent phase data, (c) 1-D depth model, and (d) convergence response of present algorithms. We have estimated the layer parameters for synthetic data corrupted with 20 % random noise for comparative analysis and found that PSO, GSA, and wPSOGSA converge at iterations 96, 556, and 187, with associated errors of 3.69, 4.04, and 3.69, respectively.

The inverted MT response by PSO (blue color), GSA (green color),
and hybrid wPSOGSA (red color) with a true model (black color) over
four-layer synthetic data as shown in

Additionally, the synthetic data corrupted with 10 % random noise are also used for the execution of inversion, keeping the search range, a number of particles, and iterations the same as before; it was observed that PSO, GSA, and wPSOGSA converge at iterations 151, 2, and 250, with associated errors of 1.7609, 1.95, and 1.76, respectively. The posterior Bayesian PDF for threshold data with 68.27 % CI is calculated similarly to a three-layer case to minimize the uncertainty in inverted results.

Comparison of the result obtained from improved differential evolution (IDE) and inverted results of PSO, GSA, and hybrid wPSOGSA obtained using the posterior PDF for four-layer synthetic apparent resistivity data with different Gauss noise levels (0 %, 10 %, and 20 %) and the true model.

For the stability evaluation of presented algorithms over four layers of synthetic MT data, similar to the three-layer case, we used two different search ranges and executed the algorithms for 1000 iterations. The method exhibits good results with four layers of synthetic data and reveals minimal variation for noise-free data. For 30 % contaminated data, the variation is approximately 10 % and 12 % in Case 1 and Case 2, respectively. The outputs do not change much across runs and provide consistent results, as shown in Table 4.

Stability analysis of a hybrid algorithm for four layers of MT synthetic data.

In one-dimensional MT data for site G5 near borehole M2
(Hutton et al., 1989), as shown
in Fig. 12, the apparent resistivity and phase values are inverted using
wPSOGSA, PSO, and GSA, keeping the same set of controlling parameters as for
noisy synthetic data, such as the swarm size, inertia weight (

The location of the MT site and geology of the island of Milos, Greece (after Stewart and McPhie, 2006).

Figure 13 shows the calculated data and model parameters as (a) a match
between observed and computed apparent resistivity data and (b) a match between
observed and computed apparent phase data; (c) the 1D inverted model; and (d) the convergence response of wPSOGSA (red color), GSA (green color), and PSO
(blue color), along with the true model (black color). Figure 13c
depicts alluvium deposits with a resistivity of 1.0

Figure 13d reveals that the algorithms converge at iterations 218, 1, and 425, with corresponding errors of 0.0494, 0.0518, and 0.0493 for PSO, GSA, and wPSOGSA, respectively. The hybrid algorithm has the least error between observed and computed data. The algorithms are executed for 1000 iterations and 10 000 models, and findings are compared with available stratigraphy. The result is derived using the Monte Carlo technique by Hutton et al. (1989). After examining our optimized effects from Fig. 13 and Table 5, hybrid wPSOGSA outperformed PSO and GSA.

The inverted MT response by PSO (blue color), GSA (green color),
and hybrid wPSOGSA (red color) with a true model (black color) over the
geothermal area, island of Milos, Greece, as shown in

A posterior Bayesian method determines the global model and related
uncertainty. Figure 14 shows another uncertainty study that examined the
six-layered resistivity model over the geothermal field, island of Milos,
Greece, and found that the peak values of the posterior PDF for all model
parameters are much closer to the actual value of the layer parameters,
providing less uncertainty. We have analyzed the wPSOGSA inverted results
from Fig. 14 and Table 5 and found that the first, second, third,
fourth, fifth, and sixth layers' resistivity with uncertainty in associated
layer parameters is

The posterior Bayesian probability density function (PDF) with 68.27 % CI for wPSOGSA over a geothermal field, island of Milos, Greece.

Table 5 compares optimized results obtained from all three presented algorithms based on the posterior Bayesian PDF under 68.27 % CI condition. However, the 1D depth model inverted from wPSOGSA shows good agreement with the available borehole M-2 (Hutton et al., 1989). As a result, the hybrid algorithm is functioning better, and the findings are encouraging.

Search range and inverted results by the posterior PDF (68.27 % CI) and PSO, GSA, and hybrid wPSOGSA for six-layered field data.

Here, a similar study of the correlation matrix is carried out for field example from the island of Milos, Greece, using all accepted models, which have posterior PDF values within 68.27 % CI. The correlation matrix of PSO, GSA, and wPSOGSA was examined over the field MT data as shown in Figs. 15, 16, and 17, demonstrating the sensitivity among inverted model parameters, and found an almost similar correlation among the layer parameters for three-layer synthetic study. From correlation analyses, we noticed that the values are showing moderate and weak correlation among parameters in the wPSOGSA case, indicating that wPSOGSA is linearly independent of layer parameters. This indicates that the parameter is less affected by other layer parameters and resistivity curves, whereas the correlation among layer parameters for field data using GSA and PSO is either strongly positive or strongly negative, which shows that the parameters are dependent on each other. Thus a change in one parameter affects the other, and the apparent resistivity curve is also very much involved.

Correlation matrix of field data taken from the geothermally rich area of the island of Milos, Greece, for PSO.

Correlation matrix of field data taken from the geothermally rich area of the island of Milos, Greece, for GSA.

Correlation matrix of field data taken from the geothermally rich area of the island of Milos, Greece, for hybrid wPSOGSA.

Another field example of MT data was picked to illustrate our technique from
Newcastleton (55.196

During 9 d, in the frequency range of 0.1 to 0.0001 Hz, the
variations of the magnetic and telluric fields concerning the time at four
sites along a line perpendicular to the anomaly's strike were recorded,
keeping a high signal-to-noise ratio, where the anisotropy ratios were very near
to 1, and the skew factor is less than 0.1 for the majority of periods. Due
to low anisotropy ratios and the skew factor, the resistivity distribution under
this location is one-dimensional (Jones and
Hutton, 1979). Here one set of MT data is inverted using PSO, GSA, and
wPSOGSA to obtain the best fitting apparent resistivity curve, apparent
phase curve, and 1D depth model, as shown in Fig. 19a, b, and c, respectively. Figure 19 shows a realistic
one-dimensional resistivity variation, with a phase response ranging from
60

The execution time for wPSOGSA (33 s) is the lowest as compared to GSA (34 s) and PSO (53 s). The convergence iterations of PSO, GSA, and wPSOGSA are 79, 101, and 65, and its associated misfit errors are 3.79, 4.72, and 3.70, respectively.

The inverted MT model is illustrated in Fig. 19c, which depicts
two low conductive zones at a depth of 21 and 400 km. The first
conductive layer (70

The location of the MT site and the geology of the Southern Uplands, Scotland (after Leslie et al., 2015).

The inverted MT response by PSO (blue color), GSA (green color),
and hybrid wPSOGSA (red color) with a true model (black color) over
Newcastleton, southern Scotland, as shown in

Search range, inverted results by the posterior PDF (68.27 % CI) using PSO, GSA, and wPSOGSA for field data.

The analysis on the two synthetic MT datasets shows that the proposed algorithm works very well and provides encouraging well-fitted calculated apparent resistivity and apparent phase data with the associated observed data. Also from the study, it is noted that the proposed algorithm is less sensitive to the search range and the constraints used in this algorithm. And the comparison of the model parameters shows that the output from the proposed algorithm is very precise in comparison to the true model and faster than the individual algorithms and the other algorithms used in the previous papers.

Inverted results by the Dawes method, Jupp–Vozoff method, the
Fischer inversion method (resistivity in

A research study conducted by Hutton et
al. (1989) compared various techniques, including Parker H+, Dawes (a
combined Monte Carlo–Hedgehog approach developed by Dawes), Jupp–Vozoff,
Fischer, and Parker D+ inversions, to interpret the subsurface geology of
the island of Milos, Greece. Among these techniques, the Dawes algorithm
proved to be the most effective in identifying a reservoir interface at the
borehole's depth. Consequently, these methods were utilized to compile 1-D
models along several traverses. All three models, Parker H+, Parker D+,
and Dawes, indicate the presence of a resistivity boundary at a depth of
approximately 1000 m (Hutton et
al., 1989), where the geothermal reservoir has been detected. The research
findings through the proposed algorithm reveal that the resistivity of the
crystalline basement beneath the geothermal site is abnormally low
(

Figure 19c presents the inverted MT model of the Southern
Uplands, Scotland, displaying two zones of low conductivity at depths of 21
and 400 km. Contrary to the findings of
Jain and Wilson (1967), there
is strong evidence suggesting that the conducting zone (70

The study presented wPSOGSA along with PSO and GSA and evaluated their efficacy and applicability to the MT data. These algorithms narrate the appraisal of 1D resistivity models from apparent resistivity, apparent phase, and frequency datasets. So, synthetic and field MT data from various geological terrains were used to demonstrate the relevance of these methods, which are further carried out by applying multiple runs, generating a large number of models that fit the apparent resistivity and apparent phase curves. Furthermore, these best-fitting models within a specified range are then chosen for statistical analysis. The statistical analysis includes the posterior PDF based on the Bayesian approach with a 68.27 % CI, the correlation matrix, and stability analysis in order to understand the accuracy of the mean model and its uncertainty. However, the solution from the posterior PDF based on the Bayesian of wPSOGSA is better than GSA and PSO, explaining the reliability of the proposed inversion algorithm. In general, conventional techniques can effectively resolve the model in random noise, but they can miscarry in methodical error or inappropriate models. Also, the performance of the proposed algorithms on field datasets has been analyzed based on the mean model, uncertainty, correlation, and stability of layered Earth models, and it was found that the results obtained from wPSOGSA are reliable, stable, and more accurate than the available results, which are well adapted to borehole lithology.

The datasets used for the present study and analysis have been taken from published papers, e.g. Shaw and Srivastava (2007), Xiong et al. (2018), Hutton et al. (1989), Jones and Hutton (1979).

M: conceptualization of the study, methodology, computer code, analysis, and drafting of the manuscript. KS: methodology, computer code, analysis, and drafting of the manuscript. UKS: supervision, suggestions, and editing.

The contact author has declared that none of the authors has any competing interests.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The authors would like to express their gratitude to the IIT(ISM), Dhanbad, for providing a pleasant environment to pursue this study and support for the research. We also express our gratitude to the chief editor, associate editor, and anonymous reviewers, whose suggestions and comments enabled us to better understand the issues and considerably improve our manuscript.

This paper was edited by Norbert Marwan and reviewed by two anonymous referees.