Next Article in Journal
Global Warming Can Lead to Depletion of Oxygen by Disrupting Phytoplankton Photosynthesis: A Mathematical Modelling Approach
Next Article in Special Issue
Gravity Data Inversion with Method of Local Corrections for Finite Elements Models
Previous Article in Journal
Matrix of Priorities for the Management of Visitation Impacts on the Geosites of Araripe UNESCO Global Geopark (NE Brazil)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Moho Depth and Crustal Architecture Beneath the Levant Basin from Global Gravity Field Model

by
Daniele Sampietro
1,*,†,
Ahmed Hamdi Mansi
2,3,† and
Martina Capponi
2,†
1
Geomatics Research & Development s.r.l., c/o ComoNExT, Via Cavour 2, 22074 Lomazzo (CO), Italy
2
Department of Civil and Environmental Engineering (DICA), Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano (MI), Italy
3
Istituto Nazionale di Geofisica e Vulcanologia, sezione di Pisa, Via U. della Faggiola 32, 56126 Pisa (PI), Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Geosciences 2018, 8(6), 200; https://doi.org/10.3390/geosciences8060200
Submission received: 17 April 2018 / Revised: 25 May 2018 / Accepted: 1 June 2018 / Published: 2 June 2018
(This article belongs to the Special Issue Numerical Methods of Geophysical Fields Inversion)

Abstract

:
The study of the discontinuity between the Earth crust and upper mantle, the so-called Moho, and of the lithospheric architecture in general, has several important applications in exploration geophysics. For instance, it is used to facilitate the inversion of seismic-related data, in order to obtain important information on the sedimentary layers or to study the Earth’s heat flux. In this paper, the Levant crustal structure is being investigated starting from the inversion of gravity disturbances coming from a global geopotential field model based on ESA GOCE satellite mission integrated with seismic derived information. In the considered area, which is of particular interest because of its richness from the resources point of view, the deep crustal structure is still a matter of study due to the presence of a thick sequence of sedimentary layers, deposited within geological eras by the Nile River. Within the current work, the shape of the Oceanic domain in correspondence to the Herodotus Basin and the Cyprus Arc has been clearly defined. Moreover the nature of the Levantine Basin and of the Eratosthenes crust has been investigated by a set of ad hoc tests, finding the presence of continental crust. Finally, the Moho depth and the crustal density distribution have been retrieved. Several localized anomalies, in the Cyprus area, have been identified and modelled too, thus confirming the presence of heavy material, with a thickness up to 10 km, in the sedimentary layer and shallower part of the crust.

Graphical Abstract

1. Introduction

In 2012, the GOCE Exploitation for Moho Modeling and Applications (GEMMA) project [1] funded by Politecnico di Milano and the European Space Agency (ESA) has demonstrated the possibility to exploit GOCE (Gravity field and steady-state Ocean Circulation Explorer) gravity gradients to model the main features of the crust on a global scale [2]. After the GEMMA project, the topic has been furthermore investigated by several ESA projects such as the GOCE + GeoExplorer [3] and GOCE + GeoExplorer II [4]. In fact, although satellite gravity observations, being acquired at high altitudes, recover only the long-wavelength of the signal, they can be successfully used in exploration activities for the so-called “Assets and Prospects” operation, which takes place before the upstream activities. The main advantage of using satellite data for such an activity consists of yielding a regional outline of the crustal architecture that cannot be achieved either by local seismic profiles nor from airborne/shipborne gravity observations.
Within this framework, in 2015, ESA, through the Value Adding Element (VAE) program, funded a new project called Gravity Inversion to Understand LIthospheric Architecture (GIULIA) in order to investigate the possibilities to exploit GOCE-based gravity models for oil and gas explorations activities. Basically, the GIULIA project was divided into two phases: analyzing the characteristics of freely available datasets and conducting a case study of crustal structure estimation on an area of interest. The Levantine Basin has been chosen for such a purpose.
In Section 2, a brief introduction to the implemented algorithms is given, while, in Section 3, the main characteristics of the study area are outlined. In Section 4, the data reduction procedure is described as well as the procedure used to define the geological provinces and to prepare the data for the following gravimetric inversion. In Section 5, the results of the inversion, their predicted accuracy, and a brief analysis of the main sources of errors are discussed. Finally, some relevant conclusions and essential remarks are drawn in Section 6.

2. Inversion Algorithm

The inversion algorithm is based on the solution presented in [1,5]. Here, only the main work-flow (schematically represented in Figure 1), which is essentially required to understand the numerical investigations presented later on, is summarized, leaving the interested reader to consult the above literature for more details.
The first step of the procedure is the classification of the study area into geological provinces, i.e., the definition of regions where the crustal mass density can be considered homogeneous. This information can be retrieved from external sources e.g., geological provinces maps, or it can be directly inferred from the gravity field itself [5]. Each single geological province (i) is then associated with a function describing an a priori relation between the depth D and the crystalline crust mass density ( ρ C ( D , i ) ), thus creating an a priori crust density model.
Once the geological provinces map and the crystalline crust density model are generated, the gravity field observations can be reduced for the effects of the topography, bathymetry, and sediments (from a priori information) as well as for the effects of the density variations inside the crust (up to the mean Moho depth, D ¯ M ) and in the upper mantle (from D ¯ M up to a certain given depth, e.g., fixed to the expected maximum depth of the Moho obtained for instance from seismic refraction data or from global Moho models). A comment is required here: within the data reduction in fact, all the masses above the sea level are removed, while the densities from the sea level down to the maximum depth of the model are brought to a constant reference density fixed to 2670 kg/m 3 for details on the data reduction, and, of the possible problematic related to this operation, the interested reader can refer to [1].
At this point, it is possible to define the density contrast at the Moho Δ ρ and apply the inversion operator I. Here, I is thought as an operator that inverts the gravitational field in order to retrieve Δ ω = Δ ρ Δ D M , with Δ D M the unknown Moho undulation with respect to D ¯ M .
The inversion operator I can be derived in different ways: it can be obtained from the well-known Oldenburg–Parker solution of the two-layer problem [6,7], or from a least squares collocation solution [8] as well as from a Wiener filter [9]. In the present study, this last technique has been exploited. With the application of the I operator, Δ ω is retrieved and Δ ρ can be adjusted (e.g., for a scale factor for each geological province) in such a way that the Δ D M = Δ ω Δ ρ is close, in a least squares sense, to the available local seismic derived information. Thereupon, D ¯ M and ρ C D , i are updated and the data reduction is applied again in order to be coherent with the new model. This procedure is iterated until convergence is met.
Basically, the presented processing scheme allows the user to estimate a final crustal model, which, on one hand, generates the expected gravitational field and, on the other hand, is as close as possible to the available seismic-derived Moho depths.

3. The Levantine Area

The study area, namely the Levantine area, is located in the Eastern region of the Mediterranean Sea and extends from 28 E to 36 E and from 30 N to 36 N. It is mainly characterized by the presence of the Levantine Basin, a small basin surrounded by continental masses in the Easternmost part of the Mediterranean Sea and the Herodotus Basin in the Western part (see Figure 2).
The Levant Basin evolved through long and complicated processes, which started in the Late Palaeozoic, when the Eastern Mediterranean Sea was characterized by the presence of Palaeotethys, a wide Ocean separating the Africa and Eurasia tectonic plates. During the Late Palaeozoic and Early Mesozoic, deep-water basins were formed by rifting along the northern margin of the African plate [10,11,12,13]. The nature and formation of the Levant Basin have been the subject of many studies since the beginning of the century [14]. Extensive reviews of previous works and models are given, for instance, by [11,15,16].
The deep crustal structure of the Eastern Mediterranean Sea has been a matter of debate for a long time. In spite of the numerous geophysical works done in this region, e.g., [17,18,19,20,21,22], there is still no agreement on the nature of the crust nor on its age. One of the main problems in obtaining information about the deep structure in this region is the existence of a thick sedimentary sequence, which, in some areas, reaches a thickness of more than 10 km [23]. A prominent unit within the sedimentary sequence was deposited during the Messinian age, at the end of Miocene, when a desiccation of the Mediterranean Sea occurred. Although this event lasted less than 1 Ma, it led to a major erosional phase on the margins of the Mediterranean Sea and to a deposition of layers of salts and evaporites in the basinal areas [24]. In the recent years, several commercial surveys have acquired numerous multichannel seismic profiles (e.g., [19]), thus providing detailed information about most of the sedimentary sequence. However, the lower part of this sequence was not penetrated by the seismic reflection data; therefore, information about the basement and Moho depths and about the nature of the crystalline crust are still scarce. In 1983, a long-range seismic refraction acquisition was conducted with the aim to address some of the aforementioned questions [25]. The results indicated a 35 km thick continental crust, under Cyprus, and an oceanic crust of 8 km thickness, within the Levant Basin between the continental margin of southern Israel and the Eratosthenes Sea-mount. However, other aspects, such as, for example, the nature of the crust under the Eratosthenes Sea-mount and in the surrounding area still remains a matter of study [26,27].

4. The Levant Case Study: Data Reduction

As presented in Section 2, before applying the gravity inversion procedure, it is necessary to reduce the gravitational field for the effects of topography, bathymetry and sediments. In this section, the results of this reduction are presented. The gravity field data have been synthesized from the GECO (GOCE and EGM2008 COmbined) model [28], which is basically a combination of GOCE and EGM2008 [29] data. In detail, the spherical harmonics of the GECO model have been synthesized in terms of gravity anomalies, on a grid of a constant altitude of 3500 m with a spatial resolution of 1 arc-min, which coincides to about 1.5 km in the West–East direction and 1.8 km in the South–North one. Within the data reduction and the inversion phases, an extension of 2 with respect to the actual extension of the study area has been considered, in order to mitigate any expected errors related to border effects [30]. The altitude of 3500 m has been selected in such a way to be as close as possible to the Earth surface, but always outside the topographic masses.
In order to assess the Global Gravity field Model (GGM) accuracy on the considered area, its synthesized signal has been compared with a set of marine gravity observations provided by the International Gravimetric Bureau (BGI) [31]. Starting from the differences and assuming that the terrestrial observations have a 1 mGal white error, an empirical localized GECO error power spectrum has been computed, thus fully representing (at least from the stochastic point of view) the main characteristics of the GGM error itself.
In Figure 3, this empirical error power spectrum is compared to the one of the gravitational signal due to the Moho undulation. The latter has been simply obtained by computing the gravitational effects of the European Moho model [32] by means of GTE software [33,34].
As it can be seen, for wavelengths smaller than about 20 km, the observation error (blue line) is larger than the Moho signal (red line), which implies that useful information is not retrievable, from GECO, below such a threshold.
To remove the gravitational effects of the topographic layers, namely the bathymetry and topography, the ETOPO1 model [35] has been used. In detail, a complete Bouguer anomaly (see Figure 4a) has been computed by estimating the gravitational effects of topography and bathymetry with the GTE package (a FFT-based software for the computation of the terrain effect) [33,34], using a reference density value of 2670 kg/m 3 . Note that the used gravitational field has a formal resolution of about 10 km. A simple test (i.e., a residual terrain correction computed as the difference of a high-resolution terrain model and a smoothed terrain model with a smoothing average size of 10 km) show that the expected magnitude of the omission error in the area is of the order of 2 mGal (std), which drops to only 0.8 mGal if the offshore area only is considered. As a consequence, the possible error due to the different spatial resolutions of the used gravitational field and of the field generated by the data reduction is smaller than the expected observation error and therefore negligible.
Regardless that the sedimentary layer, which is known to have a drastic effects on this specific area, has not been considered yet, it can be noticed that the Bouguer anomaly map already reflects the main geological characteristics of the region. In particular, the highly positive regions (reddish area in Figure 4a) well delineate the Oceanic crust of the Herodotus Basin, crossed by the Mediterranean Ridge (an accretionary wedge caused by the African Plate that subducts under the Eurasian and Anatolian plates), as well as of the Cyprus arc. Moreover, the gravitational effects of the structures lying beneath the Eratosthenes Sea-mount and the boundaries of the Levantine Basin can be identified too. The maximum value of the signal is found in correspondence of the Cyprus arc and on the boundary between the Levantine and the Herodotus basins. The former signal is due to the presence of a crustal plate subducting beneath Cyprus [36], while the latter can be attributed to the presence of localize ophiolites structures [27,37].
In order to reduce the subduction plate effects, a simple model of the plate itself was built up. Both the geometry up to a depth of 40 km and the density contrast, fixed at 380 kg/m 3 , of such a model are taken from [36], and the corresponding gravitational effect has been obtained by means of a point mass approximation, i.e., by discretizing the plate with a set of about 6000 point masses, and summing up the effect of each point mass. The reduced signal is shown in Figure 4b. It can be seen that the applied reduction leads to a much more homogeneous signal, in which the effects of the subduction plate basically disappeared, thus proving the adequacy of the, even simple, considered model.
About the others’ localized maxima, they have been removed by exploiting the stochastic characteristics of the gravitational field itself with a Kriging procedure [38]. Basically, the effects of these local anomalies are disentangled by considering a proper mask in correspondence with the anomalies, then by estimating an empirical and theoretical variogram from the data located outside the mask, and finally by interpolating on the whole study area (applying a Kriging technique using only the data outside the mask). This procedure gives a new field in which the over-mentioned localized maxima anomalies are removed. Analyzing the difference between the new field and the original one (see Figure 5a), some information on the causative bodies can be derived. Using a brute-force method, it is possible, for instance, to estimate the parameters (depth and size in the x-, y-, and z-directions) of a prism centered in correspondence of the maximum value of the anomalous field. A constant density contrast of 650 kg/m 3 (corresponding to the difference between ophiolites and sediments density, 2850 kg/m 3 and 2200 kg/m 3 , respectively) was considered here. Once the best prism parameters have been estimated, one can reduce the observed signal for the effect of the prism itself and iterate the procedure. In Figure 5c, the inversion-based estimated thickness, positions and size are presented. It can be noticed from Figure 5b that this simple modeling allows for reducing the signal standard deviation from 10 mGal to about 2 mGal. A more detailed modeling of these anomalies, at this point, seems to therefore not be required, since the residual effects are smaller than the global gravity observation error. In any case, the results are highly-consistent with previous studies [27], with local thickness of the ophiolites up to 10 km.
Once such local anomalies have been explained, a further data reduction was performed to remove the gravitational effects of the sedimentary layer. For this purpose, the geometry and depth of the basement have been retrieved from [23]. About the density distribution, the full sedimentary succession has been divided into three parts of equal thickness, in which a constant density value was assigned to each part, specifically 2000 kg/m 3 , 2300 kg/m 3 and 2600 kg/m 3 , respectively, according to [39]. The reduced gravity field is shown in Figure 4c. It can be observed that the positive regions, reflecting the effects of the superposition of the heavy oceanic crust and the light sedimentary layers, are locally emphasized giving, after the reduction, a more demarcated field over the Herodotus Basin.
The following step is the definition of the main geological provinces. We have to recall that this step is basically required to identify areas, inside the studied region, that can be considered homogeneous (from the mass density distribution point of view). For this reason, the conceptual model developed in [39,40] was analyzed, attempting to identify the boundaries of the different regions from the reduced gravitational field plotted in Figure 4b. This objective has been fulfilled by searching for the different domain boundaries by means of the tilt angle, theta angle, and their correlation angle. Briefly, the tilt angle is defined in terms of the arctangent of the ratio of the first vertical derivative to the horizontal gradient of the potential field [41], the theta angle is the arccosine of the ratio between total horizontal derivative and the analytical signal of the gravity potential field [42], and the correlation angle is computed exploiting the amplitude and gradients of both angles [43]. For more details on these angles, the interested reader can consult, for instance, [44,45].
From the results reported in Figure 6, it can be noticed how all three of the methods are suitable to detect the edges between different geological provinces, i.e., the Herodotus oceanic crust and the Mediterranean Ridge, as well as the Cyprus Arc. Moreover, the tilt angle and theta angle maps can clearly delineate the boundaries of the Levantine Basin. The correlation angle, however, seems to suggest the presence of an Oceanic crust in correspondence with the Eratosthenes Sea-mount and the Levantine Basin (values of the correlation angle close to −180 are typical of oceanic crust, while values close to 0 are found in correspondence with the continental crust [45]); these two areas will be in any case thoroughly discussed in the following section. From the analysis of these results, in order to define the geological provinces, we updated the boundaries of the model done by [40] in accordance with the observed offshore gravity field, thus obtaining the map shown in Figure 7. With respect to the densities, the initial value of 2850 kg/m 3 is assumed for the Oceanic crust [46], while values starting from 2670 kg/m 3 with a linear gradient of 9 kg/m 3 km 1 up to 40 km in the continental domain are adopted [47]. In case of larger depths, the density of the continental crust is fixed to 3030 kg/m 3 , according to [47]. A note must be taken that the continental domain, in yellow color located to the South of Cyprus, coincides to the Eratosthenes Sea-mount.
Before applying the inversion procedure, the lateral density variation inside the mantle should also be considered. In this study, it has been taken from a global model [48] directly.
Note that the proposed data reduction scheme implies having precise information on the shape and density of topography, bathymetry, and sediments as well as the mantle. While this is well-justified for the topography and bathymetry, which are the most exterior layers of the crust and therefore probably the best known, the situation is completely different for the involved sedimentary successions and the mantle mass densities. Starting from the former, namely the sedimentary layer, the one should underline that, in the study area, its gravitational effects cannot be properly disentangled from the one of the Moho. In fact, the forward computed signal of the sediments, due to the large extension and thickness of the sedimentary layer of the Nile Basin, has a complete overlap (in terms of wavelengths) with the Moho one. Therefore, it should be retrieved from external sources. In addition, the situation of the mantle density is so similar. For this reason, in the current research, we prefer to evaluate their errors and to propagate the uncertainty in terms of accuracy in the final results.

5. The Levant Crustal Structure from Gravity Inversion

In this section, the main results obtained by applying the proposed inversion algorithm on the Levant case-study are presented.
Apart from the gravitational field observations and the geological information on the possible densities of homogeneous regions, other information can be exploited in order to apply the inversion algorithm. A set of available seismic profiles (see red lines in Figure 2) has been exploited to constrain the model. In detail, seismic profiles from [39,49] bounded by the information retrieved from the receiver functions [50,51] have been used.
The numerical assessment starts by attempting to understand the actual amount of information as well as the accuracy retrievable from solely exploiting the available seismic acquisitions. In order to achieve that a simple Moho depths map was estimated by applying the Kriging procedure to the seismic-derived dataset. The results of the Kriging interpolation and its corresponding estimation errors are reported in Figure 8.
The estimation error clearly shows that the Moho is well constrained with estimates characterized by a standard deviation (STD) smaller than 2 km, in the whole Levantine Basin with important uncertainties whose STD reaches 6 km, only in correspondence to the Cyprus Arc and the Mediterranean Ridge. The above maps will be considered as a useful benchmark to validate the inversion result. The expected Moho depth in fact should be consistent with the results of the Kriging interpolation. From the statistical point of view, we can test the hypothesis D s D g (with D s and D g the Moho depth from seismic and from gravity inversion, respectively) for each pixel, supposing a confidence level of 95% and properly considering the respective predicted accuracy σ s 2 and σ g 2 [52].
As pointed out in Section 3, among the interesting features of the study area still not completely understood is whether the Eratosthenes block and the Levantine Basin crust belong to the oceanic domain or not. For this reason, the experiment described in the following was conducted. It consisted of performing several inversions assuming different crustal density functions for the two zones and then comparing the obtained results with the information provided by the available seismic profiles. In particular, for each of the two regions, a density of 2670 kg/m 3 with a linear gradient of 9 kg/m 3 /km (typical of the continental crust [47]), a density of 2850 kg/m 3 (typical of the Oceanic crust [46]), and a constant density of 2780 kg/m 3 (i.e., of a heavy continental crust [47]) have been tested. The results are reported in Table 1 together with the statistics of the differences of the seismic-based Moho depths and those obtained from the gravity inversion for the two zones. The results of this simplified test, due to the low number of the seismic profiles and the small size of the two regions, have a low significance level from the statistical point of view. In any case, they seem to suggest the presence of a continental crust in correspondence with both the Eratosthenes Block and the Levantine Basin with a constant density value of about 2780 kg/m 3 for the former and of 2670 kg/m 3 with a density gradient of 9 kg/m 3 per km for the latter.
As for the Eratosthenes Block, the above conclusion seems to be coherent with the result obtained by [26], which reveals a complex velocity structure for Eratosthenes made by a thicker middle crust and a thinner lower crust.
Once the initial density models of these two regions are fixed, it has been possible to perform the inversion of the gravitational field data. The Moho depth estimates are presented in Figure 9, where the presence of a shallower Moho beneath the Oceanic crust of the Herodotus Basin as well as in correspondence to the Cyprus Arc are evident. As for the Levantine Basin, it presents a Moho depth ranging between 20 km in the Northernmost area to about 30 km close to the Sinai coastline. The Eratosthenes Block instead shows a general crust thickening with a Moho depth reaching a value of 28 km.
Considering now the density distribution, the inversion basically decreases the density of the Eratosthenes Block with a factor of 1.3 while leaving unchanged all the others’ density functions (variations are smaller than 0.5 %).
Comparing the obtained Moho model with the seismic-derived information, Figure 8a, a good consistency between the two models where seismic profiles are available can be seen. Considering the statistics presented above, it turns out that 99% of the two maps are consistent with a confidence level of 95%.
In Figure 10, the European Moho model [51] and the GEMMA1.0 global model [1] are chosen for comparison purposes. Neglecting the different spatial resolution of the models, it can be observed that the GEMMA1.0 model is able to retrieve the major Moho variation in the considered area, even if different long-wavelengths can be easily observed. This is probably due to the fact that, in the GEMMA1.0 global model, the Mediterranean Sea region of the study area has been modeled with a unique geological province, therefore simplifying the real geological setting. Moreover, the GEMMA1.0 model does not incorporate a detailed local sediment thickness model. As for the Moho depth of the European plate, Figure 10b, it shows a completely different shape. It should be observed that this model, which has a good consistency with the CRUST1.0 model [53], in the considered area is not properly constrained with seismic profiles. This is confirmed also by the predicted map of the Moho depth uncertainty, which ranges between 4 and more than 8 km (STD), in the study area [51].
In order to evaluate the accuracy of the obtained results, three factors should at least be considered: namely the accuracy of the inversion algorithm, the accuracy of the gravitational field data, and the accuracy of the different models used within the data reduction step. The first factor, namely the accuracy of the inversion operator, can be evaluated by means of a simple closed-loop test, where a realistic Moho model in the region of interest, e.g., obtained from the European Moho, is considered. By means of the GTE software [33], the effects, in terms of gravity disturbances, of a synthetic model made of only two layers (i.e., crust and upper Mantle) were computed considering constant densities equivalent to 2850 kg/m 3 and 3250 kg/m 3 , respectively, and separated by the Moho interface. Consequently, the inversion algorithm on the synthetic field is executed. It turns out, from this synthetic experiment, that the error related to the inversion operator is of the order of 1.7 km in terms of STD, which drops to solely 0.7 km if a 2 border (where the boundary effects are expected) is considered.
As for the gravitational field related error, it can be simply propagated by means of the inversion operator I. Assuming an observation error of 5 mGal (STD), an accuracy of 0.6 km is retrievable, in terms of Moho depth.
Considering now the last factor, namely the accuracy related to the data reduction, the effects of each component is evaluated by means of a Monte Carlo analysis, in which, for each component, a certain error variance and a correlation length are assumed. Therefore, a set of Monte Carlo samples are randomly generated. The inversion algorithm is then applied to each sample and the final accuracy is evaluated in terms of the Moho depth. The statistics are reported in Table 2, where it can be seen that the effect of the lateral variation of the Mantle density is quite small with a STD of 0.7 km, while those related to the crystalline crust are more considerable with a STD of 2.1 km.
The total accuracy of the estimated Moho depth is displayed in Figure 9b and it is evaluated to be 2.3 km.

6. Conclusions

In the present work, the main results of the Levant Basin case-study, investigated within the GIULIA project, are presented. In detail, once properly reduced the data, the GECO GGM has been primarily used to cluster the considered region into a set of homogeneous (from the mass density point of view) geological provinces and afterwards to retrieve the Moho depth together with information on the density of each province.
The results clearly define the shape of the Oceanic domain in correspondence to the Herodotus Basin and the Cyprus Arc with a very thin crust and a mass density of about 2850 kg/m 3 .
A set of ad hoc tests seems to evince the presence of a heavy continental crust in correspondence of the Eratosthenes Block characterized with a Moho depth up to about 28 km and a density of 2740 kg/m 3 , while, for the Levantine Basin, the presence of Moho depths ranges from about 20 km to 30 km. Several localized anomalies have been identified and modeled in the Cyprus area, thus confirming the presence of heavy material, with a thickness up to 10 km, in the sedimentary layer and a shallower part of the crust.
In order to assess the accuracy of the Moho depth estimates, several factors such as the uncertainties of the sediment depth model and/or the mantle density have been considered. On the one hand, the most important factor is found to be the crystalline crust density uncertainty. On the other hand, the propagation of the gravitational observation errors implies an accuracy of the order of 0.6 km in the final model, thus fully proving the goodness of the GECO GGM for these kinds of activities.

Author Contributions

D.S. supervised the study; M.C., A.H.M. and D.S. conceived and performed the study; D.S. and M.C. conducted the bibliographic research; A.H.M., M.C. and D.S. wrote the paper.

Acknowledgments

The study has been partially founded by ESA in the framework of the VAE program, contract number:4000117591. Funds do not cover the costs to publish in open access.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reguzzoni, M.; Sampietro, D. GEMMA: An Earth crustal model based on GOCE satellite data. Int. J. Appl. Earth Obs. Geoinform. 2015, 35, 31–43. [Google Scholar] [CrossRef]
  2. Sampietro, D. Crustal modelling and Moho estimation with GOCE gravity data. In Remote Sensing Advances for Earth System Science; Springer: Cham, Switzerland, 2016; pp. 127–144. [Google Scholar]
  3. Bouman, J.; Ebbing, J.; Meekes, S.; Fattah, R.A.; Fuchs, M.; Gradmann, S.; Haagmans, R.; Lieb, V.; Schmidt, M.; Dettmering, D.; et al. GOCE gravity gradient data for lithospheric modeling. Int. J. Appl. Earth Obs. Geoinform. 2015, 35, 16–30. [Google Scholar] [CrossRef]
  4. Novák, P.; Tenzer, R. Gravitational gradients at satellite altitudes in global geophysical studies. Surv. Geophys. 2013, 34, 653–673. [Google Scholar] [CrossRef]
  5. Sampietro, D. Geological units and Moho depth determination in the Western Balkans exploiting GOCE data. Geophys. J. Int. 2015, 202, 1054–1063. [Google Scholar] [CrossRef]
  6. Oldenburg, D.W. The inversion and interpretation of gravity anomalies. Geophysics 1974, 39, 526–536. [Google Scholar] [CrossRef]
  7. Parker, R.; Huestis, S. The inversion of magnetic anomalies in the presence of topography. J. Geophys. Res. 1974, 79, 1587–1593. [Google Scholar] [CrossRef]
  8. Barzaghi, R.; Reguzzoni, M.; Borghi, A.; De Gaetani, C.; Sampietro, D.; Marotta, A. Global to local Moho estimate based on GOCE geopotential model and local gravity data. In VIII Hotine-Marussi Symposium on Mathematical Geodesy; Springer: Cham, Switzerland, 2015; pp. 275–282. [Google Scholar]
  9. Sampietro, D.; Reguzzoni, M.; Braitenberg, C. The GOCE estimated Moho beneath the Tibetan Plateau and Himalaya. In Earth on the Edge: Science for a Sustainable Planet; Springer: Berlin/Heidelberg, Germany, 2014; pp. 391–397. [Google Scholar]
  10. Garfunkel, Z.; Derin, B. Permian-early Mesozoic tectonism and continental margin formation in Israel and its implications for the history of the Eastern Mediterranean. Geol. Soc. Lond. Spec. Publ. 1984, 17, 187–201. [Google Scholar] [CrossRef]
  11. Robertson, A.; Dixon, J. Introduction: Aspects of the geological evolution of the Eastern Mediterranean. Geol. Soc. Lond. Spec. Publ. 1984, 17, 1–74. [Google Scholar] [CrossRef]
  12. Ben-Avraham, Z. Multiple opening and closing of the Eastern Mediterranean and South China Basins. Tectonics 1989, 8, 351–362. [Google Scholar] [CrossRef]
  13. Robertson, A.; Dixon, J.; Brown, S.; Collins, A.; Morris, A.; Pickett, E.; Sharp, I.; Ustaömer, T. Alternative tectonic models for the Late Palaeozoic-Early Tertiary development of Tethys in the Eastern Mediterranean region. Geol. Soc. Lond. Spec. Publ. 1996, 105, 239–263. [Google Scholar] [CrossRef]
  14. Argand, E. Sur l’arc des Alpes Occidentales. Eclogae Geol. Helv. 1916, 14, 145–204. [Google Scholar]
  15. Şengör, A.C. The North Anatolian transform fault: its age, offset and tectonic significance. J. Geol. Soc. 1979, 136, 269–282. [Google Scholar] [CrossRef]
  16. Dercourt, J.; Zonenshain, L.; Ricou, L.E.; Kazmin, V.; Pichon, X.L.; Knipper, A.; Grandjacquet, C.; Sbortshikov, I.; Geyssant, J.; Lepvrier, C.; et al. Geological evolution of the tethys belt from the atlantic to the pamirs since the LIAS. Tectonophysics 1986, 123, 241–315. [Google Scholar] [CrossRef]
  17. Lort, J.; Limond, W.; Gray, F. Preliminary seismic studies in the eastern Mediterranean. Earth Planet. Sci. Lett. 1974, 21, 355–366. [Google Scholar] [CrossRef]
  18. Woodside, J. Disruption of the African Plate margin in the eastern Mediterranean. Geol. Libya Symp. Geol. Libya 1991, 3, 2319–2329. [Google Scholar]
  19. Sage, L.; Letouzey, J. Convergence of the African and Eurasian plate in the eastern Mediterranean. In Petroleum and Tectonics in Mobile Belts; Imprimerie Nouvelle: Saint-Jean-de-Braye, France, 1990; pp. 49–68. [Google Scholar]
  20. Khair, K.; Tsokas, G.N.; Sawaf, T. Crustal structure of the northern Levant region: Multiple source Werner deconvolution estimates for Bouguer gravity anomalies. Geophys. J. Int. 1997, 128, 605–616. [Google Scholar] [CrossRef]
  21. Cowie, L.; Kusznir, N. Mapping crustal thickness and oceanic lithosphere distribution in the Eastern Mediterranean using gravity inversion. Pet. Geosci. 2012, 18, 373–380. [Google Scholar] [CrossRef]
  22. Cowie, L.; Kusznir, N. Gravity inversion mapping of crustal thickness and lithosphere thinning for the eastern Mediterranean. Lead. Edge 2012, 31, 810–814. [Google Scholar] [CrossRef]
  23. Rybakov, M.; Segev, A. Top of the crystalline basement in the Levant. Geochem. Geophys. Geosyst. 2004, 5. [Google Scholar] [CrossRef] [Green Version]
  24. Hsü, K.J.; Montadert, L.; Bernoulli, D.; Cita, M.B.; Erickson, A.; Garrison, R.E.; Kidd, R.B.; Mèlierés, F.; Müller, C.; Wright, R. History of the Mediterranean salinity crisis. Nature 1977, 267, 399–403. [Google Scholar] [CrossRef]
  25. Makris, J.; Abraham, Z.B.; Behle, A.; Ginzburg, A.; Giese, P.; Steinmetz, L.; Whitmarsh, R.B.; Eleftheriou, S. Seismic refraction profiles between Cyprus and Israel and their interpretation. Geophys. J. Int. 1983, 75, 575–591. [Google Scholar] [CrossRef] [Green Version]
  26. Welford, J.K.; Hall, J.; Hübscher, C.; Reiche, S.; Louden, K. Crustal seismic velocity structure from Eratosthenes Seamount to Hecataeus Rise across the Cyprus Arc, eastern Mediterranean. Geophys. J. Int. 2014, 200, 935–953. [Google Scholar] [CrossRef]
  27. Feld, C.; Mechie, J.; Hübscher, C.; Hall, J.; Nicolaides, S.; Gurbuz, C.; Bauer, K.; Louden, K.; Weber, M. Crustal structure of the Eratosthenes Seamount, Cyprus and S. Turkey from an amphibian wide-angle seismic profile. Tectonophysics 2017, 700, 32–59. [Google Scholar] [CrossRef]
  28. Gilardoni, M.; Reguzzoni, M.; Sampietro, D. GECO: A global gravity model by locally combining GOCE data and EGM2008. Stud. Geophys. Geod. 2016, 60, 228–247. [Google Scholar] [CrossRef]
  29. Pavlis, N.K.; Holmes, S.A.; Kenyon, S.C.; Factor, J.K. An earth gravitational model to degree 2160: EGM2008. EGU Gen. Assem. 2008. 2008, 13–18. [Google Scholar]
  30. Mansi, A.H.; Capponi, M.; Sampietro, D. Downward continuation of airborne gravity data by means of the change of boundary approach. Pure Appl. Geophys. 2018, 175, 977–988. [Google Scholar] [CrossRef]
  31. BGI. The International Gravimetric Bureau for access to open worldwide gravity data IAG Geodesist’s Handbook. J. Geodesy 2012, 86. [Google Scholar]
  32. Grad, M.; Tiira, T.; ESC Working Group. The Moho depth map of the European Plate. Geophys. J. Int. 2009, 176, 279–292. [Google Scholar] [CrossRef] [Green Version]
  33. Sampietro, D.; Capponi, M.; Triglione, D.; Mansi, A.; Marchetti, P.; Sansò, F. GTE: A new software for gravitational terrain effect computation: theory and performances. Pure Appl. Geophys. 2016, 173, 2435–2453. [Google Scholar] [CrossRef]
  34. Capponi, M.; Mansi, A.H.; Sampietro, D. Improving the computation of the gravitational terrain effect close to ground stations in the GTE software. Stud. Geophys. Geod. 2018, 62, 206–222. [Google Scholar] [CrossRef]
  35. Amante, C.; Eakins, B.W. ETOPO1 1 Arc-Minute Global Relief Model: Procedures, Data Sources and Analysis; US Department of Commerce, National Oceanic and Atmospheric Administration, National Environmental Satellite, Data, and Information Service, National Geophysical Data Center, Marine Geology and Geophysics Division: Colorado, CO, USA.
  36. Ergün, M.; Okay, S.; Sari, C.; Oral, E.Z.; Ash, M.; Hall, J.; Miller, H. Gravity anomalies of the Cyprus Arc and their tectonic implications. Mar. Geol. 2005, 221, 349–358. [Google Scholar] [CrossRef]
  37. Rybakov, M.; Voznesensky, V.; Ben-Avraham, Z.; Lazar, M. The Niklas anomaly southwest of Cyprus: New insights from combined gravity and magnetic data. Isr. J. Earth Sci. 2009, 57, 125–138. [Google Scholar] [CrossRef]
  38. Wackernagel, H. Multivariate Geostatistics: An Introduction with Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  39. Longacre, M.; Bentham, P.; Hanbal, I.; Cotton, J.; Edwards, R. New crustal structure of the Eastern Mediterranean basin: Detailed integration and modeling of gravity, magnetic, seismic refraction, and seismic reflection data. In Proceedings of the EGM 2007 International Workshop, Capri, Italy, 15–18 April 2007. [Google Scholar]
  40. USGS, U. Geological Survey World Petroleum Assessment 2000-Description and Results; USGS Digital Data Series DDS-60; Geological Survey (U.S.): Reston, VA, USA, 2000.
  41. Miller, H.G.; Singh, V. Potential field tilt—A new concept for location of potential field sources. J. Appl. Geophys. 1994, 32, 213–217. [Google Scholar] [CrossRef]
  42. Wijns, C.; Perez, C.; Kowalczyk, P. Theta map: Edge detection in magnetic data. Geophysics 2005, 70, L39–L43. [Google Scholar] [CrossRef]
  43. Jilinski, P.; Meju, M.A.; Fontes, S.L. Demarcation of continental-oceanic transition zone using angular differences between gradients of geophysical fields. Geophys. J. Int. 2013, 195, 276–281. [Google Scholar] [CrossRef] [Green Version]
  44. Ma, G.; Li, L. Edge detection in potential fields with the normalized total horizontal derivative. Comput. Geosci. 2012, 41, 83–87. [Google Scholar] [CrossRef]
  45. Maurya, V.; Fontes, S.; La Terra, E. An Enhanced and Effective Technique for Demarcation of Continental—Oceanic Transition-Application for South Atlantic Conjugate Margins. In Proceedings of the 79th EAGE Conference and Exhibition 2017, Paris, France, 12–15 June 2017. [Google Scholar]
  46. Carlson, R.; Raskin, G. Density of the ocean crust. Nature 1984, 311, 555–558. [Google Scholar] [CrossRef]
  47. Christensen, N.I.; Mooney, W.D. Seismic velocity structure and composition of the continental crust: A global view. J. Geophys. Res. Solid Earth 1995, 100, 9761–9788. [Google Scholar] [CrossRef]
  48. Guerri, M.; Cammarano, F.; Connolly, J.A. Effects of chemical composition, water and temperature on physical properties of continental crust. Geochem. Geophys. Geosyst. 2015, 16, 2431–2449. [Google Scholar] [CrossRef] [Green Version]
  49. Ben-Avraham, Z.; Ginzburg, A.; Makris, J.; Eppelbaum, L. Crustal structure of the Levant Basin, eastern Mediterranean. Tectonophysics 2002, 346, 23–43. [Google Scholar] [CrossRef]
  50. Tezel, T.; Shibutani, T.; Kaypak, B. Crustal thickness of Turkey determined by receiver function. J. Asian Earth Sci. 2013, 75, 36–45. [Google Scholar] [CrossRef]
  51. Grad, M.; Tiira, T. Moho depth of the European Plate from teleseismic receiver functions. J. Seismol. 2012, 16, 95–105. [Google Scholar] [CrossRef]
  52. Sampietro, D.; Reguzzoni, M.; Negretti, M. The GEMMA crustal model: First validation and data distribution. ESA Spec. Publ. 2013, 722, 30. [Google Scholar]
  53. Laske, G.; Masters, G.; Ma, Z.; Pasyanos, M. Update on CRUST1. 0—A 1-degree global model of Earth’s crust. Geophys. Res. Abstr. 2013, 15, EGU2013-2658. [Google Scholar]
Figure 1. Computational core flow chart.
Figure 1. Computational core flow chart.
Geosciences 08 00200 g001
Figure 2. The Levant Basin with the boundaries of its main geological provinces (orange lines) and the available seismic profiles (red lines).
Figure 2. The Levant Basin with the boundaries of its main geological provinces (orange lines) and the available seismic profiles (red lines).
Geosciences 08 00200 g002
Figure 3. The Power spectrum of the Moho signal (red line) and GGM errors (blue line).
Figure 3. The Power spectrum of the Moho signal (red line) and GGM errors (blue line).
Geosciences 08 00200 g003
Figure 4. The Bouguer anomaly synthesized from the GECO global gravity model and reduced for ETOPO1 the topographic layer (a). The Bouguer anomaly corrected for the effect of the subduction plate beneath the Cyprus arc (b). The final reduced field once the effects of sediments and of the ophiolites have been removed (c). Dashed blue lines represent the boundary of the main geological provinces.
Figure 4. The Bouguer anomaly synthesized from the GECO global gravity model and reduced for ETOPO1 the topographic layer (a). The Bouguer anomaly corrected for the effect of the subduction plate beneath the Cyprus arc (b). The final reduced field once the effects of sediments and of the ophiolites have been removed (c). Dashed blue lines represent the boundary of the main geological provinces.
Geosciences 08 00200 g004
Figure 5. The isolated signal due to the presence of ophiolites (a), residuals after the modeling (b) thickness and position of the main ophiolite structures as estimated in the inversion (c).
Figure 5. The isolated signal due to the presence of ophiolites (a), residuals after the modeling (b) thickness and position of the main ophiolite structures as estimated in the inversion (c).
Geosciences 08 00200 g005
Figure 6. The tilt angle (a), theta angle (b) and their correlation angle (c). A mask has been applied onshore; unit (degree). Dashed blue lines represent the boundary of the main geological provinces.
Figure 6. The tilt angle (a), theta angle (b) and their correlation angle (c). A mask has been applied onshore; unit (degree). Dashed blue lines represent the boundary of the main geological provinces.
Geosciences 08 00200 g006
Figure 7. The boundaries of the main geological provinces within the Levant Basin. Colors from yellow to red represent the continental domains, while colors from blue to azure represent the oceanic domains.
Figure 7. The boundaries of the main geological provinces within the Levant Basin. Colors from yellow to red represent the continental domains, while colors from blue to azure represent the oceanic domains.
Geosciences 08 00200 g007
Figure 8. The Kriging interpolation of the seismic data (a) and its corresponding estimation error STD (b). Dashed blue lines represent the boundary of the main geological provinces.
Figure 8. The Kriging interpolation of the seismic data (a) and its corresponding estimation error STD (b). Dashed blue lines represent the boundary of the main geological provinces.
Geosciences 08 00200 g008
Figure 9. The estimated Moho depths (a) and its accuracy in terms of standard deviation (b). Transparent regions identify the 2 border areas, dashed blue lines represent the boundary of the main geological provinces.
Figure 9. The estimated Moho depths (a) and its accuracy in terms of standard deviation (b). Transparent regions identify the 2 border areas, dashed blue lines represent the boundary of the main geological provinces.
Geosciences 08 00200 g009
Figure 10. Moho depth interpolated from the GEMMA1.0 model (a); Moho depth from the European model (b). Dashed blue lines represent the boundary of the main geological provinces.
Figure 10. Moho depth interpolated from the GEMMA1.0 model (a); Moho depth from the European model (b). Dashed blue lines represent the boundary of the main geological provinces.
Geosciences 08 00200 g010
Table 1. Differences between seismic profiles and Moho depths derived from gravity inversion executed for different density values for the Eratosthenes Block and Levantine Basin.
Table 1. Differences between seismic profiles and Moho depths derived from gravity inversion executed for different density values for the Eratosthenes Block and Levantine Basin.
FeatureDensity (kg/m 3 )Mean (km)STD (km)
Eratosthenes Block2670 + 9D−1.30.7
2780−0.10.3
28500.80.6
Levantine Basin2670 + 9D0.20.8
27800.80.8
28500.71.4
Table 2. Effect of the considered input uncertainties on the final Moho depth.
Table 2. Effect of the considered input uncertainties on the final Moho depth.
ModelError STDCorrelation LengthEffect on the Estimated Moho Depth
mantle density40 kg/m 3 60 km0.7 km
crystalline crust density120 kg/m 3 30 km2.1 km
sediment depth1 km10 km0.6 km

Share and Cite

MDPI and ACS Style

Sampietro, D.; Mansi, A.H.; Capponi, M. Moho Depth and Crustal Architecture Beneath the Levant Basin from Global Gravity Field Model. Geosciences 2018, 8, 200. https://doi.org/10.3390/geosciences8060200

AMA Style

Sampietro D, Mansi AH, Capponi M. Moho Depth and Crustal Architecture Beneath the Levant Basin from Global Gravity Field Model. Geosciences. 2018; 8(6):200. https://doi.org/10.3390/geosciences8060200

Chicago/Turabian Style

Sampietro, Daniele, Ahmed Hamdi Mansi, and Martina Capponi. 2018. "Moho Depth and Crustal Architecture Beneath the Levant Basin from Global Gravity Field Model" Geosciences 8, no. 6: 200. https://doi.org/10.3390/geosciences8060200

APA Style

Sampietro, D., Mansi, A. H., & Capponi, M. (2018). Moho Depth and Crustal Architecture Beneath the Levant Basin from Global Gravity Field Model. Geosciences, 8(6), 200. https://doi.org/10.3390/geosciences8060200

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