ALBERT

Description: In coastal marine environments, physical and biological forces can cause dynamic pH fluctuations from microscale (diffusive boundary layer [DBL]) up to ecosystem‐scale (benthic boundary layer [BBL]). In the face of ocean acidification (OA), such natural pH variations may modulate an organism's response to OA by providing temporal refugia. We investigated the effect of pH fluctuations, generated by the brown alga Fucus serratus' biological activity, on the calcifying epibionts Balanus improvisus and Electra pilosa under OA. For this, both epibionts were grown on inactive and biologically active surfaces and exposed to (1) constant pH scenarios under ambient (pH 8.1) or OA conditions (pH 7.7), or (2) oscillating pH scenarios mimicking BBL conditions at ambient (pH 7.7–8.6) or OA scenarios (pH 7.4–8.2). Furthermore, all treatment combinations were tested at 10°C and 15°C. Against our expectations, OA treatments did not affect epibiont growth under constant or fluctuating (BBL) pH conditions, indicating rather high robustness against predicted OA scenarios. Furthermore, epibiont growth was hampered and not fostered on active surfaces (fluctuating DBL conditions), indicating that fluctuating pH conditions of the DBL with elevated daytime pH do not necessarily provide temporal refugia from OA. In contrast, results indicate that factors other than pH may play larger roles for epibiont growth on macrophytes (e.g., surface characteristics, macrophyte antifouling defense, or dynamics of oxygen and nutrient concentrations). Warming enhanced epibiont growth rates significantly, independently of OA, indicating no synergistic effects of pH treatments and temperature within their natural temperature range.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: The growth and decay mechanisms of barrier layers in the northwestern tropical Atlantic are studied by investigating small‐scale processes embedded in the regional circulation of the tropical Atlantic using output from an eddy‐resolving numerical simulation at 4 km resolution forced by an atmospheric reanalysis. The simulation reproduces well the temporal and spatial patterns of barrier layer thickness (BLT) estimated with Argo and CTD in situ profiles. As seen from an analysis of the salinity and temperature vertical gradient balances, localized large barrier layers form inside North Brazil Current rings during late‐June to July because of a thickening of the isothermal layer in the rings due to horizontal temperature advection, stretching of isotherms and tilting of temperature fronts. These barrier layers decay when the isothermal layer reduces again due to the above mechanisms. Further to the north, along the North Equatorial Current, the seasonal variability of BLT is highly pronounced. Thick winter (January to early March) barrier layers locally grow as the base of the mixed layer shoals mainly due to a tilting of the salinity fronts and partly due to stretching of the isohalines, horizontal salt advection and vertical turbulent mixing. The short‐term barrier layers in this case decay due to a deepening of the mixed layer, whereas they get completely eroded in spring by a shoaling of the isothermal layer due to surface temperature stratification. This work highlights that barrier layers are localized phenomena at times growing solely due to ocean dynamics, without a surface freshwater influx.

Description: Plain Language Summary: Oceanic barrier layers exist in regions where salinity is more dominant than temperature in determining upper ocean density. Those layers lie between the bases of a constant‐density layer and a constant‐temperature layer. Barrier layers prevent vertical exchange of energy and mass between the near‐surface and the deep ocean, thus influencing air‐sea interaction. In the western tropical Atlantic, a warmer sea surface due to the presence of barrier layers can fuel hurricanes. Freshwater from the Amazon River and rainfall facilitate the growth of barrier layers, but the dynamics of their evolution are unclear. In this work, we identify/quantify the growth and decay mechanisms of barrier layers using a 4 km resolution simulation. Barrier layers grow/decay inside North Brazil Current eddies in summer because of deepening/shoaling of the constant‐temperature layer inside the eddies due to horizontal heat transport. Further north, barrier layers grow in winter as the constant‐density layer shoals mainly due to northwestward surface freshwater flow and equatorward subsurface salty water flow. Those barrier layers decay when the constant‐density layer deepens, whereas are destroyed when the constant‐temperature layer shoals in spring due to surface heating. These novel results improve the knowledge on barrier layers and help representing them in climate models.

Description: Key Points: The North Brazil and North Equatorial Currents are two regions with quasi‐permanent barrier layers in the northwestern tropical Atlantic Large barrier layer thickness (BLT) within the rings occurs during June‐July due to a thickening of the isothermal layer within the eddies Large winter BLT in North Equatorial Current is due to tilting of salinity fronts, stretching of isohalines, advection and turbulent mixing

Description: Deutsche Forschungsgemeinschaft (DFG) http://dx.doi.org/10.13039/501100001659

Description: Understanding groundwater–surface water (GW–SW) interactions is vital for water management in karstic catchments due to its impact on water quality. The objective of this study was to evaluate and compare the applicability of seven environmental tracers to quantify and localize groundwater exfiltration into a small, human‐impacted karstic river system. Tracers were selected based on their emission source to the surface water either as (a) dissolved, predominantly geogenic compounds (radon‐222, sulphate and electrical conductivity) or (b) anthropogenic compounds (predominantly) originating from wastewater treatment plant (WWTP) effluents (carbamazepine, tramadol, sodium, chloride). Two contrasting sampling approaches were compared (a) assuming steady‐state flow conditions and (b) considering the travel time of the water parcels (Lagrangian sampling) through the catchment to account for diurnal changes in inflow from the WWTP. Spatial variability of the concentrations of all tracers indicated sections of preferential groundwater inflow. Lagrangian sampling techniques seem highly relevant for capturing dynamic concentration patterns of WWTP‐derived compounds. Quantification of GW inflow with the finite element model FINIFLUX, based on observed in‐stream Rn activities led to plausible fluxes along the investigated river reaches (0.265 m3 s−1), while observations of other natural or anthropogenic environmental tracers produced less plausible water fluxes. Important point sources of groundwater exfiltration can be ascribed to locations where the river crosses geological fault lines. This indicates that commonly applied concepts describing groundwater–surface water interactions assuming diffuse flow in porous media are difficult to transfer to karstic river systems whereas concepts from fractured aquifers may be more applicable. In general, this study helps selecting the best suited hydrological tracer for GW exfiltration and leads to a better understanding of processes controlling groundwater inflow into karstic river systems.

Description: Karst aquifers represent an increased complexity when aiming to measure the interaction between groundwater and river water. Combining field‐based measurements on catchment scale and modelling, the applicability of ‘classical’ environmental groundwater tracers was compared to selected organic (micro)pollutants often considered as conservative and originally arising from a wastewater treatment plant. This study demonstrates that the choice of an appropriate tracer is crucial when either aiming to quantify groundwater exfiltration into karstic river systems, or indicating hydrological processes, applying (globally) omnipresent pollutants.

Description: German Research Foundation (DFG) http://dx.doi.org/10.13039/501100001659

Description: The general task of image classification seems to be solved due to the development of modern convolutional neural networks (CNNs). However, the high intraclass variability and interclass similarity of plankton images still prevents the practical identification of morphologically similar organisms. This prevails especially for rare organisms. Every CNN requires a vast amount of manually validated training images which renders it inefficient to train study‐specific classifiers. In most follow‐up studies, the plankton community is different from before and this data set shift (DSS) reduces the correct classification rates. A common solution is to discard all uncertain images and hope that the remains still resemble the true field situation. The intention of this North Sea Video Plankton Recorder (VPR) study is to assess if a combination of a Capsule Neural Network (CapsNet) with probability filters can improve the classification success in applications with DSS. Second, to provide a guideline how to customize automated CNN and CapsNet deep learning image analysis methods according to specific research objectives. In community analyses, our approach achieved a discard of uncertain predictions of only 5%. CapsNet and CNN reach similar precision scores, but the CapsNet has lower recall scores despite similar discard ratios. This is due to a higher discard ratio in rare classes. The recall advantage of the CNN decreases with increasing DSS. We present an alternative method to handle rare classes with a CNN achieving a mean recall of 96% by manually validating an average of 6.5% of the original images.

Description: The partial pressure of carbon dioxide (pCO2) in surface seawater is an important biogeochemical variable because, together with the pCO2 in the atmosphere, it determines the direction of air–sea carbon dioxide exchange. Large‐scale observations of pCO2 are facilitated by Ships‐of‐Opportunity (SOOP‐CO2) equipped with underway measuring instruments. The need for expanding the observation capacity and the challenges involving the sustainability and maintenance of traditional equilibrator systems led the community toward developing simpler and more autonomous systems. Here we performed a comparison between a membrane‐based sensor and a showerhead equilibration sensor installed on two SOOP‐CO2 between 2013 and 2018. We identified time‐ and space‐adequate crossovers in the Skagerrak Strait, where the two ship routes often crossed. We found a mean total difference of 1.5 ± 10.6 μatm and a root mean square error of 11 μatm. The pCO2 values recorded by the two instruments showed a strong linear correlation with a coefficient of 0.91 and a slope of 1.07 (± 0.14), despite the dynamic nature of the environment and the difficulty of comparing measurements from two different vessels. The membrane‐based sensor was integrated with a FerryBox system on a ship with a high sampling frequency in the study area. We showed the strength of having a sensor‐based network with a high spatial coverage that can be validated against conventional SOOP‐CO2 methods. Proving the validity of membrane‐based sensors in coastal and continental shelf seas and using the higher frequency measurements they provide can enable a thorough characterization of pCO2 variability in these dynamic environments.

Description: The addition of offshore wind farms (OWFs) to stratified regions of shelf seas poses an anthropogenic source of turbulence, in which the foundation structures remove power from the oceanic flow that is fed into turbulent mixing in the wake downstream. The loss of stratification within the wake of a single OWF structure is observed for the first time by means of field observations, which enable a qualitative characterization of the disturbed flow downstream. These results are complemented with high‐resolution large eddy simulations of four different stratification strengths that allow for a quantification of turbulence and mixing quantities in the wake of a foundation structure. The turbulent wake of a structure is narrow and highly energetic within the first 100 m, with the dissipation of turbulent kinetic energy well above background levels downstream of the structure. A single monopile is responsible for 7–10% additional mixing to that of the bottom mixed layer, whereby ∼10% of the turbulent kinetic energy generated by the structure is used in mixing. Although the effect of a single turbine on stratification is relatively low, large‐scale OWFs could significantly affect the vertical structure of a weakly stratified water column. Further, rough estimates show that the rate of formation of stratification in the study area is of the same order of magnitude as the additional mixing promoted by the structures, thus OWFs could modify the stratification regime and water column dynamics on a seasonal scale, depending on local conditions and farm geometries.

Description: Plain Language Summary: Advances in the renewable energy sector have enabled the construction and operation of wind farms in bodies of water deep enough to present vertical temperature differences across the water column or thermal stratification. In coastal regions dominated by tidal motion, the presence of offshore wind farm (OWF) structures brings about additional turbulence and mixing of stratification. The present study combines field measurements and numerical simulations to characterize the wake of single OWF structures and quantify the amount of turbulence and mixing generated by them. Our results suggest that the effect of OWF structures is small compared to other naturally occurring mixing mechanisms, however can be comparable to the rate of stratification buildup. Stratification in certain regions of shelf seas could be impacted by OWFs if these are built over a large area.

Description: Key Points: Enhanced mixing and disturbed stratification in the wake of monopiles is traceable in field and turbulence‐resolving numerical experiments Elevated turbulent dissipation and mixing are found in a narrow region downstream of monopiles, with a bulk mixing efficiency of 8–14% The enhanced mixing generated by the offshore wind farm structures could contribute to significant changes in stratification in shelf seas

Description: Helmholtz Foundation

Description: The intraseasonal variability of the tropical eastern boundary upwelling region in the Atlantic Ocean is investigated using multiyear mooring and satellite data. Pronounced oscillations of alongshore velocity and sea level off Angola at periods of about 90 and 120 days are observed. Similar spectral peaks are detected along the equator suggesting an equatorial forcing via equatorial and coastally trapped waves. Equatorial variability at 90 days is enhanced only in the eastern Atlantic likely forced by local zonal wind fluctuations. Variability at 120 days is generally stronger and linked to a second equatorial basin mode covering the whole equatorial basin. Besides forcing of the 120‐day variability by equatorial zonal winds, additional forcing of the resonant basin mode likely originates in the central and western tropical North Atlantic. The coastally trapped waves generated at the eastern boundary by the impinging equatorial Kelvin waves that are detected through their variations in sea level anomaly are associated with corresponding sea surface temperature anomalies delayed by about 14 days. Off Angola, those intraseasonal waves interfere with major coastal warm and cold events that occur every few years by either enhancing them as for the Benguela Niño in 1995 or damping them as for the warm event in 2001.

Description: Plain Language Summary: The tropical Angolan upwelling system hosts a highly productive ecosystem which plays a key socioeconomic role for societal development and fisheries in Angola. The eastern boundary circulation off Angola is dominated by the warm poleward‐flowing Angola Current. During austral summer, the Angola Current transports warm tropical waters into the Benguela upwelling system. Such a transport is often linked to extreme coastal warm events the so‐called Benguela Niños. The opposite of Benguela Niños are Benguela Niñas, both affecting the marine ecosystem and climate on multiyear time scale. At intraseasonal time scale, the Angola Current variability is dominated at periods of 90 and 120 days emanating from equatorial forcing. The 120‐day variability in the equatorial basin resembles a resonance of east‐ and westward‐propagating waves. This resonant basin mode transmits part of its energy poleward as coastally trapped waves forcing the variability along the Angolan coast and at the northern boundary of the Gulf of Guinea. The impact of these intraseasonal waves on the development of the extreme coastal warm or cold events can be shown by the relation between sea level and sea surface temperature anomalies in Southern Angola: maximum sea level is leading maximum sea surface temperature by about 14 days.

Description: Key Points: Intraseasonal variability of the Angola Current is linked to equatorial ocean dynamics and interfere with Benguela Niños and Niñas. Coastally trapped waves off Angola at 120‐day period are associated with equatorial basin‐mode resonance. Intraseasonal coastally trapped waves impact sea surface temperature off Angola and in the Gulf of Guinea via thermocline feedback.

Description: EU H2020 TRIATLAS project

Description: Bundesministerium für Bildung und Forschung (SACUS II)

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: The Total Exchange Flow analysis framework computes consistent bulk values quantifying the estuarine exchange flow using salinity coordinates since salinity is the main contributor to density in estuaries and the salinity budget is entirely controlled by the exchange flow. For deeper and larger estuaries temperature may contribute equally or even more to the density. That is why we included potential temperature as a second coordinate to the Total Exchange Flow analysis framework, which allows gaining insights in the potential temperature-salinity structure of the exchange flow as well as to compute consistent bulk potential temperature and therefore heat exchange values with the ocean. We applied this theory to the exchange flow of the Persian Gulf, a shallow, semienclosed marginal sea, where dominant evaporation leads to the formation of hypersaline and dense Gulf water. This drives an inverse estuarine circulation which is analyzed with special interest on the seasonal cycle of the exchange flow. The exchange flow of the Persian Gulf is numerically simulated with the General Estuarine Transport Model from 1993 to 2016 and validated against observations. Results show that a clear seasonal cycle exists with stronger exchange flow rates in the first half of the year. Furthermore, the composition of the outflowing water is investigated using passive tracers, which mark different surface waters. The results show that in the first half of the year, most outflowing water comes from the southern coast, while in the second half most water originates from the northwestern region.

Description: Substantial changes have occurred in the Arctic Ocean in the last decades. Not only sea ice has retreated significantly, but also the ocean at middepth showed a warming tendency. By using simulations we identified a mechanism that intensifies the upward trend in ocean heat supply to the Arctic Ocean through Fram Strait. The reduction in sea ice export through Fram Strait induced by Arctic sea ice decline increases the salinity in the Greenland Sea, which lowers the sea surface height and strengthens the cyclonic gyre circulation in the Nordic Seas. The Atlantic Water volume transport to the Nordic Seas and Arctic Ocean is consequently strengthened. This enhances the warming trend of the Arctic Atlantic Water layer, potentially contributing to the Arctic “Atlantification.” Our study suggests that the Nordic Seas can play the role of a switchyard to influence the heat budget of the Arctic Ocean.

Description: The hydrodynamics in estuaries is mainly governed by the competition between a horizontal density gradient, friction, and wind stress. The sensitivity of the estuarine exchange flow to the wind stress increases in the absence of tides, which is investigated here using the example of the weakly tidal Warnow river estuary in the southwestern Baltic Sea—the mouth of which is characterized by strongly varying salinities of 8 to 20 g kg−1. The interaction between a volatile salinity gradient and along-estuary wind forcing is found to cause temporary inversions of the estuarine circulation. Despite the highly dynamic conditions, the applicability of recent theories for isohaline mixing, using the framework of Total Exchange Flow, and the strength of the exchange flow, using a non-dimensional parameter space, could be confirmed. By analyzing salinity fluxes at the mouth of the estuary, a mixing completeness of 84% was calculated for the estuary. Furthermore, inversion of estuarine circulation was typically found for a local Wedderburn number (ratio of non-dimensional wind stress to non-dimensional horizontal density gradient) exceeding 0.33, indicating a high sensitivity to along-estuary wind.

Description: The effects of temporal changes in the marine geoid on estimates of the ocean dynamic topography are being investigated. Influences from mass redistribution due to changes of land hydrology, ice sheets, glacial isostatic adjustment (GIA), and ocean and atmospheric dynamics are considered, and the associated crustal deformation is included. The strongest signals are associated with the seasonal cycle caused by changes in terrestrial water storage and ice sheets as well as the redistribution of atmospheric mass. Second to this is the importance of an overall trend caused by GIA and decreasing ice sheets over Greenland and Antarctica. On long spatial scales, the amplitude of regional trends estimated for the geoid height has a sizable fraction of those from sea level anomaly (SLA) for the period 1993–2016, also after subtraction of steric height of the upper 1,000 m to analyze trends in deep ocean geostrophic currents. The estimated strong negative geoid height trend south of Greenland for the period 1993–2016 opposes changes in dynamic sea level for the same period thereby affecting past studies on changes of both the strength of the subpolar gyre based on SLA and the meridional overturning circulation on a section between Cape Farewell and Portugal applying ocean dynamic topography and hydrography. We conclude that temporal geoid height trends should be considered in studies of (multi)decadal trends in sea level and circulation on large spatial scales based on altimetry data referenced to a geoid field.

Description: Plain Language Summary: Changes in ocean surface currents are routinely obtained from satellite altimetry data. A correction for changes in the geoid, the equipotential surface of gravity closest to sea level, is considered small and thus usually neglected. We investigate temporal geoid height changes and potential implications on ocean circulation studies using space‐borne gravity data and results from ocean and atmosphere models to discover the individual processes of mass redistribution in the climate system causing thereby changes in the geoid height. We found the largest signals in the seasonal cycle for terrestrial hydrology in the Amazone basin and in negative trends for the Greenland and West Antarctic Ice sheets. For the period 1993–2016 and on spatial scale larger than 1,000 km or so the magnitude of the negative marine geoid height trend south of Greenland is similar to the strength of the negative trend in geocentric sea level from altimetry. This outcome affects past studies on changes in the strength of the subpolar gyre and the Atlantic meridional overturning circulation that neglect geoid height variations. We conclude that temporal geoid height trends should be considered in studies of (multi)decadal trends in sea level and circulation on large spatial scales based on altimetry data.

Description: Key Points: The strongest geoid height changes are associated with regionally pronounced seasonal signals and secular trends. In the Subpolar North Atlantic the geoid height trend biases circulation trend estimates based on altimetry. Altimetry data need correction for geoid height change when long‐term variations in ocean dynamics are studied.

Description: German Science Foundation (DFG)

Description: BMBF funded Verbundprojekt RACE (Regional Atlantic Circulation and Global Change)

Description: DFG http://dx.doi.org/10.13039/100004807

Description: Here, we revisit the existing concepts of the vertical structure of deep layers in the Black Sea using data from sensors deployed on profiling floats. The deep transition layer (DTL) between 700 and 1,700 m acts as an interface between the baroclinic layer and the largest bottom convective layer (BCL) of the world oceans. On top of DTL are the warm intermediate layer and deep cold intermediate layer. They both showed strong trends in the last 15 years due to warmer climate and intensification of warmer intrusions from Bosporus. A “salinity wave” was detected in 2005–2009 below ∼1,700 m, which evidenced for the first time the penetration of gravity flow from Bosporus down to the bottom. The layering of water masses was explained as resulting from the different distribution of sources of heat and salt, double diffusion, and balances between the geothermal and salinity flows in the BCL.

Description: Plain Language Summary: A total of 5,927 profiles from 31 profiling floats deployed in the Black Sea after 2005 were used to study the thermohaline dynamic. The density ratio, which sets the conditions under which double diffusion could occur qualifies most of the depths below 600 m as a double‐diffusive environment, which was confirmed by the observed staircase structure of vertical profiles. The qualitative differences between the temperature and salinity profiles in DTL are explained by the different distributions of sources and sinks of heat and salt. This layer acts as a transition zone between the bottom convective layer and the baroclinic ocean. A deep cold intermediate layer (DCIL) occupies the top of the DTL. The former became pronounced in recent years as a consequence of the warming of upper water layers. The large salinity variations from 2005 to 2009 observed at the bottom of the DTL evidenced traces of the penetration of the Bosporus plume down to the bottom for the first time. Combining known values of geothermal flow and concepts from the theory of dynamics of salinity‐stratified tanks heated from below made it possible to quantify the mean annual volume of Bosporus water reaching the bottom.

Description: Key Points: Profiling floats registered for the first time traces of Bosporus intrusions down to the bottom. The deep warm and cold intermediate layers are responsive to climate warming. Different sources of heat and salt dominate the transition between the baroclinic ocean and bottom convection layer.

Description: Processes of stratification and destratification in the German Bight region of fresh water influence (ROFI) are investigated following an extreme river discharge event in June 2013. For this purpose, a high‐resolution baroclinic ocean model is set up and validated against field data. The model results are used to study the temporal and spatial variability of stratification and the duration of persistent stratification in 2013. The relevant processes affecting stratification are investigated by analyzing the potential energy anomaly budget, with a focus on mixing and tidal straining. It is shown that the stratification in the German Bight is highly affected by the spring‐neap tidal cycle, with generally less stratification at spring tides due to dominant tidal mixing. It is also shown that the location of the river plume can modify this pattern. During spring tides, if the river plume is confined to the eastern region, stratification decreases significantly, as expected, due to the dominance of mixing over tidal straining. On the other hand, if the river plume moves toward deeper regions at spring tides, strong tidal straining becomes present. In this condition, mixing is weak, and the dominant tidal straining results in persistent stratification.

Description: Key Points: Processes impacting the German Bight stratification are investigated using a high‐resolution baroclinic model. The position of the river plume highly affects the contribution of tidal straining and mixing to changes in stratification. Strong tidal straining can result in persistent stratification even during spring tides.

Description: Bundesministerium für Bildung und Forschung (BMBF) http://dx.doi.org/10.13039/501100002347

Description: German Research Foundation http://dx.doi.org/10.13039/501100001659

Description: German Environment Agency http://dx.doi.org/10.13039/501100010809

Description: Ocean deoxygenation is a threat to marine ecosystems. We evaluated the potential of two ocean intervention technologies, that is, “artificial downwelling (AD)” and “artificial upwelling (AU),” for remedying the expansion of Oxygen Deficient Zones (ODZs). The model‐based assessment simulated AD and AU implementations for 80 years along the eastern Pacific ODZ. When AD was simulated by pumping surface seawater to the 178–457 m‐depth range of the ODZ, vertically integrated oxygen increased by up to 4.5% in the deployment region. Pumping water from 457 m depth to the surface (i.e., AU), where it can equilibrate with the atmosphere, increased the vertically integrated oxygen by 1.03%. However, both simulated AD and AU increased biological production via enhanced nutrient supply to the sea surface, resulting in enhanced export production and subsequent aerobic remineralization also outside of the actual implementation region, and an ultimate net decline of global oceanic oxygen.

Description: Key Points: Artificial downwelling (AD) and upwelling (AU) in the eastern Pacific oxygen deficient region are simulated in a global model Both technologies can effectively mitigate local expansion of intermediate‐depth oceanic oxygen deficient zones under climate change Global deoxygenation is however enhanced due to increased export production and aerobic respiration resulting from AD and AU

Description: Deutsche Forschungsgemeinschaft (DFG) http://dx.doi.org/10.13039/501100001659

Description: German Research Foundation (DFG)

Description: We present an application of generative adversarial networks (GANs) to reconstruct the sea level of the North Sea using a limited amount of data from tidal gauges (TGs). The application of this technique, which learns how to generate datasets with the same statistics as the training set, is explained in detail to ensure that interested scientists can implement it in similar or different oceanographic cases. Training is performed for all of 2016, and the model is validated on data from 3 months in 2017 and compared against reconstructions using the Kalman filter approach. Tests with datasets generated by an operational model (“true data”) demonstrated that using data from only 19 locations where TGs permanently operate is sufficient to generate an adequate reconstruction of the sea surface height (SSH) in the entire North Sea. The machine learning approach appeared successful when learning from different sources, which enabled us to feed the network with real observations from TGs and produce high‐quality reconstructions of the basin‐wide SSH. Individual reconstruction experiments using different combinations of training and target data during the training and validation process demonstrated similarities with data assimilation when errors in the data and model were not handled appropriately. The proposed method demonstrated good skill when analyzing both the full signal and the low‐frequency variability only. It was demonstrated that GANs are also skillful at learning and replicating processes with multiple time scales. The different skills in different areas of the North Sea are explained by the different signal‐to‐noise ratios associated with differences in regional dynamics.

Description: Plain Language Summary: The variability of sea level is one of the most important elements of the ocean dynamics. Basin‐wide observations are due to satellite altimeters, observations in coastal stations are provided by tidal gauges. The first are not very accurate in the coastal areas, the second do not provide basin‐wide coverage. The task in the present work is to use machine learning to reconstruct the sea‐level variability in the North Sea, which is an almost enclosed ocean region, using observations only. Using data from 19 coastal stations and data from numerical models as a representation of the true ocean (synthetic observations), we demonstrated that the generative adversarial networks reconstruct almost perfectly the sea level of the North Sea. The application of this technique, which learns how to generate datasets with the same statistics as the training set, is explained in detail to ensure that interested scientists can implement it in similar or different oceanographic cases.

Description: Key Points: Generative Adversarial Networks successfully reconstruct basin‐wide sea level in the North Sea using data from tidal gauges. Machine learning appeared successful when learning from different data sources. The proposed method is skillful at learning and replicating processes with multiple time scales.

Description: Over the last years, the number of studies that investigate or utilize the electromagnetic (EM) signals generated by ocean tides is steadily growing. However, the majority of these studies focuses on the amplitudes of EM tidal signals. This study investigates the phases of EM tidal signals and their changes. Twenty‐six years of monthly observation‐based datasets of tidal velocities, geomagnetic field, and oceanic conductivity are fed into an EM induction solver to generate varying EM tidal signals. The sensitivities of the resulting EM signals are analyzed by forbidding or allowing the input datasets to vary in time. We report on the phase's sensitivities with respect to changes in the EM properties, that is, secular variation of the geomagnetic field and changes in oceanic conductivity. Distinct temporal behavior and distinct geographic pattern for the two sensitivities can be reported. In general, apart from global phase shifts of 3–5 degrees, concentrated areas with phase shifts of up to 45 degrees occur all over the globe, over the oceans, for example, Arctic and Atlantic Ocean, as well as on coastal land regions, for example, Southwest Greenland and Japan. Very locally, phase shifts of 90 degree or higher occur.

Description: Key Points: Electromagnetic tidal signals show significant spatiotemporal phase changes. Annual and monthly phase anomalies are found to be of oceanic origin. Decadal transient phase anomalies are generated by secular variation and changing oceanic conductivity.

Description: This study presents recent observations to quantify oceanic heat fluxes along the continental slope of the Eurasian part of the Arctic Ocean, in order to understand the dominant processes leading to the observed along‐track heat loss of the Arctic Boundary Current (ABC). We investigate the fate of warm Atlantic Water (AW) along the Arctic Ocean continental margin of the Siberian Seas based on 11 cross‐slope conductivity, temperature, depth transects and direct heat flux estimates from microstructure profiles obtained in summer 2018. The ABC loses on average O(108) J m−2 per 100 km during its propagation along the Siberian shelves, corresponding to an average heat flux of 47 W m−2 out of the AW layer. The measured vertical heat flux on the upper AW interface of on average 10 W m−2 in the deep basin, and 3.7 W m−2 above the continental slope is larger than previously reported values. Still, these heat fluxes explain less than 20% of the observed heat loss within the boundary current. Heat fluxes are significantly increased in the turbulent near‐bottom layer, where AW intersects the continental slope, and at the lee side of a topographic irregularity. This indicates that mixing with ambient colder water along the continental margins is an important contribution to AW heat loss. Furthermore, the cold halocline layer receives approximately the same amount of heat due to upward mixing from the AW, compared to heat input from the summer‐warmed surface layer above. This underlines the importance of both surface warming and increased vertical mixing in a future ice‐free Arctic Ocean in summer.

Description: Plain Language Summary: Warm water from the Atlantic Ocean enters the Arctic Ocean through the Barents Sea and the Fram Strait, between Greenland and Norway, and directly influences the formation of sea ice: When the Atlantic Water (AW) is located close to the ocean's surface, as is the case shortly after its inflow in the Barents Sea, sea ice melts and new sea ice formation is hindered. This is why the Barents Sea is often ice free, even in winter. Further along the pathway, in the Laptev and East Siberian Sea study region, the AW gradually cools and dives down to deeper layers. In order to quantify the cooling and to understand how and where it happens, we measured vertical profiles of temperature and heat fluxes along a 2,500 km long part of the AW pathway. Based on these measurements, we found that the heat loss mainly occurs by mixing of warm AW with ambient cold water above the continental slope, in particular in the highly energetic region near the sea floor.

Description: Key Points: The Atlantic Water (AW) transported in the Arctic Boundary Current loses O(108)  J m−2 per 100 km during its translation along the Siberian shelves Heat fluxes are larger than previously reported values, but too small to account for this heat loss, indicating the importance of boundary mixing The heat input from the underlying AW layer to the cold halocline is of similar magnitude to the heat input from the warm surface layer above

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: NSF | GEO | Division of Ocean Sciences http://dx.doi.org/10.13039/100000141

Description: Semi-enclosed marginal seas like the Baltic Sea are often characterized by permanently anoxic deep layers, and may therefore serve as important model systems to study the causes and consequences of the predicted global expansion of oxygen minimum zones. Here, we focus on the role of lateral intrusions in maintaining the “hypoxic transition zone” (HTZ) of the Baltic Sea, which characterizes the quasi-permanent hypoxic region located between the oxygenized surface layer and the sulfidic deep-water region. Based on long-term deployments of an autonomous profiling system in the central Baltic Sea, we show that oxic mid-water intrusions are ubiquitous features, providing the most important oxygen source for the HTZ, and largely control the vertical and lateral extent of the hypoxic areas. An oxygen budget for the HTZ suggests that oxygen turnover in the HTZ is, to first order, determined by a long-term balance between sedimentary oxygen demand and oxygen import by intrusions. The downward mixing of oxygen into the HTZ is generally non-negligible but unlikely to provide a first-order contribution to the HTZ oxygen budget. On the long-term average, mid-water intrusions were shown to inject 30–60 Gmol of oxygen per year into the deep-water region below the permanent halocline. This is approximately one order of magnitude larger than the average amount of oxygen imported during the massive deep-water inflow events (Major Baltic Inflows) that occur on an approximately decadal time scale, highlighting the HTZ as a hotspot for biogeochemical turnover.

Description: Abstract Although the core velocity of the Atlantic North Equatorial Undercurrent (NEUC) is low (0.1−0.3 m s−1), it has been suggested to act as an important oxygen supply route towards the oxygen minimum zone in the eastern tropical North Atlantic. For the first time, the intraseasonal to interannual NEUC variability and its impact on oxygen are investigated based on shipboard and moored velocity observations around 5°N, 23°W. In contrast to previous studies that were mainly based on models or hydrographic data, we find hardly any seasonal cycle of NEUC transports in the central Atlantic. The NEUC transport variability is instead dominated by sporadic intraseasonal events. Only some of these events are associated with high oxygen levels suggesting an occasional eastward oxygen supply by NEUC transport events. Nevertheless, they likely contribute to the local oxygen maximum in the mean shipboard section along 23°W at the NEUC core position.

Description: In this study we demonstrate the potential of a kinetic energy backscatter scheme for use in global ocean simulations. Ocean models commonly employ (bi)harmonic eddy viscosities causing excessive dissipation of kinetic energy in eddy-permitting simulations. Overdissipation not only affects the smallest resolved scales but also the generation of eddies through baroclinic instabilities, impacting the entire wave number spectrum. The backscatter scheme returns part of this overdissipated energy back into the resolved flow. We employ backscatter in the FESOM2 multiresolution ocean model with a quasi-uniform 1/4° mesh. In multidecadal ocean simulations, backscatter increases eddy activity by a factor 2 or more, moving the simulation closer to observational estimates of sea surface height variability. Moreover, mean sea surface height, temperature, and salinity biases are reduced. This amounts to a globally averaged bias reduction of around 10% for each field, which is even larger in the Antarctic Circumpolar Current. However, in some regions such as the coastal Kuroshio, backscatter leads to a slight overenergizing of the flow and, in the Antarctic, to an unrealistic reduction of sea ice. Some of the bias increases can be reduced by a retuning of the model, and we suggest related adjustments to the backscatter scheme. The backscatter simulation is about 2.5 times as expensive as a simulation without backscatter. Most of the increased cost is due to a halving of the time step to accommodate higher simulated velocities.

Description: Recent evidence shows that wind-driven ocean currents, like the western boundary currents, are strongly affected by global warming. However, due to insufficient observations both on temporal and spatial scales, the impact of climate change on large-scale ocean gyres is still not clear. Here, based on satellite observations of sea surface height and sea surface temperature, we find a consistent poleward shift of the major ocean gyres. Due to strong natural variability, most of the observed ocean gyre shifts are not statistically significant, implying that natural variations may contribute to the observed trends. However, climate model simulations forced with increasing greenhouse gases suggest that the observed shift is most likely to be a response of global warming. The displacement of ocean gyres, which is coupled with the poleward shift of extratropical atmospheric circulation, has broad impacts on ocean heat transport, regional sea level rise, and coastal ocean circulation.

Description: Simulating Arctic Ocean mesoscale eddies in ocean circulation models presents a great challenge because of their small size. This study employs an unstructured-mesh ocean-sea ice model to conduct a decadal-scale global simulation with a 1-km Arctic. It provides a basinwide overview of Arctic eddy energetics. Increasing model resolution from 4 to 1 km increases Arctic eddy kinetic energy (EKE) and total kinetic energy (TKE) by about 40% and 15%, respectively. EKE is the highest along main currents over topography slopes, where strong conversion from available potential energy to EKE takes place. It is high in halocline with a maximum typically centered in the depth range of 70–110 m, and in the Atlantic Water layer of the Eurasian Basin as well. The seasonal variability of EKE along the continental slopes of southern Canada and eastern Eurasian basins is similar, stronger in fall and weaker in spring.

Description: A method using a linear shallow water model is presented for decomposing the temporal variability of the barotropic stream function in a high-resolution ocean model. The method is based on the vertically averaged momentum equations and is applied to the time series of annual mean stream function from the model configuration VIKING20 for the northern North Atlantic. An important result is the role played by the nonlinear advection terms in VIKING20 for driving transport. The method is illustrated by examining how the Gulf Stream transport in the recirculation region responds to the winter North Atlantic Oscillation (NAO). While no statistically significant response is found in the year overlapping with the winter NAO index, there is a tendency for the Gulf Stream transport to increase as the NAO becomes more positive. This becomes significant in lead years 1 and 2 when the mean flow advection and eddy momentum flux contributions, associated with nonlinear momentum advection, dominate. Only after 2 years, does the potential energy term, associated with the density field, start to play a role and it is only after 5 years that the transport dependence on the NAO ceases to be significant. It is also shown that the potential energy contribution to the transport stream function has significant memory of up to 5 years in the Labrador and Irminger Seas. However, it is only around the northern rim of these seas that VIKING20 and the transport reconstruction exhibit similar memory. This is due to masking by the mean flow advection and eddy momentum flux contributions.

Description: Phytoplankton functional-type (PFT) data are assimilated into the global coupled ocean-ecosystem model MITgcm-REcoM2 for two years using a local ensemble Kalman filter. The ecosystem model has two PFTs: small phytoplankton (SP) and diatoms. Three different sets of satellite PFT data are assimilated: Ocean-Color-Phytoplankton Functional Type (OC-PFT), Phytoplankton Differential Optical Absorption Spectroscopy (PhytoDOAS), and SYNergistic exploitation of hyper- and multi-spectral precursor SENtinel measurements to determine Phytoplankton Functional Types (SynSenPFT), which is a synergistic product combining the independent PFT products OC-PFT and PhytoDOAS. The effect of assimilating PFT data is compared with the assimilation of total chlorophyll data (TChla), which constrains both PFTs through multivariate assimilation. While the assimilation of TChla already improves both PFTs, the assimilation of PFT data further improves the representation of the phytoplankton community. The effect is particularly large for diatoms where, compared to the assimilation of TChla, the SynSenPFT assimilation results in 57% and 67% reduction of root-mean-square error and bias, respectively, while the correlation is increased from 0.45 to 0.54. For SP the assimilation of SynSenPFT data reduces the root-mean-square error and bias by 14% each and increases the correlation by 30%. The separate assimilation of the PFT data products OC-PFT, SynSenPFT, and joint assimilation of OC-PFT and PhytoDOAS data leads to similar results while the assimilation of PhytoDOAS data alone leads to deteriorated SP but improved diatoms. When both OC-PFT and PhytoDOAS data are jointly assimilated, the representation of diatoms is improved compared to the assimilation of only OC-PFT. The results show slightly lower errors than when the synergistic SynSenPFT data are assimilated, which shows that the assimilation successfully combines the separate data sources.

Description: Internal-tide generation has been quantified using both pressure work and energy conversion. When calculating the pressure work from simulated or observed data, the internal-tide pressure has to be decomposed from the full pressure, for which various options exist. We show that the conversion, that has to be derived from the depth-integrated energy equations, contains the work done by both the form drag at the bottom and that at the surface, with the latter being about 1% of the former. For calculating the pressure work, the internal-tide pressure identified as the deviation from the depth-averaged pressure perturbation has to be used. We analyzed the work done by the bottom form drag in STORMTIDE2, a concurrent simulation of circulations and tides. As expected, the identified internal-tide pressure reveals the characteristic pressure drop from the windward to the leeward side of an obstacle. The M2 internal-tide generation in STORMTIDE2 is more strongly controlled by the barotropic tide than by the topographic slope, partly because the tidal velocity can change up to one order of magnitude from the top to the foot of a high ridge within a short distance, a feature only produced by a high-resolution model. Consequently, the intense generation maps the immediate proximities of the summits of high ridges, making the global generation to be strongest near 1,200 m and decreasing drastically below 3,000 m. The depth structure of the generation differs in different basins, which could impact differently on circulations in different basins.

Description: The space-borne geodetic temporal mean dynamic topography (MDT) is obtained from the difference of altimetric mean sea surface (MSS) h and the geoid height N. With the geostrophic surface currents obtained from its gradient, the MDT is an essential parameter when describing the ocean dynamics. Spectral consistency of h and N is crucial to minimize MDT errors. Frequently, h is globalized to allow for a spherical harmonic analysis, and small scales beyond maximum degree and order (d/o) resolved in the geoid are cut off. However, common globalization causes ocean-land steps in h−N and spectral inconsistencies of N and h over land. To overcome both issues, a methodology is proposed based on globalization of the MDT. A Laplacian smoother with the coastal MDT values as boundary condition is applied, resulting in a smooth surface over land and a continuous ocean-land transition. The new methodology strongly reduces Gibbs effects and the need to work with high-resolution MDTs to minimize them. Reduction of resolution is tested to reduce MDT uncertainties caused by the commission error expected to increase with decreasing scale. Applying drifter data and a high-resolution hydrodynamic ocean model, it is shown that for the Gulf Stream and the Kuroshio, geodetic MDTs applying recent combined geoid models contain physical information up to at least d/o 420 (48-km spatial scale). For oceanic regions with strong geoid gradients, a higher-resolution MDT might be needed to prevent Gibbs effects caused by remaining inconsistencies between the geoid and the MSS.

Description: The Atlantic Subtropical Cells (STCs) are shallow wind-driven overturning circulations connecting the tropical upwelling areas to the subtropical subduction regions. In both hemispheres, they are characterized by equatorward transport at thermocline level, upwelling at the equator, and poleward Ekman transport in the surface layer. This study uses recent data from Argo floats complemented by ship sections at the western boundary as well as reanalysis products to estimate the meridional water mass transports and to investigate the vertical and horizontal structure of the STCs from an observational perspective. The seasonally varying depth of meridional velocity reversal is used as the interface between the surface poleward flow and the thermocline equatorward flow. The latter is bounded by the 26.0 kg m−3 isopycnal at depth. We find that the thermocline layer convergence is dominated by the southern hemisphere water mass transport (9.0 ± 1.1 Sv from the southern hemisphere compared to 2.9 ± 1.3 Sv from the northern hemisphere) and that this transport is mostly confined to the western boundary. Compared to the asymmetric convergence at thermocline level, the wind-driven Ekman divergence in the surface layer is more symmetric, being 20.4 ± 3.1 Sv between 10°N and 10°S. The net poleward transports (Ekman minus geostrophy) in the surface layer concur with values derived from reanalysis data (5.5 ± 0.8 Sv at 10°S and 6.4 ± 1.4 Sv at 10°N). A diapycnal transport of about 3 Sv across the 26.0 kg m−3 isopycnal is required in order to maintain the mass balance of the STC circulation.

Description: The TICON (TIdal CONstants) dataset contains harmonic constants of 40 tidal constituents computed for 1,145 tide gauges distributed globally. The tidal estimations are based on publicly available sea level records of the second version of the Global Extreme Sea Level Analysis (GESLA) project and were derived through a least squares-based harmonic analysis on the single time series. A preliminary screening was performed on all records to exclude doubtful observations. Only the records containing more than 70% of valid measurements were processed, that correspond to 89.7% of the total 1,276 original public GESLA records. The results are stored in a text file, and include additional information on the position of the stations, the starting and ending years of the analysed record, the estimated error of the fit, a code that corresponds to the source of the record and additional information on the single time series. In ocean tide models, data from in situ stations are used for validation purposes, and TICON is a useful and easy-to-handle data set that allows the users to select the records according to different criteria most suitable for their purposes. The data are provided with DOI identification in the PANGAEA repository.

Description: The Denmark Strait overflow is the major export route of dense water from the Arctic Mediterranean into the North Atlantic. At the Strait's shallow sill, the overflow is a bottom-intensified cold and dense plume, bound to the east by a thermal front formed with the warmer, northward flowing North Icelandic Irminger Current. More than two decades of observations at the sill show strong fluctuations of volume flux on daily time scales. To better understand the source of this variability, a five-mooring array was installed at the sill, capturing nearly 1 year of velocity and bottom temperature measurements at a high temporal and spatial resolution. Bottom temperature fluctuations that exceed 4 °C indicate a meandering of the front between the plume and the North Icelandic Irminger Current. Current vector rotation shows trains of alternating cyclones and anticyclones at the sill. An eddy crosses the sill every 3 to 6 days with a mean velocity of 0.4 m/s and a typical diameter of 30 to 40 km. The results suggest that anticyclones, with centers passing through the deepest part of the sill, may be responsible for periods of increased volume flux—also referred to as boluses and pulses in previous studies. Although the relationship between eddies, pulses, and boluses is still unclear, the results show that eddies are directly linked to fluctuations in the strength, thickness, and position of the overflow plume.

Description: The Filchner-Ronne Ice Shelf, fringing the southern Weddell Sea, is Antarctica's second largest ice shelf. At present, basal melt rates are low due to active dense water formation; however, model projections suggest a drastic increase in the future due to enhanced inflow of open-ocean warm water. Mooring observations from 2014 to 2016 along the eastern flank of the Filchner Trough (76°S) revealed a distinct seasonal cycle with inflow if Warm Deep Water during summer and autumn. Here we present extended time series showing an exceptionally warm and long inflow in 2017, with maximum temperatures exceeding 0.5°C. Warm temperatures persisted throughout winter, associated with a fresh anomaly, which lead to a change in stratification over the shelf, favoring an earlier inflow in the following summer. We suggest that the fresh anomaly developed upstream after anomalous summer sea ice melting and contributed to a shoaling of the shelf break thermocline.

Description: The Weddell Sea is of global importance in the formation of dense bottom waters associated with sea ice formation and ocean-ice sheet interaction occurring on the shelf areas. In this context, the Weddell Sea boundary current system (BCS) presents a major conduit for transporting relatively warm water to the Weddell Sea ice shelves and for exporting some modified form of Wedell Sea deep and bottom waters into the open ocean. This study investigates the downstream evolution of the structure and the seasonality of the BCS along the Weddell Sea continental slope, combining ocean data collected for the past two decades at three study locations. The interannual-mean geostrophic flow, which follows planetary potential vorticity contours, shifts from being surface intensified to bottom intensified along stream. The shift occurs due to the densification of water masses and the decreasing surface stress that occurs westward, toward the Antarctic Peninsula. A coherent along-slope seasonal acceleration of the barotropic flow exists, with maximum speed in austral autumn and minimum speed in austral summer. The barotropic flow significantly contributes to the seasonal variability in bottom velocity along the tip of the Antarctic Peninsula. Our analysis suggests that the winds on the eastern/northeastern side of the gyre determines the seasonal acceleration of the barotropic flow. In turn, they might control the export of Weddell Sea Bottom Water on seasonal time scales. The processes controlling the baroclinic seasonality of the flow need further investigation.

Description: We analyzed velocity and hydrographic data from 23 moorings in the northeast Chukchi Sea from 2011 to 2014. In most years the eastern side of Hanna Shoal was strongly stratified year-round, while weakly stratified regions prevailed on the shelf south and west of the Shoal. Stratification differences cause differential vertical mixing rates, which in conjunction with advection of different bottom water properties resulted in seasonally varying along-isobath density gradients. In agreement with numerical models, we find that bottom waters flow anticyclonically around the Shoal. Whereas most of the shelf responded barotropically to wind-forcing, there was a strong baroclinic component to the flow field northeast of Hanna Shoal, resulting in no net vertically integrated transport on average. In contrast there is a net eastward transport from west of the Shoal, which implies convergence north of the Shoal. Convergence and along-isobath density gradients may foster cross-shelf exchange north of Hanna Shoal. Modal analyses indicate that the shelf south of the Shoal and Barrow Canyon responded coherently to local and remote winds, whereas the wind-current response around Hanna Shoal was less coherent. Barotropic topographic waves, of ~3-day period, were generated episodically northeast of the Shoal and propagate clockwise around Hanna Shoal, but are blocked from entering Barrow Canyon and are possibly scattered by the horizontally sheared flow and converging isobaths on the western side of the Shoal. Analysis of water properties on the western side of Hanna Shoal suggests that these include contributions from the western and southern portions of the Chukchi Sea.

Description: Ocean bottom pressure (OBP) variability serves as a proxy of ocean mass variability, the knowledge of which is needed in geophysical applications. The question of how well it can be modeled by the present general ocean circulation models on time scales in excess of 1 day is addressed here by comparing the simulated OBP variability with the observed one. To this end, a new multiyear data set is used, obtained with an array of bottom pressure gauges deployed deeply along a transect across the Southern Ocean. We present a brief description of OBP data and show large-scale correlations over several thousand kilometers at all time scales using daily and monthly averaged data. Annual and semiannual cycles are weak. Close to the Agulhas Retroflection, signals of up to 30 cm equivalent water height are detected. Further south, signals are mostly intermittent and noisy. It is shown that the models simulate consistent patterns of bottom pressure variability on monthly and longer scales except for areas with high mesoscale eddy activity, where high resolution is needed to capture the variability due to eddies. Furthermore, despite good agreement in the amplitude of variability, the in situ and simulated OBP show only modest correlation.

Description: The influence of large-scale oceanic circulation on salinity in the northern North Sea has lead to the hypothesis that nutrient concentrations in this region are also driven by remote oceanic anomalies. Here, using a newly established biogeochemical data set of the North Sea, we show that interannual to decadal variability in winter nutrient concentrations exhibits distinct phase deviations from salinity. The variability in salinity is explained by zonal shifts in the position of the subpolar front (SPF) in the eastern North Atlantic and the associated advective delay. However, the high correlation and absence of advective delay between the position of the SPF and winter nutrient concentrations in the Shetland region (59–61°N, 1°W to 3°E) point to the role of atmospheric variability in driving concurrent changes in winter nutrient concentrations and the SPF position. Our analysis suggests that the prevailing wind direction and local distribution of winter nutrient concentrations together determine the interannual to decadal variability in winter nutrient concentrations in this region. In the analyzed observations, we find a strong spatial gradient in mean winter nutrient concentrations northwest of the Shetland region, which is absent in salinity. The horizontal shift of this spatial gradient, forced by changes in wind direction, has a larger influence on winter nutrient concentration in the Shetland region than the nutrient signal in oceanic anomalies originating from the eastern subpolar North Atlantic. Overall, we conclude that interannual to decadal variability in the observed nutrient concentrations is mainly driven by atmospheric variability here expressed as wind direction.

Description: Deep convection over tropical oceans often appears intensified at the edge of convectively active regions, both in idealized studies and in observations. This edge intensification of convection is studied in detail here, using the steady state of a radiative-convective equilibrium study, marked by a single convective cluster with deep convection intensified at the edge of this cluster. The cause for edge intensification and its dependence on the cluster area is investigated by comparing the spatial distribution of deep convection to different variables known to be important for convection. Analysis of the simulation suggests that the edge is marked by an increased probability for the triggering of convection rather than by stronger updrafts. In particular, while the edge of the moist region is not thermodynamically more favorable, we find strong surface convergence and therefore dynamical lifting at this edge. The surface convergence is shown to result from two opposing flows. On the one hand, there is, as expected from previous radiative-convective equilibrium simulations, a low-level inflow directed toward the moist region. On the other hand, there is a positive density anomaly at the surface which is the result of continuously forming cold pools within the convectively active region, creating a super-cold-pool. As the velocity of the low-level inflow approximately matches the potential propagation speed of the super-cold-pool boundary, these opposing flows explain the presence of strong convergence at the edge of this region. Whether the resulting lifting induces the formation of deep convection is shown to depend on the large-scale instability.

Description: We investigate the relative importance of ecosystem complexity and phytoplankton light absorption for climate studies. While the complexity of Earth System models (ESMs) with respect to marine biota has increased over the past years, the relative importance of biological processes in driving climate‐relevant mechanisms such as the biological carbon pump and phytoplankton light absorption is still unknown. The climate effects of these mechanisms have been studied separately, but not together. To shed light on the role of biologically mediated feedbacks, we performed different model experiments with the EcoGENIE ESM. The model experiments have been conducted with and without phytoplankton light absorption and with two or 12 plankton functional types. For a robust comparison, all simulations are tuned to have the same primary production. Our model experiments show that phytoplankton light absorption changes ocean physics and biogeochemistry. Higher sea surface temperature decreases the solubility of CO2 which in turn increases the atmospheric CO2 concentration, and finally the atmospheric temperature rises by 0.45°C. An increase in ecosystem complexity increases the export production of particulate organic carbon but decreases the amount of dissolved organic matter. These changes in the marine carbon cycling, however, hardly reduces the atmospheric CO2 concentrations and slightly decreases the atmospheric temperature by 0.034°C. Overall we show that phytoplankton light absorption has a higher impact on the carbon cycle and on the climate system than a more detailed representation of the marine biota.

Description: Plain Language Summary: Climate models are used to evaluate how climate will continue to change on timescales relevant for human beings. However, they are complex and need to be simplified to run them even on supercomputers. Usually physical processes are very well resolved while biological processes gain less attention. We focus on biological processes and investigate which of them turn out to be particularly important for the climate system so that they should be included in climate models. We are aware that phytoplankton changes the optical properties, redistributing the heat in the ocean. We have also known for a long time that the marine organisms remove carbon from the system by CO2‐uptake and downward transport of organically bound carbon by migration or sinking. Yet, the relative importance of these two processes are unknown. Our study suggests that phytoplankton light absorption has a higher impact on the climate system than increasing ecosystem complexity. Phytoplankton light absorption should thus be included by default in models carrying out climate change scenarios.

Description: Key Points: We investigate the relative importance of ecosystem complexity and phytoplankton light absorption with the Earth System model (ESM) Eco Grid‐ENabled Integrated ESM (EcoGENIE). Phytoplankton light absorption increases the sea surface temperature (SST) while a higher ecosystem complexity leads to a slightly deeper downward flux of organic matter. Our results suggest that phytoplankton light absorption has a higher impact on the climate system than increasing ecosystem complexity.

Description: Universität Hamburg http://dx.doi.org/10.13039/501100005711

Description: The Filchner‐Ronne Ice Shelf (FRIS) is characterized by moderate basal melt rates due to the near‐freezing waters that dominate the wide southern Weddell Sea continental shelf. We revisited the region in austral summer 2018 with detailed hydrographic and noble gas surveys along FRIS. The FRIS front was characterized by High Salinity Shelf Water (HSSW) in Ronne Depression, Ice Shelf Water (ISW) on its eastern flank, and an inflow of modified Warm Deep Water (mWDW) entering through Central Trough. Filchner Trough was dominated by Ronne HSSW‐sourced ISW, likely forced by a recently intensified circulation beneath FRIS due to enhanced sea ice production in the Ronne polynya since 2015. Glacial meltwater fractions and tracer‐based water mass dating indicate two separate ISW outflow cores, one hugging the Berkner slope after a two‐year travel time, and the other located in the central Filchner Trough following a ∼six year‐long transit through the FRIS cavity. Historical measurements indicate the presence of two distinct modes, in which water masses in Filchner Trough were dominated by either Ronne HSSW‐derived ISW (Ronne‐mode) or more locally derived Berkner‐HSSW (Berkner‐mode). While the dominance of these modes has alternated on interannual time scales, ocean densities in Filchner Trough have remained remarkably stable since the first surveys in 1980. Indeed, geostrophic velocities indicated outflowing ISW‐cores along the trough's western flank and onto Berkner Bank, which suggests that Ronne‐ISW preconditions Berkner‐HSSW production. The negligible density difference between Berkner‐ and Ronne‐mode waters indicates that each contributes cold dense shelf waters to protect FRIS against inflowing mWDW.

Description: Plain Language Summary: We visited the largest floating Antarctic ice shelf in the southern Weddell Sea in 2018 with an icebreaker expedition, and measured ocean temperature, salinity, meltwater content, and other parameters in front of the FRIS. We found that the ocean conditions were still dominated by the very cold and dense waters needed to protect the ice shelf from inflowing warm waters from the deep ocean. We compared the 2018 conditions with earlier surveys since the 1980s and concluded that, in spite of climate change and in contrast to other Antarctic regions, the water masses on the southern Weddell Sea shelf remained relatively stable overall. We found that most of the stations we visited near the Filchner Ice Shelf edge were dominated by cold ISW, which forms when water masses interact with the underside of the shelf ice. Our measurements helped improve our understanding regarding the currents and water masses on the southern Weddell Sea continental shelf.

Description: Key Points: Hydrographic status update with the first comprehensive CTD survey along the entire FRIS front since 1995. Strong and stable presence of High Salinity Shelf Water in Ronne Depression over decades. Dominance of Ronne‐sourced Ice Shelf Water in Filchner Trough in 2018 points to intensified sub‐FRIS circulation.

Description: Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI) http://dx.doi.org/10.13039/501100003207

Description: Atlantic Water (AW), which is transported northward by the West Spitsbergen Current (WSC), partly recirculates (i.e., turns westward) in Fram Strait. This determines how much heat and salt reaches the Arctic Ocean, and how much joins the East Greenland Current on its southward path. We describe the AW recirculation's location, seasonality, and mesoscale variability by analyzing the first observations from moored instruments at five latitudes in central Fram Strait, spanning a period from August 2016 to July 2018. We observe recirculation on the prime meridian at 78°50’N and 80°10’N, respectively south and north of the Molloy Hole, and no recirculation further south at 78°10’N and further north at 80°50’N. At a fifth mooring location northwest of the Molloy Hole at 79°30’N, we observe some influence of the two recirculation pathways. The southern recirculation is observed as a continuous westward flow that carries AW throughout the year, though it may be subject to broadening and narrowing. It is affected by eddies in spring, likely due to the seasonality of mesoscale instability in the WSC. The northern recirculation is observed solely as passing eddies on the prime meridian, which are strongest during late autumn and winter, and absent during summer. This seasonality is likely affected both by the conditions set by the WSC and by the sea ice. Open ocean eddies originating from the WSC interact with the sea ice edge when they subduct below the fresher, colder water.

Description: Plain Language Summary: West of Spitsbergen, the so‐called West Spitsbergen Current (WSC) carries relatively warm water towards the Arctic Ocean. Part of this water turns westward before it reaches the Arctic Ocean. Another current transports it back south along the Greenland shelf. We look at observations of temperature, salinity, and flow velocity from two years in the region. We find that the warm water flows westwards at two latitudes. At 78°50’N it flows westwards throughout the year, but is affected by eddies in the ocean during spring. These eddies also carry warm water westwards and likely originate from the WSC. At 80°10’N, the warm water passes only over periods of a few days, and only during late autumn, winter, and spring. That means that we do not observe the warm water flowing westwards at this latitude, but only eddies that carry it westwards. These eddies also originate from the WSC, and their properties likely change, when they meet the ice edge and have to push underneath.

Description: Key Points: We present the first multi‐year mooring‐based observations of the Atlantic Water recirculation at the prime meridian. At 78°50’N we observe a continuous westward flow, present throughout the year and most affected by eddies in spring. At 80°10′N we only observe passing eddies, which are strongest during late autumn and winter, and absent during summer.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Helmholtz Infrastructure Initiative FRAM

Description: The efficiency of sediment routing from land to the ocean depends on the position of submarine canyon heads with regard to terrestrial sediment sources. We aim to identify the main controls on whether a submarine canyon head remains connected to terrestrial sediment input during Holocene sea‐level rise. Globally, we identified 798 canyon heads that are currently located at the 120m‐depth contour (the Last Glacial Maximum shoreline) and 183 canyon heads that are connected to the shore (within a distance of 6 km) during the present‐day highstand. Regional hotspots of shore‐connected canyons are the Mediterranean active margin and the Pacific coast of Central and South America. We used 34 terrestrial and marine predictor variables to predict shore‐connected canyon occurrence using Bayesian regression. Our analysis shows that steep and narrow shelves facilitate canyon‐head connectivity to the shore. Moreover, shore‐connected canyons occur preferentially along active margins characterized by resistant bedrock and high river‐water discharge.

Description: Plain Language Summary: Since the last glaciation about 20,000 years ago, sea level has risen by about 120 m. As a consequence, most coastlines have migrated landward, inundating large shelf areas. Some of these areas are now dissected by submarine canyons. However, with only 4% of the world’s submarine canyons reaching today's coastline, these canyons remain the exception. Here, we aim to identify the environmental factors and processes that control whether rates of headward canyon incision can keep pace with landward migration of the coastline during the Holocene. We determine 34 variables that potentially predict whether a canyon remains connected to the coastline. We find that shore‐connected canyons preferentially occur along continental margins with narrow and steep shelves, such as the Mediterranean active margin and the Pacific coast of Central and South America. Moreover, our analysis supports the occurrence of such canyons offshore river basins, that are characterized by resistant bedrock and high water discharge. Such rivers deliver coarse‐grained sediment to submarine canyons, which can erode the canyon head and floor. To this end, our analysis offers new insights into the formation and maintenance of submarine canyons that are required to efficiently transport sediments, pollutants, and organic carbon from rivers to the ocean floor.

Description: Key Points: Presently, 183 submarine canyons are connected to the shoreline (within 6 km) along the world’s major continents. Narrow shelves and high shelf gradients precondition the maintenance of canyon‐head connectivity to the shore. Canyon heads preferentially remain connected to the shore offshore river catchments with resistant bedrock and high water discharge.

Description: Mechanisms controlling the heat budget of the North Sea are investigated based on a combination of satellite sea surface temperature measurements and numerical model simulations. Lateral heat fluxes across the shelf edge and into the Baltic Sea as well as vertical ocean‐atmosphere heat exchange are considered. A 3‐D variational (3DVAR) data assimilation (DA) scheme is applied, which contains assumed model error correlations that depend on the mixed layer depth derived from a coupled circulation/ocean wave model. The analysis balances pressure gradients introduced by temperature modifications. Significant hydrodynamic model response to DA was found, which should be considered in the heat budget estimations. The observed change of the current velocity field decreases the lateral advective volume/heat exchanges between the North Sea and the Atlantic, yielding an increased heat flux from the Atlantic into the North Sea and more heat flux from the sea to the atmosphere. The largest DA impact on volume/heat transport is in the Norwegian Channel, where the dominant process is Eulerian transport, followed by tidal pumping and wind pumping. Further analysis reveals an acceleration of the along‐shelf current at the northern edge of the North Sea, a decrease in the horizontal pressure gradient from the Atlantic to the North Sea, and a reduction of the Eulerian transport of volume/heat outward the North Sea. Furthermore, the coupling between the circulation model and the wave model has significant impacts on lateral heat advection in the DA run, which is due to the wave impact on the mixed layer depth.

Description: Plain Language Summary: Seawater temperature simulations are important for climate change research, fishery management, coastline protection, ecological balance maintenance, and weather predictions. To improve the seawater temperature prediction capability, a data assimilation (DA) scheme is often applied to combine data from measurements, such as from satellites, buoys, and ships, with data provided by climate models that consider circulation, wave, atmosphere, and ice components. For decades, various DA methods have been developed with a focus on implementing sophisticated mathematical techniques. However, few studies have focused on the impacts of DA on physical processes and the secondary effects of DA. We used a model and satellite data to investigate the impacts of sea surface temperature (SST) assimilations on the volume and heat budgets over the North Sea. We found that DA improved SST modeling, thereby modifying the volume and heat budgets between the North Sea and the Atlantic. The largest change occurs at the Norwegian Channel, where the total water/heat transport from the North Sea outward is reduced. Moreover, SST assimilation also changes the air‐sea heat exchange. This study improves our understanding of the relations between model physics and DA.

Description: Key Points: Even with balanced pressure gradients there is a significant response of the three‐dimensional model circulation to sea surface temperature (SST) assimilation in the North Sea. Volume transports through the Norwegian Channel are modified by SST assimilation leading to changes in lateral heat advection. Mechanisms controlling the heat budget of the North Sea are different at west, north, and east boundaries.

Description: the German Helmholtz Association (HGF)

Description: the Federal Ministry of Education and Research (BMBF)

Description: Helmholtz Association http://dx.doi.org/10.13039/501100009318

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Microbial degradation of dissolved organic matter (DOM) contributes to the formation and preservation of oxygen minimum zones (OMZs) in the ocean, but information on the spatial distribution and molecular composition of DOM in OMZ regions is scarce. We quantified molecular components of DOM that is, dissolved amino acids (DAA) and dissolved combined carbohydrates (DCCHO), in the upwelling region off Peru. We found the highest concentrations of DCCHO in fully oxygenated surface waters steeply declining at shallow depth. The highest DAA concentrations were observed close to the surface also, but attenuation of DAA concentration over depth was less pronounced. Compositional changes of DCCHO were strongest within more oxygenated waters. Compositional changes of DAA were also evident under suboxic conditions (〈5 µmol O2 kg−1) and indicated bacterial peptide degradation. Moreover, specific free amino acids (alanine and threonine) were enhanced within suboxic waters, pointing to a potential production of dissolved organic nitrogen under suboxic conditions. Our results therewith suggest that deoxygenation supports a spatial decoupling of DCCHO and DAA production and degradation dynamics and give new insights to carbon and nitrogen cycling in the OMZ off Peru.

Description: Plain Language Summary: The molecular composition of dissolved organic matter (DOM) in seawater mainly depends on production and consumption processes by plankton organisms. The biological consumption of DOM reduces the amount of oxygen in seawater, but little is known about the composition of DOM in so‐called oxygen minimum zones (OMZs). We studied the distribution of DOM, specifically of dissolved sugars and amino acids in the OMZ off Peru. We observed that the reworking and degradation of sugars and amino acids differed between more and less oxygenated waters. We also found indications for the production of specific amino acids within the OMZ. Our study therewith suggests, that oxygen concentrations affect dissolved amino acid and sugar dynamics and may help to better understand carbon and nitrogen cycling within OMZs.

Description: Key Points: Decoupling of dissolved combined carbohydrates and dissolved amino acids dynamics, in the oxygen minimum zone off Peru. Free amino acids indicate production of dissolved organic nitrogen within suboxic waters.

Description: Hemholtz Association

Description: German Science Foundation

Description: Climate change is expected to enhance the hydrological cycle in northern latitudes reducing the salinity in the Baltic Sea, a land-locked marginal sea with a large catchment area located in northern Europe. With the help of ocean simulations forced by historical atmospheric and hydrological reconstructions and local observations, we analyzed long-term changes in the sea surface salinity of the Baltic Sea as well as its latitudinal gradient. The variability of both is dominated by multidecadal oscillations with a period of about 30 years, while both atmospheric variables, wind and river runoff, contribute to this variability. Centennial changes show a statistically significant positive trend in the North-South gradient of sea surface salinity for 1900–2008. This change is mainly attributed to increased river runoff from the northernmost catchment indicating a footprint of the anthropogenic impact on salinity with consequences for the marine ecosystem and species distributions.

Description: An approach to analyze high-end sea level rise is presented to provide a conceptual framework for high-end estimates as a function of time scale, thereby linking robust sea level science with stakeholder needs. Instead of developing and agreeing on a set of high-end sea level rise numbers or using an expert consultation, our effort is focused on the essential task of providing a generic conceptual framework for such discussions and demonstrating its feasibility to address this problem. In contrast, information about high-end sea level rise projections was derived previously either from a likely range emerging from the highest view of emissions in the Intergovernmental Panel on Climate Change assessment (currently the Representative Concentration Pathway 8.5 scenario) or from independent ad hoc studies and expert solicitations. Ideally, users need high-end sea level information representing the upper tail of a single joint sea level frequency distribution, which considers all plausible yet unknown emission scenarios as well as involved physical mechanisms and natural variability of sea level, but this is not possible. In the absence of such information we propose a framework that would infer the required information from explicit conditional statements (lines of evidence) in combination with upper (plausible) physical bounds. This approach acknowledges the growing uncertainty in respective estimates with increasing time scale. It also allows consideration of the various levels of risk aversion of the diverse stakeholders who make coastal policy and adaptation decisions, while maintaining scientific rigor.

Description: C-band microwave sea surface radar backscatter observations from the FINO-2 tower in the western Baltic are analyzed with respect to their relevance for air-sea CO2 transfer velocity parameterizations. The scatterometer measurements observed from a height of 25 m above the sea surface using a multifrequency scatterometer instrument of the University of Hamburg were obtained quasi-simultaneously with eddy covariance CO2 flux measurements. Both data sets are merged here to derive a gas transfer velocity parameterization based on radar backscatter measurements. At the location of the FINO-2 tower, the resulting time-averaged gas transfer velocity amounts to 26.95 cm/hr. In combination with ΔPCO2 measurements available from the vicinity of the FINO-2 platform, a time-mean CO2 flux of 0.23 μmol·m−2·s−1 into the Baltic was estimated. Applied to monthly mean satellite-based C-band ASCAT scatterometer data, the newly derived gas transfer velocity parameterization provides estimates of seasonal and annual mean global maps of air-sea transfer velocities. The new results agree in their general pattern with previous estimates using wind speed parameterization. However, the backscatter-based transfer velocities appear smaller at higher latitudes. Globally averaged air-sea CO2 fluxes would thereby be reduced by 20%. To what extent this is a robust result, or if it depends on the fact that the training data set did not represent conditions, has to be investigated in the future.

Description: A regional seismic survey on the southeastern Lomonosov Ridge (LR) and adjacent basins provides constraints on the coupled evolution of ocean circulations, depositional regime, and tectonic processes. First, Mesozoic strata on the LR, its faulted flanks and the initial Amundsen Basin were covered with syn-rift sediments of Paleocene to early Eocene age. Numerous vertical faults indicate differential compaction of possibly anoxic sediments deposited in the young, still isolated Eurasian Basin. The second stage, as indicated by a prominent high-amplitude-reflector sequence covering the ridge, was a time of widespread changes in deposition conditions, likely controlled by the ongoing subsidence of the LR and gradual opening of the Fram Strait. Episodic incursions of water masses from the North Atlantic probably were the consequences and led to the deposition of thin sedimentary layers of different lithology. The third stage is marked by continuous deposition since the early Miocene (20 Ma). At that time, the ridge no longer posed an obstacle between the Amerasia and Eurasia Basins and pelagic sedimentation was established. Drift bodies, sediment waves, and erosional structures indicate the onset of circulation. Faulting on the ridge slope has led to a series of terraces where sediment drifts have accumulated since the early Miocene. It is suggested that ongoing sagging of the ridge and currents may have shaped the steep sediment free flanks of the terraces. Lastly, a sequence of high-amplitude reflectors marks the transition to the early Pliocene large-scale Northern Hemisphere glaciations.

Description: The paper describes an update of the GECCO (German contribution to the Estimating the Circulation and Climate of the Ocean project) ocean synthesis, now in its version 3, and provides an evaluation of the results with assimilated and independent data. GECCO3 covers the 71-year period 1948–2018 and differs from its predecessor by returning to a single assimilation window instead of partitioning the period in 5-year-long overlapping windows which was previously necessary to yield convergence. A solution to the convergence problem is presented. GECCO3 is intended to be used for the initialization of coupled climate models and is configured for the higher-resolution version of the earth system model (MPI-ESM) developed at the Max Planck Institute for Meteorology. It uses the bathymetry and grid of the MPI-ESM with quasi-uniform resolution of 0.4°, thereby providing the first global eddy-permitting synthesis based on the adjoint method. The synthesis additionally features the estimation of various mixing parameters and can regionally choose between explicit or parametrized eddy fluxes. Except for the altimeter data in tropical regions, GECCO3 is in better agreement with the assimilated data than GECCO2. The improvements relative to the in situ data partly result from the much larger amount of Argo data, which show lower model–data differences. Global heat content changes are in good agreement with recent estimates, but show uptake almost exclusively in the top 700 m. An alternative version of GECCO3, created by starting from different first-guess control parameters, was used to evaluate the uncertainty of the estimated parameters and state due to lack of convergence. This estimate suggests a large uncertainty related to the uptake of heat into the lower layers, while estimates of mean meridional transport of heat and freshwater are not affected.

Description: Low-mode internal waves propagate over large distances and provide energy for turbulent mixing when they break far from their generation sites. A realistic representation of the oceanic energy cycle in ocean and climate models requires a consistent implementation of their generation, propagation, and dissipation. Here we combine the long-term mean energy flux from satellite altimetry with results from a 1/10° global ocean general circulation model that resolves the low modes of internal waves and in situ observations of stratification and horizontal currents to study energy flux and dissipation along a 1000 km internal tide beam in the eastern North Atlantic. Internal wave fluxes were estimated from twelve 36- to 48-hr stations in along- and across-beam direction to resolve both the inertial period and tidal cycle. The observed internal tide energy fluxes range from 5.9 kW m−1 near the generation sites to 0.5 kW m−1 at distant stations. Estimates of energy dissipation come from both finestructure and upper ocean microstructure profiles and range, vertically integrated, from 0.5 to 3.3 mW m−2 along the beam. Overall, the in situ observations confirm the internal tide pattern derived from satellite altimetry, but the in situ energy fluxes are more variable and decrease less monotonically along the beam. Internal tides in the model propagate over shorter distances compared to results from altimetry and in situ measurements, but more spatial details close the main generation sites are resolved.

Description: We examine the relative dispersion and the contribution of tides on the relative diffusivities of surface drifters in the North Sea. The drifters are released in two clusters, yielding 43 pairs, in the vicinity of a tidal mixing front in the German Bight, which is located in the southeastern area of the North Sea. Both clusters indicate decreasing dispersion when crossing the tidal mixing front, followed by exponentially increasing dispersion with e-folding times of 0.5 days for Cluster 1 and 0.3 days for Cluster 2. A transition of the dispersion regimes is observed at scales of the order of the Rossby radius of deformation (10 km). After that, the relative dispersion grows with a power-law dependency with a short period of ballistic dispersion (quadratic growth), followed by a Richardson regime (cubic growth) in the final phase. Scale-dependent metrics such as the relative diffusivities are consistent with these findings, while the analysis of the finite-scale Lyapunov exponents (FSLEs) shows contradictory results for the submesoscales. In summary, the analysis of various statistical Lagrangian metrics suggests that tracer stirring at the submesoscales is nonlocal and becomes local at separation scales larger than 10 km. The analysis of meridional and zonal dispersion components indicates anisotropic dispersion at the submesoscales, which changes into isotropic dispersion on the mesoscales. Spectral analysis of the relative diffusivity gives evidence that semidiurnal and shallow-water tides influence relative diffusivity at the mesoscales, especially for drifter separations above 50 km.

Description: Storms are infrequent, intense, physical forcing events that represent a potentially significant driver of ocean ecosystems. The objective of this study was to assess changes in water column structure and turbulent fluxes caused by storms using an autonomous underwater glider, as well as the chlorophyll a (Chl a) response to the altered physical environment. The glider was able to measure throughout the complete life cycle of Storm Bertha as it passed over the North Sea in August 2014, from its arrival to dissipation. Storm Bertha triggered rapid mixing of the thermocline through shear instability, increasing vertical fluxes by nearly an order of magnitude, and promoting increases in surface layer Chl a. The results demonstrate that storms represent a significant fraction of seasonal vertical turbulent fluxes, with potentially important consequences for biological production in shelf seas.

Description: Ocean heat transport is often thought to play a secondary role for Arctic surface warming in part because warm water which flows northward is prevented from reaching the surface by a cold and stable halocline layer. However, recent observations in various regions indicate that occasionally, warm water is found directly below the surface mixed layer. Here we investigate Arctic Ocean surface energy fluxes and the cold halocline layer in climate model simulations from the Coupled Model Intercomparison Project Phase 5. An ensemble of 15 models shows decreased sea ice formation and increased ocean energy release during fall, winter, and spring for a high-emission future scenario. Along the main pathways for warm water advection, this increased energy release is not locally balanced by increased Arctic Ocean energy uptake in summer. Because during Arctic winter, the ocean mixed layer is mainly heated from below, we analyze changes of the cold halocline layer in the monthly mean Coupled Model Intercomparison Project Phase 5 data. Fresh water acts to stabilize the upper ocean as expected based on previous studies. We find that in spite of this stabilizing effect, periods in which warm water is found directly or almost directly below the mixed layer and which occur mainly in winter and spring become more frequent in high-emission future scenario simulations, especially along the main pathways for warm water advection. This could reduce sea ice formation and surface albedo.

Description: Changes in ocean gateway configuration can induce basin-scale rearrangements in ocean current characteristics. However, there is large uncertainty in the relative timing of the Oligocene/Miocene subsidence histories of the Greenland-Scotland Ridge (GSR) and the Fram Strait (FS). By using a climate model, we investigate the temperature and salinity changes in response to the subsidence of these two key ocean gateways during early to middle Miocene. For a singular subsidence of the GSR, we detect warming and a salinity increase in the Nordic Seas and the Arctic Ocean. As convection sites shift to the north of Iceland, North Atlantic Deep Water (NADW) is formed at cooler temperatures. The associated deep ocean cooling and upwelling of deep waters to the Southern Ocean surface can cause a cooling in the southern high latitudes. These characteristic responses to the GSR deepening are independent of the FS being shallow or deep. An isolated subsidence of the FS gateway for a deep GSR shows less pronounced warming and salinity increase in the Nordic Seas. Arctic temperatures remain unaltered, but a stronger salinity increase is detected, which further increases the density of NADW. The increase in salinity enhances the contribution of NADW to the abyssal ocean at the expense of the colder southern source water component. These relative changes largely counteract each other and cause a negligible warming in the upwelling regions of the Southern Ocean.

Description: Various alternative formulations of the LES equations have been explored in which additional evolution equations for variables such as the acceleration, the subgrid-scale stress tensor, or the subgrid-scale force are explicitly carried. Statistics of the velocity field obtained from the equation for the acceleration are shown to depend strongly on the initial conditions. This feature, which is independent of LES modeling issues, seems to prove that the velocity-acceleration formulation of the Navier-Stokes is not useful for numerical simulation. Equations for the subgrid-scale quantities appear to be much more stable. However, models required by this formulation of the LES problem still require additional study.

Description: Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It is found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It is found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.

Description: The bubble coalescence and interfacial instabilities that are important to modeling critical heat flux (CHF) in reduced-gravity systems can be sensitive to even minute body forces. Understanding these complex phenomena is vital to the design and safe implementation of two-phase thermal management loops proposed for space and planetary-based thermal systems. While reduced gravity conditions cannot be accurately simulated in 1g ground-based experiments, such experiments can help isolate the effects of the various forces (body force, surface tension force and inertia) which influence flow boiling CHF. In this project, the effects of the component of body force perpendicular to a heated wall were examined by conducting 1g flow boiling experiments at different orientations. FC-72 liquid was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface at conditions approaching CHF. High-speed video imaging was employed to capture dominant CHF mechanisms. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed great sensitivity to orientation for flow velocities below 0.2 m/s, where very small CHF values where measured, especially with downflow and downward-facing heated wall orientations. High flow velocities dampened the effects of orientation considerably. Figure I shows representative images for the different CHF regimes. The Wavy Vapor Layer regime was dominant for all high velocities and most orientations, while all other regimes were encountered at low velocities, in the downflow and/or downward-facing heated wall orientations. The Interfacial Lift-off model was modified to predict the effects of orientation on CHF for the dominant Wavy Vapor Layer regime. The photographic study captured a fairly continuous wavy vapor layer travelling along the heated wall while permitting liquid contact only in wetting fronts, located in the troughs of the interfacial waves. CHF commenced when wetting fronts near the outlet were lifted off the wall. The Interfacial Lift-off model is shown to be an effective tool for predicting the effects of body force on CHF at high velocities.

Description: Probabilistic CFD design is needed because we are asked to do more with less. To cost effectively accomplish the design task, we need to formally quantify the effect of uncertainties (variables) in the design. Probabilistic design is one effective method to formally quantify the effect of uncertainties. Our objective is to establish a revolutionary new early design process, by developing non-deterministic physics-based probabilistic design tools, which will include all the life cycle processes. Breakthroughs will be sought in speed, accuracy, intelligence, and usability of the system. This paper is concerned with the usefulness of parametric optimization method coupled with a Navier-Stokes analysis code for the aero-thermodynamic design of turbomachinery combustor liner. The interconnection between the CFD code and NESSUS codes facilitated the coupling between the thermal profiles and structural design. We have developed new concepts for reducing the computational cost of unsteady, three-dimensional, compressible aerodynamic analyses for multistage turbomachinery flows. The flow was modeled by the three-dimensional Favre-Reynolds-averaged Navier-Stokes equations using the k-epsilon turbulence closure, which was integrated using an implicit third-order upwind solver. The methodology developed in this paper is expected to lead to the design optimization of turbomachinery blades.

Description: Parallel Plate Plastometer (PPP) is a device commonly used for measuring the viscosity of high polymers at low rates of shear in the range 10(exp 4) to 10(exp 9) poises. This device is being validated for use in measuring the viscosity of liquid glasses at high temperatures having similar ranges for the viscosity values. PPP instrument consists of two similar parallel plates, both in the range of 1 inch in diameter with the upper plate being movable while the lower one is kept stationary. Load is applied to the upper plate by means of a beam connected to shaft attached to the upper plate. The viscosity of the fluid is deduced from measuring the variation of the plate separation, h, as a function of time when a specified fixed load is applied on the beam. Operating plate speeds measured with the PPP is usually in the range of 10.3 cm/s or lower. The flow field within the PPP can be simulated using the equations of motion of fluid flow for this configuration. With flow speeds in the range quoted above the flow field between the two plates is certainly incompressible and laminar. Such flows can be easily simulated using numerical modeling with computational fluid dynamics (CFD) codes. We present below the mathematical model used to simulate this flow field and also the solutions obtained for the flow using a commercially available finite element CFD code.

Description: The Physics of Colloids in Space (PCS) experiment was accommodated within International Space Station (ISS) EXpedite the PRocessing of Experiments to Space Station (EXPRESS) Rack 2 and was remotely operated from early June 2001 until February 2002 from NASA Glenn Research Center's Telescience Support Center in Cleveland, Ohio and from a remote site at Harvard University in Cambridge, Massachusetts. PCS is an experiment conceived by Professor David A. Weitz of Harvard University (the Principal Investigator), focusing on the behavior of three different classes of colloid mixtures. The sophisticated light scattering instrumentation comprising PCS is capable of color imaging, and dynamic and static light scattering from 11 to 169 degrees, Bragg scattering over the range from 10 to 60 degrees, and laser light scattering at low angles from 0.3 to 6.0 degrees. The PCS instrumentation performed remarkably well, demonstrating a flexibility that enabled experiments to be performed that had not been envisioned prior to launch. While on-orbit, PCS accomplished 2400 hours of science operations, and was declared a resounding success. Each of the eight sample cells worked well and produced interesting and important results. Crystal nucleation and growth and the resulting structures of two binary colloidal crystal alloys were studied, with the long duration microgravity environment of the ISS facilitating extended studies on the growth and coarsening characteristics of the crystals. In another experiment run, the de-mixing of the colloid-polymer critical-point sample was studied as it phase-separates into two phases, one that resembles a gas and one that resembles a liquid. This process was studied over four decades of length scale, from 1 micron to 1 centimeter, behavior that cannot be observed in this sample on Earth because sedimentation would cause the colloids to fall to the bottom of the cell faster than the de-mixing process could occur. Similarly, the study of gelation and aging of another colloid-polymer sample, the colloid-polymer gel, also provided valuable information on gelation mechanisms, as did investigations on the extremely the low concentration silica and polystyrene fractal gel samples.

Description: An explicit time filter is applied to the Navier-Stokes equation prior to a space filter. The time filter is supposed to be smooth, and an exact expansion depending on the time derivatives of the velocity is derived for the associated stress tensor. On the contrary, the effect of the space filter is treated as usual and an eddy viscosity model is introduced in the LES equation. The total stress is thus represented using a new class of mixed models combining time and space derivatives of the LES field.

Description: As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Description: This paper compares the filtering used in Coherent Vortex Simulation (CVS) decomposition with an orthogonal wavelet basis, with the Proper Orthogonal Decomposition (POD) or Fourier filtering. Both methods are applied to a field of Direct Numerical Simulation (DNS) data of 3D forced homogeneous isotropic turbulence at microscale Reynolds number R(sub lambda) = 168. We show that, with only 3%N retained modes, CVS filtering separates the coherent vortex tubes from the incoherent background flow. The latter is structureless, has an equipartition energy spectrum, and has a Gaussian velocity probability distribution function (PDF) and an exponential vorticity PDF. On the other hand, the Fourier basis does not extract the coherent vortex tubes cleanly and leaves organized structures in the residual high wavenumber modes whose PDFs are stretched exponentials for both the velocity and the vorticity.

Description: The objective of this study is to investigate the use of the second generation bi-orthogonal wavelet transform for the field decomposition in the Coherent Vortex Simulation of turbulent flows. The performances of the bi-orthogonal second generation wavelet transform and the orthogonal wavelet transform using Daubechies wavelets with the same number of vanishing moments are compared in a priori tests using a spectral direct numerical simulation (DNS) database of isotropic turbulence fields: 256(exp 3) and 512(exp 3) DNS of forced homogeneous turbulence (Re(sub lambda) = 168) and 256(exp 3) and 512(exp 3) DNS of decaying homogeneous turbulence (Re(sub lambda) = 55). It is found that bi-orthogonal second generation wavelets can be used for coherent vortex extraction. The results of a priori tests indicate that second generation wavelets have better compression and the residual field is closer to Gaussian. However, it was found that the use of second generation wavelets results in an integral length scale for the incoherent part that is larger than that derived from orthogonal wavelets. A way of dealing with this difficulty is suggested.

Description: The main goal of this presentation is to give some of the objectives of the testing program. This includes: develop jet noise data base for separate flow nozzles with bypass ratio's 5 to 14; evaluate effect of pylon on noise; develop low performance impact noise suppression concepts; and evaluate potential for active control of jet noise.

Description: Coherent Vortex Simulation (CVS) filtering has been applied to Direct Numerical Simulation (DNS) data of forced and unforced time-developing turbulent mixing layers. CVS filtering splits the turbulent flow into two orthogonal parts, one corresponding to coherent vortices and the other to incoherent background flow. We have shown that the coherent vortices can be represented by few wavelet modes and that these modes are sufficient to reproduce the vorticity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of CVS filtering, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

Description: This presentation will discuss Ir for online diagnostics with acoustics; Ptot, Ttot, Pstat rake surveys for mean flow measurements; focused Schleiren for density and some turbulence structure; Laser sheet visualization for near-nozzle diagnostics; and Two-point hotwire measurements for turbulence models.

Description: The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Description: The overall objective of this work is to study nucleate boiling heat transfer under microgravity conditions in such a way that while providing basic knowledge of the phenomena, it also leads to development of simulation models and correlations that can be used as design tools for a wide range of gravity levels. In the study a building block type of approach is used and both pool and low velocity flow boiling are investigated. Starting with experiments using a single bubble, the complexity of the experiments is increased to two or three inline bubbles, to five bubbles placed on a two-dimensional grid. Finally, experiments are conducted where a large number of prescribed cavities nucleate on the heater and when a commercial surface is used. So far experiments have been conducted at earth normal gravity and in the reduced gravity environment of the KC-135 aircraft whereas experiments on the space station are planned. Modeling/complete numerical simulation of the boiling process is an integral part of the total effort. Experiments conducted with single bubbles formed on a nucleation site microfabricated on a polished silicon wafer show that for gravity levels (g) varying from 1.5g(sub e) to 0.01g(sub e), the bubble diameter at departure varies approximately as (g(sub e)/g)(exp 1/2) and the growth period as (g(sub e)/g). When bubbles merge either inline or in a plane, the bubble diameter at departure is found to be smaller than that obtained for a single bubble and shows a weaker dependence on the level of gravity. The possible reason is that as the bubbles merge they create fluid circulation around the bubbles, which in turn induces a lift force that is responsible for the earlier departure of the bubbles. The verification of this proposition is being sought through numerical simulations. There is a merger of two inline, three inline, and several bubbles in a plane in the low gravity environment of the KC-135 aircraft. After merger and before departure, a mushroom type of bubble with several stems attached to the heater surface is clearly evident. Local heat fluxes during growth and departure of a single bubble were also measured. It was found that during most of the growth period of the bubble, generally the wall heat flux decreased with time because of the increased dry area under the bubble. However, the heat flux increased rapidly just prior to departure of the bubble because of the transient conduction into the cold liquid rushing to fill the space vacated by the bubble as the bubble base shrinks. The measured heat fluxes at various radial locations are found to be in qualitative agreement with the numerical predictions. Single bubble studies at earth normal gravity have also been performed on surfaces oriented at different angles to the gravitational acceleration with flow parallel to the surface. It is found that in all cases the bubbles slide along the surface before lift-off from the surface. The lift force generated as a result of the relative motion between the sliding bubbles and the imposed flow is found to play an important role when the normal force due to buoyancy is reduced. An experimental apparatus for the study of the bubble behavior with imposed flow under reduced gravity conditions has been developed and will soon be employed for experiments in the KC-135 aircraft.

Description: Dynamic wetting, the displacement of one fluid by another immiscible fluid on a surface, controls many natural and technological phenomena, such as coating, printing, spray painting and lubricating. Particularly in coating and spraying applications, contact lines advance across pre-existing fluid films. Most previous work has focused on contact lines advancing across films sufficiently thick that they behave as simple Newtonian fluids. Ultrathin films, where the film thickness may impinge on fundamental length scales in the fluid, have received less attention. In this talk, we will discuss the effects of ultrathin polymer films on dynamic wetting. We measure the interface shape within microns of moving contact lines advancing across preexisting films and compare the measurements to existing models of viscous bending for interfaces advancing across dry surfaces and 'thick' (in the sense that they behave as liquids) films. In the experiments, we advance a contact line of 10-poise and 1-poise polydimethylsiloxane (silicone oil) across pre-coated films of the same fluid with thickness from a single chain thickness (approx. 10 A) through a couple of radii of gyration (100-200 A) to films so thick they are likely bulk in behavior (10(exp 3) A). All films are physisorbed, i.e. they readily rinse from the surface. Thus, molecules in the film are not anchored to the surface and can move within the film if the hydrodynamics dictate such motion. For films of the thickness of a single chain (approx. 10 A), our experiments indicate that the advancing fluid behaves just as it would if it advanced over a dry surface. For the thicker films (10(exp 3) A), we find behavior indicating that the molecules in the film are acting as a fluid with the bulk properties. In this regime, results for the two different fluids are identical when the experiments are performed at the same pre-existing film thickness and advancing capillary number, Ca. For film of thickness of a few radii of gyration (approx. 100-200 A), the behavior depends on Ca of the advancing meniscus. At low Ca, the viscous bending of the interface near the contact line does not behave as it would on a dry surface. It has a lower curvature than expected. However, at higher Ca, the viscous bending is described by the model for spreading over a dry surface. These results show that the fluid flow in the film does behave differently than bulk as the film thickness becomes comparable to molecular length scale. But even more intriguing is the unusual velocity dependence of that behavior where the film behaves more solid-like at higher contact line speeds. We will discuss these results in terms of the properties of confined polymer melts.