Structure of the earth''s crust on the territory of the U.S.S.R.

The compilation of statistical data for 269 seismic crustal sections (total length: 81,000 km) which are available in the U.S.S.R. has shown that the preliminary conclusions drawn on relations between the elevation of the surface relief and Bouguer anomalies on one hand and crustal thickness (depth to the M-discontinuity) on the other hand are not fulfilled for the continental part of the U.S.S.R. The level of isostatic compensation has been found to be much deeper than the base of the earth's crust due to density inhomogeneities of the crust and upper mantle down to a depth of 150 km.

The results of seismic investigations have revealed a great diversity of relations between shallow geological and deep crustal structures:

Changes in the relief of the M-discontinuity have been found within the ancient platforms which are conformable with the Precambrian structures and which can exceed 20 km. In the North Caspian syneclise, extended areas devoid of the “granitic” layer have been discovered for the first time in continents. The crust was found to be thicker in the syneclises and anteclises of the Turanian EpiHercynian plate. In the West Siberian platforms these relations are reversed to a great extent.

Substantial differences in crustal structure and thickness were found in the crust of the Palaeo zoides and Mesozoides. Regions of substantial neotectonic activity in the Tien-Shan Palaeozoides do not greatly differ in crustal thickness if compared to the Kazakhstan Palaeozoides which were little active in Cenozoic time. The same is true for the South Siberian Palaeozoides.

The Alpides of the southern areas in the U.S.S.R. display a sharply differing surface relief and a strongly varying crustal structure. Mountains with roots (Greater Caucasus, Crimea) and without roots (Kopet-Dagh, Lesser Caucasus) were found there.

The Cenozoides of the Far East are characterized by a rugged topography of the M-discontinuity, a thinner crust and a less-pronounced “granitic” layer. A relatively small thickness of the crust was discovered in the Baikal rift zone.

The effective thickness of the magnetized domains of the crust as well as other calculations show that the temperature at the depth of the M-discontinuity (i.e., at depths of 40–50 km) is not higher than 300–400° C for most parts of the U.S.S.R.     

Crustal investigations of the U.S.S.R. by means of earthquake-generated converted waves

The method of earthquake-generated converted waves which is based on the simultaneous recording of longitudinal (P), transverse (S) and converted (PS and SP) waves has been used in the U.S.S.R. since 1956. Converted phases generated at crustal and upper-mantle discontinuities in the seismic focal area are carefully analyzed. In this paper some dynamic characteristics of transmitted PS- and SP-waves arising at different types of boundaries are described.

From a comparison of the properties of the recorded converted waves with the results of theoretical computations conclusions can be drawn with regard to the possible structure of the “exchange” (conversion) boundaries. Seismic cross-sections are presented which illustrate the high effectiveness of the method in regional investigations of crustal structure. Owing to its significance and the multitude of observations (over 22,000 km of observation lines) the method is one of the principal seismic techniques used in the U.S.S.R. which give quantitative information on crustal layering.     

Crustal structure of the Rhinegraben area

Numerous ge ological and geophysical investigations within the past decades have shown that the Rhinegraben is the most pronounced segment of an extended continental rift system in Europe. The structure of the upper and lower crust is significantly different from the structure of the adjacent “normal” continental crust.

Two crustal cross-sections across the central and southern part of the Rhinegraben have been constructed based on a new evaluation of seismic refraction and reflection measurements. The most striking features of the structure derived are the existence of a well-developed velocity reversal in the upper crust and of a characteristic cushion-like layer with a compressional velocity of 7.6–7.7 km/sec in the lower crust above a normal mantle with 8.2 km/sec. Immediately below the sialic low-velocity zone in the middle part of the crust, an intermediate layer with lamellar structure and of presumably basic composition could be mapped.

It is interesting to note that the asymmetry of the sedimentary fill in the central Rhinegraben seems to extend down deeper into the upper crust as indicated by the focal depths of earthquakes. The top of the rift “cushion” shows a marked relief which has no obvious relation to the crustal structure above it or the visible rift at the surface.     

On certain regularities of crustal structure associated with the major geologic features of Southeastern Europe

The analysis of numerous seismic studies from various geological provinces has demonstrated that variations in crustal thickness depend primarily on the thickness of the “basaltic” layer. In some areas two M discontinuities can be found — the present one and an ancient one. The lower crust, formed in Proterozoic time is apparently still preserved. Roots exist under the former Proterozoic orogens, in spite of the complete denudation of the orogenic mountains. Younger (Paleozoic-Mesozoic) subsurface structures are not so clearly pronounced in the crustal structure. More active reconstruction of the crust seems to have taken place in the course of Alpine orogenesis.     

The deep-seated crustal structure in the western part of the Black Sea and adjacent areas: Seismic reflection measurement

In recent years the northwestern Black Sea has been investigated by a great number of geophysical methods. Charts of the M discontinuity and (isopachous) charts of the “granitic”, the “basaltic”, the Paleozoic, the Jurassic-Triassic, the Upper and Lower Cretaceous and the Eocene layers were plotted based on the results of the combined data of these investigations together with associated drilling data. The data for different velocity levels confirms the concept of layered-block structure of the crust, where large blocks are divided by deep faults penetrating to the upper mantle. Sedimentation within each block is continuous while reverse fault zones, dividing the East European Platform with a crustal thickness of more than 40 km and the Scythian Platform with a crust of about 30 km thick, and the latter from the Black Sea depression with crust of about 20 km, are discontinuous. Therefore, one can speak of a continuous-discontinuous nature of the sedimentation.

An inverse relationship in thicknesses of the “granitic” and sedimentary layers has been established. In places of intensive sedimentation the thickness of the “granitic” layer is less than that within the stable unsagging blocks. On the whole the greater the thickness of “basaltic” layer, the greater is the crustal thickness.

The relationship between the main geological structures of the area should be sought in the nature of structure of these “granitic” and “basaltic” layers.     

Deep structure under Yellowstone National Park U.S.A.: A continental “hot spot”

In order to understand the origin of long-lived loci of volcanism (sometimes called “hot spots”) and their possible role in global tectonic processes, it is essential to know their deep structure. Even though some work has been done on the crustal, upper-mantle, and deep-mantle structure under some of these “hot spots”, the picture is far from clear. In an attempt to study the structure under the Yellowstone National Park U.S.A., which is considered to be such a “hot spot”, we recorded teleseisms using 26 telemetered seismic stations and three groups of portable stations. The network was operated within a 150 km radius centered on the Yellowstone caldera, the major, Quaternary volcanic feature of the Yellowstone region. Teleseismic delays of about 1.5 sec are found inside the caldera, and the delays remain high over a 100 km wide area around the caldera. The spatial distribution and magnitude of the delays indicate the presence of a large body of low-velocity material with horizontal dimensions corresponding approximately to the caldera size (40 km × 80 km) near the surface and extending to a depth of 200–250 km under the caldera. Using ray-tracing and inversion techniques, it is estimated that the compressional velocity inside the anomalous body is lower than in the surrounding rock by about 15% in the upper crust and by 5% in the lower crust and upper mantle. It is postulated that the body is partly composed of molten rock with a high degree of partial melting at shallow depths and is responsible for the observed Yellowstone volcanism. The large size of the partially molten body, taken together with its location at the head of a 350 km zone of volcanic propagation along the axis of the Snake River Plain, indicates that the volcanism associated with Yellowstone has its origin below the lithosphere and is relatively stationary with respect to plate motion. Using our techniques, we are unable to detect any measurable velocity contrast in the mantle beneath the low-velocity body, and, hence, we are unable to determine whether the Yellowstone melting anomaly is associated with a deep heat source or with any deep phenomenon such as a convection plume, chemical plume, or gravitational anchor.     

The Oslo Rift: P-T relations and lithospheric structure

Extrusion temperatures for basaltic lavas in the Permo-Carboniferous Oslo Rift, estimated from whole rock major element compositions, are estimated to be 1270 to 1340°C. This means that magmatism during the Oslo rifting event was not associated with a large temperature anomaly in the underlying upper mantle. Partial melting is believed to be caused by a combination of crustal extension, a weak temperature anomaly in the underlying asthenosphere, and/or high fluid-contents in the mantle source region (“wet-spot”). Petrological and gcochemical data imply that large masses of cumulate rocks were deposited in the deep crust during the Oslo rifting event. The densities and seismic velocities (V_p) of these cumulate rocks are estimated to be 2.8–3.5 g/cm³ and 7.5–8.0 km/s. A rough estimate suggests that cumulus minerals alone account for a net transfer of at least 2 × 10¹⁷ kg of magmatic material from the mantle into the deep crust. In addition comes material representing

1. (a) cumulate minerals corresponding to eroded magmatic surface and subsurface rocks

2. (b) intercumulus material, and

3. (c) magmas crystallized to completion in the deep crust.

Estimates based exclusively on geophysical data tend to underestimate the true transfer of mass into the lower crust as gabbroic cumulate rocks, and melts crystallizing to completion in the lower crust have densities and seismic velocities similar to those of lower crustal wallrocks.     

Main features of crustal structure in western and southern Europe based on data of explosion seismology

During the past two decades deep seismic sounding measurements have been carried out in western and southern Europe, mainly using the refraction method. These investigations were performed partly on a national basis but as well within international cooperative programs under the sponsorship of the European Seismological Commission.

In France, a systematic study has been executed to determine the main feature of deep structures under the Central Massif and the Paris Basin. In the Forez and Margeride regions, the sub-crustal velocity is lower (7.2 km/sec) than the normal value (8.0 km/sec) observed in the adjacent areas.

The central and southern part of Western Germany is covered by an extensive network of refraction profiles. The crustal thickness varies, similarly to France, from 25 to 35 km. A great amount of deep reflection data was obtained by commercial and special reflection work. The crust beneath the Rhinegraben area shows the typical “rift system” structure with a low subcrustal velocity (7.4–7.7 km/sec).

Very intensive refraction work has been carried out in the Alpine area. The maximum crustal thickness found near the axis of the negative gravity anomaly is about 55–60 km. Furthermore, a clear lowvelocity layer at a depth between 10 and 30 km has been detected. A key position with regard to the geotectonic structure of the Alps is held by the zone of Ivrea characterized by a pronounced gravity high. From the refraction work it may be concluded that there material of the lower crust and the upper mantle (7.2–7.5 km/sec) is overlying a layer of extremely low velocity (5.0 km/sec) which is interpreted as sialic crust.

Three years ago, a systematic study of crustal structure of the Italian peninsula has been started. Reversed profiles were observed on Sicily, in Calabria, and in Puglia. On Sicily, the structure is very complicated; the crust of the western part looks like a transition between a continental and oceanic structure whereas the eastern side shows a continental-type crust. In Calabria and Puglia, the crustal thickness has been determined to be about 25–35 km.     

Near-vertical seismic reflection image using a novel acquisition technique across the Vrancea Zone and Foscani Basin, south-eastern Carpathians (Romania)

The DACIA PLAN (Danube and Carpathian Integrated Action on Process in the Lithosphere and Neotectonics) deep seismic sounding survey was performed in August–September 2001 in south-eastern Romania, at the same time as the regional deep refraction seismic survey VRANCEA 2001. The main goal of the experiment was to obtain new information on the deep structure of the external Carpathians nappes and the architecture of Tertiary/Quaternary basins developed within and adjacent to the seismically-active Vrancea zone, including the Focsani Basin. The seismic reflection line had a WNW–ESE orientation, running from internal East Carpathians units, across the mountainous south-eastern Carpathians, and the foreland Focsani Basin towards the Danube Delta. There were 131 shot points along the profile, with about 1 km spacing, and data were recorded with stand-alone RefTek-125s (also known as “Texans”), supplied by the University Texas at El Paso and the PASSCAL Institute. The entire line was recorded in three deployments, using about 340 receivers in the first deployment and 640 receivers in each of the other two deployments. The resulting deep seismic reflection stacks, processed to 20 s along the entire profile and to 10 s in the eastern Focsani Basin, are presented here. The regional architecture of the latter, interpreted in the context of abundant independent constraint from exploration seismic and subsurface data, is well imaged. Image quality within and beneath the thrust belt is of much poorer quality. Nevertheless, there is good evidence to suggest that a thick (10 km) sedimentary basin having the structure of a graben and of indeterminate age underlies the westernmost part of the Focsani Basin, in the depth range 10–25 km. Most of the crustal depth seismicity observed in the Vrancea zone (as opposed to the more intense upper mantle seismicity) appears to be associated with this sedimentary basin. The sedimentary successions within this basin and other horizons visible further to the west, beneath the Carpathian nappes, suggest that the geometry of the Neogene and recent uplift observed in the Vrancea zone, likely coupled with contemporaneous rapid subsidence in the foreland, is detached from deeper levels of the crust at about 10 km depth. The Moho lies at a depth of about 40 km along the profile, its poor expression in the reflection stack being strengthened by independent estimates from the refraction data. Given the apparent thickness of the (meta)sedimentary supracrustal units, the crystalline crust beneath this area is quite thin (< 20 km) supporting the hypothesis that there may have been delamination of (lower) continental crust in this area involved in the evolution of the seismic Vrancea zone.     

Structure of the crust and upper mantle beneath the central european rift system

The crustal structure of the central European rift system has been investigated by seismic methods with varying success. Only a few investigations deal with the upper-mantle structure. Beneath the Rhinegraben the Moho is elevated, with a minimum depth of 25 km. Below the flanks it is a first-order discontinuity, while within the graben it is replaced by a transition zone with the strongest velocity gradient at 20–22 km depth. An anomalously high velocity of up to 8.6 km/s seems to exist within the underlying upper mantle at 40–50 km depth. A similar structure is also found beneath the Limagnegraben and the young volcanic zones within the Massif Central of France, but the velocity within the upper mantle at 40–50 km depth seems to be slightly lower. Here, the total crustal thickness reaches only 25 km. The crystalline crust becomes extremely thin beneath the southern Rhônegraben, where the sediments reach a thickness of about 10 km while the Moho is found at 24 km depth. The pronounced crustal thinning does not continue along the entire graben system. North of the Rhinegraben in particular the typical graben structure is interrupted by the Rhenohercynian zone with a “normal” West-European crust of 30 km thickness evident beneath the north-trending Hessische Senke. A single-ended profile again indicates a graben-like crustal structure west of the Leinegraben north of the Rhenohercynian zone. No details are available for the North German Plain where the central European rift system disappears beneath a sedimentary sequence of more than 10 km thickness.     

Crustal structure in the region of the British Isles

Progress in the application of seismic refraction methods to the determination of crustal structure for the British Isles and surrounding sea areas, is reviewed for the period which follows the publications of 1965. The work has been strongly oriented towards the application of “Time Term” interpretation to land-based observations of explosions in the English Channel, to the south of Ireland, in the Irish Sea, in the North Sea, and off the west coast of Scotland.

All of the surveys have included determinations of velocities for P_n and P_g, with some indication of an increase in P_n-velocity with range. In part of the area, indications of lower-crust velocities ranging from 6.9 to 7.3 km/sec have been found.     

Reviewing pre-TRANSALP DSS models

The interpretation of DSS (deep seismic soundings) profiles in Central and Eastern Alps is recalled in the paper and the models of the lower crust and Moho proposed several years ago are compared to the results of the TRANSALP seismic reflection profile. This evaluation highlights a good agreement as far as the geometry of the deep crustal structure is concerned. Therefore, the reliability of the interpretative models, previously exclusively based on DSS profiles, becomes improved. The deep structure beneath the whole Alpine range is examined reconsidering the map of the Moho boundary and the structural model already proposed for the central-eastern sector. Five main interpretative transects are put side by side, starting from the Western Alps and moving eastwards to the Swiss–Lombardian Central Alps (“European Geotraverse”), to the cross section from southern Bavaria to the Euganei Hills, to the TRANSALP profile, and finally to the easternmost profile available so far (southern Bavaria–Trieste). The comparison outlines lateral variations of the deep crustal structure as well as a sharp contrast between the Adria and the European lower crust and Moho. The transition from the Adria plate to the Dinaric domain remains, up to now, undefined.     

Geophysical Study of the Valle Fértil Lineament Between 28°45′S and 31°30′S: Boundary Between the Cuyania and Pampia Terranes

A gravimetric and magnetometric study was carried out in the north-eastern portion of the Cuyania terrane and adjacent Pampia terrane. Gravimetric models permitted to interpret the occurrence of dense materials at the suture zone between the latter terranes. Magnetometric models led to propose the existence of different susceptibilities on either side of the suture. The Curie temperature point depth, representing the lower boundary of the magnetised crust, was found to be located at 25 km, consistent with the lower limit of the brittle crust delineated by seismic data; this unusually thick portion of the crust is thought to release stress producing significant seismicity.

Moho depths determined from seismic studies near western Sierras Pampeanas are significantly greater than those obtained from gravimetric crustal models.

Considering mass and gravity changes originated by the flat-slab Nazca plate along Cuyania and western Pampia terranes, it is possible to reconcile Moho thickness obtained either by seismic or by gravity data. Thus, topography and crustal thickness are controlled not only by erosion and shortening but by upper mantle heterogeneities produced by: (a) the oceanic subducted Nazca plate with “normal slope” also including asthenospheric materials between both continental and oceanic lithospheres; (b) flat-slab subducted Nazca plate (as shown in this work) without significant asthenospheric materials between both lithospheres. These changes influence the relationship between topographic altitudes and crustal thickness in different ways, differing from the simple Airy system relationship and modifying the crustal scale shortening calculation. These changes are significantly enlarged in the study area. Future changes in Nazca Plate slope will produce changes in the isostatic balance.     

Comments on “The crustal structure of the Balearic Sea” — in light of deep-sea drilling in the Mediterranean

Supplementary to the paper by K. Hinz on The crustal structure of the Balearic Sea, some results are reported from the deep-sea drilling cruise of the “Glomar Challenger”.     

Palaeozoic structural development along the Tornquist Zone, Kattegat area, Denmark

The interpretation of newly released commercial 2D reflection seismic data in the Kattegat area, Denmark, has provided us with a better understanding of the Palaeozoic tectonic processes along the Tornquist Fault Zone. A Base Palaeozoic time structure map, a Lower Palaeozoic TWT isopach map, a “true” Lower Palaeozoic TWT isopach map, an Upper Carboniferous/Lower Permian syn-rift TWT isopach map, a Top pre-Zechstein time structure map and a Zechstein combined TWT isopach and Palaeogeography map have been generated. The uniform Lower Palaeozoic sequence thickness in the Kattegat, both inside and outside the Tornquist Zone indicates only minor lateral movements if any, whereas the extensive Upper Silurian sequence, increasing in thickness to the north, indicates a relatively fast regional subsidence. The Base Palaeozoic time structure map and the Late Palaeozoic syn-rift isopach map show a clear Late Palaeozoic extension in the area. The syn-rift isopach map, in combination with the time-equivalent opening of the Skagerrak graben at right angles to the Tornquist Zone in the Kattegat, indicates that this extensional tectonic event had a dextral slip component. Measurements on internal extensional faults in the Tornquist Zone, give a minimum right-lateral displacement of 10.4 km. The footwall blocks were deeply eroded during the Early Permian rifting, and at Zechstein times the area became a peneplane. The Tornquist Zone was later exposed to several tectonic phases, where dextral slip played a role, indicated by the “push up” and “pull down” structures caused by restraining and releasing bends of the Børglum Fault. The dextral displacement along the Børglum Fault since the beginning of the Permian is in the order of 5–7 km based on the displacement of a Lower Palaeozoic local depocentre. Early Permian depocentres and faults, which gives a total amount of right-lateral displacement since the Early Palaeozoic in the order of 15–20 km. The continuously repeated tectonic episodes along the Tornquist Zone throughout most of the Phanerozoic, show that the zone was easily reactivated, implying deep-seated basement faults. The Tornquist Zone can be seen as a “buffer zone”, between continental blocks, whenever changes in the regional stress field are induced.     

Lower lithospheric structure beneath the Trans-European Suture Zone from POLONAISE''97 seismic profiles

The large-scale POLONAISE'97 seismic experiment investigated the velocity structure of the lithosphere in the Trans-European Suture Zone (TESZ) region between the Precambrian East European Craton (EEC) and Palaeozoic Platform (PP). In the area of the Polish Basin, the P-wave velocity is very low (Vp <6.1 km/s) down to depths of 15–20 km, and the consolidated basement (Vp5.7–5.8 km/s) is 5–12 km deep. The thickness of the crust is 30 km beneath the Palaeozoic Platform, 40–45 km beneath the TESZ, and 40–50 km beneath the EEC. The compressional wave velocity of the sub-Moho mantle is >8.25 km/s in the Palaeozoic Platform and 8.1 km/s in the Precambrian Platform. Good quality record sections were obtained to the longest offsets of about 600 km from the shot points, with clear first arrivals and later phases of waves reflected/refracted in the lower lithosphere. Two-dimensional interpretation of the reversed system of travel times constrains a series of reflectors in the depth range of 50–90 km. A seismic reflector appears as a general feature at around 10 km depth below Moho in the area, independent of the actual depth to the Moho and sub-Moho seismic velocity. “Ringing reflections” are explained by relatively small-scale heterogeneities beneath the depth interval from 90 to 110 km. Qualitative interpretation of the observed wave field shows a differentiation of the reflectivity in the lower lithosphere. The seismic reflectivity of the uppermost mantle is stronger beneath the Palaeozoic Platform and TESZ than the East European Platform. The deepest interpreted seismic reflector with zone of high reflectivity may mark a change in upper mantle structure from an upper zone characterised by seismic scatterers of small vertical dimension to a lower zone with vertically larger seismic scatterers, possible caused by inclusions of partial melt.     

Structures in the Canyon Mountain Ophiolite indicate an island-arc intrusion

Lithosphere created in an interarc basin is expected to be characterized by features distinguishing it from “normal” oceanic lithosphere. Apart from island-arc geochemical affinities and from the occurrence of hydrous high-T parageneses in the mantle and deep crustal sequences, it is expected that due to a low rate of spreading, vertical transport prevails over lateral drifting.

The Canyon Mountain complex located in an island-arc environment of Permo-Triassic age offers a remarkable illustration of these expected geological characteristics. In particular, mantle diapirism is deduced from the structural study. Smaller diapirs are formed in crustal formations. The intrusions took place at variable temperatures (1300°–800°C) and were accompanied by multistage melting in hydrous conditions.     

Double-sided onshore–offshore seismic imaging of a plate boundary: “super-gathers” across South Island, New Zealand

Onshore–offshore seismic refraction profiling allows for the determination of crustal and mantle structures in the transition between continental and oceanic environments. Islands and narrow landmasses have the unique geometry of allowing for double-sided onshore–offshore experiments that favor the construction of composite “super-gathers” using the acquisition of onshore–offshore and ocean-bottom seismometer receiver gathers, land explosion shot gathers, and near-vertical incidence multichannel seismic (MCS) profiling. A number of sites at plate boundaries are amenable to the application of double-sided onshore–offshore imaging, including the Indo-Australian/Pacific transform boundary on South Island, New Zealand. By comparing the ratio of island width to mantle refraction (Pn) “maximum” crossover distance, using nondimensional distances, we provide an indicator of raypath “coverage” for crustal illumination. Islands or narrow land masses whose widths are less than twice their maximum crossover distance are candidates for double-sided onshore–offshore experiments. The SIGHT (South Island GeopHysical invesTigation) experiment in New Zealand is located where the width of South Island is sufficiently narrow with respect to its crustal thickness that a double-sided onshore–offshore experiment allows for complete crustal imaging of the associated plate boundary.     

The Earth''s crust in the northwestern part of the Pacific mobile belt

Results of seismic investigations in the transition zone from the Asian Continent to the Pacific Ocean are reported in detail. At the bottom of the sedimentary sequence presumably Cretaceous rocks are found in depressions of the sea floor. The “granitic” layer in the transition zone consists of igneous-sedimentary rocks in different stages of granitization. The “basaltic” layer is developed irregularly in thickness and seismic velocities; its origin is obscure. Apparently the earth's crust in the transition zone is still under formation.     

Earthquakes and soil liquefaction in flood stories of the ancient Near East

Flood stories in the Hebrew Bible and the Koran appear to be derived from earlier flood stories like those in the Gilgamesh Epic and still earlier in the Atrahasis. All would have their source from floods of the Tigris and Euphrates rivers.

The Gilgamesh Epic magnifies the catastrophe by having the flood begin with winds, lightning, and a shattering of the earth, or earthquake. Elsewhere in Gilgamesh, an earthquake can be shown to have produced pits and chasms along with gushing of water. It is commonly observed that earthquake shaking causes water to gush from the ground and leaves pits and open fissures. The process is known as soil liquefaction. Earthquake is also a possible explanation for the verse “all the fountains of the great deep (were) broken up” that began the Flood in Genesis. Traditionally, the “great deep” was the ocean bottom. A more recent translation substitutes “burst” for “broken up” in describing the fountains, suggesting that they erupted at the ground surface and were caused by an earthquake with soil liquefaction. Another relation between soil liquefaction and the Flood is found in the Koran where the Flood starts when “water gushed forth from the oven”. Soil liquefaction observed erupting preferentially into houses during an earthquake provides a logical interpretation if the oven is seen as a tiny house. A case can be made that earthquakes with soil liquefaction are embedded in all of these flood stories.