On Late Cretaceous Tethyan upwelling sediments from the Mishash/Ghareb Formation (Negev, Israel), bulk geochemical and biomarker analyses were performed to explain the high proportion of phosphates in the lower part and of organic matter (OM) preserved in upper parts of the studied section. The profile is composed of three facies types; the underlying Phosphate Member (PM), the Oil Shale Member (OSM) and the overlying Marl Member (MM). Total organic carbon (TOC) contents are highly variable over the whole profile reaching from 0.6% in the MM, to 24.5% in the OSM. Total iron (TFe) varies from 0.1% in the PM to 3.3% in the OSM. Total sulfur (TS) ranges between 0.1% in the MM and 3.4% in the OSM, resulting in a high C/S ratio of 6.5 in the OSM section. A mean proportion of 11.5% total phosphorus (TP) in the PM changed abruptly with the facies to a mean value of only 0.9% in the OSM and the MM. The TOC/TOCOR ratios argue for a high bacterial sulfate reduction activity and in addition, results from fatty acid analyses indicate that the activity of sulfide-oxidizing activity of bacteria was high during deposition of the PM, while decreasing during the deposition of the OSM. The upwelling conditions effected a high primary productivity and consequently the presence of abundant OM. This, in combination with high sulfate availability in the sediments of the PM resulted in a higher sulfide production due to the activity of sulfate-reducing bacteria. Iron availability was a limiting factor during the deposition of the whole section, affecting the incorporation of S into OM. This resulted in the preservation of a substantial part of OM against microbial degradation due to naturally-occurring sulfurization processes expressed by the high C/S ratio of 6.5 in the OSM. Further, the abundant sulfide in the pore water supported the growth of sulfide-oxidizing bacteria promoting the deposition of P, which amounted to as much as 15% in the PM. These conditions changed drastically from the PM to the OSM, resulting in a significant reduction of the apatite precipitation and a high concentration of reactive S species reacting with the OM.

To estimate the potential inventory of natural gas hydrates (NGH) in the Levant Basin, southeastern Mediterranean Sea, we correlated the gas hydrate stability zone (GHSZ), modeled with local thermodynamic parameters, with seismic indicators of gas. A compilation of the oceanographic measurements defines the >1 km deep water temperature and salinity to 13.8 °C and 38.8‰ respectively, predicting the top GHSZ at a water depth of ~1250 m. Assuming sub-seafloor hydrostatic pore-pressure, water-body salinity, and geothermal gradients ranging between 20 to 28.5 °C/km, yields a useful first-order GHSZ approximation. Our model predicts that the entire northwestern half of the Levant seafloor lies within the GHSZ, with a median sub-seafloor thickness of ~150 m. High amplitude seismic reflectivity (HASR), correlates with the active seafloor gas seepage and is distributed across the deep-sea fan of the Nile within the Levant Basin. Trends observed in the distribution of the HASR are suggested to represent: (1) Shallow gas and possibly hydrates within buried channel-lobe systems 25 to 100 mbsf; and (2) a regionally discontinuous bottom simulating reflection (BSR) broadly matching the modeled base of GHSZ. We therefore estimate the potential methane hydrates resources within the Levant Basin at ~100 trillion cubic feet (Tcf) and its carbon content at ~1.5 gigatonnes.