This project seeks to study the water depth anomalies in the Gulf of Aden and its rifted continental Margins. The oceanic crust adjacent to continental rifted margins often shows anomalous water depth (bathymetry) and subsidence history compared to the predictions of ocean lithosphere plate theory. These water depth anomalies are important both for their implications for geodynamic theories of lithosphere temperature structure and also for deep-water oil and gas exploration at rifted continental margins.
A broad region of positive anomalous bathymetry is observed on the young oceanic crust against the Yemen rifted margin in the Gulf of Aden. One proposed hypothesis for explaining this anomalous bathymetry is that the young oceanic lithosphere is underline at depth by relict thicker continental lithosphere with a deep temperature structure different to that of oceanic lithosphere. More mundane explanations for these bathymetric anomalies are that the oceanic crust could simply be thicker than usual or could be underlain by very thin continental crust so giving shallower bathymetry though isostasy. This study therefore uses satellite gravity inversion to map crustal thickness and ocean-continent transition location in the Gulf of Aden and its rifted continental margins. Bathymetric anomalies corrected for sediment loading and lithosphere age is compared with crustal thickness determined from gravity inversion.
1.0 Introduction
The Gulf of Aden within the Red Sea divides the gulf and the horn of Africa. It presents a unique ecosystem and resources that deserves scientific attention. It is characterized by dense, salty water formed by net evaporation with rates up to 1.4 - 2.0 m yr-1 (Hastenrath Lamb 1979). Approximately three million barrels of oil is being transported daily through the Gulf of Aden. In the olden days, and even now, this Gulf used and still provide a considerable amount of sea food for the inhabitants of the surrounding arid lands (Al Saafani 2008). This Gulf is on the other hand important in transport and livelihoods of the people in the region. It serves as a highway for international trade between east and west. Its importance is fundamental to the fact that the present and future generations of peoples dependence on fishing is at stake.
The Gulf of Aden is a small oceanic basin bounded by young conjugate passive margins well preserved beneath a thin postrift sedimentary cover (Leroy et al 2004). It is an important characterized by the processes of rifting, and break-up of continental lithosphere and evolution of young oceanic basin. According to Marcos (1970), the Gulf of Aden is characterized by a deep convection in the northern section that leads to the formation of a deep water mass flowing out into the Gulf of Aden underneath a layer of less saline inflowing water.
The existence of the Arabian monsoon is a phenomenon that dominates the region and affects the general oceanography and meteorology of the region. Heileman and Mistafa (2008) confirm that a northeast monsoon winds extend well into the Gulf of Aden and the southern Red Sea during winters, causing a seasonal reversal in the winds over this entire region. The prevalent seasonal monsoon reversal and the local coastal configuration combined during summer season forces a radically different circulation pattern composed of a thin surface outflow and an intermediate inflowing layer of Gulf of Aden thermocline water and a drastically reduced outflowing deep layer (Patzert 1974). The general surface circulation within the basin is cyclonic.
2.0 Geological Background of the area
The Gulf of Aden comprises the southern limit of the Arabian plate that started to drift away from Africa around 17.6 million years ago. It is a young ocean basin formed by the rifting of Asia from Africa. Its opening is usually considered as the westward propagation of a lithosphere crack from the Carlsberg ridge to the Afar volcanic area, which now represents the centre of a broad anomalous region with low S-wave mantle velocity, high elevation and much more evidence of mantle plume activity Lucazeau et al (2008). It has a well-defined continental margin, small oceanic basin, and an active mid-ocean ridge (Sheba) in the center characterized by a rift valley and fracture zone. The geophysical survey of the basin, reveals continental, oceanic domains and ocean continent transition (OCT) domains with distinct morphological and sedimentological characteristics (Leroy et al 2004)
According to Al Saafani (2008), the Gulf of Aden is approximately 900 km in length and varies in width from 26 km at Bab el Mandab to about 320 km at Ras Asir. The middle region is the deepest at approximately between 20002500m and the shallow area is less than 1000 m. It has an average depth of 1800m, which increases from west to east. The Gulf of Aden is therefore a modern analogue for the early stages of mature margins (Lucazeau et al 2008), such as the Atlantic margins. The conjugate margins of the Encens Sheba survey area were formed during the last period of rifting. The thickness of the oceanic crust of the Gulf of Aden varies from 4.8 to 8.4 km, as interpreted from a seismic refraction survey (Cochran, 1982). In these margins, the actual thermal regimes can be analyzed directly by means of surface heat flow, seismic tomography or regional isostasy, and compared with the observed tectonic and magmatic styles.
Through the ArabiaSomalia Plate, successive Jurassic and Cretaceous rifting episodes have created major basins (Leroy et al 2004 Cochran, 1982). These basins have a predominant eastwest to northwestsoutheast orientation. Seafloor spreading is more recent in the western part of the Gulf of Aden than in its eastern part according to previous magnetic anomaly studies, and that a magnetic quiet zone corresponds to an area of thinned crust (Leroy et al 2004). The structures and evolution of the Gulf of Aden margins are related to successive development of the formation of the Gulf of Aden during the Oligo-Miocene and the continental margin of the Indian Ocean during the Early Cretaceous up to the Paleocene (Leroy et al 2004 Besse and Courtillot, 1988).
According to dAcremont et al (2004), during the Oligocene times, the Afro Arabian Plate began to separate due to the creation of two divergent basins, which both evolved in the development of two narrow oceanic basins that comprised of the Gulf of Aden to the south, between Arabia and Somalia and the Red Sea to the west, between Arabia and Africa (Nubia). The East African rift forms the third branch of the Afar ridgeridgeridge type (RRR) triple junction (dAcremont et al 2004 Wolfenden et al. 2004) that is still in the rifting stage. This later Oligo-Miocene stretching episode of the ArabiaSomalia Plate reactivated inherited structures of these Mesozoic basins (dAcremont et al 2004 Ellis et al.1996 Granath 2001 Bellahsen 2002).
3.0 Water depth anomaly in the Gulf of Aden and why it is shallower than it should be.
The water depth anomaly in the Gulf of Aden is shallower than it should be because of eddies, outflows and evaporation due to high and recent rise in sea surface temperatures (SST) during summers seasons. According to Bower et al (2005) the Saline, dense Red Sea Water (RSW) originates in the northern Red Sea because of an excess of evaporation over precipitation. It enters the Gulf of Aden (GOA) in the northwestern Indian Ocean as a dense overflow through the shallow Bab el Mandeb. The Red Sea Outflow Water (RSOW) is the entrained mixed product water that descends from the INCLUDEPICTURE httpjournals.ametsoc.orgna101homeliteratumpublisheramsjournalsentities223C.gif MERGEFORMATINET 150 m deep Hanish Sill in the northern BAM Strait, less dense fresher water. However, during winter, the Red Sea outflow transport is typically about two times the annual mean because of monsoon winds and seasonal fluctuations in buoyancy forcing. Outflow transport reaches a maximum ( INCLUDEPICTURE httpjournals.ametsoc.orgna101homeliteratumpublisheramsjournalsentities223C.gif MERGEFORMATINET 0.6 Sv) in winter (OctoberMay), when prevailing monsoon winds over the region are from the south-southeast (Bower et al 2002).
There have been recently numerous and new oceanographic observations in the Gulf of Aden along the northwestern Indian Ocean. Some studies have shown large, energetic, deep-reaching mesoscale eddies (Bower et al 2002) that contribute and influence the spreading rates and pathways of intermediate-depth Red Sea Water (RSW). The other important pieces of the thermal puzzle are the present-day topography, gravity, seismic velocities and all the geological aspects contribute to the shallowness. Comparison of salinity and direct velocity measurements (Bower et al 2002) indicates that the eddies advect and stir RSW through the Gulf of Aden and the anomalous water properties in the center of the anticyclonic eddy point to a possible formation site in the Somali Current System.
The exchange flow has a two-layer structure, with dense, saline RSW flowing out at the bottom, and less dense GOA surface water flowing toward the Red Sea in the upper layer. During summer (JuneSeptember), prevailing winds from the north-northwest drive a surface flow out of the Red Sea, and GOA water flows in via an intermediate layer sandwiched between the surface layer and a thin layer of outflowing dense RSW, producing a three-layer exchange flow. The overflow results from an excess of evaporation over precipitation in the entire Red Sea is estimated to be about 2 m yr_1.
The existence of mesoscale eddies and currents contribute to the water depth. Low oxygen water perists throughout a wide depth range in the intermediate waters of the Arabian Sea. The extensive exchange of water between the Red Sea, the Gulf of Aden and the Arabian Sea, the strong evaporation and the monsoonal winds that blow over the region, all assist in the formation of complex vertical structures in the water column of the Gulf of Aden (Al Saafani 2008). As argued by Yamanaka et al (2008), the resultant outflow from the Red Sea through the Gulf of Aden is considered to play a strong role in determining the properties of these intermediate waters
Several studies and research have been done to ascertain the reasons for shallowness in relation to water depth anomalies in this Gulf. Bower et al (2005) conducted a study dubbed the Red Sea Outflow Experiment (REDSOX), which was the first comprehensive field study of the hydrography and circulation of both the descending Red Sea outflow plume, and the equilibrated RSOW in the GOA. The study aimed at achieving the following
describing the pathways and downstream evolution of the descending outflow plumes in the western Gulf of Aden, quantifying the processes that control the final depth of the equilibrated RSOW, and
Identifying the transport processes and mechanisms that advect RSOW and its properties through the GOA and into the Indian Ocean.
This study have revealed new large, energetic, deep-reaching eddies in the Gulf of Aden that fundamentally influence the spreading of RSW. With the use of Shipboat Acoustic Doppler Current Profiler (SADCP), measurements at 100 and 300 m reveal that currents in the Gulf of Aden were strong and organized into coherent eddy structures in February 2001 (Bower et al 2002). The most conspicuous feature is an energetic cyclonic eddy adjacent to the Somali coast in the southwestern part of the gulf.
From these studies we can deduce that the outflow waters reached neutral buoyancy where the bathymetric channels empty into the deep Tadjura Rift during winters and the northern channel was the source for the most saline, deepest and densest salinity maximum. Outflow currents are lower during summer due to evaporation. The outflow water reached neutral buoyancy somewhat upstream of the channel exits, with lower salinity, temperature, and potential density in comparison with winter (Bower et al 2002). It is also evident that the two deeper salinity maxima at nominal depths of 600 and 800 m, are initially more confined by the walls of the Tadjura Rift. Waters associated with both maxima are generally swept southward out of the rift through gaps in the southern rift wall and also along the continental slope.,
Other theories suggest that magma accumulation, heat and changing atmospheric pressure. Lucazeau et al (2008) alleges that accumulation of magma below the OCT crust could produce the expected thermal and density anomaly by advection of hotter and lighter material and release of latent heat. He further explains that another alternative is serpentinization of the continental mantle, which might be exposed and fractured in the OCT and then percolated by water flows, as observed at several locations on mid-ocean ridges. Al Saafani (2008) noted that the atmospheric pressure at mean sea-level is highest during January and lowest during July at Aden with a range of 10 m bar.
How to get the data to produce gravity inversion maps
There are several ways of getting data and deriving gravity inversion maps. The methods are presented here.
Gravity
An appropriate method for producing gravity inversion maps is by modeling geological zones having anomalous density. It basically involves determining the top of the anomalous zone from non-potential fields data. For this case, for example, the westward continuation into the Gulf of the axial trough and linear magnetic anomalies of the Sheba Ridge is often used. A potential fields data is then used to derive the lower boundary of the geologic anomalous zones.
Figure 1 Gravity map
A lower boundary to an anomalous zone is formulated by predicting parameters representing the lower boundary within predetermined limits. Gravity and bathymetry profiles across the Carlsberg Ridge in the Indian Ocean are analyzed to investigate the isostatic compensation of this part of the Indian Ocean spreading ridge system. This is done by using a gravity inversion process on the potential fields data, such as measurements of gravity data andor magnetic data. These may be in both vector and tensor form. The potential fields data is compared to the predicted fields from the results of the inversion process to obtain a difference between the two. If the difference exceeds a predetermined value, the parameters representing the anomalous zone are adjusted to improve the fit. When the lower boundary limits are reached or the difference between the model and the data is less than the predetermined value or convergence is attained, the anomalous zone has been determined.
Bathymetry
The Gulf of Aden comprise two bathymetric channels that originate seaward of Perim Narrows and direct the dense RSW from the strait to the open ocean. One can use 3-D seismic array to map crustal stretching across and along the margins. User receiver functions can be used to retrieve crustal thickness and crustal composition along the margins. Another option is utilizing the ocean wave-shoaling photographic imagery. The instrument used for the acquisition of sea-surface image sequences is a ground-based nautical X-band radar with horizontal polarization. The device utilized during the experiments is a softwarehardware combination consisting of a commercial, navigational Furuno X-band radar antenna and radar device. In the case of datum determination, a tidal-gauge-measurement is used.
SHAPE MERGEFORMAT
Figure 2 Bathymetry map
Other cross-spectral techniques are employed to obtain an admittance function for isostatic studies. Synthetic topography and gravity profiles are computed for a cooling plate model are subtracted from the observed, to estimate the median valley signatures. The residual gravity anomaly over the ridge axis is mostly explained by the median valley topography with a uniform crustal thickness. Spatial profiles are developed using contours, seismic data or satellite imagery to create 3D maps. Then the profiles are created with a known parameter, a prerequisite for checking depth and depth inversions. Depth inversion occurs when an observation has a shallower depth than the observation directly preceding it.
The sonar instrument uses a transducer that is usually mounted on the bottom of a ship. The sound pulses are sent from the transducer straight down into the water where the sound reflects off the seafloor and returns to the transducer. It is claimed that acoustic penetration into the sea bed increases with decreasing frequency. The distance to the seafloor is calculated based on the time the sound takes to travel to the bottom and back to the surface. Water depth is estimated by using the speed of sound through the water. The sound pulses are sent out regularly as the ship of opportunity moves along the surface, which produces a line showing the depth of the ocean beneath the ship. This continuous depth data is used to create bathymetry maps of the survey area.
4.3 Age
The age is mapped by shipborne magnetometers which allow delineation of zones of normal and reversed magnetic polarity. The magnetic zones form distinct stripes on the map as the oceanic plates grow. Heat flow measurements are another method of age determination and inversion map production. The average value of heat flow measurements calculated from measurements and other collected data indicates that the age of basins and margins.
SHAPE MERGEFORMAT
Figure 3 A map showing age
The thermal gravity anomaly can be conditioned using ocean isochrons from plate reconstruction models to provide the age and location of oceanic lithosphere (Chappell and Kusznir 2008). There is a relationship of water depth with age. The agreement in age from both heat flow and water depth data favors aids in age determination. The map production using all these characteristics helps in determining geologic events with respect to their dates.
4.4 Crust thickness
According to Lucazeau et al (2008), the high resolution 3-D forward modeling approaches reveal a possible crustal thickness and density distribution beneath. The use of satellite remotely sensed imagery give precise information about the sediments and sea bed. The seabed, basement and mantle boundaries are defined by a series of triangular facets, whose size varies as the amount of constraining data changes. The sediment and basement boundaries and the base of the crust are defined by larger facets than those defining the bathymetry.
SHAPE MERGEFORMAT
Figure 4 Crust thickness map
Another method is the use of a mechanical bed level detection in combination with DGPS. The bed level soundings are often performed by use of a vehicle moving through the surf zone. Alternatively the water depth can be measure by using the single-beam echo sounder. The sonar instrument is used to measure the crust thickness. The Instrument measures the vertical distribution of the Turbidity turbidity levels in the water column, transition from water column to bed based on the scattering of light from the suspended particles and the bed material particles and transition from water column to air.
5.0 How the data is collected for Satellite gravity inversion.
Analysis of satellite gravity data offers an opportunity to rapidly evaluate the sedimentary structure of basins. The accurate and evenly distributed measurements of the gravity field determined from satellite orbits contain information about the bathymetry, age, sediment thickness and crustal structure of the worlds oceans. The data is collected from satellite based methods.
Many organizations have carried out seismic surveys, opportunity vessels and satellite sensors to collect data for gravity inversion maps. These institutions have archived them and available to the general public. For instance, the sea level data from satellite based altimeters is used and obtained from, for example, Archiving, Validation, and Interpretation of Satellite Oceanography Data (AVISO) operations canter. AVISO distributes sea surface heights (SSH) and sea level anomalies (SLA) measured by Jason and ERS12 satellites (Al Saafani 2008) Jason and ERS12 altimeter data helped in better resolving the mesoscale variability and the data provides more homogeneous and reduced mapping errors than the individual data set. The National Oceanic and Atmospheric Administration (NOAA) as a federal US agency focused on the condition of the oceans have a lot of satellite data.
The measurements made by the sensors mounted on board satellites are used, especially, Seas Surface Temperatures (SST) and the sea level heights are measured by the satellite based altimeters and the sea winds measured by the scatterometer. Other profiling instrument such as the ConductivityTemperatureDepth (CTD), Expendable Bathythermograph (XBT) and Mechanical Bathythermograph (MBT) and the Ocean Station Data (OSD), obtained using reversing bottles attached with reversing thermometers are used (Al Saafani 2008).
The use of hyperspectral sensing as one of the Remote sensing remote sensing techniques can be used to generate maps of the seafloor. A model technique is the Hyperspectral Mapping (HyMap) that is widely used in bathymetry and sea floor mapping. This process involves transformation of subsurface reflectance to the bottom albedo. According to Heege et al (2003), the unknown input value of depth is calculated iteratively in combination with the spectral unmixing of the respective bottom reflectance. The unmixing procedure produces the sea floor coverage of three main bottom components and the residual error between the model bottom reflectance and the calculated reflectance. The final depth, bottom reflectance and bottom coverage is achieved at the minimum value of the residual error. The final step of the thematic processing classifies the bottom reflectance due to the spectral signature of different bottom types and species using a Fuzzy Logic method and assignment of individual probability functions for each defined sea floor component (Heege et al 2003).
The gravity inversion maps
Thermal Gravity anomaly
Thermal gravity anomaly is generated by the gravitational admittance that modifies the topography. The thermal gravity anomaly is attributed the early stages in the formation of divergent margins when the lithosphere experiences large changes in temperature which play a key role in the anomaly of the Gulf of Aden. Evidence suggest that in this regions, thermal anomaly in the upper mantle that has persisted after continental break-up. According to Chappell and Kusznir (2008), the oceanic lithosphere thermal model is always used to predict the lithospheres thermal gravity anomaly. In the deeper parts of the margin the heat flow is high and constant, but it decreases abruptly near the shelf-slope.
Thermal gravity anomaly parameters contain information on the state of isostacy for a surface topography feature. The sensitivity of the lithosphere thermal gravity anomaly and the predicted ocean depth from gravity inversion maps are always correlated. Variations of the sea floor bathymetry constitute a load distribution on the oceanic lithosphere. The presence of shallow-water sediments deposited after the opening of the Atlantic Ocean hints at lower subsidence than would have occurred in the absence of persistent thermalanomalies.
Figure 5 Thermal gravity anomaly maps (mgal) a) with sed and b) with sed and volcanic corrections.
As shown in figure , the curve along the gulf is more shallower than towards the ocean, hence reduction in water discharge, the latent heat due to high temperatures and eddies contribute to the water depth anomaly in the shores of Yemen.
Thinning Factor
The lateral density changes caused by the elevated geotherm in thinned continental margins such as the Gulf of Aden and adjacent ocean basin lithosphere yield important thermal gravity anomaly. According to Leroy et al (2004), magnetic quiet zone corresponds to an area of thinned crust. Lucazeau et al (2008) argues that the lithosphere in the deep margin should be locally hotter and more buoyant than any homogeneous margin. There are several methods to determine crustal depth, lithosphere thinning and the location of the ocean-continent transition at rifted continental margins using 3-D gravity inversion which includes a correction for the large negative lithosphere thermal gravity anomaly within continental margin lithosphere.
Figure 6 Thinning Factor maps a) with sed correction and b) with sed and volcanic corrections
Sediment Thickness
The sediment thickness can be obtained from geophysical studies and associated gravity inversion maps. For instance, a three-dimensional Bouguer anomaly map is produced for dimensional inverse approach to gravity data interpretation. The Bouguer anomalies along the axial portion of the rift floor, as deduced from the results of the regional and residual separation, are mainly caused by deep-seated structures.
Figure 7 Sediment Thickness map
The use of satellite remotely sensed imagery and hyperspectral remote sensing gives precise information about the sediments and sea bed. The seabed, basement and crust boundaries are defined by a series of triangular facets with their sizes varying as the amount of constraining data changes. In this case, the sediment and basement boundaries and the base of the crust are defined by larger facets than those defining the bathymetry.
Figure 8 Volcanic addition map (m)
Mantle Residual Anomaly
There are prominent slow anomalies within the Gulf of Aden and the entire Indian Ocean unlike the fast the central Atlantic and the older parts of the Pacific. Since most of the large slow anomalies define geoid highs, Leroy et al (2004) indicates that there is poor overall correlation between velocities in the upper mantle and the geoid because subduction zones in general are associated with geoid highs and regions of fast velocity below.
Figure 9 Mantle Residual Anomaly (mgal) a) with sed correction and b) with sed and volcanic corrections
These maps identify target area of study using satellite gravity inversion. The maps are compiled from seismic data, surface and wave data and contoured to produce the exact heights of the crust at various points. Embedded thermal correction and prediction of Gulf of Aden crust thickness is used to map crustal thickness. The map shows thickening of the crust from Red Sea and Gulf of Aden eastwards and northwards. This crust is inherent in the middle of the Arabian plate.
7.0 Summary
The water depth anomaly has been analyzed to describe the vertical and horizontal structure of Red Sea. The outflow water in the western Gulf of Aden, at the location where it is first injected into the open ocean is one contributor to shallowness despite the high atmospheric temperatures and eddies among the contributing factors. An important region for oil exploration and shipping, the application of 1km resolution imagery aids in the studies and development of this region.
This project is focused on how to create gravity inversion maps and the general application of satellite data in bathymetry, SST estimation, thinning and age determination. It is evident that the Gulf of Aden existing in the slow Indian Ocean margin is experiencing water depth anomalies especially during summers. Remote Sensing and Earth observation does not only cover the global and regional survey of geophysical parameters by satellite- or airborne radars, it also includes local observations by ground based radar techniques.