Physical oceanography during Maria S. Merian cruise MSM21/2

Flemish Pass, located at the western subpolar margin, is a passage (sill depth 1200 m) that is constrained by the Grand Banks and the underwater plateau Flemish Cap. In addition to the Deep Western Boundary Current (DWBC) pathway offshore of Flemish Cap, Flemish Pass represents another southward transport pathway for two modes of Labrador Sea Water (LSW), the lightest component of North Atlantic Deep Water carried with the DWBC. This pathway avoids potential stirring regions east of Flemish Cap and deflection into the interior North Atlantic. Ship-based velocity measurements between 2009 and 2013 at 47°N in Flemish Pass and in the DWBC east of Flemish Cap revealed a considerable southward transport of Upper LSW through Flemish Pass (15-27%, -1.0 to -1.5 Sv). About 98% of the denser Deep LSW were carried around Flemish Cap as Flemish Pass is too shallow for considerable transport of Deep LSW. Hydrographic time series from ship-based measurements show a significant warming of 0.3°C/decade and a salinification of 0.03/decade of the Upper LSW in Flemish Pass between 1993 and 2013. Almost identical trends were found for the evolution in the Labrador Sea and in the DWBC east of Flemish Cap. This indicates that the long-term hydrographic variability of Upper LSW in Flemish Pass as well as in the DWBC at 47°N is dominated by changes in the Labrador Sea, which are advected southward. Fifty years of numerical ocean model simulations in Flemish Pass suggest that these trends are part of a multidecadal cycle.

Organic geochemistry of sediment cores from the Arbian Sea

Two deep-sea sediment cores from the northeastern and the southeastern Arabian Sea were studied in order to reconstruct the palaeoenvironments of the past glacial cycles. Core 136KL was recovered from the high-productivity area off Pakistan within the modern oxygen-minimum zone (OMZ). By contrast, modern primary productivity at the site of MD900963 close to Maldives is moderate and bottom waters are today well oxygenated. For both cores, we reconstructed the changes in palaeoproductivity using a set of biomarkers (alkenones, dinosterol and brassicasterol); the main result is that primary productivity is enhanced during glacial stages and lowered during interstadials. The proxies associated with productivity show a 23 kyr cyclicity corresponding to the precession-related insolation cycle. Palaeoredox conditions were studied in both cores using a new organic geochemical parameter (C35/C31-n-alkane ratio) developed by analysing surface sediments from a transect across the OMZ off Pakistan. The value of this ratio in core 136KL shows many variations during the last 65 kyr, indicating that the OMZ was not stable during this time: it disappeared completely during Heinrich- and the Younger Dryas events, pointing to a connection between global oceanic circulation and the stability of the OMZ. The C35/C31 ratio determined in sediments of core MD900963 shows that bottom waters remained rather well oxygenated over the last 330 kyr, which is confirmed by comparison with authigenic metal concentrations in the same sediments. A zonally averaged, circulation-biogeochemical ocean model was used to explore how the intermediate Indian Ocean responds to a freshwater flux anomaly at the surface of the North Atlantic. As suggested by the geochemical time series, both the abundance of Southern Ocean Water and the oxygen concentration are significantly increased in response to this freshwater perturbation.

Simulated temperature changes as a result of Antarctic ice sheet growth and CO2 reduction

During the Middle Miocene climate transition about 14 million years ago, the Antarctic ice sheet expanded to near-modern volume. Surprisingly, this ice sheet growth was accompanied by a warming in the surface waters of the Southern Ocean, whereas a slight deep-water temperature increase was delayed by more than 200 thousand years. Here we use a coupled atmosphere-ocean model to assess the relative effects of changes in atmospheric CO2 concentration and ice sheet growth on regional and global temperatures. In the simulations, changes in the wind field associated with the growth of the ice sheet induce changes in ocean circulation, deep-water formation and sea-ice cover that result in sea surface warming and deep-water cooling in large swaths of the Atlantic and Indian ocean sectors of the Southern Ocean. We interpret these changes as the dominant ocean surface response to a 100-thousand-year phase of massive ice growth in Antarctica. A rise in global annual mean temperatures is also seen in response to increased Antarctic ice surface elevation. In contrast, the longer-term surface and deep-water temperature trends are dominated by changes in atmospheric CO2 concentration. We therefore conclude that the climatic and oceanographic impacts of the Miocene expansion of the Antarctic ice sheet are governed by a complex interplay between wind field, ocean circulation and the sea-ice system.

Holocene sediment transport from the glaciated margin of East/Northeast Greenland to the N Iceland shelves

We examine variations in the ice-rafted sources for sediments in the Iceland/East Greenland offshore marine archives by utilizing a sediment unmixing model and link the results to a coupled iceberg-ocean model. Surface samples from around Iceland and along the E/NE Greenland shelf are used to define potential sediment sources, and these are examined within the context of the down-core variations in mineralogy in the <2 mm sediment fraction from a transect of cores across Denmark Strait. A sediment unmixing model is used to estimate the fraction of sediment <2 mm off NW and N Iceland exported across Denmark Strait; this averaged between 10 and 20%. Both the sediment unmixing model and the coupled iceberg-ocean model are consistent in finding that the fraction of "far-travelled" sediments in the Denmark Strait environs is overwhelmingly of local, mid-East Greenland, provenance, and therefore with a significant cross-channel component to their travel. The Holocene record of ice-rafted sediments denotes a three-part division of the Holocene in terms of iceberg sediment transport with a notable increase in the process starting ca 4000 cal yr BP. This latter increase may represent the re-advance during the Neoglacial period of land-terminating glaciers on the Geikie Plateau to become marine-terminating. The contrast in spectral signals between these cores and the 1500-yr cycle at VM28-14, just south of the Denmark Strait, combined with the coupled iceberg-model results, leads us to speculate that the signal at VM28-14 reflects pulses in overflow waters, rather than an ice-rafted signal.

Benthic foraminifera across the Cretaceous/Paleogene boundary in ODP Hole 198-1210A

Deep-sea benthic foraminifera show important but transient assemblage changes at the Cretaceous/Paleogene (K/Pg) boundary, when many biota suffered severe extinction. We quantitatively analyzed benthic foraminiferal assemblages from lower bathyal-upper abyssal (1500-2000 m) northwest Pacific ODP Site 1210 (Shatsky Rise) and compared the results with published data on assemblages at lower bathyal (~ 1500 m) Pacific DSDP Site 465 (Hess Rise) to gain insight in paleoecological and paleoenvironmental changes at that time. At both sites, diversity and heterogeneity rapidly decreased across the K/Pg boundary, then recovered. Species assemblages at both sites show a similar pattern of turnover from the uppermost Maastrichtian into the lowermost Danian: 1) The relative abundance of buliminids (indicative of a generally high food supply) increases towards the uppermost Cretaceous, and peaks rapidly just above the K/Pg boundary, coeval with a peak in benthic foraminiferal accumulation rate (BFAR), a proxy for food supply. 2) A peak in relative abundance of Stensioeina beccariiformis, a cosmopolitan form generally more common at the middle than at the lower bathyal sites, occurs just above the buliminid peak. 3) The relative abundance of Nuttallides truempyi, a more oligotrophic form, decreases at the boundary, then increases above the peak in Stensioeina beccariiformis. The food supply to the deep sea in the Pacific Ocean thus apparently increased rather than decreased in the earliest Danian. The low benthic diversity during a time of high food supply indicates a stressed environment. This stress might have been caused by reorganization of the planktic ecosystem: primary producer niches vacated by the mass extinction of calcifying nannoplankton may have been rapidly (<10 kyr) filled by other, possibly opportunistic, primary producers, leading to delivery of another type of food, and/or irregular food delivery through a succession of opportunistic blooms. The deep-sea benthic foraminiferal data thus are in strong disagreement with the widely accepted hypothesis that the global deep-sea floor became severely food-depleted following the K/Pg extinction due to the mass extinction of primary producers ("Strangelove Ocean Model") or to the collapse of the biotic pump ("Living Ocean Model").

On the steric and mass-induced contributions to the annual sea level variations in the Mediterranean Sea

The sea level variation (SLVtotal) is the sum of two major contributions: steric and mass-induced. The steric SLVsteric is that resulting from the thermal and salinity changes in a given water column. It only involves volume change, hence has no gravitational effect. The mass-induced SLVmass, on the other hand, arises from adding or subtracting water mass to or from the water column and has direct gravitational signature. We examine the closure of the seasonal SLV budget and estimate the relative importance of the two contributions in the Mediterranean Sea as a function of time. We use ocean altimetry data (from TOPEX/Poseidon, Jason 1, ERS, and ENVISAT missions) to estimate SLVtotal, temperature, and salinity data (from the Estimating the Circulation and Climate of the Ocean ocean model) to estimate SLVsteric, and time variable gravity data (from Gravity Recovery and Climate Experiment (GRACE) Project, April 2002 to July 2004) to estimate SLVmass. We find that the annual cycle of SLVtotal in the Mediterranean is mainly driven by SLVsteric but moderately offset by SLVmass. The agreement between the seasonal SLVmass estimations from SLVtotal – SLVsteric and from GRACE is quite remarkable; the annual cycle reaches the maximum value in mid-February, almost half a cycle later than SLVtotal or SLVsteric, which peak by mid-October and mid-September, respectively. Thus, when sea level is rising (falling), the Mediterranean Sea is actually losing (gaining) mass. Furthermore, as SLVmass is balanced by vertical (precipitation minus evaporation, P–E) and horizontal (exchange of water with the Atlantic, Black Sea, and river runoff) mass fluxes, we compared it with the P–E determined from meteorological data to estimate the annual cycle of the horizontal flux.

Physical oceanography during three Maria S. Merian cruises (MSM21/2, MSM27, MSM28)

Flemish Pass, located at the western subpolar margin, is a passage (sill depth 1200 m) that is constrained by the Grand Banks and the underwater plateau Flemish Cap. In addition to the Deep Western Boundary Current (DWBC) pathway offshore of Flemish Cap, Flemish Pass represents another southward transport pathway for two modes of Labrador Sea Water (LSW), the lightest component of North Atlantic Deep Water carried with the DWBC. This pathway avoids potential stirring regions east of Flemish Cap and deflection into the interior North Atlantic. Ship-based velocity measurements between 2009 and 2013 at 47°N in Flemish Pass and in the DWBC east of Flemish Cap revealed a considerable southward transport of Upper LSW through Flemish Pass (15-27%, -1.0 to -1.5 Sv). About 98% of the denser Deep LSW were carried around Flemish Cap as Flemish Pass is too shallow for considerable transport of Deep LSW. Hydrographic time series from ship-based measurements show a significant warming of 0.3°C/decade and a salinification of 0.03/decade of the Upper LSW in Flemish Pass between 1993 and 2013. Almost identical trends were found for the evolution in the Labrador Sea and in the DWBC east of Flemish Cap. This indicates that the long-term hydrographic variability of Upper LSW in Flemish Pass as well as in the DWBC at 47°N is dominated by changes in the Labrador Sea, which are advected southward. Fifty years of numerical ocean model simulations in Flemish Pass suggest that these trends are part of a multidecadal cycle.

Predictable hydrodynamic conditions explain temporal variations in the density of benthic foraging seabirds in a tidal stream environment

Tidal stream turbines could have several direct impacts upon pursuit-diving seabirds foraging within tidal stream environments (mean horizontal current speeds > 2 ms−1), including collisions and displacement. Understanding how foraging seabirds respond to temporally variable but predictable hydrodynamic conditions immediately around devices could identify when interactions between seabirds and devices are most likely to occur; information which would quantify the magnitude of potential impacts, and also facilitate the development of suitable mitigation measures. This study uses shore-based observational surveys and Finite Volume Community Ocean Model outputs to test whether temporally predictable hydrodynamic conditions (horizontal current speeds, water elevation, turbulence) influenced the density of foraging black guillemots Cepphus grylle and European shags Phalacrocorax aristotelis in a tidal stream environment in Orkney, United Kingdom, during the breeding season. These species are particularly vulnerable to interactions with devices due to their tendency to exploit benthic and epi-benthic prey on or near the seabed. The density of both species decreased as a function of horizontal current speeds, whereas the density of black guillemots also decreased as a function of water elevation. These relationships could be linked to higher energetic costs of dives in particularly fast horizontal current speeds (>3 ms−1) and deeper water. Therefore, interactions between these species and moving components seem unlikely at particularly high horizontal current speeds. Combining this information, with that on the rotation rates of moving components at lower horizontal current speeds, could be used to assess collision risk in this site during breeding seasons. It is also likely that moderating any device operation during both lowest water elevation and lowest horizontal current speeds could reduce the risk of collisions for these species in this site during this season. The approaches used in this study could have useful applications within Environmental Impact Assessments, and should be considered when assessing and mitigating negative impacts from specific devices within development sites.

Stable carbon isotope ratios and accumulation rates of organic carbon and nutrients of DSDP Hole 79-545 (Appendix A)

Future warming is predicted to shift the Earth system into a mode with progressive increase and vigour of extreme climate events possibly stimulating other mechanisms that invigorate global warming. This study provides new data and modelling investigating climatic consequences and biogeochemical feedbacks that happened in a warmer world ~112 Myr ago. Our study focuses on the Cretaceous Oceanic Anoxic Event (OAE) 1b and explores how the Earth system responded to a moderate ~25,000 yr lasting climate perturbation that is modelled to be less than 1 °C in global average temperature. Using a new chronological model for OAE 1b we present high-resolution elemental and bulk carbon isotope records from DSDP Site 545 from Mazagan Plateau off NW Africa and combine this information with a coupled atmosphere-land-ocean model. The simulations suggest that a perturbation at the onset of OAE 1b caused almost instantaneous warming of the atmosphere on the order of 0.3 °C followed by a longer (~45,000 yr) period of ~0.8 °C cooling. The marine records from DSDP Site 545 support that these moderate swings in global climate had immediate consequences for African continental supply of mineral matter and nutrients (phosphorous), subsequent oxygen availability, and organic carbon burial in the eastern subtropical Atlantic, however, without turning the ocean anoxic. The match between modelling results and stratigraphic isotopic data support previous studies [summarized in Jenkyns 2003, doi:10.1098/rsta.2003.1240] in that methane emission from marine hydrates, albeit moderate in dimension, may have been the trigger for OAE 1b, though we can not finally rule out alternative mechanisms. Following the hydrate mechanism a total of 1.15 * 10**18 g methane carbon (delta13C=-60 ?), equivalent to about 10% to the total modern gas hydrate inventory, generated the delta13Ccarb profile recorded in the section. Modelling suggests a combination of moderate-scale methane pulses supplemented by continuous methane emission at elevated levels over ~25,000 yr. The proposed mechanism, though difficult to finally confirm in the geological past, is arguably more likely to occur in a warmer world and apparently perturbs global climate and ocean chemistry almost instantaneously. This study shows that, once set-off, this mechanism can maintain Earth's climate in a perturbed mode over geological time leading to pronounced changes in regional climate.

Rotina computacional de obtenção, tratamento e visualização de dados oceanográficos

Atualmente, sensores remotos e computadores de alto desempenho estão sendo utilizados como instrumentos principais na coleta e produção de dados oceanográficos. De posse destes dados, é possível realizar estudos que permitem simular e prever o comportamento do oceano por meio de modelos numéricos regionais. Dentre os fatores importantes no estudo da oceanografia, podem ser destacados àqueles referentes aos impactos ambientais, de contaminação antrópica, utilização de energias renováveis, operações portuárias e etc. Contudo, devido ao grande volume de dados gerados por instituições ambientais, na forma de resultados de modelos globais como o HYCOM (Hybrid Coordinate Ocean Model) e dos programas de Reanalysis da NOAA (National Oceanic and Atmospheric Administration), torna-se necessária a criação de rotinas computacionais para realizar o tratamento de condições iniciais e de contorno, de modo que possam ser aplicadas a modelos regionais como o TELEMAC3D (www.opentelemac.org). Problemas relacionados a baixa resolução, ausência de dados e a necessidade de interpolação para diferentes malhas ou sistemas de coordenadas verticais, tornam necessária a criação de um mecanismo computacional que realize este tratamento adequadamente. Com isto, foram desenvolvidas rotinas na linguagem de programação Python, empregando interpoladores de vizinho mais próximo, de modo que, a partir de dados brutos dos modelos HYCOM e do programa de Reanalysis da NOAA, foram preparadas condições iniciais e de contorno para a realização de uma simulação numérica teste. Estes resultados foram confrontados com outro resultado numérico onde, as condições foram construídas a partir de um método de interpolação mais sofisticado, escrita em outra linguagem, e que já vem sendo utilizada no laboratório. A análise dos resultados permitiu concluir que, a rotina desenvolvida no âmbito deste trabalho, funciona adequadamente para a geração de condições iniciais e de contorno do modelo TELEMAC3D. Entretanto, um interpolador mais sofisticado deve ser desenvolvido de forma a aumentar a qualidade nas interpolações, otimizar o custo computacional, e produzir condições que sejam mais realísticas para a utilização do modelo TELEMAC3D.

Study of the development and evolution of seasonal thermocline and turbulent processes in the northern parts of Persian Gulf using numerical modeling

The Persian Gulf (PG) is a semi-enclosed shallow sea which is connected to open ocean through the Strait of Hormuz. Thermocline as a suddenly decrease of temperature in subsurface layer in water column leading to stratification happens in the PG seasonally. The forcing comprise tide, river inflow, solar radiation, evaporation, northwesterly wind and water exchange with the Oman Sea that influence on this process. In this research, analysis of the field data and a numerical (Princeton Ocean Model, POM) study on the summer thermocline development in the PG are presented. The Mt. Mitchell cruise 1992 salinity and temperature observations show that the thermocline is effectively removed due to strong wind mixing and lower solar radiation in winter but is gradually formed and developed during spring and summer; in fact as a result of an increase in vertical convection through the water in winter, vertical gradient of temperature is decreased and thermocline is effectively removed. Thermocline development that evolves from east to west is studied using numerical simulation and some existing observations. Results show that as the northwesterly wind in winter, at summer transition period, weakens the fresher inflow from Oman Sea, solar radiation increases in this time interval; such these factors have been caused the thermocline to be formed and developed from winter to summer even over the northwestern part of the PG. The model results show that for the more realistic monthly averaged wind experiments the thermocline develops as is indicated by summer observations. The formation of thermocline also seems to decrease the dissolved oxygen in water column due to lack of mixing as a result of induced stratification. Over most of PG the temperature difference between surface and subsurface increases exponentially from March until May. Similar variations for salinity differences are also predicted, although with smaller values than observed. Indeed thermocline development happens more rapidly in the Persian Gulf from spring to summer. Vertical difference of temperature increases to 9 centigrade degrees in some parts of the case study zone from surface to bottom in summer. Correlation coefficients of temperature and salinity between the model results and measurements have been obtained 0.85 and 0.8 respectively. The rate of thermcline development was found to be between 0.1 to 0.2 meter per day in the Persian Gulf during the 6 months from winter to early summer. Also it is resulted from the used model that turbulence kinetic energy increases in the northwestern part of the PG from winter to early summer that could be due to increase in internal waves activities and stability intensified through water column during this time.

Tide-topography interaction along the eastern Brazilian shelf

In the Tropics, continental shelves governed by western boundary currents are considered to be among the least productive ocean margins in the world, unless eddy-induced shelf-edge upwelling becomes significant. The eastern Brazilian shelf in the Southwest Atlantic is one of these, and since the slight nutrient input from continental sources is extremely oligotrophic. It is characterized by complex bathymetry with the presence of shallow banks and seamounts. In this work, a full three-dimensional nonlinear primitive equation ocean model is used to demonstrate that the interaction of tidal currents and the bottom topography of the east Brazil continental shelf is capable of producing local upwelling of South Atlantic Central Water, bringing nutrients up from deep waters to the surface layer. Such upper layer enrichment is found to be of significance in increasing local primary productivity. (c) 2005 Elsevier Ltd. All rights reserved.

Coastal modelling studies for forecasting and remediation solutions

The coastal area along the Emilia-Romagna (ER), in the Italian side of the northern Adriatic Sea, is considered to implement an unstructured numerical ocean model with the aim to develop innovative tools for the coastal management and a forecasting system for the storm surge risk reduction. The Adriatic Sea has been the focus of several studies because of its peculiar dynamics driven by many forcings acting at basin and local scales. The ER coast is particularly exposed to storm surge events. In particular conditions, winds, tides and seicehs may combine and contribute to the flooding of the coastal area. The global sea level rise expected in the next decades will increase even more the hazard along the ER and Adriatic coast. Reliable Adriatic and Mediterranean scale numerical ocean models are now available to allow the dynamical downscaling of very high-resolution models in limited coastal areas. In this work the numerical ocean model SHYFEM is implemented in the Goro lagoon (named GOLFEM) and along the ER coast (ShyfER) to test innovative solutions against sea related coastal hazards. GOLFEM was succesfully applied to analyze the Goro lagoon dynamics and to assess the dynamical effects of human interventions through the analysis of what-if scenarios. The assessment of storm surge hazard in the Goro lagoon was carried out through the development of an ensemble storm surge forecasting system with GOLFEM using forcing from different operational meteorological and ocean models showing the fundamental importance of the boundary conditions. The ShyfER domain is used to investigate innovative solutions against storm surge related hazard along the ER coast. The seagrass is assessed as a nature-based solution (NBS) for coastal protection under present and future climate conditions. The results show negligible effects on sea level but sensible effects in reducing bottom current velocity.

Study of the equatorial Atlantic Ocean mixing layer using a one-dimensional turbulence model

The General Ocean Turbulence Model (GOTM) is applied to the diagnostic turbulence field of the mixing layer (ML) over the equatorial region of the Atlantic Ocean. Two situations were investigated: rainy and dry seasons, defined, respectively, by the presence of the intertropical convergence zone and by its northward displacement. Simulations were carried out using data from a PIRATA buoy located on the equator at 23º W to compute surface turbulent fluxes and from the NASA/GEWEX Surface Radiation Budget Project to close the surface radiation balance. A data assimilation scheme was used as a surrogate for the physical effects not present in the one-dimensional model. In the rainy season, results show that the ML is shallower due to the weaker surface stress and stronger stable stratification; the maximum ML depth reached during this season is around 15 m, with an averaged diurnal variation of 7 m depth. In the dry season, the stronger surface stress and the enhanced surface heat balance components enable higher mechanical production of turbulent kinetic energy and, at night, the buoyancy acts also enhancing turbulence in the first meters of depth, characterizing a deeper ML, reaching around 60 m and presenting an average diurnal variation of 30 m.