Gravity and magnetic study of the subsurface geology

in Mount Carmel and the Yizre’el Valley

M. Rybakov, A. Segev and V. Lyakhovsky

 STATE OF ISRAEL ‫מדינת ישראל‬
THE MINISTRY OF NATIONAL INFRASTRUCTURES ‫משרד התשתיות הלאומיות‬
EARTH SCIENCE RESEARCH ADMINISTRATION ‫מינהל המחקר למדעי האדמה‬

Gravity and magnetic study of the subsurface geology

in Mount Carmel and the Yizre’el Valley

Michael Rybakov, Amit Segev and Vladimir Lyakhovsky

Geological Survey of Israel

December 2009

ES-34-2009 GSI/37/2009
Table of contents

Abstract
Introduction 1
Analysis of previous gravity and magnetic surveys 2
Magnetic ground survey 5
Discussion 8
References 10
Figures 11
Hebrew abstract 27
Publication documentation page 29

List of Figures

Fig. 1. Magnetic ground coverage of the study area.

Fig. 2. Gravity coverage map.
Fig. 3. Gravity maps compiled for various searching radii (see text for details).
Fig. 4. Bouguer gravity maps.
Fig. 5. Aeromagnetic anomaly maps.
Fig. 6. Structural map of the top Judea Group.
Fig. 7. Regional Bouguer and replaced gravity maps.
Fig. 8. Replaced residual gravity map
Fig. 9. Curvature derivative of the Bouguer gravity.
Fig. 10. Curvature derivative of the replaced gravity.
Fig. 11. Curvature derivative of the replaced residual gravity.
Fig. 12. GSM-19 magnetometer at the time of magnetic profiling.
Fig. 13. Graph of the total magnetic intensity (TMI) measurements along the southern
Yizre'el Valley.
Fig. 14. Total magnetic intensity (TMI) measurements.
Fig. 15. Total magnetic intensity (TMI) anomaly map.
Fig. 16. Reduced to pole (RTP) magnetic map using the results of the ground magnetic
measurements.
Fig. 17. Upward continuation (900 m) of the ground magnetic data (color map) in compare
with the aeromagnetic data (~3000 feet above MSL; black hatched contours).
ABSTRACT

The present report concentrates on two items:

a. Analysis of the previous gravity and magnetic surveys and a recompilation of
gravity and magnetic maps. The Bouguer and the replaced gravity maps are
presented, as well as the aeromagnetic anomaly map and the curvature maps. These
regional maps help to better understand the Carmel-Gilboa fracture system and to
delineate subsurface magmatic bodies.
b. A magnetic ground survey that includes: (1) data acquisition (during 2004-2009);
(2) processing; (3) compilation of the magnetic anomaly map; (4) reduction to pole
(RTP) of the land magnetic data; (5) comparison between ground and aeromagnetic
data.

The main result of this research is the detailed (1:50,000 scale) RTP magnetic
map of the study area compiled by ~200,000 magnetic stations. The data were filtered,
gridded and presented firstly as a total magnetic intensity (TMI) map that was later
reduced to a pole (RTP) magnetic map. The pattern of the previous detected
aeromagnetic anomalies was well confirmed by upward continuation (elevation of 900
m) of the present magnetic ground survey.

The new maps provide important tools for detailed study of the subsurface
magmatic bodies and structural features.
INTRODUCTION

The most prominent Carmel structure is an elongated, tilted block extending

from the onshore into the shelf of the Haifa Bay area. This structure comprises a NNE
trending folding system, traversed by several NW trending faults. The folding is
probably part of the Late Cretaceous to Tertiary "Syrian Arc" compressional phase. The
Carmel structure onshore is separated from the Qishon graben to the north by the
Carmel fault, which has a vertical displacement of more than 1000 m. A number of
multidisciplinary studies have been carried out in this area (e.g., Ben-Gai and Ben-
Avraham, 1995, and references therein, Schattner et al., 2006, Schattner et al., 2007).
The regional gravity and magnetic data were recently reinterpreted by Rybakov et al.
(2000), using 3-D forward gravity and magnetic modeling, inversion schemes and
updated petrophysical data of the Jurassic, Cretaceous and Tertiary sedimentary layers
and volcanic rock units. These authors suggested that volcanic bodies are the main
cause of magnetic anomalies onshore and offshore northern Israel, and the magmatic
root of the Asher volcanics is most probably located close to the Carmel fault.

Active tectonics in the Yizre’el valley was previously studied by Rotstein et al.
(2004), using high-resolution seismic reflection and gravity data. They suggested that
the Yizre'el valley is a complex structure where the faults cut it at different locations
and directions and the block pattern is more complex than that of the Dead Sea rift.
They also suggested that the Yizre'el valley is a NW-SE boundary between the
Shomeron and the Lower Galilee areas.

Magnetic data acquisition (Fig. 1) and interpretation of gravity and magnetic

data from the Yizre'el Valley were commissioned by the Geological Survey of Israel
within the framework of a regional study of the Carmel-Gilboa fracture system (Segev
et al., 2006). Acquisition of gravity and magnetic data in the Yizre’el Valley has been
continued in 2009. The main purpose of this survey was to provide a detailed magnetic
map (at 1:50,000 scale) of the study area. The present report concentrates on the two
main goals:
1. Analysis of the previous gravity and magnetic surveys and a recompilation of the
gravity and magnetic maps.

1
2. A magnetic ground survey, processing, compilation of magnetic anomaly maps,
reduction to pole of the magnetic data, and comparison between ground and
aeromagnetic data.

ANALYSIS OF THE PREVIOUS GRAVITY AND MAGNETIC SURVEYS

The regional gravity and magnetic data used in this study are part of the Israeli
national database. The general gravity data coverage of the study area is shown in
Figure 2. The data coverage in the Menashe Hills (western part of the study area) and
Mount Gilboa is poor, but is better in regions covered by alluvium, where the spacing
of the measuring stations is on the order of 400–500 m. Compilation of gravity maps
using such inhomogeneous distribution of the gravity stations causes variations in
reliability at the different parts of the area. Figure 3 illustrates the maps compiled for
various searching radii. Radius options control which data points are considered by the
gridding operation when interpolating grid nodes. And although the gravity map at the
bottom right panel (Fig. 3) seems to be full, its reliability outside of well-covered
valleys is less than that inside the valleys. This should be taken into account in
interpreting gravity.

All the gravity data (reference ellipsoid 1967) were reduced to a Bouguer
density of 2670 kgr/m3. The terrain corrections for all land gravity stations were
calculated up to a 20 km distance, using a model with a 25 m grid adopted from the
Digital Terrain Model DTM (Hall, 1993). The gravity data have a reliability and
accuracy that enables interpolation to a 2 mGal contour interval. The data set was
gridded (spacing - 0.25 km) using the inverse distance method and gently smoothed.
The gravity map (Fig. 2) presents the data that can be used also for regional–residual
gravity separation, horizontal and vertical gravity derivatives calculation and gravity–
magnetic correlation. The 3-D gravity modeling and resulting map compilation are
based also on this grid.

The magnetic field was measured in 1967 and 1970 by the aeromagnetic
surveys having traverse line spacing of about 2 km, along a N55°W direction, and the
control line spacing was about 10 km in a N35°E direction, with a flying altitude of

2
1,000 m above sea level (Domzalski, 1967, Folkman, 1970). Due to analog recordings,
the 1:100,000 magnetic contour maps were compiled without the use of a computer
database. Values of the isomagnetic lines were 250, 50, 10 and 5 nanoTesla (nTl).
These magnetic-intensity maps were later digitized to create a digital magnetic
database, and a 1:500,000 magnetic anomaly map of Israel was published by Rybakov
et al. (1994B). This magnetic data set was checked for erroneous values, gridded and
slightly smoothed using the inverse distance method. The International Geomagnetic
Reference Field (IGRF) was removed and, therefore, the magnetic data represent the
local magnetic anomalies (Fig. 4). The resulting grid can be used for reduction to the
pole and upward magnetic continuation, pseudo-gravity calculation and gravity–
magnetic correlations. 3-D magnetic modeling and inversion were also based on this
grid.

The detailed Bouguer gravity map of the study area (Fig. 5) presents an
interesting pattern of the gravity anomalies with varying values from -20 to 60 mGal. A
prominent, mainly NE–SW curved and narrow gravity ridge found is associated with
the Gideon and Afula faults. Comparing gravity and geological maps (Fig. 5) one can
see relatively negative anomalies appear to be above the areas covered by a thick
succession of the Mt. Scopus and Avedat groups (Jenin area), as well as above thick
Neogene sediments (Afula graben). Contrarily, the outcropped Umm el Fahm anticline
and its partial subsurface continuation below the HaYogev–Mizra high are
characterized by relatively positive gravity values. Such a pattern is typical not only for
this area but for all of Israel (Ginzburg, 1960).

In order to use the gravity data for studying the deep subsurface structure we
had to remove the gravity effect of the post-Turonian sediments from the Bouguer
gravity. This effect was calculated in 3-D mode using three rectangular grids which
define the distribution of mass: the top surface (topography); the bottom surface (top
Judea Gr., Fig. 6) and the density. In every grid cell the effective density contrast was
calculated using the thickness-density correlation. Removing the gravity effect of the
post-Turonian sediments from the Bouguer gravity resulted in the replaced gravity map
for most parts of Israel (Fig. 7). The pattern of the local gravity anomalies in the latter
map seems to be essentially different from the pattern of the local anomalies in the

3
Bouguer gravity. For example, the prominent negative Bouguer anomaly, located in
nearest offshore Israel, is completely absent on the replaced gravity map (Rybakov et
al, 1998). Also, the Carmel gravity high and Iron local maximum has vanished
(Rybakov, 2008). However, the regional trend of the replaced gravity map is similar to
the Bouguer trend. Therefore, in order to further eliminate the residual gravity
anomalies we calculated the crustal gravity effect associated with the boundary
between the consolidated crust and mantle (Moho boundary). The geometry of the
Moho was adopted from Segev et al. (2006). This gravity effect was removed to
compile the replaced residual gravity map (Fig. 8), which is analyzed below together
with the magnetic data.

For additional enhancement of local features in the gravity field we calculated

the curvature of the gravity field using the program developed by Dr. S. Schmidt (Free
University of Berlin). The curvature is an attribute of a curve or a surface in the 2-D or
3-D space, respectively, that describes how much the curve deviates from a straight line
(planar surface), and how bent a curve/surface is at a particular point. Therefore, the
curvature is a surface-related attribute which provides insight into particular aspects or
properties of a surface which are otherwise difficult or impossible to observe. From the
methodological point of view, it is related to the second derivative of the surface, which
enhances linear features (gravity lineaments) and may help to identify general strike
directions.

This approach has been tested and applied on gravity data in various regions
(e.g., the Hula valley in northern Israel; Rybakov et al., 2003) where the curvature
attributes showed the main fault pattern. Figure 9 depicts the curvature derivative of the
Bouguer gravity in the study area. We applied the curvature approach also to the
replaced and replaced residual gravity (Figs. 10 and 11, respectively). Figure 10 shows
clearly that the curvature pattern appears to be different and simplified relative to the
previous one. Detailed analysis of the curvature maps will be carried out next stage of
the study together with the geological evidence from the surface (geological maps) and
the subsurface (high resolution seismic reflection lines).

4
MAGNETIC GROUND SURVEY

Acquisition of gravity and magnetic data in the Yizre’el Valley has been started
out in 2004 by the Geological Survey of Israel within the framework of regional studies
of the Carmel-Gilboa fracture system (Segev et al., 2006). The main purpose of this
survey was to provide a first-time compilation of a detailed magnetic map (on a
1:50,000 scale) of the study area. This section of this report concentrates on the
following items:
1. Processing of the magnetic ground measurements collected during 2004-2009.
2. Compilation of the magnetic anomaly map.
3. Reduction to pole magnetic data as a starting point for geological interpretation.
4. Comparison between ground and aeromagnetic data.

Magnetic ground surveys before 2009 were described in detail in our previous
report in the framework of BSF project 2004046 (The subsurface geology of the Lower
Galilee using magnetic and gravity data, Rybakov et al., 2008). Generally, the new
2009 survey was carried out using the same field design as well as previous surveys.
Measurements were made with a Gem Systems GSM-19 Magnetometer with automatic
total intensity magnetic field recording at 1 sec time increment, which corresponds to
about 1-m magnetic station spacing. Figure 1 illustrates overall magnetic ground
coverage of the study area. Gravity measurements were made during one survey day at
a limited number of locations (a few gravity stations on the Carmel ridge). These data
will be included in the next report.

A total of about of two-hundred thousand magnetic stations were collected to

date in the study area. Positioning of the stations in 2009 was accomplished by using
built-in GPS antenna attached above a magnetic sensor (Fig. 12). Lat/Long data were
recorded using a UTM easting/northing coordinate system (WGS84 zone 36)
simultaneously with the magnetic data. The raw data were transferred through the RS-
232 connection to a Note-Book computer and the coordinates of magnetic
measurements were converted to the Israel grid.

Figure 13 illustrates one of the recently measured magnetic lines. It should be

noted that the survey area consists of numerous man-made objects (fences, bridges,

5
metallic water pipes, electrical lines and others) which created the strong local
disturbances of the magnetic field. Daily magnetic variations have been measured in
the Bar-Gyora observatory (Survey of Israel) and used for correction of our survey
data. During the time of our magnetic measurement, magnetic variations were no more
than 15nT without the heavy spikes.

Processing of the magnetic field ground measurements included: (a) removing

drop-out points (zero readings); (b) identifying and removing spikes (non zero
anomalous measurements over a short duration); and (c) removing the measurements
with low signal quality. Values of the IGRF (International Geomagnetic Reference
Field), which corresponded to the location and time of the survey, were derived from
the Internet and removed from the values measured. The IGRF values increase in the
northern direction by about 4 nT/km.

The magnetic TMI ground data were gridded (spacing 50 m) using the inverse
distance-to-a-power technique with small and large radii that define the local
neighborhood of points to consider when interpolating each grid node. This defines the
distance in data units from the grid node that the gridding program looks to find data
points when calculating grid nodes. Data points outside the searching radius are not
considered during grid node interpolation. For example, for a searching radius set to be
smaller than half the distance between data points, a significant number of grid nodes
may be blanked in the grid file. Blanked grid nodes cause truncation of contour lines
within contour maps, as one can see in Figure 14. Based on the TMI grids, the total
magnetic intensity (TMI) magnetic anomaly maps (Figs. 14 and 15, respectively) were
compiled. And although below we only use the map based on a grid calculated with a
large searching radius (Fig. 15), the map on Figure 14 reminds us not to forget that the
magnetic coverage of the area is very uneven.

Because of the low inclination of the Earth’s magnetic field vector in the study
area, the total magnetic intensity map is difficult to interpret for location and strike of
magnetic bodies. The theoretical magnetic anomaly, caused by a body magnetized by
this normal geomagnetic field, contains the conjugated maximum and minimum, the

6
latter lying to the north. This signature is important for geological understanding of the
pattern of the observed magnetic anomalies. Magnetic maxima are displaced from their
magnetic source rocks. North-south trending magnetic anomalies are attenuated, and
east-west trending anomalies are exaggerated. Randomly distributed sources may cause
an east-west anomaly pattern. To remove the distortion caused by the low latitude, the
magnetic grid was processed using a reduction-to-pole (RTP) program. The RTP
removes anomaly asymmetry caused by inclination and locates anomalies above the
causative bodies, thus converting measurements which were recorded in an inclined
magnetic field of Earth to what would have resulted if the magnetic field were vertical.
The RTP produces a magnetic map that more directly represents the geology of the
study region. It should be noted that the RTP anomalies are better correlated with the
gravity anomalies in the case of magnetic and dense magmatics.

It is worth pointing out that no RTP results may be accurate enough in this
particular case because of remanent magnetization. Volcanic rocks are often
accompanied by strong remanent magnetization that varies in both magnitude and
direction. Unfortunately, we had to ignore the remanent magnetization because of a
lack of relevant information. Regarding the use of magnetic data for quantitative
interpretation and modeling, we have to stress that RTP is a transformation and is
designed to improve our ability to interpret magnetic data through visual inspection.
Therefore, it is strongly recommended using RTP for qualitative interpretation and not
for quantitative interpretation, even when there is no remanent magnetization and at
medium-to-high latitudes (inclinations). And the RTP grid should not be used for
quantitative interpretation as this grid is a doubly smoothed presentation of the original
magnetic data. The first smoothing was done in converting the observed magnetic
values to the TMI grid. Such double smoothing necessarily resulted in a decrease in the
magnitude and gradients of the magnetic anomalies. Therefore for quantitative
interpretation we preferred to use the TMI data directly observed along the surveying
lines. This enables better delineation of the magnetic gradients and gives more reliable
and accurate depth estimations.

7
DISCUSSION

Concluding the present report, we compared the historical aeromagnetic data

(Domzalski, 1967) with magnetic ground measurements. An aeromagnetic survey of
Israel was carried out by Spartan Air Services Ltd. (Domzalski, 1967). The surveying
parameters were: nominal traverse line spacing - about 2 km in a N55°W direction;
nominal control line spacing - about 10 km in a N35°E direction and the nominal flying
altitude – 3000 ft above m.s.l. Original flight-line data from the survey are no longer
available, and only the hardcopy 1:100,000 magnetic contour maps remain. These maps
were compiled without the use of a computer database. Values of the isomagnetic lines
were 250, 50, 10 and 5nT. To recapture the data, these maps have been digitized,
primarily along flight lines where actual data were recorded, but sometimes between
flight lines wherever necessary. Digitizing reproduces the map, but reproducing the
measured magnetic field depends on the proficiency of the mapmaker. The randomly
digitized points were gridded, converting them to a representation of the magnetic field
at equally spaced locations, which further modifies the original data. A magnetic digital
database (grid spacing 500 m) was created and a new, 1:500,000 magnetic anomaly
map of Israel was compiled and published in 1994.
For a long time the aeromagnetic data was the only magnetic source providing insights
for a better understanding of the subsurface geologic, tectonic and magmatic events in
Israel.

In a few locations (Makhtesh Ramon, Arava, Atlit, etc.) these data were
compared with separate magnetic ground profiles, however, only now do we have a
large enough area covered by ground survey and can compare the maps (Fig. 17). For
comparison the magnetic ground data was continued upward 900m, corresponding to
flight elevation. One can see that the pattern of aeromagnetic anomalies is well
confirmed by the magnetic ground survey. Some discrepancies between these maps
appear to be at the locations where there is a lack of magnetic ground data. Another
reason the upward continued and measured data is not completely similar appears to be
the difference in the spacing of the ground and aero networks: flight spacing was 2km
while the set of ground measurements was essentially denser. In these cases the
magnetic gradients and locations of magnetic extremes will be defined more accurately
by the magnetic ground data that was even continued upward to 900 m elevation.

8
 We should note that in the next stage of our study the magnetic anomaly maps
(TMI and RTP – Figs. 15 and 16) will be analyzed together with all available
geological, gravity and seismic reflection information in order to better understand the
Carmel-Gilboa fracture system and subsurface magmatic events.

9
REFERENCES

Ben-Gai, Y., Ben-Avraham, Z., 1995. Tectonic processes in offshore northern Israel
and the evolution of the Carmel structure. Mar. Pet. Geol. 12 _5., 533–548.
Domzalski, W., 1967. Aeromagnetic survey of Israel: interpretation, IPRG Report
SMA/482/67, 62 p.
Folkman, J., 1970. Aeromagnetic survey in the areas of the Golan, Shomron,, Judea
and Gaza offshore. IPRG Report SMA/760/70, 34 p.
Ginzburg, A., 1960. Geophysical studies in the central and northern coastal plain and
the western Emeq. Ph.D. Thesis. Hebrew University, Jerusalem.
Hall, J.K., 1993. The GSI Digital Terrain Model_DTM completed. In: Bogoch, R.,
Eshet, Y._Eds.., GSI Curr. Res.. pp. 47–50, Jerusalem.
Rotstein,Y., Shaliv, G. and Rybakov, M., 2004. The use of high resolution seismics in
an active tectonics study: the Yizre’el graben, Israel. Tectonophysics 382 (1-2), 31-
50
Rybakov, M., 2008. Gravity and magnetic study of the Zion's Joseph license and Asher
permit areas. Prepared for Zion Oil & Gas, Inc. Job 904/324/08.
Rybakov, M., Fleischer, L. and ten Brink, U. 2003. Structure of base fill in the Hula
basin inferred from gravity data. , Isr. J. Earth Sci.; 52:113–122
Rybakov, M., Goldshmidt V., Fleischer L. and Rotstein Y., 1998. Geological
"stripping" of the gravity field of Israel. Isr. J. Earth Sci., 47, 69-74.
Rybakov, M., Goldshmidt, V., Fleischer, L. and Ben-Gai, Y. 2000. 3-D gravity and
magnetic interpretation for the Haifa Bay area (Israel), J. Appl. Geophys. 44 (4),
353-367
Rybakov, M., Goldshmidt, V., Folkman, Y., Rotstein, Y., Ben-Avraham, Z. and Hall,
J., 1994B. Magnetic anomaly map of Israel, scale 1:500,000, Tel Aviv. GII and
Survey of Israel.
Rybakov, M. and Segev, A., 2004. The top of the crystalline basement in the Levant.
Geochem. Geophys. Geosyst., 5, Q09001, doi:10.1029/2004GC000690.
Rybakov, M., Segev, A., and Lyakhovsky V., 2008. The subsurface geology of the
Lower Galilee using magnetic and gravity data. BSF report 2004046
Schattner, U., Ben-Avraham, Z., Reshef, M., Bar-Am, G.,and Lazar, M., 2006.
Oligocene-Miocene formation of the Haifa basin: Qishon-Sirhan rifting coeval with
the Red Sea-Suez rift system. Tectonophysics 419(1-4), 1-12.
Schattner, U., Ben-Avraham, and Z., Lazar, M., 2007. Transform margin of the
northern Levant, eastern Mediterranean: From formation to reactivation. Tectonics
26, 1-17.
Segev, A., Rybakov, M., Lyakhovsky, V., Hofstetter, A., Tibor, G., Goldshmidt, V.,and
Ben-Avraham, Z., 2006. The structure, isostasy and gravity field of the Levant
continental margin and the southeast Mediterranean area. Tectonophysics 425, 137-
157.

10
FIGURES

Fig. 1. Magnetic ground coverage of the study area. White dots denote the magnetic
stations measured in the framework of magnetic surveys of the Geological
Survey. The red dots correspond to the magnetic survey carried out in Ramot
Menashe by Shefa Yamim Ltd.

11
Fig. 2. Gravity coverage map.

12
Fig. 3. Gravity maps compiled for various searching radii (see text for details).

13
Fig. 4. Bouguer gravity map in color (left panel) and the same in contours
superimposed on the geological map (right panel).

14
Fig. 5. Aeromagnetic anomaly map in color (left panel) and the same in contours
superimposed on the geological map (right panel; see text for survey design).

15
Fig. 6. Structural map of the top Judea Group.

16
A) B)

Fig. 7. A) Regional Bouguer gravity map superimposed by typical density logs. See
correlation between the main gravity anomalies and the density data.
B) Replaced regional gravity map after Rybakov et al, 1998. Insert- Tectonic
position of the study area.

17
Fig. 8. Residual replaced gravity map.

18
Fig. 9. Curvature derivative of the Bouguer gravity.

19
Fig. 10. Curvature derivative of the replaced gravity.

20
Fig. 11. Curvature derivative of the residual replaced gravity.

21
Fig. 12. GSM-19 magnetometer at the time of magnetic profiling by M. Rybakov.

Fig. 13. Graph of the total magnetic intensity measurements(TMI) (red dots) along the
southern Yizre'el Valley. Note the strong noise which doesn’t interfere with the large
positive magnetic anomaly, named Gilboa.

22
Fig. 14. Total magnetic intensity (TMI) measurements. Results of the gridding with
small searching radius.

23
Fig. 15. Total magnetic intensity (TMI) anomaly map. White dots denote location of
the magnetic measuring stations.

24
Fig. 16. Reduced to pole (RTP) magnetic map.

25
Fig. 17. Upward continuation (900 m) of the ground magnetic data (color map) in
comparison with the aeromagnetic data (~3000 feet above m.s.l.; black hatched
contours). Contour interval is 25 nT for both maps (see Figs. 1 and 15 for
location of the magnetic ground measurements).

26
 ‫תקציר‬
‫העבודה הנוכחית עוסקת בשני נושאים עיקריים‪:‬‬
‫‪ .1‬עיבוד וניתוח של סקרים גרבימטריים ומגנטיים קודמים ובניית מפות מסכמות‪ .‬מוצגות מפות‬
‫גרביטי ‪ Bouguer‬ו‪ Replaced -‬וכן מפת האנומליות המגנטיות האוויריות ומפת העקמומיות‬
‫)‪ (curvature‬של אנומליות הגרביטי‪ .‬מפות אזוריות אלה מאפשרות להבין את מערכת השבירה‬
‫של הכרמל‪-‬גלבוע וכן מצביעות על נוכחותם של גופים מגמטיים בתת‪-‬הקרקע‪.‬‬
‫‪ .2‬סקר מגנטי רגלי אשר כלל‪) :‬א( ביצוע מדידות )במהלך ‪ 2004‬עד ‪) ;(2009‬ב( עיבוד נתונים;‬
‫)ג( עריכת מפת האנומליות המגנטיות; )ד( חישוב בשיטת ‪ (RTP) reduction to pole‬של‬
‫האנומליות המגנטיות הרגליות; )ה( השוואה בין הסקר המגנטי הרגלי לבין הסקר האווירי‪.‬‬

‫התוצאה העיקרית של הסקר הנוכחי היא המפה המגנטית ‪ RTP‬המפורטת )בקנ"מ‬

‫‪ (1:50,000‬של אזור המחקר אשר חושבה מכ‪ 200,000-‬נקודות מדידה‪ .‬הנתונים עברו סינון ויצירת‬
‫שריג )‪ (grid‬של נתוני השדה המגנטי הכללי )‪ (total magnetic intensity -TMI‬ורק לאחר מכן‬
‫חושבו כ‪. RTP-‬‬
‫המארג של המפה המגנטית האווירית הישנה נמצא מתאים למדי לזה שחושב ע"י ‪upward‬‬
‫‪) continuation‬לגובה של ‪ 900‬מ'( לנתוני הסקר המגנטי הרגלי‪.‬‬
‫המפות המוצגות בדו"ח זה מהוות כלי חשוב למחקר מעמיק של גופים מגמטיים בתת‪-‬‬
‫הקרקע וכן של מאפייני המבנה של תת‪-‬הקרקע באזור עמק יזרעאל והר הכרמל‪.‬‬

‫‪27‬‬
 PUBLICATION DOCUMENTATION PAGE

1. Publication No. 2. 3. Recipient Accession No.

ES-34-2009
4.Title and Subtitle: 5.Publication Date: December 2009
Gravity and magnetic study of the subsurface
geology in Mount Carmel and the Yizre’el 6. Performing Organization Code
Valley
7. Authors: Rybakov, Segev and Lyakhovsky 8. Performing Organization Rep. No.
GSI/37/2009
9. Performing Organization Name and 10. Project/ Task / Work Unit No.
Address: GSI, 30 Malkhei Israel St.,
Jerusalem 11. Contract No. 28-17-47

12.Sponsoring Organization Name and 13. Type of report and period covered
Address:
Final report
(a)The Ministry of National Infrastructures
P.O.B. Box 13106, 91130 Jerusalem 14. Sponsoring Organiation Code

15. Supplementary Notes

16. Abstract. The present report concentrates on two items: (A) Analysis of the previous gravity
and magnetic surveys and a recompilation of gravity and magnetic maps. The Bouguer and the
replaced gravity maps are presented, as well as the aeromagnetic anomaly map and the
curvature maps. These regional maps help to better understand the Carmel-Gilboa fracture
system and to delineate subsurface magmatic bodies. (B) A magnetic ground survey that
includes: (1) data acquisition (during 2004-2009); (2) processing; (3) compilation of the
magnetic anomaly map; (4) reduction to pole (RTP) the land magnetic data; (5) comparison
between ground and aeromagnetic data.
The main result of this research is the detailed (1:50,000 scale) RTP magnetic map of
the study area compiled by ~200,000 magnetic stations. The data were filtered, gridded and
presented firstly as a total magnetic intensity (TMI) map that was later reduced to pole (RTP).
The pattern of the previously detected aeromagnetic anomalies was well confirmed by upward
continuation (elevation of 900 m) of the present magnetic ground survey.
The new maps provide important tools for detailed study of the subsurface magmatic bodies
and structural features.

17. Keywords: Mount Carmel, Yizre’el Valley, magnetic ground survey.

18. Avaibility Statement 19. Security Class 20. Security Class

(This Report) (This Page)
21. No. of the Pages 30 22. Price