Development of a regional model for the North Indian Ocean

We have developed a one-way nested Indian Ocean regional model. The model combines the National Oceanic and Atmospheric Administration (NOAA) Geophysical Fluid Dynamics Laboratory's (GFDL) Modular Ocean Model (MOM4p1) at global climate model resolution (nominally one degree), and a regional Indian Ocean MOM4p1 configuration with 25 km horizontal resolution and 1 m vertical resolution near the surface. Inter-annual global simulations with Coordinated Ocean-Ice Reference Experiments (CORE-II) surface forcing over years 1992-2005 provide surface boundary conditions. We show that relative to the global simulation, (i) biases in upper ocean temperature, salinity and mixed layer depth are reduced, (ii) sea surface height and upper ocean circulation are closer to observations, and (iii) improvements in model simulation can be attributed to refined resolution, more realistic topography and inclusion of seasonal river runoff. Notably, the surface salinity bias is reduced to less than 0.1 psu over the Bay of Bengal using relatively weak restoring to observations, and the model simulates the strong, shallow halocline often observed in the North Bay of Bengal. There is marked improvement in subsurface salinity and temperature, as well as mixed layer depth in the Bay of Bengal. Major seasonal signatures in observed sea surface height anomaly in the tropical Indian Ocean, including the coastal waveguide around the Indian peninsula, are simulated with great fidelity. The use of realistic topography and seasonal river runoff brings the three dimensional structure of the East India Coastal Current and West India Coastal Current much closer to observations. As a result, the incursion of low salinity Bay of Bengal water into the southeastern Arabian Sea is more realistic. (C) 2013 Elsevier Ltd. All rights reserved.

The deep ocean density structure at the Last Glacial Maximum : What was it and why?

The search for reliable proxies of past deep ocean temperature and salinity has proved difficult, thereby limiting our ability to understand the coupling of ocean circulation and climate over glacial-interglacial timescales. Previous inferences of deep ocean temperature and salinity from sediment pore fluid oxygen isotopes and chlorinity indicate that the deep ocean density structure at the Last Glacial Maximum (LGM, approximately 20,000 years BP) was set by salinity, and that the density contrast between northern and southern sourced deep waters was markedly greater than in the modern ocean. High density stratification could help explain the marked contrast in carbon isotope distribution recorded in the LGM ocean relative to that we observe today, but what made the ocean's density structure so different at the LGM? How did it evolve from one state to another? Further, given the sparsity of the LGM temperature and salinity data set, what else can we learn by increasing the spatial density of proxy records?

We investigate the cause and feasibility of a highly and salinity stratified deep ocean at the LGM and we work to increase the amount of information we can glean about the past ocean from pore fluid profiles of oxygen isotopes and chloride. Using a coupled ocean--sea ice--ice shelf cavity model we test whether the deep ocean density structure at the LGM can be explained by ice--ocean interactions over the Antarctic continental shelves, and show that a large contribution of the LGM salinity stratification can be explained through lower ocean temperature. In order to extract the maximum information from pore fluid profiles of oxygen isotopes and chloride we evaluate several inverse methods for ill-posed problems and their ability to recover bottom water histories from sediment pore fluid profiles. We demonstrate that Bayesian Markov Chain Monte Carlo parameter estimation techniques enable us to robustly recover the full solution space of bottom water histories, not only at the LGM, but through the most recent deglaciation and the Holocene up to the present. Finally, we evaluate a non-destructive pore fluid sampling technique, Rhizon samplers, in comparison to traditional squeezing methods and show that despite their promise, Rhizons are unlikely to be a good sampling tool for pore fluid measurements of oxygen isotopes and chloride.

Year-to-year changes of vertical temperature distribution in the California Current region: 1954 to 1986

Studies by Enfield and Allen (1980), McLain et al (1985), and others have shown that anomalously warm years in the northern coastal California Current correspond to El Niño conditions in the equatorial Pacific Ocean. Ocean model studies suggest a mechanical link between the northern coastal California Current and the equatorial ocean through long waves that propagate cyclonically along the ocean boundary (McCreary 1976; Clarke 1983; Shriver et al 1991). However, distinct observational evidence of such an oceanic connection is not extensive. Much of the supposed El Niño variation in temperature and sea level data from the coastal California Current region can be associated with the effects of anomalously intense north Pacific atmospheric cyclogenesis, which is frequently augmented during El Niño years (Wallace and Gutzler 1981; Simpson 1983; Emery and Hamilton 1984). This study uses time series of ocean temperature data to distinguish between locally forced effects, initiated by north Pacific atmospheric changes, and remotely forced effects, initiated by equatorial Pacific atmospheric changes related to El Niño events.

A description of water masses in Nootka Sound and estimation of diffusivity and upward velocity in Muchalat Inlet using temperature, salinity and dissolved oxygen

Senior thesis written for Oceanography 445

Isolating the signal of ocean global warming

Identifying the signature of global warming in the world's oceans is challenging because low frequency circulation changes can dominate local temperature changes. The IPCC fourth assessment reported an average ocean heating rate of 0.21 ± 0.04 Wm−2 over the period 1961–2003, with considerable spatial, interannual and inter-decadal variability. We present a new analysis of millions of ocean temperature profiles designed to filter out local dynamical changes to give a more consistent view of the underlying warming. Time series of temperature anomaly for all waters warmer than 14°C show large reductions in interannual to inter-decadal variability and a more spatially uniform upper ocean warming trend (0.12 Wm−2 on average) than previous results. This new measure of ocean warming is also more robust to some sources of error in the ocean observing system. Our new analysis provides a useful addition for evaluation of coupled climate models, to the traditional fixed depth analyses.

Idealized model for changes in equilibrium temperature, mixed layer depth, and boundary layer cloud over land in a doubled CO2 climate

An idealized equilibrium model for the undisturbed partly cloudy boundary layer (BL) is used as a framework to explore the coupling of the energy, water, and carbon cycles over land in midlatitudes and show the sensitivity to the clear‐sky shortwave flux, the midtropospheric temperature, moisture, CO2, and subsidence. The changes in the surface fluxes, the BL equilibrium, and cloud cover are shown for a warmer, doubled CO2 climate. Reduced stomatal conductance in a simple vegetation model amplifies the background 2 K ocean temperature rise to an (unrealistically large) 6 K increase in near‐surface temperature over land, with a corresponding drop of near‐surface relative humidity of about 19%, and a rise of cloud base of about 70 hPa. Cloud changes depend strongly on changes of mean subsidence; but evaporative fraction (EF) decreases. EF is almost uniquely related to mixed layer (ML) depth, independent of background forcing climate. This suggests that it might be possible to infer EF for heterogeneous landscapes from ML depth. The asymmetry of increased evaporation over the oceans and reduced transpiration over land increases in a warmer doubled CO2 climate.

Attribution of the spatial pattern of CO2-forced sea level change to ocean surface flux changes

Climate models taking part in the coupled model intercomparison project phase 5 (CMIP5) all predict a global mean sea level rise for the 21st century. Yet the sea level change is not spatially uniform and differs among models. Here we evaluate the role of air–sea fluxes of heat, water and momentum (windstress) to find the spatial pattern associated to each of them as well as the spread they can account for. Using one AOGCM to which we apply the surface flux changes from other AOGCMs, we show that the heat flux and windstress changes dominate both the pattern and the spread, but taking the freshwater flux into account as well yields a sea level change pattern in better agreement with the CMIP5 ensemble mean. Differences among the CMIP5 control ocean temperature fields have a smaller impact on the sea level change pattern.

Heat tolerance, behavioural temperature selection and temperature-dependent respiration in larval Octopus huttoni

This study reports temperature effects on paralarvae from a benthic octopus species, Octopus huttoni, found throughout New Zealand and temperate Australia. We quantified the thermal tolerance, thermal preference and temperature-dependent respiration rates in 1-5 days old paralarvae. Thermal stress (1°C increase h^-1) and thermal selection (~10-24°C vertical gradient) experiments were conducted with paralarvae reared for 4 days at 16°C. In addition, measurement of oxygen consumption at 10, 15, 20 and 25°C was made for paralarvae aged 1, 4 and 5 days using microrespirometry. Onset of spasms, rigour (CT_max) and mortality (upper lethal limit) occurred for 50% of experimental animals at, respectively, 26.0±0.2°C, 27.8±0.2°C and 31.4±0.1°C. The upper, 23.1±0.2°C, and lower, 15.0±1.7°C, temperatures actively avoided by paralarvae correspond with the temperature range over which normal behaviours were observed in the thermal stress experiments. Over the temperature range of 10°C-25°C, respiration rates, standardized for an individual larva, increased with age, from 54.0 to 165.2nmol larvae^-1h^-1 in one-day old larvae to 40.1-99.4nmol h^-1 at five days. Older larvae showed a lesser response to increased temperature: the effect of increasing temperature from 20 to 25°C (Q₁₀) on 5 days old larvae (Q₁₀=1.35) was lower when compared with the 1 day old larvae (Q₁₀=1.68). The lower Q₁₀ in older larvae may reflect age-related changes in metabolic processes or a greater scope of older larvae to respond to thermal stress such as by reducing activity. Collectively, our data indicate that temperatures >25°C may be a critical temperature. Further studies on the population-level variation in thermal tolerance in this species are warranted to predict how continued increases in ocean temperature will limit O. huttoni at early larval stages across the range of this species.

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.

Biogenic habitat structure provided by temperate macroalgae change along a latitudinal gradient in ocean climate

[EN] Global warming is affecting all major ecosystems, including temperate reefs where canopy-forming seaweeds provide biogenic habitat. In contrast to the rapidly growing recognition of how climate affects the performance and distribution of individuals and populations, relatively little is known about possible links between climate and biogenic habitat structure. We examined the relationship between several ocean temperature characteristics, expressed on time-scales of days, months and years, on habitat patch characteristics on 24 subtidal temperate reefs along a latitudinal gradient (Western Australia; ca 34 to 27º S). Significant climate related variation in habitat structure was observed, even though the landscape cover of kelp and fucalean canopies did not change across the climate gradient: monospecific patches of kelp became increasingly dominant in warmer climates, at the expense of mixed kelp-fucalean canopies. The decline in mixed canopies was associated with an increase in the abundance of Sargassum spp., replacing a more diverse canopy assemblage of Scytothalia doryocarpa and several other large fucoids. There were no observed differences in the proportion of open gaps or gap characteristics. These habitat changes were closely related to patterns in minimum temperatures and temperature thresholds (days > 20 °C), presumably because temperate algae require cool periods for successful reproduction and recruitment (even if the adults can survive warmer temperatures). Although the observed habitat variation may appear subtle, similar structural differences have been linked to a range of effects on canopy-associated organisms through the provision of habitat and ecosystem engineering. Consequently, our study suggests that the magnitude of projected temperature increase is likely to cause changes in habitat structure and thereby indirectly affect numerous habitat-dependent plants and animals

Temporal dynamics of benthic assemblages along a gradient of ocean acidification at a CO2 vent’s system of the Ischia Island

Global warming and ocean acidification, due to rising atmospheric levels of CO2, represent an actual threat to terrestrial and marine environments. Since Industrial Revolution, in less of 250 years, pH of surface seawater decreased on average of 0.1 unit, and is expected to further decreases of approximately 0.3-0.4 units by the end of this century. Naturally acidified marine areas, such as CO2 vent systems at the Ischia Island, allow to study acclimatation and adaptation of individual species as well as the structure of communities, and ecosystems to OA. The main aim of this thesis was to study how hard bottom sublittoral benthic assemblages changed trough time along a pH gradient. For this purpose, the temporal dynamics of mature assemblages established on artificial substrates (volcanic tiles) over a 3 year- period were analysed. Our results revealed how composition and dynamics of the community were altered and highly simplified at different level of seawater acidification. In fact, extreme low values of pH (approximately 6.9), affected strongly the assemblages, reducing diversity both in terms of taxa and functional groups, respect to lower acidification levels (mean pH 7.8) and ambient conditions (8.1 unit). Temporal variation was observed in terms of species composition but not in functional groups. Variability was related to species belonging to the same functional group, suggesting the occurrence of functional redundancy. Therefore, the analysis of functional groups kept information on the structure, but lost information on species diversity and dynamics. Decreasing in ocean pH is only one of many future global changes that will occur at the end of this century (increase of ocean temperature, sea level rise, eutrophication etc.). The interaction between these factors and OA could exacerbate the community and ecosystem effects showed by this thesis.

Impact of a Reduced Arctic Sea Ice Cover on Ocean and Atmospheric Properties

The Arctic sea ice cover declined over the last few decades and reached a record minimum in 2007, with a slight recovery thereafter. Inspired by this the authors investigate the response of atmospheric and oceanic properties to a 1-yr period of reduced sea ice cover. Two ensembles of equilibrium and transient simulations are produced with the Community Climate System Model. A sea ice change is induced through an albedo change of 1 yr. The sea ice area and thickness recover in both ensembles after 3 and 5 yr, respectively. The sea ice anomaly leads to changes in ocean temperature and salinity to a depth of about 200 m in the Arctic Basin. Further, the salinity and temperature changes in the surface layer trigger a “Great Salinity Anomaly” in the North Atlantic that takes roughly 8 yr to travel across the North Atlantic back to high latitudes. In the atmosphere the changes induced by the sea ice anomaly do not last as long as in the ocean. The response in the transient and equilibrium simulations, while similar overall, differs in specific regional and temporal details. The surface air temperature increases over the Arctic Basin and the anomaly extends through the whole atmospheric column, changing the geopotential height fields and thus the storm tracks. The patterns of warming and thus the position of the geopotential height changes vary in the two ensembles. While the equilibrium simulation shifts the storm tracks to the south over the eastern North Atlantic and Europe, the transient simulation shifts the storm tracks south over the western North Atlantic and North America. The authors propose that the overall reduction in sea ice cover is important for producing ocean anomalies; however, for atmospheric anomalies the regional location of the sea ice anomalies is more important. While observed trends in Arctic sea ice are large and exceed those simulated by comprehensive climate models, there is little evidence based on this particular model that the seasonal loss of sea ice (e.g., as occurred in 2007) would constitute a threshold after which the Arctic would exhibit nonlinear, irreversible, or strongly accelerated sea ice loss. Caution should be exerted when extrapolating short-term trends to future sea ice behavior.