Authigenic carbonate mineral formation in the Pagassitikos palaeolake during the latest Pleistocene, central Greece

The Pagassitikos Gulf in Greece is a semi-enclosed bay with a maximum depth of 102 m. According to the present-day bathymetric configuration and the sea level during the latest Pleistocene, the gulf would have been isolated from the open sea, forming a palaeolake since ~32 cal. ka b.p. Sediment core B-4 was recovered from the deepest sector of the gulf and revealed evidence of a totally different depositional environment in the lowest part of the core: this contained light grey-coloured sediments, contrasting strongly with overlying olive grey muds. Multi-proxy analyses showed the predominance of carbonate minerals (aragonite, dolomite and calcite) and gypsum in the lowest part of the core. Carbonate mineral deposition can be attributed to autochthonous precipitation that took place in a saline palaeolake with high evaporation rates during the last glacial–early deglacial period; the lowest core sample to be AMS ¹⁴C dated provided an age of 19.53 cal. ka b.p. The palaeolake was presumably reconnected to the open sea at ~13.2 cal. ka b.p. during the last sea-level rise, marking the commencement of marine sedimentation characterised by the predominance of terrigenous aluminosilicates and fairly constant depositional conditions lasting up to the present day.     

Evidence of NW extension of the North Anatolian Fault Zone in the Marmara Sea: a new interpretation of the Marmara Sea (İzmit) earthquake on 17 August 1999

Active faults aligning in NW–SE direction and forming flower structures of strike-slip faults were observed in shallow seismic data from the shelf offshore of Avc?lar in the northern Marmara Sea. By following the parallel drainage pattern and scarps, these faults were traced as NW–SE-directed lineaments in the morphology of the northern onshore sector of the Marmara Sea (eastern Thrace Peninsula). Right-lateral displacements in two watersheds of drainage and on the coast of the Marmara Sea and Black Sea are associated with these lineaments. This right-lateral displacement along the course of these faults suggests a new, active strike-slip fault zone located at the NW extension of the northern boundary fault of the Ç?narc?k Basin in the Marmara Sea. This new fault zone is interpreted as the NW extension of the northern branch of the North Anatolian Fault Zone (NAFZ), extending from the Ç?narc?k Basin of the Marmara Sea to the Black Sea coast of the Thrace Peninsula, and passing through B üy ük çekmece and K ü ç ük çekmece lagoons. These data suggest that the rupture of the 17 August 1999 earthquake in the NAFZ may have extended through Avc?lar. Indeed, Avc?lar and ?zmit, both located on the Marmara Sea coast along the rupture route, were strongly struck by the earthquake whereas the settlements between Avc?lar and ?zmit were much less affected. Therefore, this interpretation can explain the extraordinary damage in Avc?lar, based on the newly discovered rupture of the NAFZ in the Marmara Sea. However, this suggestion needs to be confirmed by further seismological studies.     

Sensitive determination of enzymatically labile dissolved organic phosphorus and its vertical profiles in the oligotrophic western North Pacific and East China Sea

Trace concentrations of labile dissolved organic phosphorus (LDOP) in oligotrophic seawater were measured by use of an enzymatic procedure and a nanomolar phosphate analytical system consisting of a gas-segmented continuous flow analyzer with a liquid waveguide capillary cell. LDOP, defined as DOP hydrolyzed by alkaline phosphatase (AP) from Escherichia coli, was quantified as the difference between the phosphate concentrations of the seawater sample with and without AP treatment. For sensitive measurement of LDOP, we found that phosphate contamination derived from commercially available AP must be corrected, and azide treatment before AP treatment proved effective in removing biological effect that occurs during DOP hydrolysis. Field observations at six stations of the western North Pacific and the East China Sea during the boreal summer revealed that, in the upper 200 m of the water column, LDOP concentrations ranged from the detection limit (3 nM) to 243 nM, and phosphate concentrations ranged from 5 to 374 nM. The spatial distribution patterns of LDOP were similar to those of phosphate. Most of the depth profiles for LDOP and phosphate showed concentrations were extremely low, <25 nM, between the surface and the deep chlorophyll maximum layer (DCML) and increased below the DCML. Strongly depleted LDOP and phosphate above the DCML suggest that LDOP is actively hydrolyzed under phosphate-depleted conditions and utilized by microbes.     

Sessile animals on tar globules in the waters around the Ryukyu Islands

A large number of tar globules with sessile organisms were collected from the surface tows taken with larval nets in the waters around the Ryukyu Islands during November and December, 1973. Bryozoans (one species), tubeworms (Serpulidae,Janua (Dexiospira) foraminosa (Moore &Bush)) and goose barnacles (Lepas pectinata Spengler,L. anatifera Linné) were the most important species of sessile animals found on these tar globules. Sinking of tar globules byLepas of middle or large size was suggested from the differences in their specific gravity.     

Discussion: a critique of Possible waterways between the Marmara Sea and the Black Sea in the late Quaternary: evidence from ostracod and foraminifer assemblages in lakes ?znik and Sapanca, Turkey, Geo-Marine Letters, 2011

The identification of past connection routes between the Black Sea and the Sea of Marmara, other than the traditional one through to the Bosphorus Strait, would be of considerable interest to the international scientific community. Nazik et al. (Geo-Mar Lett 31:75?C86 (2011) doi:10.1007/s00367-010-0216-9) suggest the possibility of two alternative waterway connections via lakes Sapanca and ?znik. Their Black Sea to Sea of Marmara multi-connection hypothesis, which is based on undated marine fossils collected in both lakes from surficial grab samples, conflicts with many earlier studies. In this contribution, the hypothesis and the underlying data are discussed in the light of previous tectonic, sedimentological and limnological findings showing that it is impossible to have had marine connections through lakes Sapanca and ?znik during the last 11.5?ka. Global sea-level trends and tectonic uplift rates would accommodate a connection between the Sea of Marmara and Lake ?znik in the middle Pleistocene. Uplift rates for the northern block of the North Anatolian Fault, when compared with the global sea-level curve, clearly indicate that there cannot have been a connection through the ?zmit Gulf?CLake Sapanca?CSakarya Valley for at least the past 500?ka. Moreover, borehole sediments along the western shores of Lake Sapanca, which reach down to the bedrock, do not contain any marine fossils.     

The volcanic debris avalanche on the SE submarine slope of Nisyros volcano, Greece: geophysical exploration and implications for subaerial eruption history

A spectacular hummocky topography was discovered offshore of the south-eastern slope of the Nisyros island volcano in the eastern sector of the Aegean volcanic arc in 2000–2001, using multibeam bathymetric mapping and seismic profiling, and interpreted as part of a volcanic debris avalanche originating onland. During E/V Nautilus cruise NA011 in 2010, a detailed side-scan sonar and ROV exploration aimed at evaluating the surface morphology of this avalanche field. Combining the new data with selected older datasets reveals that the debris avalanche is characterized by numerous (at least 78) variously sized and shaped hummocks. Some of these are distinctly round, either scattered or aligned in groups, whereas others are elongated in the form of ridges. This is consistent with existing models accounting for variations in the longitudinal and lateral velocity ratio of landslides. Maximum dimensions reach 60 m in height above the sea bottom, 220 m in length and 230 m in width. The structures outline a large tongue-shaped, submarine hummock field of about 22.2 km², approx. 4.8 km wide and 4.6 km long and with an estimated volume of 0.277 km³. Due to its characteristic shape, the collapsed volcanic flank is interpreted to represent a singular failing event, involving a rapid and virtually instantaneous downslope movement of the slide mass into the sea. Indeed, the H/L (height of 280 m vs. run-out of 7 km) ratio for the Nisyros slide is 0.04; plotted against volume, this falls within the theoretical bounds as well as measured values typical of submarine landslides. The timing of the event is probably related to the extrusion of Nikia lavas and their subsequent failure and formation of a main scarp observed at about 120 m depth on an 8-km-long seismic profile and a map of slope angle distribution, at the depth where the palaeo-coastline was located 40 ka ago. An inferred age of ca. 40 ka for the avalanche awaits confirmation based on dating of core material.     

How does multiscale modelling and inclusion of realistic palaeobathymetry affect numerical simulation of the Storegga Slide tsunami?

The ∼8.15 ka Storegga submarine slide was a large (∼3000 km³), tsunamigenic slide off the coast of Norway. The resulting tsunami had run-up heights of around 10–20 m on the Norwegian coast, over 12 m in Shetland, 3–6 m on the Scottish mainland coast and reached as far as Greenland. Accurate numerical simulations of Storegga require high spatial resolution near the coasts, particularly near tsunami run-up observations, and also in the slide region. However, as the computational domain must span the whole of the Norwegian-Greenland sea, employing uniformly high spatial resolution is computationally prohibitive. To overcome this problem, we present a multiscale numerical model of the Storegga slide-generated tsunami where spatial resolution varies from 500 m to 50 km across the entire Norwegian-Greenland sea domain to optimally resolve the slide region, important coastlines and bathymetric changes. We compare results from our multiscale model to previous results using constant-resolution models and show that accounting for changes in bathymetry since 8.15 ka, neglected in previous numerical studies of the Storegga slide-tsunami, improves the agreement between the model and inferred run-up heights in specific locations, especially in the Shetlands, where maximum run-up height increased from 8 m (modern bathymetry) to 13 m (palaeobathymetry). By tracking the Storegga tsunami as far south as the southern North sea, we also found that wave heights were high enough to inundate Doggerland, an island in the southern North Sea prior to sea level rise over the last 8 ka.     

Revisiting the February 6th 1783 Scilla (Calabria, Italy) landslide and tsunami by numerical simulation

On February 6th, 1783, a landslide of about 5 × 10⁶ m³ triggered by a 5.8 M earthquake occurred near the village of Scilla (Southern Calabria, Italy). The rock mass fell into the sea as a rock avalanche, producing a tsunami with a run-up as high as 16 m. The tsunami killed about 1,500 people, making it one of the most catastrophic tsunamis in Italian history. A combined landslide-tsunami simulation is proposed in this paper. It is based on an already performed reconstruction of the landslide, derived from subaerial and submarine investigation by means of geomorphological, geological and geomechanical surveys. The DAN3D model is used to simulate the landslide propagation both in the subaerial and in the submerged parts of the slope, while a simple linear shallow water model is applied for both tsunami generation and propagation. A satisfying back-analysis of the landslide propagation has been achieved in terms of run-out, areal distribution and thickness of the final deposit. Moreover, landslide velocities comparable to similar events reported in the literature are achieved. Output data from numerical simulation of the landslide are used as input parameters for tsunami modelling. It is worth noting that locations affected by recordable waves according to the simulation correspond to those ones recorded by historical documents. With regard to run-up heights a good agreement is achieved at some locations (Messina, Catona, Punta del Faro) between computed and real values, while in other places modelled heights are overestimated. The discrepancies, which were most significant at locations characterized by a very low slope gradient in the vicinity of the landslide, were probably caused by effects such as wave breaking, for which the adopted tsunami model does not account, as well as by uncertainties in the historical data.     

Factors controlling the morphological evolution of the Çanakkale Strait (Dardanelles, Turkey)

Seismic profiling, bathymetric and physical oceanographic data collected from the Çanakkale Strait revealed that the morphological evolution of the strait has been controlled by tectonic activity, and sediment erosion and deposition. Sediments in the strait have been sourced mostly by rivers draining the Biga Peninsula during lowstand periods. In highstand periods, by contrast, deposits in the strait were reworked by currents. The seafloor morphology of the Çanakkale Strait is also controlled by a sequence of factors ranging from tectonics to current erosion and deposition. Channel deposits overlying the basement are being eroded at the narrower, meandering central section of the strait (the Nara Passage) due to high current velocities. The eroded sediments are deposited in the relatively linear and wider, northern and southern sectors of the strait exposed to low current velocities. As a result, the high-energy areas are more deeply incised due to the erosion, whereas deposition elevates the seafloor in the areas exposed to lower current energy. Three strike-slip faults, which possibly relate to the activity of the North Anatolian Fault Zone, are responsible for the irregular shape of the strait and this, in turn, controls the current velocity along the strait. The high-energy conditions probably commenced with the latest invasion of Mediterranean waters some 12 ka b.p., and have continued as a two-layered current system to the present day.     

Contributions to the dynamics of the Antarctic Circumpolar Current (A. C. C.) II. Integral relationship

Non-dimensional equations of motion are derived for the A.C.C. of the barotropic mode, including the bottom friction and the horizontal eddy viscosity. Integration of the vorticity equation along a streamline leads to the zeroth order stream function which is dependent only on depth divided by Coriolis parameter. Integration of the momentum equation along a streamline yields the relation between the momentum input by wind stress and its dissipation by the bottom friction and by the horizontal eddy viscosity. This relation determines the magnitude of the stream function. It explains differences in the total transport of the A.C.C. obtained byBryan andCox (1972), though it gives only one third of the total transport obtained byKamenkovich (1972) with his vertical eddy viscosity of 10²cm² s^?1. With 1 cm² s^?1 of this viscosity,Bryan andCox obtained the transport of about 650 or less than 32×10⁶m³s^?1 for constant or variable depth models, respectively. The higher transport is mainly due to broadening of the width of the A.C.C., whereas the lower value is due to its narrowing and meandering which in turn make the horizontal eddy viscosity more effective (by exercising friction on both sides of the A.C.C.) and the wind stress input smaller than the almost zonal streamlines for constant depth. In the Appendix dynamics of the bottom boundary layer is treated to give rational estimates of the bottom stress in terms of the geostrophic flow and is compared to the recent observations of the benthic boundary current in the Straits of Florida and off San Diego.     

Evidence for active faults in Küçükçekmece Lagoon (Marmara Sea,Turkey), inferred from high-resolution seismic data

A total of 42 km of high-resolution seismic reflection and bathymetric data were collected for the first time to document stratigraphic and structural features of the uppermost 5 m of the Holocene sedimentary infill of Küçükçekmece Lagoon along the Marmara Sea coast of Turkey. The lagoon gradually deepens from 1 m off the northern coast to a maximum of 20 m in the southern basin. Stratigraphically, the uppermost seismic unit is characterized by a generally parallel reflection configuration, indicating deposition under low-energy conditions. In the southern basin of the lagoon, the sub-bottom is locally characterized by frequency attenuated and chaotic reflections interpreted as gas-charged sediments. Structurally, the soft sediment of the first 5 m below the lagoon floor is locally deformed by active strike-slip fault zones, here named FZ1, FZ2, and FZ3. These fault zones are NW–SE oriented and follow the long axis of the lagoon, compatible with the geographic alignment of the lagoon, the onland drainage pattern, and the scarps of the surrounding terrain. Moreover, the fault zones in Küçükçekmece Lagoon are well correlated with active offshore faults mapped during previous studies. This suggests that the FZ1, FZ2, and FZ3 fault zones are not merely local fault systems deforming the Küçükçekmece Lagoon bottom, but that they may be part of a regional fault zone extending both north and southward to merge with the northern branch of the North Anatolian Fault Zone (NAFZ) in the Ç?narc?k Basin. This, however, needs to be confirmed by further structural and seismological studies around Küçükçekmece Lagoon in order to more firmly establish its link with the NAFZ in the Marmara Sea, and to highlight potential seismic risks for the densely populated Istanbul metropolitan area.     

Morpho-tectonic evolution of the Marmara Sea inferred from multi-beam bathymetric and seismic data

In an initial stage, the Sea of Marmara developed as a graben and, in due course, considerable volumes of sediments were deposited in this basin. Before 200 ka, a new fault (New Marmara Fault) cutting through the whole basin developed, which postdated large sub-marine land sliding in the western part of the basin. This mass movement created the Western Ridge. The initiation of this strike-slip fault indicates that the extensional stress regime was replaced by a new, shearing stress field. In the eastern part of the Marmara Basin, the New Marmara Fault consists of two branches. The northern one replaces the normal faulting at the bottom of the northeastern slope of the basin. As a result, this slope has been rejuvenated. The southern branch is located along the central axis of the basin, forming the major extension of the North Anatolian Fault Zone within the region. Two restraining bends were formed because of the counterclockwise rotation of that part of the Anatolian Block. This resulted the uplifting of the Eastern Ridge and the formation of the positive flower structure within the Tekirdag Basin. The establishment of the compressional regime around the Sea of Marmara also resulted in the northwest–southeast shortening of the initial Marmara Basin.     

Postglacial floodings of the Marmara Sea: molluscs and sediments tell the story

The early Holocene marine flooding of the Black Sea has been the subject of intense scientific debate since the “Noah’s Flood” hypothesis was proposed in the late 1990s. The chronology of the flooding is not straightforward because the connection between the Black Sea and the Mediterranean Sea involves the intermediate Marmara Sea Basin via two sills (Dardanelles and Bosphorus). This study explores the chronology of late Pleistocene–Holocene flooding by examining sedimentary facies and molluscs from 24 gravity cores spanning shelf to slope settings in the southern Marmara Sea Basin. A late Pleistocene Ponto-Caspian (Neoeuxinian) mollusc association is found in 12 of the cores, comprising 14 mollusc species and dominated by brackish (oligohaline–lower mesohaline) endemic taxa (dreissenids, hydrobiids). The Neoeuxinian association is replaced by a Turritella–Corbula association at the onset of the Holocene. The latter is dominated by marine species, several of which are known to thrive under dysoxic conditions in muddy bottoms. This association is common in early Holocene intervals as well as sapropel intervals in younger Holocene strata. It is an indicator of low-salinity outflows from the Black Sea into the Marmara Sea that drive stratification. A marine Mediterranean association (87 species) represents both soft bottom and hard substrate faunas that lived in well-ventilated conditions and upper mesohaline–polyhaline salinities (ca. 25 psu). Shallower areas were occupied by hard substrate taxa and phytopdetritic communities, whereas deeper areas had soft bottom faunas. The middle shelf part of the northern Gemlik Gulf has intervals with irregular and discontinuous sedimentary structures admixed with worn Neoeuxinian and euryhaline Mediterranean faunas. These intervals represent reworking events (slumping) likely related to seismic activity rooted in the North Anatolian Fault system. The core data and faunas indicate an oscillating postglacial sea-level rise and phases of increased/decreased ventilation in the Marmara Sea during the Holocene, as well as palaeobiogeographic reorganisations of Ponto-Caspian and Mediterranean water bodies since the latest Pleistocene (<30 ka). The findings contribute to arguments against a single catastrophic flooding of the Black Sea at about 7.5 ka (Noah’s Flood).     

Acoustic monitoring of gas emissions from the seafloor. Part II: a case study from the Sea of Marmara

A rotating, acoustic gas bubble detector, BOB (Bubble OBservatory) module was deployed during two surveys, conducted in 2009 and 2011 respectively, to study the temporal variations of gas emissions from the Marmara seafloor, along the North Anatolian Fault zone. The echosounder mounted on the instrument insonifies an angular sector of 7° during a given duration (of about 1 h). Then it rotates to the next, near-by angular sector and so forth. When the full angular domain is insonified, the “pan and tilt system” rotates back to its initial position, in order to start a new cycle (of about 1 day). The acoustic data reveal that gas emission is not a steady process, with observed temporal variations ranging between a few minutes and 24 h (from one cycle to the other). Echo-integration and inversion performed on the acoustic data as described in the companion paper of Leblond et al. (Mar Geophys Res, 2014), also indicate important variations in, respectively, the target strength and the volumetric flow rates of individual sources. However, the observed temporal variations may not be related to the properties of the gas source only, but reflect possible variations in sea-bottom currents, which could deviate the bubble train towards the neighboring sector. During the 2011 survey, a 4-component ocean bottom seismometer (OBS) was co-located at the seafloor, 59 m away from the BOB module. The acoustic data from our rotating, monitoring system support, but do not provide undisputable evidence to confirm, the hypothesis formulated by Tary et al. (2012), that the short-duration, non-seismic micro-events recorded by the OBS are likely produced by gas-related processes within the near seabed sediments. Hence, the use of a multibeam echosounder, or of several split beam echosounders should be preferred to rotating systems, for future experiments.     

Temporal variability of tsunami arrival detection distance revealed by virtual tsunami observation experiments using numerical simulation and 1-month HF radar observation

The combination of a high-frequency ocean surface radar and a tsunami detection method should be assessed as the onshore-offshore distribution of tsunami detection probability, because the probability will vary in accordance with the signal-to-noise ratio (SNR) and the tsunami magnitude in addition to the radar system specifications. Here, we statistically examine the tsunami detection distance based on virtual tsunami observation experiments by using signals received by a high-frequency radar in February 2014 installed on the southern coast of Japan and numerically simulated velocities induced by a Nankai Trough earthquake. In the experiments, the Doppler frequencies associated with the simulated velocities were superimposed on the receiving signals of the radar, and the radial velocities were calculated from the synthesized signals by the fast Fourier transform. Tsunami arrival was then detected based on the temporal change in the cross-correlation of the velocities, before and after tsunami arrival, between two points 3 km apart along a radar beam. We found that the possibility of tsunami detection primarily depends on the kinetic energy ratio between tsunami current and background current velocities. The monthly average detection probability is over 90% when the energy ratio exceeds 5 (offshore distance: 9 km ≤ L ≤ 36 km) and reduces to 50% when the energy ratio is approximately 1 (L = 42 km) over the shelf slope. The ratio varied with the background current physics and SNR, which was mainly affected by ocean surface wave heights and ionospheric electron density.     

Evidence of NW extension of the North Anatolian Fault Zone in the Marmara Sea: a new interpretation of the Marmara Sea (Izmit) earthquake on 17 August 1999

Active faults aligning in NW–SE direction and forming flower structures of strike-slip faults were observed in shallow seismic data from the shelf offshore of Avcılar in the northern Marmara Sea. By following the parallel drainage pattern and scarps, these faults were traced as NW–SE-directed lineaments in the morphology of the northern onshore sector of the Marmara Sea (eastern Thrace Peninsula). Right-lateral displacements in two watersheds of drainage and on the coast of the Marmara Sea and Black Sea are associated with these lineaments. This right-lateral displacement along the course of these faults suggests a new, active strike-slip fault zone located at the NW extension of the northern boundary fault of the ?ınarcık Basin in the Marmara Sea. This new fault zone is interpreted as the NW extension of the northern branch of the North Anatolian Fault Zone (NAFZ), extending from the ?ınarcık Basin of the Marmara Sea to the Black Sea coast of the Thrace Peninsula, and passing through B üy ük ?ekmece and K ü ? ük ?ekmece lagoons. These data suggest that the rupture of the 17 August 1999 earthquake in the NAFZ may have extended through Avcılar. Indeed, Avcılar and İzmit, both located on the Marmara Sea coast along the rupture route, were strongly struck by the earthquake whereas the settlements between Avcılar and İzmit were much less affected. Therefore, this interpretation can explain the extraordinary damage in Avcılar, based on the newly discovered rupture of the NAFZ in the Marmara Sea. However, this suggestion needs to be confirmed by further seismological studies.     

Late Pleistocene–Holocene evolution of the southern Marmara shelf and sub-basins: middle strand of the North Anatolian fault,southern Marmara Sea,Turkey

Although there are many research studies on the northern and southern branches of the North Anatolian fault, cutting through the deep basins of the Sea of Marmara in the north and creating a series of pull-apart basins on the southern mainland, little data is available about the geometrical and kinematical characteristics of the middle strand of the North Anatolian fault. The first detailed geometry of the middle strand of the North Anatolian fault along the southern Marmara shelf, including the Gemlik and Band?rma Bay, will be given in this study, by a combined interpretation of different seismic data sets. The characteristic features of its segments and their importance on the paleogeographic evolution of the southern shelf sub-basins were defined. The longest one of these faults, the Armutlu-Band?rma segment, is a 75-km long dextral strike-slip fault which connects the W–E trending Gençali segment in the east and NE–SW trending Kap?da?-Edincik segment in the west. In this context, the Gemlik Bay opened as a pull-apart basin under the control of the middle strand whilst a new fault segment developed during the late Pleistocene, cutting through the eastern rim of the bay. In this region, a delta front forming the paleoshoreline of the Gemlik paleolake was cut and shifted approximately 60 ± 5 m by the new segment. The same offset on this fault was also measured on a natural scarp of acoustic basement to the west and integrated with this paleoshoreline forming the slightly descending topset–foreset reflections of the delta front. Therefore the new segment is believed to be active at least for the last 30,000 years. The annual lateral slip rate representing this period of time will be 2 mm, which is quite consistent with modern GPS measurements. Towards the west, the Band?rma Bay is a rectangular transpressional basin whilst the Erdek Bay is a passive basin under the control of NW–SE trending faults. When the water level of the paleo-Marmara lake dropped down to ?90 m, the water levels of the suspended paleolakes of Band?rma and Gemlik on the southern shelf were ?50.3 (?3.3 Global Isostatic Adjustment—GIA) and ?60.5 (?3.3 GIA) m below the present mean sea level, respectively. As of today a similar example can be seen between the Sea of Marmara and the shallow freshwater lakes of Manyas and Uluabat. Similarly, the paleolakes of Gemlik and Bandirma were affected by the water level fluctuations at different time periods, even though both lakes were isolated from the Sea of Marmara during the glacial periods.     

Characterizing air–sea CO2 exchange dynamics during winter in the coastal water off the Hugli-Matla estuarine system in the northern Bay of Bengal, India

The distribution of the fugacity of CO₂ ( $ f_{{{\text{CO}}_{ 2} }} $ ) and air–sea CO₂ exchange were comprehensively investigated in the outer estuary to offshore shallow water region (lying adjacent to the Sundarban mangrove forest) covering an area of ~2,000 km² in the northern Bay of Bengal during the winter. A total of ten sampling surveys were conducted between 1 December, 2011 and 21 February, 2012. Physico-chemical variables like sea surface temperature (SST), salinity, pH, total alkalinity (TAlk), dissolved inorganic carbon (DIC) and in vivo chlorophyll-a along with atmospheric variables were measured in order to study their role in controlling the CO₂ flux. Surface water $ f_{{{\text{CO}}_{ 2} }} $ ranged between 111 and 459 μatm which correlated significantly with the SST (r = 0.71, p < 0.001, n = 62). Neither DIC nor TAlk showed any linear relationship with varying salinity in the estuarine mixing zone, demonstrating the significant presence of non-carbonate alkalinity. An overall net biological control on the surface $ f_{{{\text{CO}}_{ 2} }} $ distribution was established during the study, although no significant correlation was found between chlorophyll-a and $ f_{{{\text{CO}}_{ 2} }} $ (water). The shallow water region studied was mostly under-saturated with CO₂ and acted as a sink for atmospheric CO₂. The difference between surface water and atmospheric $ f_{{{\text{CO}}_{ 2} }} $ ( $ \Updelta f_{{{\text{CO}}_{ 2} }} $ ) ranged from ?274 to 69 μatm, with an average seaward flux of ?10.5 ± 12.6 μmol m^?2 h^?1. The $ \Updelta f_{{{\text{CO}}_{ 2} }} $ and hence the air–sea CO₂ exchange was primarily regulated by the variation in sea surface $ f_{{{\text{CO}}_{ 2} }} $ , since atmospheric $ f_{{{\text{CO}}_{ 2} }} $ varied over a comparatively narrow range of 361.23–399.05 μatm.     

Reply to Discussion: a critique of Possible waterways between the Marmara Sea and the Black Sea in the late Quaternary: evidence from ostracod and foraminifer assemblages in lakes ?znik and Sapanca, Turkey, Geo-Marine Letters, 2011

In their discussion of our 2011 paper dealing with possible waterways between the Marmara Sea and the Black Sea in the ??late?? Quaternary, based on data from ostracod and foraminifer assemblages in lakes ?znik and Sapanca, Turkey, Yalt?rak et al. (Geo-Mar Lett 32:267?C274, 2012) essentially reject the idea of any links whatsoever, be they between the Marmara Sea and the lakes ?znik and Sapanca, or further to the Black Sea via the valley of the Sakarya River. The evidence they provide in support of their view, however, is essentially circumstantial, in part conjectural, and also inconclusive considering the findings in favour of linkage between the Marmara Sea and the lakes at the very least, while the proposed connection with the Sakarya River valley remains speculative because of the lack of unambiguous data. On the other hand, Yalt?rak et al. (Geo-Mar Lett 32:267?C274, 2012) do raise valid points of concern which deserve careful future investigation, the most important being the possibility of sample contamination from dumped marine sediment used for construction purposes along some parts of the shore of Lake ?znik. We agree that a concerted multidisciplinary effort is required to address the many unresolved issues in connection with the potential waterways proposed by us and others before us.     

Geophysical investigations of crust-scale structural model of the Qiongdongnan Basin, Northern South China Sea

Rifting of the Qiongdongnan Basin was initiated in the Cenozoic above a pre-Cenozoic basement, which was overprinted by extensional tectonics and soon after the basin became part of the rifted passive continental margin of the South China Sea. We have integrated available grids of sedimentary horizons, wells, seismic reflection data, and the observed gravity field into the first crust-scale structural model of the Qiongdongnan Basin. Many characteristics of this model reflect the tectonostratigraphic history of the basin. The structure and isopach maps of the basin allow us to reconstruct the history of the basin comprising: (a) The sediments of central depression are about 10 km thicker than on the northern and southern sides; (b) The sediments in the western part of the basin are about 6 km thicker than that in the eastern part; (c) a dominant structural trend of gradually shifting depocentres from the Paleogene sequence (45–23.3 Ma) to the Neogene to Quaternary sequence (23.3 Ma–present) towards the west or southwest. The present-day configuration of the basin reveals that the Cenozoic sediments are thinner towards the east. By integrating several reflection seismic profiles, interval velocity and performing gravity modeling, we model the sub-sedimentary basement of the Qiongdongnan Basin. There are about 2–4 km thick high-velocity bodies horizontal extended for a about 40–70 km in the lower crust (v > 7.0 km/s) and most probably these are underplated to the lower stretched continental crust during the final rifting and early spreading phase. The crystalline continental crust spans from the weakly stretched domains (about 25 km thick) near the continental shelf to the extremely thinned domains (<2.8 km) in the central depression, representing the continental margin rifting process in the Qiongdongnan Basin. Our crust-scale structural model shows that the thinnest crystalline crust (<3 km) is found in the Changchang Sag located in the east of the basin, and the relatively thinner crystalline crust (<3.5 km) is in the Ledong Lingshui Sag in the west of the basin. The distribution of crustal extension factor β show that β in central depression is higher (>7.0), while that on northern and southern sides is lower (<3.0). This model can illuminate future numerical simulations, including the reconstruction of the evolutionary processes from the rifted basin to the passive margin and the evolution of the thermal field of the basin.