Ocean heat content variability and change in an ensemble of ocean reanalyses

Accurate knowledge of the location and magnitude of ocean heat content (OHC) variability and change is essential for understanding the processes that govern decadal variations in surface temperature, quantifying changes in the planetary energy budget, and developing constraints on the transient climate response to external forcings. We present an overview of the temporal and spatial characteristics of OHC variability and change as represented by an ensemble of dynamical and statistical ocean reanalyses (ORAs). Spatial maps of the 0–300 m layer show large regions of the Pacific and Indian Oceans where the interannual variability of the ensemble mean exceeds ensemble spread, indicating that OHC variations are well-constrained by the available observations over the period 1993–2009. At deeper levels, the ORAs are less well-constrained by observations with the largest differences across the ensemble mostly associated with areas of high eddy kinetic energy, such as the Southern Ocean and boundary current regions. Spatial patterns of OHC change for the period 1997–2009 show good agreement in the upper 300 m and are characterized by a strong dipole pattern in the Pacific Ocean. There is less agreement in the patterns of change at deeper levels, potentially linked to differences in the representation of ocean dynamics, such as water mass formation processes. However, the Atlantic and Southern Oceans are regions in which many ORAs show widespread warming below 700 m over the period 1997–2009. Annual time series of global and hemispheric OHC change for 0–700 m show the largest spread for the data sparse Southern Hemisphere and a number of ORAs seem to be subject to large initialization ‘shock’ over the first few years. In agreement with previous studies, a number of ORAs exhibit enhanced ocean heat uptake below 300 and 700 m during the mid-1990s or early 2000s. The ORA ensemble mean (±1 standard deviation) of rolling 5-year trends in full-depth OHC shows a relatively steady heat uptake of approximately 0.9 ± 0.8 W m−2 (expressed relative to Earth’s surface area) between 1995 and 2002, which reduces to about 0.2 ± 0.6 W m−2 between 2004 and 2006, in qualitative agreement with recent analysis of Earth’s energy imbalance. There is a marked reduction in the ensemble spread of OHC trends below 300 m as the Argo profiling float observations become available in the early 2000s. In general, we suggest that ORAs should be treated with caution when employed to understand past ocean warming trends—especially when considering the deeper ocean where there is little in the way of observational constraints. The current work emphasizes the need to better observe the deep ocean, both for providing observational constraints for future ocean state estimation efforts and also to develop improved models and data assimilation methods.

An ensemble of eddy-permitting global ocean reanalyses from the MyOcean project

A set of four eddy-permitting global ocean reanalyses produced in the framework of the MyOcean project have been compared over the altimetry period 1993–2011. The main differences among the reanalyses used here come from the data assimilation scheme implemented to control the ocean state by inserting reprocessed observations of sea surface temperature (SST), in situ temperature and salinity profiles, sea level anomaly and sea-ice concentration. A first objective of this work includes assessing the interannual variability and trends for a series of parameters, usually considered in the community as essential ocean variables: SST, sea surface salinity, temperature and salinity averaged over meaningful layers of the water column, sea level, transports across pre-defined sections, and sea ice parameters. The eddy-permitting nature of the global reanalyses allows also to estimate eddy kinetic energy. The results show that in general there is a good consistency between the different reanalyses. An intercomparison against experiments without data assimilation was done during the MyOcean project and we conclude that data assimilation is crucial for correctly simulating some quantities such as regional trends of sea level as well as the eddy kinetic energy. A second objective is to show that the ensemble mean of reanalyses can be evaluated as one single system regarding its reliability in reproducing the climate signals, where both variability and uncertainties are assessed through the ensemble spread and signal-to-noise ratio. The main advantage of having access to several reanalyses differing in the way data assimilation is performed is that it becomes possible to assess part of the total uncertainty. Given the fact that we use very similar ocean models and atmospheric forcing, we can conclude that the spread of the ensemble of reanalyses is mainly representative of our ability to gauge uncertainty in the assimilation methods. This uncertainty changes a lot from one ocean parameter to another, especially in global indices. However, despite several caveats in the design of the multi-system ensemble, the main conclusion from this study is that an eddy-permitting multi-system ensemble approach has become mature and our results provide a first step towards a systematic comparison of eddy-permitting global ocean reanalyses aimed at providing robust conclusions on the recent evolution of the oceanic state.

Southern Ocean isopycnal mixing and ventilation changes driven by winds

Observed and predicted changes in the strength of the westerly winds blowing over the Southern Ocean have motivated a number of studies of the response of the Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation (MOC) to wind perturbations and led to the discovery of the``eddy-compensation" regime, wherein the MOC becomes insensitive to wind changes. In addition to the MOC, tracer transport also depends on mixing processes. Here we show, in a high-resolution process model, that isopycnal mixing by mesoscale eddies is strongly dependent on the wind strength. This dependence can be explained by mixing-length theory and is driven by increases in eddy kinetic energy; the mixing length does not change strongly in our simulation. Simulation of a passive ventilation tracer (analogous to CFCs or anthropogenic CO$_2$) demonstrates that variations in tracer uptake across experiments are dominated by changes in isopycnal mixing, rather than changes in the MOC. We argue that, to properly understand tracer uptake under different wind-forcing scenarios, the sensitivity of isopycnal mixing to winds must be accounted for.

Decadal predictability of regional scale wind speed and wind energy potentials over Central Europe

Decadal predictions on timescales from one year to one decade are gaining importance since this time frame falls within the planning horizon of politics, economy and society. The present study examines the decadal predictability of regional wind speed and wind energy potentials in three generations of the MiKlip (‘Mittelfristige Klimaprognosen’) decadal prediction system. The system is based on the global Max-Planck-Institute Earth System Model (MPI-ESM), and the three generations differ primarily in the ocean initialisation. Ensembles of uninitialised historical and yearly initialised hindcast experiments are used to assess the forecast skill for 10 m wind speeds and wind energy output (Eout) over Central Europe with lead times from one year to one decade. With this aim, a statistical-dynamical downscaling (SDD) approach is used for the regionalisation. Its added value is evaluated by comparison of skill scores for MPI-ESM large-scale wind speeds and SDD-simulated regional wind speeds. All three MPI-ESM ensemble generations show some forecast skill for annual mean wind speed and Eout over Central Europe on yearly and multi-yearly time scales. This forecast skill is mostly limited to the first years after initialisation. Differences between the three ensemble generations are generally small. The regionalisation preserves and sometimes increases the forecast skills of the global runs but results depend on lead time and ensemble generation. Moreover, regionalisation often improves the ensemble spread. Seasonal Eout skills are generally lower than for annual means. Skill scores are lowest during summer and persist longest in autumn. A large-scale westerly weather type with strong pressure gradients over Central Europe is identified as potential source of the skill for wind energy potentials, showing a similar forecast skill and a high correlation with Eout anomalies. These results are promising towards the establishment of a decadal prediction system for wind energy applications over Central Europe.

South Atlantic Ocean cyclogenesis climatology simulated by regional climate model (RegCM3)

A detailed climatology of the cyclogenesis over the Southern Atlantic Ocean (SAO) from 1990 to 1999 and how it is simulated by the RegCM3 (Regional Climate Model) is presented here. The simulation used as initial and boundary conditions the National Centers for Environmental Prediction-Department of Energy (NCEP/DOE) reanalysis. The cyclones were identified with an automatic scheme that searches for cyclonic relative vorticity (zeta(10)) obtained from a 10-m height wind field. All the systems with zeta(10) a parts per thousand currency sign -1.5 x 10(-5) s(-1) and lifetime equal or larger than 24 h were considered in the climatology. Over SAO, in 10 years were detected 2,760 and 2,787 cyclogeneses in the simulation and NCEP, respectively, with an annual mean of 276.0 +/- A 11.2 and 278.7 +/- A 11.1. This result suggests that the RegCM3 has a good skill to simulate the cyclogenesis climatology. However, the larger model underestimations (-9.8%) are found for the initially stronger systems (zeta(10) a parts per thousand currency sign -2.5 x 10(-5) s(-1)). It was noted that over the SAO the annual cycle of the cyclogenesis depends of its initial intensity. Considering the systems initiate with zeta(10) a parts per thousand currency sign -1.5 x 10(-5) s(-1), the annual cycle is not well defined and the higher frequency occurs in the autumn (summer) in the NCEP (RegCM3). The stronger systems (zeta(10) a parts per thousand currency sign -2.5 x 10(-5) s(-1)) have a well-characterized high frequency of cyclogenesis during the winter in both NCEP and RegCM3. This work confirms the existence of three cyclogenetic regions in the west sector of the SAO, near the South America east coast and shows that RegCM3 is able to reproduce the main features of these cyclogenetic areas.

STUDY OF THE EQUATORIAL ATLANTIC OCEAN MIXING LAYER USING A ONE-DIMENSIONAL TURBULENCE MODEL

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

South Atlantic Ocean cyclogenesis climatology simulated by regional climate model (RegCM3)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelo OLAM (ocean-land-atmosphere-model): descrição, aplicações, e perspectivas

O modelo OLAM foi desenvolvido com objetivo de estender a capacidade de representar os fenômenos de escala global e regional simultaneamente. Este modelo apresenta inovações quanto aos processos dinâmicos, configuração de grade, estrutura de memória e técnicas de solução numérica das equações prognósticas. As equações de Navier-Stokes são resolvidas através da técnica de volumes finitos que conservam massa, momento e energia. No presente trabalho, apresenta-se uma descrição sucinta do OLAM e alguns resultados de sua aplicação em simulações climáticas da precipitação mensal para a região norte da América do Sul, bem como em rodadas para previsão numérica de tempo regional. Os resultados mostram que o modelo consegue representar bem os aspectos meteorológicos de grande escala. Em geral, seu desempenho melhora quando são adotadas grades de maior resolução espacial, nas quais se verificam melhorias significativas tanto na estimativa da precipitação mensal regional, quanto na previsão numérica de tempo.

Daily morphological changes determined by high-energy events on an embayed beach: a qualitative model

This paper examines the daily morphological responses of Sununga Beach, an embayed beach located on the south-eastern Brazilian coast, to storms in the South Atlantic Ocean. The main mechanisms and timing of beach erosion and accretion, the relationship between wave height and direction, and beach volume changes are considered, to establish a qualitative model for short-term embayed beach morphological changes. The methodology consisted of daily topographic surveys during the month of May in 2001, 2002, and 2003, using an RTK-GPS (real-time kinematics global positioning system). Weather and wave model results were used to correlate hydrodynamics and beach morphology. The results indicate that the morphodynamics of Sununga Beach are characterized by a process of beach rotation, which occurred more or less clearly during all three surveys. Unlike what has been commonly described in the literature for longer time intervals and alternations of fair and stormy weather, the beach rotation processes on Sununga Beach occurred under conditions of moderate-to-high wave energy change (wave heights greater than 2 m). An integrated evaluation of the behaviour of the meteorological aspects, together with beach morphology, enabled us to recognize that extra-tropical cyclones were the most important agent in remobilizing the beach planform, whether in beach rotation or in cross-shore erosion. Copyright (C) 2010 John Wiley & Sons, Ltd.

Mean and Variability of the Tropical Atlantic Ocean in the CCSM4

This study analyzes important aspects of the tropical Atlantic Ocean from simulations of the fourth version of the Community Climate System Model (CCSM4): the mean sea surface temperature (SST) and wind stress, the Atlantic warm pools, the principal modes of SST variability, and the heat budget in the Benguela region. The main goal was to assess the similarities and differences between the CCSM4 simulations and observations. The results indicate that the tropical Atlantic overall is realistic in CCSM4. However, there are still significant biases in the CCSM4 Atlantic SSTs, with a colder tropical North Atlantic and a hotter tropical South Atlantic, that are related to biases in the wind stress. These are also reflected in the Atlantic warm pools in April and September, with its volume greater than in observations in April and smaller than in observations in September. The variability of SSTs in the tropical Atlantic is well represented in CCSM4. However, in the equatorial and tropical South Atlantic regions, CCSM4 has two distinct modes of variability, in contrast to observed behavior. A model heat budget analysis of the Benguela region indicates that the variability of the upper-ocean temperature is dominated by vertical advection, followed by meridional advection.

The Southern Ocean and Its Climate in CCSM4

The new Community Climate System Model, version 4 (CCSM4), provides a powerful tool to understand and predict the earth's climate system. Several aspects of the Southern Ocean in the CCSM4 are explored, including the surface climatology and interannual variability, simulation of key climate water masses (Antarctic Bottom Water, Subantarctic Mode Water, and Antarctic Intermediate Water), the transport and structure of the Antarctic Circumpolar Current, and interbasin exchange via the Agulhas and Tasman leakages and at the Brazil-Malvinas Confluence. It is found that the CCSM4 has varying degrees of accuracy in the simulation of the climate of the Southern Ocean when compared with observations. This study has identified aspects of the model that warrant further analysis that will result in a more comprehensive understanding of ocean-atmosphere-ice dynamics and interactions that control the earth's climate and its variability.

Turbulence as a driver for vertical plankton distribution in the subsurface upper ocean

[EN] Vertical distributions of turbulent energy dissipation rates and fluorescence were measured simultaneously with a high-resolution micro-profiler in four different oceanographic regions, from temperate to polar and from coastal to open waters settings. High fluorescence values, forming a deep chlorophyll maximum (DCM), were often located in weakly stratified portions of the upper water column, just below layers with maximum levels of turbulent energy dissipation rate. In the vicinity of the DCM, a significant negative relationship between fluorescence and turbulent energy dissipation rate was found. We discuss the mechanisms that may explain the observed patterns of planktonic biomass distribution within the ocean mixed layer, including a vertically variable diffusion coefficient and the alteration of the cells sinking velocity by turbulent motion. These findings provide further insight into the processes controlling the vertical distribution of the pelagic community and position of the DCM.

A parallel future for Ocean research?

[EN] Today, science is difficult to pursue because funding is so tenuous. In such a financial climate, researchers need to consider parallel alternatives to ensure that scientific research can continue. Based on this thinking, we created BIOCEANSolutions, a company born of a research group. A great variety of environmental regulations and standards have emerged over recent years with the purpose of protecting natural ecosystems. These have enabled us to link our research to the market of environmental management. Marine activities can alter environmental conditions, resulting in changes in physiological states, species diversity, abundance, and biomass in the local biological communities. In this way, we can apply our knowledge, to plankton ecophysiology and biochemical oceanography. We measure enzyme activities as bio-indicators of energy metabolism and other physiological rates and biologic-oceanographic processes in marine organisms. This information provides insight into the health of marine communities, the stress levels of individual organisms, and potential anomalies that may be affecting them. In the process of verifying standards and complying with regulations, we can apply our analytic capability and knowledge. The main analyses that we offer are: (1) the activity of the electron transport system (ETS) or potential respiration (Φ), (2) the physiological measurement of respiration (oxygen consumption), (3) the activity of Isocitrate dehydrogenase (IDH), (4) the respiratory CO2 production, and (5) the activity of Glutamate dehydrogenase (GDH) and (6) the physiological measurement of ammonium excretion. In addition, our experience in a productive research group allows us to pursue and develop technical-experimental activities such as marine and freshwater aquaculture, oceanographic field sampling, as well as providing guidance, counseling, and academic services. In summary, this new company will permit us to create a symbiosis between public and private sectors that serve clients and will allow us to grow and expand as a research team.

Study of a wind-wave numerical model and its integration with ocean and oil-spill numerical models

The ability to represent the transport and fate of an oil slick at the sea surface is a formidable task. By using an accurate numerical representation of oil evolution and movement in seawater, the possibility to asses and reduce the oil-spill pollution risk can be greatly improved. The blowing of the wind on the sea surface generates ocean waves, which give rise to transport of pollutants by wave-induced velocities that are known as Stokes’ Drift velocities. The Stokes’ Drift transport associated to a random gravity wave field is a function of the wave Energy Spectra that statistically fully describe it and that can be provided by a wave numerical model. Therefore, in order to perform an accurate numerical simulation of the oil motion in seawater, a coupling of the oil-spill model with a wave forecasting model is needed. In this Thesis work, the coupling of the MEDSLIK-II oil-spill numerical model with the SWAN wind-wave numerical model has been performed and tested. In order to improve the knowledge of the wind-wave model and its numerical performances, a preliminary sensitivity study to different SWAN model configuration has been carried out. The SWAN model results have been compared with the ISPRA directional buoys located at Venezia, Ancona and Monopoli and the best model settings have been detected. Then, high resolution currents provided by a relocatable model (SURF) have been used to force both the wave and the oil-spill models and its coupling with the SWAN model has been tested. The trajectories of four drifters have been simulated by using JONSWAP parametric spectra or SWAN directional-frequency energy output spectra and results have been compared with the real paths traveled by the drifters.

Hydraulic evaluation of the gyro electric wave energy converter

Motivation Thanks for a scholarship offered by ALma Mater Studiorum I could stay in Denmark for six months during which I could do physical tests on the device Gyro PTO at the Departmet of Civil Engineering of Aalborg University. Aim The goal of my thesis is an hydraulic evaluation of the device: Gyro PTO, a gyroscopic device for conversion of mechanical energy in ocean surface waves to electrical energy. The principle of the system is the application of the gyroscopic moment of flywheels equipped on a swing float excited by waves. The laboratory activities were carried out by: Morten Kramer, Jan Olsen, Irene Guaraldi, Morten Thøtt, Nikolaj Holk. The main purpose of the tests was to investigate the power absorption performance in irregular waves, but testing also included performance measures in regular waves and simple tests to get knowledge about characteristics of the device, which could facilitate the possibility of performing numerical simulations and optimizations. Methodology To generate the waves and measure the performance of the device a workstation was created in the laboratory. The workstation consist of four computers in each of wich there was a different program. Programs have been used : Awasys6, LabView, Wave lab, Motive optitrack, Matlab, Autocad Main Results Thanks to the obtained data with the tank testing was possible to make the process of wave analisys. We obtained significant wave height and period through a script Matlab and then the values of power produced, and energy efficiency of the device for two types of waves: regular and irregular. We also got results as: physical size, weight, inertia moments, hydrostatics, eigen periods, mooring stiffness, friction, hydrodynamic coefficients etc. We obtained significant parameters related to the prototype in the laboratory after which we scale up the results obtained for two future applications: one in Nissun Brending and in the North Sea. Conclusions The main conclusion on the testing is that more focus should be put into ensuring a stable and positive power output in a variety of wave conditions. In the irregular waves the power production was negative and therefore it does not make sense to scale up the results directly. The average measured capture width in the regular waves was 0.21 m. As the device width is 0.63 m this corresponds to a capture width ratio of: 0.21/0.63 * 100 = 33 %. Let’s assume that it is possible to get the device to produce as well in irregular waves under any wave conditions, and lets further assume that the yearly absorbed energy can be converted into electricity at a PTO-efficiency of 90 %. Under all those assumptions the results in table are found, i.e. a Nissum Bredning would produce 0.87 MWh/year and a North Sea device 85 MWh/year.