Consistent and contrasting decadal Arctic sea ice thickness predictions from a highly optimized sea ice model

[1] Decadal hindcast simulations of Arctic Ocean sea ice thickness made by a modern dynamic-thermodynamic sea ice model and forced independently by both the ERA-40 and NCEP/NCAR reanalysis data sets are compared for the first time. Using comprehensive data sets of observations made between 1979 and 2001 of sea ice thickness, draft, extent, and speeds, we find that it is possible to tune model parameters to give satisfactory agreement with observed data, thereby highlighting the skill of modern sea ice models, though the parameter values chosen differ according to the model forcing used. We find a consistent decreasing trend in Arctic Ocean sea ice thickness since 1979, and a steady decline in the Eastern Arctic Ocean over the full 40-year period of comparison that accelerated after 1980, but the predictions of Western Arctic Ocean sea ice thickness between 1962 and 1980 differ substantially. The origins of differing thickness trends and variability were isolated not to parameter differences but to differences in the forcing fields applied, and in how they are applied. It is argued that uncertainty, differences and errors in sea ice model forcing sets complicate the use of models to determine the exact causes of the recently reported decline in Arctic sea ice thickness, but help in the determination of robust features if the models are tuned appropriately against observations.

Optimization of a sea ice model using basinwide observations of Arctic sea ice thickness, extent, and velocity

A stand-alone sea ice model is tuned and validated using satellite-derived, basinwide observations of sea ice thickness, extent, and velocity from the years 1993 to 2001. This is the first time that basin-scale measurements of sea ice thickness have been used for this purpose. The model is based on the CICE sea ice model code developed at the Los Alamos National Laboratory, with some minor modifications, and forcing consists of 40-yr ECMWF Re-Analysis (ERA-40) and Polar Exchange at the Sea Surface (POLES) data. Three parameters are varied in the tuning process: Ca, the air–ice drag coefficient; P*, the ice strength parameter; and α, the broadband albedo of cold bare ice, with the aim being to determine the subset of this three-dimensional parameter space that gives the best simultaneous agreement with observations with this forcing set. It is found that observations of sea ice extent and velocity alone are not sufficient to unambiguously tune the model, and that sea ice thickness measurements are necessary to locate a unique subset of parameter space in which simultaneous agreement is achieved with all three observational datasets.

Improving the spatial distribution of modeled Arctic sea ice thickness

The spatial distribution of ice thickness/draft in the Arctic Ocean is examined using a sea ice model. A comparison of model predictions with submarine observations of sea ice draft made during cruises between 1987 and 1997 reveals that the model has the same deficiencies found in previous studies, namely ice that is too thick in the Beaufort Sea and too thin near the North Pole. We find that increasing the large scale shear strength of the sea ice leads to substantial improvements in the model's spatial distribution of sea ice thickness, and simultaneously improves the agreement between modeled and ERS-derived 1993–2001 mean winter ice thickness.

Dependence of sea ice yield-curve shape on ice thickness

In this note, the authors discuss the contribution that frictional sliding of ice floes (or floe aggregates) past each other and pressure ridging make to the plastic yield curve of sea ice. Using results from a previous study that explicitly modeled the amount of sliding and ridging that occurs for a given global strain rate, it is noted that the relative contribution of sliding and ridging to ice stress depends upon ice thickness. The implication is that the shape and size of the plastic yield curve is dependent upon ice thickness. The yield-curve shape dependence is in addition to plastic hardening/weakening that relates the size of the yield curve to ice thickness. In most sea ice dynamics models the yield-curve shape is taken to be independent of ice thickness. The authors show that the change of the yield curve due to a change in the ice thickness can be taken into account by a weighted sum of two thickness-independent rheologies describing ridging and sliding effects separately. It would be straightforward to implement the thickness-dependent yield-curve shape described here into sea ice models used for global or regional ice prediction.

Automated identification and tracking of polar-cap plasma patches at solar minimum

A method of automatically identifying and tracking polar-cap plasma patches, utilising data inversion and feature-tracking methods, is presented. A well-established and widely used 4-D ionospheric imaging algorithm, the Multi-Instrument Data Assimilation System (MIDAS), inverts slant total electron content (TEC) data from ground-based Global Navigation Satellite System (GNSS) receivers to produce images of the free electron distribution in the polar-cap ionosphere. These are integrated to form vertical TEC maps. A flexible feature-tracking algorithm, TRACK, previously used extensively in meteorological storm-tracking studies is used to identify and track maxima in the resulting 2-D data fields. Various criteria are used to discriminate between genuine patches and "false-positive" maxima such as the continuously moving day-side maximum, which results from the Earth's rotation rather than plasma motion. Results for a 12-month period at solar minimum, when extensive validation data are available, are presented. The method identifies 71 separate structures consistent with patch motion during this time. The limitations of solar minimum and the consequent small number of patches make climatological inferences difficult, but the feasibility of the method for patches larger than approximately 500 km in scale is demonstrated and a larger study incorporating other parts of the solar cycle is warranted. Possible further optimisation of discrimination criteria, particularly regarding the definition of a patch in terms of its plasma concentration enhancement over the surrounding background, may improve results.

If the British left: agricultural policy outside the CAP?

If the Conservative party wins Britain’s General Election in 2015 the incoming Government intends to negotiate a new settlement with the EU, and then hold an in-out referendum. Opinion polls suggest a narrow majority would vote to leave. To exercise an informed choice the electorate would need to know what feasible alternatives are available, so that these can be assessed and debated. Furthermore, in the event of a no vote the Government must be ready to implement its alternative strategy: Plan B. This implies that both a new EU settlement, and Plan B, need to be negotiated concurrently. There has been no serious attempt to explain what Plan B would mean for UK farm policy. Would, for example, British farmers continue to receive the level of support that they currently enjoy; would support be more focussed on environmental objectives? WTO rules on most-favoured-nation treatment and regional trade agreements would apply. Could a WTO compatible agri-food trade agreement be negotiated with its former EU partners, or would Irish and Brazilian beef face the same tariff barriers on imports into the British market? For the EU, the loss of a major net contributor to EU finances could result in a re-examination of the CAP budget.

Will Arctic sea ice thickness initialization improve seasonal forecast skill?

Arctic sea ice thickness is thought to be an important predictor of Arctic sea ice extent. However, coupled seasonal forecast systems do not generally use sea ice thickness observations in their initialization and are therefore missing a potentially important source of additional skill. To investigate how large this source is, a set of ensemble potential predictability experiments with a global climate model, initialized with and without knowledge of the sea ice thickness initial state, have been run. These experiments show that accurate knowledge of the sea ice thickness field is crucially important for sea ice concentration and extent forecasts up to 8 months ahead, especially in summer. Perturbing sea ice thickness also has a significant impact on the forecast error in Arctic 2 m temperature a few months ahead. These results suggest that advancing capabilities to observe and assimilate sea ice thickness into coupled forecast systems could significantly increase skill.

Motion of the dayside polar cap boundary during substorm cycles: II. Generation of poleward-moving events and polar cap patches by pulses in the magnetopause reconnection rate

Using data from the EISCAT (European Incoherent Scatter) VHF and CUTLASS (Co-operative UK Twin- Located Auroral Sounding System) HF radars, we study the formation of ionospheric polar cap patches and their relationship to the magnetopause reconnection pulses identified in the companion paper by Lockwood et al. (2005). It is shown that the poleward-moving, high-concentration plasma patches observed in the ionosphere by EISCAT on 23 November 1999, as reported by Davies et al. (2002), were often associated with corresponding reconnection rate pulses. However, not all such pulses generated a patch and only within a limited MLT range (11:00–12:00 MLT) did a patch result from a reconnection pulse. Three proposed mechanisms for the production of patches, and of the concentration minima that separate them, are analysed and evaluated: (1) concentration enhancement within the patches by cusp/cleft precipitation; (2) plasma depletion in the minima between the patches by fast plasma flows; and (3) intermittent injection of photoionisation-enhanced plasma into the polar cap. We devise a test to distinguish between the effects of these mechanisms. Some of the events repeat too frequently to apply the test. Others have sufficiently long repeat periods and mechanism (3) is shown to be the only explanation of three of the longer-lived patches seen on this day. However, effect (2) also appears to contribute to some events. We conclude that plasma concentration gradients on the edges of the larger patches arise mainly from local time variations in the subauroral plasma, via the mechanism proposed by Lockwood et al. (2000).

Motion of the dayside polar cap boundary during substorm cycles: I. Observations of pulses in the magnetopause reconnection rate

Using data from the EISCAT (European Incoherent Scatter) VHF radar and DMSP (Defense Meteorological Satellite Program) spacecraft passes, we study the motion of the dayside open-closed field line boundary during two substorm cycles. The satellite data show that the motions of ion and electron temperature boundaries in EISCAT data, as reported by Moen et al. (2004), are not localised around the radar; rather, they reflect motions of the open-closed field line boundary at all MLT throughout the dayside auroral ionosphere. The boundary is shown to erode equatorward when the IMF points southward, consistent with the effect of magnetopause reconnection. During the substorm expansion and recovery phases, the dayside boundary returns poleward, whether the IMF points northward or southward. However, the poleward retreat was much faster during the substorm for which the IMF had returned to northward than for the substorm for which the IMF remained southward – even though the former substorm is much the weaker of the two. These poleward retreats are consistent with the destruction of open flux at the tail current sheet. Application of a new analysis of the peak ion energies at the equatorward edge of the cleft/cusp/mantle dispersion seen by the DMSP satellites identifies the dayside reconnection merging gap to extend in MLT from about 9.5 to 15.5 h for most of the interval. Analysis of the boundary motion, and of the convection velocities seen near the boundary by EISCAT, allows calculation of the reconnection rate (mapped down to the ionosphere) from the flow component normal to the boundary in its own rest frame. This reconnection rate is not, in general, significantly different from zero before 06:45 UT (MLT<9.5 h) – indicating that the X line footprint expands over the EISCAT field-of-view to earlier MLT only occasionally and briefly. Between 06:45 UT and 12:45UT (9.5

Ground-based observations of the auroral zone and polar cap ionospheric responses to dayside transient reconnection

Observations from the EISCAT VHF incoherent scatter radar system in northern Norway, during a run of the common programme CP-4, reveal a series of polewardpropagating F-region electron density enhancements in the pre-noon sector on 23 November 1999. These plasma density features, which are observed under conditions of a strongly southward interplanetary magnetic field, exhibit a recurrence rate of under 10 min and appear to emanate from the vicinity of the open/closed field-line boundary from where they travel into the polar cap; this is suggestive of their being an ionospheric response to transient reconnection at the dayside magnetopause (flux transfer events). Simultaneous with the density structures detected by the VHF radar, polewardmoving radar auroral forms (PMRAFs) are observed by the Finland HF coherent scatter radar. It is thought that PMRAFs, which are commonly observed near local noon by HF radars, are also related to flux transfer events, although the specific mechanism for the generation of the field-aligned irregularities within such features is not well understood. The HF observations suggest, that for much of their existence, the PMRAFs trace fossil signatures of transient reconnection rather than revealing the footprint of active reconnection itself; this is evidenced not least by the fact that the PMRAFs become narrower in spectral width as they evolve away from the region of more classical, broad cusp scatter in which they originate. Interpretation of the HF observations with reference to the plasma parameters diagnosed by the incoherent scatter radar suggests that as the PMRAFs migrate away from the reconnection site and across the polar cap, entrained in the ambient antisunward flow, the irregularities therein are generated by the presence of gradients in the electron density, with these gradients having been formed through structuring of the ionosphere in the cusp region in response to transient reconnection.

A Summary of the NATO ASI on Polar Cap Boundary Phenomena

The polar cap boundary is a subject of central importance to current magnetosphere-ionosphere research and its applications in “space weather” activities. The problems are that it has a number of definitions, and that the most physically meaningful definition (namely the open-closed field line boundary) is very difficult to identify in observations. New understanding of the importance of the structure and dynamics of the boundary region made the time right for a meeting reviewing our knowledge in this area. The Advanced Study Institute (ASI) on Svalbard in June 1997 discussed the boundary on both the dayside and the nightside, mapping magnetically to the dayside magnetopause and to tail plasma sheet/lobe interface, respectively. We held a “brainstorming” session, in which different ideas which arose from the presented papers were discussed and developed, and a summary session, in which session convenors gave a personal view of progress that has been made and problems which still need solving. Both were designed as ways of promoting further discussion. This paper attempts to distil some of the themes that emerged from these discussions.

The response of ionospheric convection in the polar cap to substorm activity

We report multi-instrument observations during an isolated substorm on 17 October 1989. The EISCAT radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71°-78°. SAMNET and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. IMP-8 magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux to the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71° by the time of the expansion phase onset. A westward electrojet, centred at 65.4°, occurred at the onset of the expansion phase. This electrojet subsequently moved poleward to a maximum of 68.1° at 2000 UT and also widened. During the expansion phase, there is evidence of bursts of plasma flow which are spatially localised at longitudes within the substorm current wedge and which occurred well poleward of the westward electrojet. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the “distant” neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase. It is not obvious whether the electrojet mapped to a near-Earth neutral line, but at its most poleward, the expanded electrojet does not reach the estimated latitude of the polar cap boundary.

Rapid changes in LLBL thickness

We employ data from two spacecraft, at the dawn flank of the magnetopause, to investigate fluctuations in the thickness of the low-latitude boundary layer (LLBL). We show the LLBL is considerably thinner shortly after the detection of a flux transfer event than it was during the event. These data are shown to be consistent with the theory of transient increases in the open LLBL thickness caused by a pulse of enhanced reconnection at the magnetopause.

Ion flows and heating at a contracting polar-cap boundary: GISMOS observations indicating viscous-like interaction on the flanks of the magnetotail

This paper complements that in this issue by Clauer et al. concerning the international GISMOS campaign of 3–5 June 1987. From a detailed study of the EISCAT data, the polar-cap boundary, as defined by an almost shear east-west convection reversal, is found to contract across the EISCAT field of view between 04 and 07 MLT. An annulus of enhanced ion temperature and non-thermal plasma is observed immediately equatorward of the contracting boundary due to the lag in the response of the neutral-wind pattern to the change in ion flows. The ion flow inside the polar cap and at the boundary is shown to be relatively smooth, compared with that in the auroral oval, at 15-second resolution. The flow at the boundary is directed poleward, with velocities which exceed that of the boundary itself. The effect of velocity shears on the beamswinging technique used to derive the ion flows has been analysed in detail and it is found that spurious flows across a moving boundary can be generated. However, these are much smaller than the observed flows into the polar cap and cannot explain the 7 kV potential difference across the observed segment of the cap boundary between 04:30–06:30 UT. The ion temperature enhancements at the two observing azimuths is used to define the boundary orientation. The results are consistent with recent observations of slow anti-sunward flow of closed field lines on the flanks of the geomagnetic tail, which appears to be generated by some form of “viscous” coupling to the magnetosheath plasma.