Harvesting the sea: the exploitation of marine resources in the Roman Mediterranean

Harvesting the Sea provides the first systematic treatment of the exploitation of various marine resources, such as large-scale fishing, fish salting, salt and purple-dye production, and oyster and fish-farming, in the Roman world and its role within the ancient economy. Bringing together literary, epigraphic, and legal sources, with a wealth of archaeological data collected in recent years, the book shows that these marine resources were an important feature of the Roman economy and, in scope and market-oriented production, paralleled phenomena taking place in the Roman agricultural economy on land. The book also examines the importance of technological innovations, the organization of labour, and the use of the existing legal framework in defence of economic interests against competitors for the same natural resource.

Climate on the Southern Black Sea coast during the Holocene: implications from the Sofular Cave record

We present the updated Holocene section of the Sofular Cave record from the southernBlackSeacoast (northern Turkey); an area with considerably different present-day climate compared to that of the neighboring Eastern Mediterranean region. Stalagmite δ13C, growth rates and initial (234U/238U) ratios provide information about hydrological changes above the cave; and prove to be more useful than δ18O for deciphering Holocene climatic variations. Between ∼9.6 and 5.4 ka BP (despite a pause from ∼8.4 to 7.8 ka BP), the Sofular record indicates a remarkable increase in rainfall amount and intensity, in line with other paleoclimate studies in the Eastern Mediterranean. During that period, enhanced summertime insolation either produced much stronger storms in the following fall and winter through high sea surface temperatures, or it invoked a regional summer monsoon circulation and rainfall. We suggest that one or both of these climatic mechanisms led to a coupling of the BlackSea and the Mediterranean rainfall regimes at that time, which can explain the observed proxy signals. However, there are discrepancies among the Eastern Mediterranean records in terms of the timing of this wet period; implying that changes were probably not always occurring through the same mechanism. Nevertheless, the Sofular Cave record does provide hints and bring about new questions about the connection between regional and large scale climates, highlighting the need for a more extensive network of high quality paleoclimate records to better understand Holoceneclimate.

Landscape, regional and global estimates of nitrogen flux from land to sea: errors and uncertainties

Regional to global scale modelling of N flux from land to ocean has progressed to date through the development of simple empirical models representing bulk N flux rates from large watersheds, regions, or continents on the basis of a limited selection of model parameters. Watershed scale N flux modelling has developed a range of physically-based approaches ranging from models where N flux rates are predicted through a physical representation of the processes involved, through to catchment scale models which provide a simplified representation of true systems behaviour. Generally, these watershed scale models describe within their structure the dominant process controls on N flux at the catchment or watershed scale, and take into account variations in the extent to which these processes control N flux rates as a function of landscape sensitivity to N cycling and export. This paper addresses the nature of the errors and uncertainties inherent in existing regional to global scale models, and the nature of error propagation associated with upscaling from small catchment to regional scale through a suite of spatial aggregation and conceptual lumping experiments conducted on a validated watershed scale model, the export coefficient model. Results from the analysis support the findings of other researchers developing macroscale models in allied research fields. Conclusions from the study confirm that reliable and accurate regional scale N flux modelling needs to take account of the heterogeneity of landscapes and the impact that this has on N cycling processes within homogenous landscape units.

Vegetation and environmental dynamics in the southern Black Sea region since 18 kyr BP derived from the marine core 22-GC3

Sediments from the Black Sea, a region historically dominated by forests and steppe landscapes, are a valuable source of detailed information on the changes in regional terrestrial and aquatic environments at decadal to millennial scales. Here we present multi-proxy environmental records (pollen, dinoflagellate cysts, Ca, Ti and oxygen isotope data) from the uppermost 305 cm of the core 22-GC3 (42°13.53′N, 36°29.55′E) collected from a water depth of 838 m in the southern part of the Black Sea in 2007. The records span the last ~ 18 kyr (all ages are given in cal kyr BP). The pollen data reveal the dominance of the Artemisia-steppe in the region, suggesting rather dry/cold environments ~ 18–14.5 kyr BP. Warming/humidity increase during melt-water pulses (~ 16.1–14.5 kyr BP), indicated by δ18O records from the 22-GC3 core sediment and from the Sofular Cave stalagmite, is expressed in more negative δ13C values from the Sofular Cave, usually interpreted as the spreading of C3 plants. The records representing the interstadial complex (~ 14.5–12.9 kyr BP) show an increase in temperature and moisture, indicated by forest development, increased primary productivity and reduced surface run-off, whereas the switch from primary terrigenous to primary authigenic Ca origin occurs ~ 500 yr later. The Younger Dryas cooling is clearly demonstrated by more negative δ13C values from the Sofular Cave and a reduction of pines. The early Holocene (11.7–8.5 kyr BP) interval reveals relatively dry conditions compared to the mostly moist and warm middle Holocene (8.5–5 kyr BP), which is characterized by the establishment of the species-rich warm mixed and temperate deciduous forests in the low elevation belt, temperate deciduous beech-hornbeam forests in the middle and cool conifer forest in upper mountain belt. The border between the early and middle Holocene in the vegetation records coincides with the opening of the Mediterranean corridor at ~ 8.3 kyr BP, as indicated by a marked change in the dinocyst assemblages and in the sediment lithology. Changes in the pollen assemblages indicate a reduction in forest cover after ~ 5 kyr BP, which was likely caused by increased anthropogenic pressure on the regional vegetation.

The link between the Barents Sea and ENSO events simulated by NEMO model

An analysis of observational data in the Barents Sea along a meridian at 33°30' E between 70°30' and 72°30' N has reported a negative correlation between El Niño/La Niña Southern Oscillation (ENSO) events and water temperature in the top 200 m: the temperature drops about 0.5 °C during warm ENSO events while during cold ENSO events the top 200 m layer of the Barents Sea is warmer. Results from 1 and 1/4-degree global NEMO models show a similar response for the whole Barents Sea. During the strong warm ENSO event in 1997–1998 an anomalous anticyclonic atmospheric circulation over the Barents Sea enhances heat loses, as well as substantially influencing the Barents Sea inflow from the North Atlantic, via changes in ocean currents. Under normal conditions along the Scandinavian peninsula there is a warm current entering the Barents Sea from the North Atlantic, however after the 1997–1998 event this current is weakened. During 1997–1998 the model annual mean temperature in the Barents Sea is decreased by about 0.8 °C, also resulting in a higher sea ice volume. In contrast during the cold ENSO events in 1999–2000 and 2007–2008, the model shows a lower sea ice volume, and higher annual mean temperatures in the upper layer of the Barents Sea of about 0.7 °C. An analysis of model data shows that the strength of the Atlantic inflow in the Barents Sea is the main cause of heat content variability, and is forced by changing pressure and winds in the North Atlantic. However, surface heat-exchange with the atmosphere provides the means by which the Barents sea heat budget relaxes to normal in the subsequent year after the ENSO events.

Pleistocene water intrusions from the Mediterranean and Caspian Seas into the Black Sea

The hydrological balance of the Black Sea is governed by riverine input and by the exchange with the Mediterranean Sea. A speleothem record from a cave in northern Turkey that tracks the isotopic signature of Black Sea surface water suggests an open connection to the Mediterranean Sea in at least twelve periods in the past 670,000 years.

The role of northern Arabian Sea surface temperature biases in CMIP5 model simulations and future projections of Indian summer monsoon rainfall

Many climate models have problems simulating Indian summer monsoon rainfall and its variability, resulting in considerable uncertainty in future projections. Problems may relate to many factors, such as local effects of the formulation of physical parametrisation schemes, while common model biases that develop elsewhere within the climate system may also be important. Here we examine the extent and impact of cold sea surface temperature (SST) biases developing in the northern Arabian Sea in the CMIP5 multi-model ensemble, where such SST biases are shown to be common. Such biases have previously been shown to reduce monsoon rainfall in the Met Office Unified Model (MetUM) by weakening moisture fluxes incident upon India. The Arabian Sea SST biases in CMIP5 models consistently develop in winter, via strengthening of the winter monsoon circulation, and persist into spring and summer. A clear relationship exists between Arabian Sea cold SST bias and weak monsoon rainfall in CMIP5 models, similar to effects in the MetUM. Part of this effect may also relate to other factors, such as forcing of the early monsoon by spring-time excessive equatorial precipitation. Atmosphere-only future time-slice experiments show that Arabian Sea cold SST biases have potential to weaken future monsoon rainfall increases by limiting moisture flux acceleration through non-linearity of the Clausius-Clapeyron relationship. Analysis of CMIP5 model future scenario simulations suggests that, while such effects are likely small compared to other sources of uncertainty, models with large Arabian Sea cold SST biases suppress the range of potential outcomes for changes to future early monsoon rainfall.

Incorporation of a physically based melt pond scheme into the sea ice component of a climate model

The extent and thickness of the Arctic sea ice cover has decreased dramatically in the past few decades with minima in sea ice extent in September 2005 and 2007. These minima have not been predicted in the IPCC AR4 report, suggesting that the sea ice component of climate models should more realistically represent the processes controlling the sea ice mass balance. One of the processes poorly represented in sea ice models is the formation and evolution of melt ponds. Melt ponds accumulate on the surface of sea ice from snow and sea ice melt and their presence reduces the albedo of the ice cover, leading to further melt. Toward the end of the melt season, melt ponds cover up to 50% of the sea ice surface. We have developed a melt pond evolution theory. Here, we have incorporated this melt pond theory into the Los Alamos CICE sea ice model, which has required us to include the refreezing of melt ponds. We present results showing that the presence, or otherwise, of a representation of melt ponds has a significant effect on the predicted sea ice thickness and extent. We also present a sensitivity study to uncertainty in the sea ice permeability, number of thickness categories in the model representation, meltwater redistribution scheme, and pond albedo. We conclude with a recommendation that our melt pond scheme is included in sea ice models, and the number of thickness categories should be increased and concentrated at lower thicknesses.

Effect of shear rupture on aggregate scale formation in sea ice

A discrete element model is used to study shear rupture of sea ice under convergent wind stresses. The model includes compressive, tensile, and shear rupture of viscous elastic joints connecting floes that move under the action of the wind stresses. The adopted shear rupture is governed by Coulomb’s criterion. The ice pack is a 400 km long square domain consisting of 4 km size floes. In the standard case with tensile strength 10 times smaller than the compressive strength, under uniaxial compression the failure regime is mainly shear rupture with the most probable scenario corresponding to that with the minimum failure work. The orientation of cracks delineating formed aggregates is bimodal with the peaks around the angles given by the wing crack theory determining diamond-shaped blocks. The ice block (floe aggregate) size decreases as the wind stress gradient increases since the elastic strain energy grows faster leading to a higher speed of crack propagation. As the tensile strength grows, shear rupture becomes harder to attain and compressive failure becomes equally important leading to elongation of blocks perpendicular to the compression direction and the blocks grow larger. In the standard case, as the wind stress confinement ratio increases the failure mode changes at a confinement ratio within 0.2–0.4, which corresponds to the analytical critical confinement ratio of 0.32. Below this value, the cracks are bimodal delineating diamond shape aggregates, while above this value failure becomes isotropic and is determined by small-scale stress anomalies due to irregularities in floe shape.

A model of the three-dimensional evolution of Arctic melt ponds on first-year and multiyear sea ice

During winter the ocean surface in polar regions freezes over to form sea ice. In the summer the upper layers of sea ice and snow melts producing meltwater that accumulates in Arctic melt ponds on the surface of sea ice. An accurate estimate of the fraction of the sea ice surface covered in melt ponds is essential for a realistic estimate of the albedo for global climate models. We present a melt-pond–sea-ice model that simulates the three-dimensional evolution of melt ponds on an Arctic sea ice surface. The advancements of this model compared to previous models are the inclusion of snow topography; meltwater transport rates are calculated from hydraulic gradients and ice permeability; and the incorporation of a detailed one-dimensional, thermodynamic radiative balance. Results of model runs simulating first-year and multiyear sea ice are presented. Model results show good agreement with observations, with duration of pond coverage, pond area, and ice ablation comparing well for both the first-year ice and multiyear ice cases. We investigate the sensitivity of the melt pond cover to changes in ice topography, snow topography, and vertical ice permeability. Snow was found to have an important impact mainly at the start of the melt season, whereas initial ice topography strongly controlled pond size and pond fraction throughout the melt season. A reduction in ice permeability allowed surface flooding of relatively flat, first-year ice but had little impact on the pond coverage of rougher, multiyear ice. We discuss our results, including model shortcomings and areas of experimental uncertainty.

Sea ice rheology

The polar oceans of Earth are covered by sea ice. On timescales much greater than a day, the motion and deformation of the sea ice cover (i.e., its dynamics) are primarily determined by atmospheric and oceanic tractions on its upper and lower surfaces and by internal ice forces that arise within the ice cover owing to its deformation. This review discusses the relationship between the internal ice forces and the deformation of the ice cover, focusing on representations suitable for inclusion within global climate models. I first draw attention to theories that treat the sea ice cover as an isotropic continuum and then to the recent development of anisotropic models that deal with the presence of oriented weaknesses in the ice cover, known as leads.

Impact of melt ponds on Arctic sea ice simulations from 1990 to 2007

The extent and thickness of the Arctic sea ice cover has decreased dramatically in the past few decades with minima in sea ice extent in September 2007 and 2011 and climate models did not predict this decline. One of the processes poorly represented in sea ice models is the formation and evolution of melt ponds. Melt ponds form on Arctic sea ice during the melting season and their presence affects the heat and mass balances of the ice cover, mainly by decreasing the value of the surface albedo by up to 20%. We have developed a melt pond model suitable for forecasting the presence of melt ponds based on sea ice conditions. This model has been incorporated into the Los Alamos CICE sea ice model, the sea ice component of several IPCC climate models. Simulations for the period 1990 to 2007 are in good agreement with observed ice concentration. In comparison to simulations without ponds, the September ice volume is nearly 40% lower. Sensitivity studies within the range of uncertainty reveal that, of the parameters pertinent to the present melt pond parameterization and for our prescribed atmospheric and oceanic forcing, variations of optical properties and the amount of snowfall have the strongest impact on sea ice extent and volume. We conclude that melt ponds will play an increasingly important role in the melting of the Arctic ice cover and their incorporation in the sea ice component of Global Circulation Models is essential for accurate future sea ice forecasts.

Impact of a new anisotropic rheology on simulations of Arctic sea ice

new rheology that explicitly accounts for the subcontinuum anisotropy of the sea ice cover is implemented into the Los Alamos sea ice model. This is in contrast to all models of sea ice included in global circulation models that use an isotropic rheology. The model contains one new prognostic variable, the local structure tensor, which quantifies the degree of anisotropy of the sea ice, and two parameters that set the time scale of the evolution of this tensor. The anisotropic rheology provides a subcontinuum description of the mechanical behavior of sea ice and accounts for a continuum scale stress with large shear to compression ratio and tensile stress component. Results over the Arctic of a stand-alone version of the model are presented and anisotropic model sensitivity runs are compared with a reference elasto-visco-plastic simulation. Under realistic forcing sea ice quickly becomes highly anisotropic over large length scales, as is observed from satellite imagery. The influence of the new rheology on the state and dynamics of the sea ice cover is discussed. Our reference anisotropic run reveals that the new rheology leads to a substantial change of the spatial distribution of ice thickness and ice drift relative to the reference standard visco-plastic isotropic run, with ice thickness regionally increased by more than 1 m, and ice speed reduced by up to 50%.

Systematic winter sea-surface temperature biases in the northern Arabian Sea in HiGEM and the CMIP3 models

Analysis of 20th century simulations of the High resolution Global Environment Model (HiGEM) and the Third Coupled Model Intercomparison Project (CMIP3) models shows that most have a cold sea-surface temperature (SST) bias in the northern Arabian Sea during boreal winter. The association between Arabian Sea SST and the South Asian monsoon has been widely studied in observations and models, with winter cold biases known to be detrimental to rainfall simulation during the subsequent monsoon in coupled general circulation models (GCMs). However, the causes of these SST biases are not well understood. Indeed this is one of the first papers to address causes of the cold biases. The models show anomalously strong north-easterly winter monsoon winds and cold air temperatures in north-west India, Pakistan and beyond. This leads to the anomalous advection of cold, dry air over the Arabian Sea. The cold land region is also associated with an anomalously strong meridional surface temperature gradient during winter, contributing to the enhanced low-level convergence and excessive precipitation over the western equatorial Indian Ocean seen in many models.