Reconnection at the high-latitude magnetopause during northward interplanetary magnetic field conditions

The Polar spacecraft had a prolonged encounter with the high-latitude dayside magnetopause on May 29, 1996. This encounter with the magnetopause occurred when the interplanetary magnetic field was directed northward. From the three-dimensional electron and ion distribution functions measured by the Hydra instrument, it has been possible to identify nearly all of the distinct boundary layer regions associated with high-latitude reconnection. The regions that have been identified are (1) the cusp; (2) the magnetopause current layer; (3) magnetosheath field lines that have interconnected in only the Northern Hemisphere; (4) magnetosheath field lines that have interconnected in only the Southern Hemisphere; (5) magnetosheath field lines that have interconnected in both the Northern and Southern Hemispheres; (6) magnetosheath that is disconnected from the terrestrial magnetic field; and (7) high-latitude plasma sheet field lines that are participating in magnetosheath reconnection. Reconnection over this time period was occurring at high latitudes over a broad local-time extent, interconnecting the magnetosheath and lobe and/or plasma sheet field lines in both the Northern and Southern Hemispheres. Newly closed boundary layer field lines were observed as reconnection occur-red first at high latitudes in one hemisphere and then later in the other. These observations establish the location of magnetopause reconnection during these northward interplanetary magnetic field conditions as being at high latitudes, poleward of the cusp, and further reinforce the general interpretation of electron and ion phase space density signatures as indicators of magnetic reconnection and boundary layer formation.

Evidence of component merging equatorward of the cusp

The Polar spacecraft passed through a region near the dayside magnetopause on May 29, 1996, at a geocentric distance of similar to 8 R-E and high, northern magnetic latitudes. The interplanetary magnetic field (IMF) was northward during the pass. Data from the Thermal Ion Dynamics Experiment revealed the existence of low-speed (similar to 50 km s(-1)) ion D-shaped distributions mixed with cold ions (similar to 2 eV) over a period of 2.5 hours. These ions were traveling parallel to the magnetic field toward the Northern Hemisphere ionosphere and were convecting primarily eastward. The D-shaped distributions are distinct from a convecting Maxwellian and, along with the magnetic field direction, are taken as evidence that the spacecraft was inside the magnetosphere and not in the magnetosheath. Furthermore, the absence of ions in the antiparallel direction is taken as evidence that low-shear merging was occurring at a location southward of the spacecraft and equatorward of the Southern Hemisphere cusp. The cold ions were of ionospheric origin, with initially slow field-aligned speeds, which were accelerated upon reflection from the magnetopause. These observations provide significant new evidence consistent with component magnetic merging sites equatorward of the cusp for northward IMF.

Ionospheric and geomagnetic responses to changes in IMF B_Z: a superposed epoch study

Superposed epoch studies have been carried out in order to determine the ionospheric response at mid-latitudes to southward turnings of the interplanetary magnetic field (IMF). This is compared with the geomagnetic response, as seen in the indices K p, AE and Dst. The solar wind, IMF and geomagnetic data used were hourly averages from the years 1967–1989 and thus cover a full 22-year cycle in the solar magnetic field. These data were divided into subsets, determined by the magnitudes of the southward turnings and the concomitant increase in solar wind pressure. The superposed epoch studies were carried out using the time of the southward turning as time zero. The response of the mid-latitude ionosphere is studied by looking at the F-layer critical frequencies, f o F2, from hourly soundings by the Slough ionosonde and their deviation from the monthly median values, δf o F2. For the southward turnings with a change in B z of δB z > 11.5 nT accompanied by a solar wind dynamic pressure P exceeding 5 nPa, the F region critical frequency, f o F2, shows a marked decrease, reaching a minimum value about 20 h after the southward turning. This recovers to pre-event values over the subsequent 24 h, on average. The Dst index shows the classic storm-time decrease to about −60 nT. Four days later, the index has still to fully recover and is at about −25 nT. Both the K p and AE indices show rises before the southward turnings, when the IMF is strongly northward but the solar wind dynamic pressure is enhanced. The average AE index does register a clear isolated pulse (averaging 650 nT for 2 h, compared with a background peak level of near 450 nT at these times) showing enhanced energy deposition at high latitudes in substorms but, like K p, remains somewhat enhanced for several days, even after the average IMF has returned to zero after 1 day. This AE background decays away over several days as the Dst index recovers, indicating that there is some contamination of the currents observed at the AE stations by the continuing enhanced equatorial ring current. For data averaged over all seasons, the critical frequencies are depressed at Slough by 1.3 MHz, which is close to the lower decile of the overall distribution of δf o Fl values. Taking 30-day periods around summer and winter solstice, the largest depression is 1.6 and 1.2 MHz, respectively. This seasonal dependence is confirmed by a similar study for a Southern Hemisphere station, Argentine Island, giving peak depressions of 1.8 MHz and 0.5 MHz for summer and winter. For the subset of turnings where δB z > 11.5 nT and P ≤ 5 nPa, the response of the geomagnetic indices is similar but smaller, while the change in δf o F2 has all but disappeared. This confirms that the energy deposited at high latitudes, which leads to the geomagnetic and ionospheric disturbances following a southward turning of the IMF, increases with the energy density (dynamic pressure) of the solar wind flow. The magnitude of all responses are shown to depend on δB z . At Slough, the peak depression always occurs when Slough rotates into the noon sector. The largest ionospheric response is for southward turnings seen between 15–21 UT.

Dependence of convective flows and particle precipitation in the high-latitude dayside ionosphere on theXandYcomponents of the interplanetary magnetic field

The asymmetries in the convective flows, current systems, and particle precipitation in the high-latitude dayside ionosphere which are related to the equatorial plane components of the interplanetary magnetic field (IMF) are discussed in relation to the results of several recent observational studies. It is argued that all of the effects reported to date which are ascribed to the y component of the IMF can be understood, at least qualitatively, in terms of a simple theoretical picture in which the effects result from the stresses exerted on the magnetosphere consequent on the interconnection of terrestrial and interplanetary fields. In particular, relaxation under the action of these stresses allows, in effect, a partial penetration of the IMF into the magnetospheric cavity, such that the sense of the expected asymmetry effects on closed field lines can be understood, to zeroth order, in terms of the “dipole plus uniform field” model. In particular, in response to IMF By, the dayside cusp should be displaced in longitude about noon in the same sense as By in the northern hemisphere, and in the opposite sense to By in the southern hemisphere, while simultaneously the auroral oval as a whole should be shifted in the dawn-dusk direction in the opposite sense with respect to By. These expected displacements are found to be consistent with recently published observations. Similar considerations lead to the suggestion that the auroral oval may also undergo displacements in the noon-midnight direction which are associated with the x component of the IMF. We show that a previously published study of the position of the auroral oval contains strong initial evidence for the existence of this effect. However, recent results on variations in the latitude of the cusp are more ambiguous. This topic therefore requires further study before definitive conclusions can be drawn.

An assessment of the impact of local processes on dust lifting in martian climate models

Simulation of the lifting of dust from the planetary surface is of substantially greater importance on Mars than on Earth, due to the fundamental role that atmospheric dust plays in the former’s climate, yet the dust emission parameterisations used to date in martian global climate models (MGCMs) lag, understandably, behind their terrestrial counterparts in terms of sophistication. Recent developments in estimating surface roughness length over all martian terrains and in modelling atmospheric circulations at regional to local scales (less than O(100 km)) presents an opportunity to formulate an improved wind stress lifting parameterisation. We have upgraded the conventional scheme by including the spatially varying roughness length in the lifting parameterisation in a fully consistent manner (thereby correcting a possible underestimation of the true threshold level for wind stress lifting), and used a modification to account for deviations from neutral stability in the surface layer. Following these improvements, it is found that wind speeds at typical MGCM resolution never reach the lifting threshold at most gridpoints: winds fall particularly short in the southern midlatitudes, where mean roughness is large. Sub-grid scale variability, manifested in both the near-surface wind field and the surface roughness, is then considered, and is found to be a crucial means of bridging the gap between model winds and thresholds. Both forms of small-scale variability contribute to the formation of dust emission ‘hotspots’: areas within the model gridbox with particularly favourable conditions for lifting, namely a smooth surface combined with strong near-surface gusts. Such small-scale emission could in fact be particularly influential on Mars, due both to the intense positive radiative feedbacks that can drive storm growth and a strong hysteresis effect on saltation. By modelling this variability, dust lifting is predicted at the locations at which dust storms are frequently observed, including the flushing storm sources of Chryse and Utopia, and southern midlatitude areas from which larger storms tend to initiate, such as Hellas and Solis Planum. The seasonal cycle of emission, which includes a double-peaked structure in northern autumn and winter, also appears realistic. Significant increases to lifting rates are produced for any sensible choices of parameters controlling the sub-grid distributions used, but results are sensitive to the smallest scale of variability considered, which high-resolution modelling suggests should be O(1 km) or less. Use of such models in future will permit the use of a diagnosed (rather than prescribed) variable gustiness intensity, which should further enhance dust lifting in the southern hemisphere in particular.

Possible impacts of a future grand solar minimum on climate: stratospheric and global circulation changes

It has been suggested that the Sun may evolve into a period of lower activity over the 21st century. This study examines the potential climate impacts of the onset of an extreme ‘Maunder Minimum like’ grand solar minimum using a comprehensive global climate model. Over the second half of the 21st century, the scenario assumes a decrease in total solar irradiance of 0.12% compared to a reference RCP8.5 experiment. The decrease in solar irradiance cools the stratopause (~1 hPa) in the annual and global mean by 1.4 K. The impact on global mean near-surface temperature is small (~−0.1 K), but larger changes in regional climate occur during the stratospheric dynamically active seasons. In Northern hemisphere (NH) winter-time, there is a weakening of the stratospheric westerly jet by up to ~3-4 m s1, with the largest changes occurring in January-February. This is accompanied by a deepening of the Aleutian low at the surface and an increase in blocking over northern Europe and the north Pacific. There is also an equatorward shift in the Southern hemisphere (SH) midlatitude eddy-driven jet in austral spring. The occurrence of an amplified regional response during winter and spring suggests a contribution from a top-down pathway for solar-climate coupling; this is tested using an experiment in which ultraviolet (200–320 nm) radiation is decreased in isolation of other changes. The results show that a large decline in solar activity over the 21st century could have important impacts on the stratosphere and regional surface climate.

Observational constraints on atmospheric and oceanic cross-equatorial heat transports: revisiting the precipitation asymmetry problem in climate models

Satellite based top-of-atmosphere (TOA) and surface radiation budget observations are combined with mass corrected vertically integrated atmospheric energy divergence and tendency from reanalysis to infer the regional distribution of the TOA, atmospheric and surface energy budget terms over the globe. Hemispheric contrasts in the energy budget terms are used to determine the radiative and combined sensible and latent heat contributions to the cross-equatorial heat transports in the atmosphere (AHT_EQ) and ocean (OHT_EQ). The contrast in net atmospheric radiation implies an AHT_EQ from the northern hemisphere (NH) to the southern hemisphere (SH) (0.75 PW), while the hemispheric difference in sensible and latent heat implies an AHT_EQ in the opposite direction (0.51 PW), resulting in a net NH to SH AHT_EQ (0.24 PW). At the surface, the hemispheric contrast in the radiative component (0.95 PW) dominates, implying a 0.44 PW SH to NH OHT_EQ. Coupled model intercomparison project phase 5 (CMIP5) models with excessive net downward surface radiation and surface-to-atmosphere sensible and latent heat transport in the SH relative to the NH exhibit anomalous northward AHT_EQ and overestimate SH tropical precipitation. The hemispheric bias in net surface radiative flux is due to too much longwave surface radiative cooling in the NH tropics in both clear and all-sky conditions and excessive shortwave surface radiation in the SH subtropics and extratropics due to an underestimation in reflection by clouds.

The solsticial pause on Mars: 2 modelling and investigation of causes

The martian solsticial pause, presented in a companion paper (Lewis et al., this issue), was investigated further through a series of model runs using the UK version of the LMD/UK Mars Global Climate Model. It was found that the pause could not be adequately reproduced if radiatively active water ice clouds were omitted from the model. When clouds were used, along with a realistic time-dependent dust opacity distribution, a substantial minimum in near-surface transient eddy activity formed around solstice in both hemispheres. The net effect of the clouds in the model is, by altering the thermal structure of the atmosphere, to decrease the vertical shear of the westerly jet near the surface around solstice, and thus reduce baroclinic growth rates. A similar effect was seen under conditions of large dust loading, implying that northern midlatitude eddy activity will tend to become suppressed after a period of intense flushing storm formation around the northern cap edge. Suppression of baroclinic eddy generation by the barotropic component of the flow and via diabatic eddy dissipation were also investigated as possible mechanisms leading to the formation of the solsticial pause but were found not to make major contributions. Zonal variations in topography were found to be important, as their presence results in weakened transient eddies around winter solstice in both hemispheres, through modification of the near-surface flow. The zonal topographic asymmetry appears to be the primary reason for the weakness of eddy activity in the southern hemisphere relative to the northern hemisphere, and the ultimate cause of the solsticial pause in both hemispheres. The meridional topographic gradient was found to exert a much weaker influence on near-surface transient eddies.

The solsticial pause on Mars: 1. A planetary wave reanalysis

Large-scale planetary waves are diagnosed from an analysis of profiles retrieved from the Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft during its scientific mapping phase. The analysis is conducted by assimilating thermal profiles and total dust opacity retrievals into a Mars global circulation model. Transient waves are largest throughout the northern hemisphere autumn, winter and spring period and almost absent during the summer. The southern hemisphere exhibits generally weaker transient wave behaviour. A striking feature of the low-altitude transient waves in the analysis is that they show a broad subsidiary minimum in amplitude centred on the winter solstice, a period when the thermal contrast between the summer hemisphere and the winter pole is strongest and baroclinic wave activity might be expected to be strong. This behaviour, here called the ‘solsticial pause,’ is present in every year of the analysis. This strong pause is under-represented in many independent model experiments, which tend to produce relatively uniform baroclinic wave activity throughout the winter. This paper documents and diagnoses the transient wave solsticial pause found in the analysis; a companion paper investigates the origin of the phenomenon in a series of model experiments.

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.

Climatology and variability of trace gases in extratropical double-tropopause regions from MLS, HIRDLS, and ACE-FTS measurements

Upper tropospheric and lower stratospheric measurements from the Aura Microwave Limb Sounder (MLS), the Aura High Resolution Dynamics Limb Sounder (HIRDLS), and the Atmospheric Chemistry Experiment-Fourier transform spectrometer (ACE-FTS) are used to present the first global climatological comparison of extratropical, nonpolar trace gas distributions in double-tropopause (DT) and single-tropopause (ST) regions. Stratospheric tracers, O3, HNO3, and HCl, have lower mixing ratios ∼2–8 km above the primary (lowermost) tropopause in DT than in ST regions in all seasons, with maximum Northern Hemisphere (NH) differences near 50% in winter and 30% in summer. Southern Hemisphere winter differences are somewhat smaller, but summer differences are similar in the two hemispheres. H2O in DT regions of both hemispheres shows strong negative anomalies in November through February and positive anomalies in July through October, reflecting the strong seasonal cycle in H2O near the tropical tropopause. CO and other tropospheric tracers examined have higher DT than ST values 2–7 km above the primary tropopause, with the largest differences in winter. Large DT-ST differences extend to high NH latitudes in fall and winter, with longitudinal maxima in regions associated with enhanced wave activity and subtropical jet variations. Results for O3 and HNO3 agree closely between MLS and HIRDLS, and differences from ACE-FTS are consistent with its sparse and irregular midlatitude sampling. Consistent signatures in climatological trace gas fields provide strong evidence that transport from the tropical upper troposphere into the layer between double tropopauses is an important pathway for stratosphere-troposphere exchange.

The influence of sea ice dynamics on the climate sensitivity and memory to increased Antarctic sea ice

The study analyzes the sensitivity and memory of the Southern Hemisphere coupled climate system to increased Antarctic sea ice (ASI), taking into account the persistence of the sea ice maxima in the current climate. The mechanisms involved in restoring the climate balance under two sets of experiments, which differ in regard to their sea ice models, are discussed. The experiments are perturbed with extremes of ASI and integrated for 10 yr in a large 30-member ensemble. The results show that an ASI maximum is able to persist for ; 4 yr in the current climate, followed by a negative sea ice phase. The sea ice insulating effect during the positive phase reduces heat fluxes south of 60 8 S, while at the same time these are intensified at the sea ice edge. The increased air stability over the sea ice field strengthens the polar cell while the baroclinicity increases at midlatitudes. The mean sea level pressure is reduced (increased) over high latitudes (midlatitudes), typical of the southern annular mode (SAM) positive phase. The Southern Ocean (SO) becomes colder and fresher as the sea ice melts mainly through sea ice lateral melting, the consequence of which is an increase in the ocean stability by buoyancy and mixing changes. The climate sensitivity is triggered by the sea ice insulating process and the resulting freshwater pulse (fast response), while the climate equilibrium is restored by the heat stored in the SO subsurface layers (long response). It is concluded that the time needed for the ASI anomaly to be dissipated and/or melted is shortened by the sea ice dynamical processes.

Glacial lake drainage in Patagonia (13-8 kyr) and response of the adjacent Pacific Ocean

Large freshwater lakes formed in North America and Europe during deglaciation following the Last Glacial Maximum. Rapid drainage of these lakes into the Oceans resulted in abrupt perturbations in climate, including the Younger Dryas and 8.2 kyr cooling events. In the mid-latitudes of the Southern Hemisphere major glacial lakes also formed and drained during deglaciation but little is known about the magnitude, organization and timing of these drainage events and their e ect on regional climate. We use 16 new single-grain optically stimulated luminescence (OSL) dates to de ne three stages of rapid glacial lake drainage in the Lago General Carrera/Lago Buenos Aires and Lago Cohrane/ Pueyrredón basins of Patagonia and provide the rst assessment of the e ects of lake drainage on the Paci c Ocean. Lake drainage occurred between 13 and 8 kyr ago and was initially gradual eastward into the Atlantic, then subsequently reorganized westward into the Paci c as new drainage routes opened up during Patagonian Ice Sheet deglaciation. Coupled ocean-atmosphere model experiments using HadCM3 with an imposed freshwater surface “hosing” to simulate glacial lake drainage suggest that a negative salinity anomaly was advected south around Cape Horn, resulting in brief but signi cant impacts on coastal ocean vertical mixing and regional climate.

ENSO related variation of equatorial MRG and Rossby waves and forcing from higher latitudes

The contrasting behaviour of westward-moving mixed Rossby-gravity (WMRG) and the first Rossby (R1) waves in El Niño (EN) and La Niña (LN) seasons is documented with a focus on the Northern Hemisphere winter. The eastward-moving variance in the upper troposphere is dominated by WMRG and R1 structures that appear to be Doppler-shifted by the flow and are referred to as WMRG-E and R1-E. In the East Pacific and Atlantic the years with stronger equatorial westerly winds have the stronger WMRG and WMRG- E. In the East Pacific, R1 is also a maximum in LN. However, R1-E exhibits an eastward-shift between LN and EN. The changes with ENSO phase provide a test-bed for the understanding of these waves. In the East Pacific and Atlantic, the stronger WMRG-E and WMRG with stronger westerlies are in accord with the dispersion relation with simple Doppler-shifting by the zonal flow. The possible existence of free waves can also explain stronger R1 in EN in the Eastern Hemisphere. 1-D free wave propagation theory based on wave activity conservation is also important for R1. However, this theory is unable to explain the amplitude maxima for other waves observed in the strong equatorial westerly regions in the Western Hemisphere, and certainly not their ENSO-related variation. The forcing of equatorial waves by higher latitude wave activity and its variation with ENSO phase is therefore examined. Propagation of extratropical eastward-moving Rossby wave activity through the westerly ducts into the equatorial region where it triggers WMRG-E is favoured in the stronger westerlies, in LN in the East Pacific and EN in the Atlantic. It is also found that WMRG is forced by Southern Hemisphere westward-moving wavetrains arching into the equatorial region where they are reflected. The most significant mechanism for both R1 and R1-E appear to be lateral forcing by subtropical wavetrains.

Resonant Wave Interactions in the Presence of a Diurnally Varying Heat Source

Resonant interactions among equatorial waves in the presence of a diurnally varying heat source are studied in the context of the diabatic version of the equatorial beta-plane primitive equations for a motionless, hydrostatic, horizontally homogeneous and stably stratified background atmosphere. The heat source is assumed to be periodic in time and of small amplitude [i.e., O(epsilon)] and is prescribed to roughly represent the typical heating associated with deep convection in the tropical atmosphere. In this context, using the asymptotic method of multiple time scales, the free linear Rossby, Kelvin, mixed Rossby-gravity, and inertio-gravity waves, as well as their vertical structures, are obtained as leading-order solutions. These waves are shown to interact resonantly in a triad configuration at the O(e) approximation, and the dynamics of these interactions have been studied in the presence of the forcing. It is shown that for the planetary-scale wave resonant triads composed of two first baroclinic equatorially trapped waves and one barotropic Rossby mode, the spectrum of the thermal forcing is such that only one of the triad components is resonant with the heat source. As a result, to illustrate the role of the diurnal forcing in these interactions in a simplified fashion, two kinds of triads have been analyzed. The first one refers to triads composed of a k = 0 first baroclinic geostrophic mode, which is resonant with the stationary component of the diurnal heat source, and two dispersive modes, namely, a mixed Rossby-gravity wave and a barotropic Rossby mode. The other class corresponds to triads composed of two first baroclinic inertio-gravity waves in which the highest-frequency wave resonates with a transient harmonic of the forcing. The integration of the asymptotic reduced equations for these selected resonant triads shows that the stationary component of the diurnal heat source acts as an ""accelerator"" for the energy exchanges between the two dispersive waves through the excitation of the catalyst geostrophic mode. On the other hand, since in the second class of triads the mode that resonates with the forcing is the most energetically active member because of the energy constraints imposed by the triad dynamics, the results show that the convective forcing in this case is responsible for a longer time scale modulation in the resonant interactions, generating a period doubling in the energy exchanges. The results suggest that the diurnal variation of tropical convection might play an important role in generating low-frequency fluctuations in the atmospheric circulation through resonant nonlinear interactions.