Evaluation of the sea ice simulation in a new coupled atmosphere-ocean climate model (HadGEM1)

[ 1] A rapid increase in the variety, quality, and quantity of observations in polar regions is leading to a significant improvement in the understanding of sea ice dynamic and thermodynamic processes and their representation in global climate models. We assess the simulation of sea ice in the new Hadley Centre Global Environmental Model (HadGEM1) against the latest available observations. The HadGEM1 sea ice component uses elastic-viscous-plastic dynamics, multiple ice thickness categories, and zero-layer thermodynamics. The model evaluation is focused on the mean state of the key variables of ice concentration, thickness, velocity, and albedo. The model shows good agreement with observational data sets. The variability of the ice forced by the North Atlantic Oscillation is also found to agree with observations.

The response of the sea ice edge to atmospheric and oceanic jet formation

The sea ice edge presents a region of many feedback processes between the atmosphere, ocean, and sea ice (Maslowski et al.). Here the authors focus on the impact of on-ice atmospheric and oceanic flows at the sea ice edge. Mesoscale jet formation due to the Coriolis effect is well understood over sharp changes in surface roughness such as coastlines (Hunt et al.). This sharp change in surface roughness is experienced by the atmosphere and ocean encountering a compacted sea ice edge. This paper presents a study of a dynamic sea ice edge responding to prescribed atmospheric and oceanic jet formation. An idealized analytical model of sea ice drift is developed and compared to a sea ice climate model [the Los Alamos Sea Ice Model (CICE)] run on an idealized domain. The response of the CICE model to jet formation is tested at various resolutions. It is found that the formation of atmospheric jets at the sea ice edge increases the wind speed parallel to the sea ice edge and results in the formation of a sea ice drift jet in agreement with an observed sea ice drift jet (Johannessen et al.). The increase in ice drift speed is dependent upon the angle between the ice edge and wind and results in up to a 40% increase in ice transport along the sea ice edge. The possibility of oceanic jet formation and the resultant effect upon the sea ice edge is less conclusive. Observations and climate model data of the polar oceans have been analyzed to show areas of likely atmospheric jet formation, with the Fram Strait being of particular interest.

Regina Maris and the command of the sea: the sixteenth century origins of modern maritime strategy

The concept of the command of the sea has its roots in medieval notions of the sovereignty of coastal waters, as claimed by several monarchs and polities of Europe. In the sixteenth century, a surge of intellectual creativity, especially in Elizabethan England, fused this notion with the Thucydidean term ‘thalassocracy’ – the rule of the sea. In the light of the explorations of the oceans, this led to a new conceptualisation of naval warfare, developed in theory and then put into practice. This falsifies the mistaken but widespread assumption that there was no significant writing on naval strategy before the nineteenth century.

Sea-breeze modification of the growth of a marine internal boundary layer

A numerical mesoscale model is used to make a high-resolution simulation of the marine boundary layer in the Persian Gulf, during conditions of offshore flow from Saudi Arabia. A marine internal boundary layer (MIBL) and a sea-breeze circulation (SBC) are found to co-exist. The sea breeze develops in the mid-afternoon, at which time its front is displaced several tens of kilometres offshore. Between the coast and the sea-breeze system, the MIBL that occurs is consistent with a picture described in the existing literature. However, the MIBL is perturbed by the SBC, the boundary layer deepening significantly seaward of the sea-breeze front. Our analysis suggests that this strong, localized deepening is not a direct consequence of frontal uplift, but rather that the immediate cause is the retardation of the prevailing, low-level offshore wind by the SBC. The simulated boundary-layer development can be accounted for by using a simple 1D Lagrangian model of growth driven by the surface heat flux. This model is obtained as a straightforward modification of an established MIBL analytic growth model.

The possibility of lowering oil content of potato chips by combining atmospheric frying with postfrying vacuum application

Four protocols involving the application of low pressures, either toward the end of frying or after frying, were investigated with the aim of lowering the oil content of potato chips. Protocol 1 involving frying at atmospheric pressure followed by a 3 min draining time constituted the control. Protocol 2 involved lowering of pressure to 13.33 kPa, 40 s before the end of frying, followed by draining for 3 min at the same pressure. Protocol 3 was the same as protocol 2, except that the pressure was lowered 3 s before the end of frying. Protocol 4 involved lowering the pressure to 13.33 kPa after the product was lifted from the oil and holding it at this value over the draining time of 3 min. Protocol 4 gave a product having the lowest oil content (37.12 g oil/100 g defatted dry matter), while protocol 2 gave the product with highest oil content (71.10 g oil/100 g defatted dry matter), followed by those obtained using protocols 1 and 3(68.48 g oil/100 g defatted dry matter and 52.50 g oil/100 g defatted dry matter, respectively). Protocol 4 was further evaluated to study the effects of draining times and vacuum applied, and compared with the control. It was noted that over the modest range of pressures investigated, there was no significant effect of the vacuum applied on the oil content of the product. This study demonstrates that the oil content of potato chips can be lowered significantly by combining atmospheric frying with draining under vacuum.

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.

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.

Intercomparison of the Arctic sea ice cover in global ocean–sea ice reanalyses from the ORA-IP project

Ocean–sea ice reanalyses are crucial for assessing the variability and recent trends in the Arctic sea ice cover. This is especially true for sea ice volume, as long-term and large scale sea ice thickness observations are inexistent. Results from the Ocean ReAnalyses Intercomparison Project (ORA-IP) are presented, with a focus on Arctic sea ice fields reconstructed by state-of-the-art global ocean reanalyses. Differences between the various reanalyses are explored in terms of the effects of data assimilation, model physics and atmospheric forcing on properties of the sea ice cover, including concentration, thickness, velocity and snow. Amongst the 14 reanalyses studied here, 9 assimilate sea ice concentration, and none assimilate sea ice thickness data. The comparison reveals an overall agreement in the reconstructed concentration fields, mainly because of the constraints in surface temperature imposed by direct assimilation of ocean observations, prescribed or assimilated atmospheric forcing and assimilation of sea ice concentration. However, some spread still exists amongst the reanalyses, due to a variety of factors. In particular, a large spread in sea ice thickness is found within the ensemble of reanalyses, partially caused by the biases inherited from their sea ice model components. Biases are also affected by the assimilation of sea ice concentration and the treatment of sea ice thickness in the data assimilation process. An important outcome of this study is that the spatial distribution of ice volume varies widely between products, with no reanalysis standing out as clearly superior as compared to altimetry estimates. The ice thickness from systems without assimilation of sea ice concentration is not worse than that from systems constrained with sea ice observations. An evaluation of the sea ice velocity fields reveals that ice drifts too fast in most systems. As an ensemble, the ORA-IP reanalyses capture trends in Arctic sea ice area and extent relatively well. However, the ensemble can not be used to get a robust estimate of recent trends in the Arctic sea ice volume. Biases in the reanalyses certainly impact the simulated air–sea fluxes in the polar regions, and questions the suitability of current sea ice reanalyses to initialize seasonal forecasts.

Numerical modeling of the propagation environment in the atmospheric boundary layer over the Persian Gulf

Strong vertical gradients at the top of the atmospheric boundary layer affect the propagation of electromagnetic waves and can produce radar ducts. A three-dimensional, time-dependent, nonhydrostatic numerical model was used to simulate the propagation environment in the atmosphere over the Persian Gulf when aircraft observations of ducting had been made. A division of the observations into high- and low-wind cases was used as a framework for the simulations. Three sets of simulations were conducted with initial conditions of varying degrees of idealization and were compared with the observations taken in the Ship Antisubmarine Warfare Readiness/Effectiveness Measuring (SHAREM-115) program. The best results occurred with the initialization based on a sounding taken over the coast modified by the inclusion of data on low-level atmospheric conditions over the Gulf waters. The development of moist, cool, stable marine internal boundary layers (MIBL) in air flowing from land over the waters of the Gulf was simulated. The MIBLs were capped by temperature inversions and associated lapses of humidity and refractivity. The low-wind MIBL was shallower and the gradients at its top were sharper than in the high-wind case, in agreement with the observations. Because it is also forced by land–sea contrasts, a sea-breeze circulation frequently occurs in association with the MIBL. The size, location, and internal structure of the sea-breeze circulation were realistically simulated. The gradients of temperature and humidity that bound the MIBL cause perturbations in the refractivity distribution that, in turn, lead to trapping layers and ducts. The existence, location, and surface character of the ducts were well captured. Horizontal variations in duct characteristics due to the sea-breeze circulation were also evident. The simulations successfully distinguished between high- and low-wind occasions, a notable feature of the SHAREM-115 observations. The modeled magnitudes of duct depth and strength, although leaving scope for improvement, were most encouraging.

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%.

A preliminary investigation of the release of oil from contaminated clays and other fine-grained sediments during consolidation

Fine-grained sediments on land, or in a freshwater or marine environment, may become contaminated with a wide range of pollutants including hydrocarbons. This paper is concerned with preliminary studies of the mobilization and transportation of hydrocarbons, during the process of consolidation, to adjacent sediments or water bodies. A modified Rowe Cell was used to measure the consolidation properties of prepared kaolinite and bentonite clay-water slurries, with and without the addition of oil, along with hydrocarbon-bearing drill-cuttings samples taken from the sea-bed adjacent to two North Sea oil-well platforms. The consolidation properties of the kaolinite and bentonite clay slurries were little altered by the addition of oil, which was present at concentrations of between 8073 and 59 572 mg kg(-1). During each consolidation stage, samples of the expelled pore-fluids were collected and analysed for oil content. These values were very low in comparison with the original oil concentration in the samples and changed little between each consolidation stage. Analysis of the slurry samples both before and after consolidation confirms that, proportionally, little oil is removed as a result of consolidation. The implication of these results is that, for the range of samples tested, the very high hydraulic gradients and particle rearrangements that occur during the process of consolidation are capable of releasing only proportionally small amounts of oil bound to the fine-grained clay and silt particles.

A preliminary investigation of the release of oil from contaminated clays and other fine-grained sediments during consolidation

Fine-grained sediments on land, or in a freshwater or marine environment, may become contaminated with a wide range of pollutants including hydrocarbons. This paper is concerned with preliminary studies of the mobilization and transportation of hydrocarbons, during the process of consolidation, to adjacent sediments or water bodies. A modified Rowe Cell was used to measure the consolidation properties of prepared kaolinite and bentonite clay-water slurries, with and without the addition of oil, along with hydrocarbon-bearing drill-cuttings samples taken from the sea-bed adjacent to two North Sea oil-well platforms. The consolidation properties of the kaolinite and bentonite clay slurries were little altered by the addition of oil, which was present at concentrations of between 8073 and 59 572 mg kg(-1). During each consolidation stage, samples of the expelled pore-fluids were collected and analysed for oil content. These values were very low in comparison with the original oil concentration in the samples and changed little between each consolidation stage. Analysis of the slurry samples both before and after consolidation confirms that, proportionally, little oil is removed as a result of consolidation. The implication of these results is that, for the range of samples tested, the very high hydraulic gradients and particle rearrangements that occur during the process of consolidation are capable of releasing only proportionally small amounts of oil bound to the fine-grained clay and silt particles.

Assessing the date of the global oil peak: the need to use 2P reserves

Combining geological knowledge with proved plus probable ('2P') oil discovery data indicates that over 60 countries are now past their resource-limited peak of conventional oil production. The data show that the global peak of conventional oil production is close. Many analysts who rely only on proved ('1P') oil reserves data draw a very different conclusion. But proved oil reserves contain no information about the true size of discoveries, being variously under-reported, over-reported and not reported. Reliance on 1P data has led to a number of misconceptions, including the notion that past oil forecasts were incorrect, that oil reserves grow very significantly due to technology gain, and that the global supply of oil is ensured provided sufficient investment is forthcoming to 'turn resources into reserves'. These misconceptions have been widely held, including within academia, governments, some oil companies, and organisations such as the IEA. In addition to conventional oil, the world contains large quantities of non-conventional oil. Most current detailed models show that past the conventional oil peak the non-conventional oils are unlikely to come on-stream fast enough to offset conventional's decline. To determine the extent of future oil supply constraints calculations are required to determine fundamental rate limits for the production of non-conventional oils, as well as oil from gas, coal and biomass, and of oil substitution. Such assessments will need to examine technological readiness and lead-times, as well as rate constraints on investment, pollution, and net-energy return. (C) 2007 Elsevier Ltd. All rights reserved.