Flux transfer events at the magnetopause and in the ionosphere

On December 1, 1986 the ISEE 1 and 2 spacecraft pair passed through the dayside magnetopause at a location which mapped approximately to ionospheric field-line foot-points near the fields of view of the EISCAT radar and photometers and an all-sky camera on Svalbard. The magnetosheath magnetic field was southward and duskward at the time, and flux transfer events (FTEs) were observed at the ISEE location. At the same time, the EISCAT radar observed ionospheric flow bursts of up to 1 km s−1. The peak of each burst followed an FTE observation at ISEE by a few minutes. The bursts, each lasting ten or fifteen minutes, were comprised of first a westward then a poleward flow. An all-sky camera at Ny Ålesund observed dayside auroral breakup forms during or shortly after the flow bursts, moving westward then poleward. While these flow bursts and associated dayside auroral forms have been previously reported in association with southward IMF orientations, this is the first observation of a direct link to FTEs at the magnetopause. On this occasion, the lower limit on the inferred potential associated with the FTEs is roughly 10 kV. Their inferred east-west extent in the ionosphere ranges between 700 and 1000 km, corresponding to a 3 – 5 RE local time extent at the average magnetopause.

The ionospheric signatures of flux transfer events and solar wind dynamic pressure changes

The generation of flow and current vortices in the dayside auroral ionosphere has been predicted for two processes ocurring at the dayside magnetopause. The first of these mechanisms is time-dependent magnetic reconnection, in “flux transfer events” (FTEs); the second is the action of solar wind dynamic pressure changes. The ionospheric flow signature of an FTE should be a twin vortex, with the mean flow velocity in the central region of the pattern equal to the velocity of the pattern as a whole. On the other hand, a pulse of enhanced or reduced dynamic pressure is also expected to produce a twin vortex, but with the central plasma flow being generally different in speed from, and almost orthogonal to, the motion of the whole pattern. In this paper, we make use of this distinction to discuss recent observations of vortical flow patterns in the dayside auroral ionosphere in terms of one or other of the proposed mechanisms. We conclude that some of the observations reported are consistent only with the predicted signature of FTEs. We then evaluate the dimensions of the open flux tubes required to explain some recent simultaneous radar and auroral observations and infer that they are typically 300 km in north–south extent but up to 2000 km in longitudinal extent (i.e., roughly 5 hours of MLT). Hence these observations suggest that recent theories of FTEs which invoke time-varying reconnection at an elongated neutral line may be correct. We also present some simultaneous observations of the interplanetary magnetic field (IMF) and solar wind dynamic pressure (observed using the IMP8 satellite) and the ionospheric flow (observed using the EISCAT radar) which are also only consistent with the FTE model. We estimate that for continuously southward IMF ( ≈ 5 nT) these FTEs contribute about 30 kV to the mean total transpolar voltage (∼30%).

Low-altitude signatures of the cusp and flux transfer events

The usual interpretation of a flux transfer event (FTE) at the magnetopause, in terms of time-dependent and possibly patchy reconnection, demands that it generate an ionospheric signature. Recent ground-based observations have revealed that auroral transients in the cusp/cleft region have all the characteristics required of FTE effects. However, signatures in the major available dataset, namely that from low-altitude polar-orbiting satellites, have not yet been identified. In this paper, we consider a cusp pass of the DE-2 spacecraft during strongly southward IMF. The particle detectors show magnetosheath ion injection signatures. However, the satellite motion and convection are opposed, and we discuss how the observed falling energy dispersion of the precipitating ions can have arisen from a static, moving or growing source. The spatial scale of the source is typical of an FTE. A simple model of the ionospheric signature of an FTE reproduces the observed electric and magnetic field perturbations. Precipitating electrons of peak energy ∼100eV are found to lie on the predicted boundary of the newly-opened tube, very similar to those found on the edges of FTEs at the magnetopause. The injected ions are within this boundary and their dispersion is consistent with its growth as reconnection proceeds. The reconnection potential and the potential of the induced ionospheric motion are found to be the same (≃25kV). The scanning imager on DE-1 shows a localised transient auroral feature around DE-2 at this time, similar to the recent optical/radar observations of FTEs.

Dayside auroral activity and magnetic flux transfer from the solar wind

Combined observations by meridian-scanning photometers and the EISCAT radar show that the "midday-auroral breakup" phenomenon is associated with major increases in ionospheric flow. A sequence of nine events is observed in the early afternoon MLT sector during a period when the IMF is strongly southward with a large positive By component. Each auroral structure is seen at both 630 and 557.7nm and initially moves westward, accompanied by an increase in potential of 30-60kV across the north-south dimension of the EISCAT field-of-view. After a few minutes the arc (or arc fragment) moves into the polar cap and fades, and the velocities observed by the radar swing from westward toward northward. We conclude that dayside auroral breakup is closely associated with momentum transfer across the magnetopause which occurs in a series of events 5-15 minutes apart. The largest of the observed events has dimensions of about 300km (in the direction of westward motion) by 700km, is bounded on its poleward edge by a 5kR arc and is associated with a potential of at least 80kV.

Ionospheric ion upwelling in the wake of flux transfer events at the dayside magnetopause

The effects of flux transfer events (FTE) on the dayside auroral ionosphere are studied, using a simple twin-vortex model of induced ionospheric plasma flow. It is shown that the predicted and observed velocities of these flows are sufficient to drive nonthermal plasma in the F region, not only within the newly opened flux tube of the FTE, but also on the closed, or "old" open, field lines around it. In fact, with the expected poleward neutral wind, the plasma is more highly nonthermal on the flanks of, but outside, the open flux tube: EISCAT observations indicate that plasma is indeed driven into nonthermal distributions in these regions. The nonthermal plasma is thereby subject to additional upforce due to the resulting ion temperature anisotropy and transient expansion due to Joule heating and also to ion accelerations associated with the FTE field aligned current system. Any upflows produced on closed field lines in the vicinity of the FTE are effectively bunched-up in the "wake" of the FTE. Observations from the AMPTE-UKS satellite at the magnetopause reveal ion upflows of energy ∼100 eV flowing out from the ionosphere on closed field lines which are only found in the wake of the FTE. Such flows are also only found shortly after two, out of all the FTEs observed by AMPTE-UKS. The outflow from the ionosphere is two orders of magnitude greater than predicted for the "classical" polar wind. It is shown that such ionospheric ion flows are only expected in association with FTEs on the magnetopause which are well removed from the sub-solar point-either towards dusk or, as in the UKS example discussed here, towards dawn. It is suggested that such ionospheric ions will only be observed if the center of the FTE open flux tube passes very close to the satellite. Consequently, we conclude the ion upflows presented here are probably driven by the second of two possible source FTEs and are observed at the satellite with a lag after the FTE which is less than their time-of-flight.

A process-based analysis of ocean heat uptake in an AOGCM with an eddy-permitting ocean component

About 90% of the anthropogenic increase in heat stored in the climate system is found the oceans. Therefore it is relevant to understand the details of ocean heat uptake. Here we present a detailed, process-based analysis of ocean heat uptake (OHU) processes in HiGEM1.2, an atmosphere-ocean general circulation model (AOGCM) with an eddy-permitting ocean component of 1/3 degree resolution. Similarly to various other models, HiGEM1.2 shows that the global heat budget is dominated by a downward advection of heat compensated by upward isopycnal diffusion. Only in the upper tropical ocean do we find the classical balance between downward diapycnal diffusion and upward advection of heat. The upward isopycnal diffusion of heat is located mostly in the Southern Ocean, which thus dominates the global heat budget. We compare the responses to a 4xCO2 forcing and an enhancement of the windstress forcing in the Southern Ocean. This highlights the importance of regional processes for the global ocean heat uptake. These are mainly surface fluxes and convection in the high latitudes, and advection in the Southern Ocean mid-latitudes. Changes in diffusion are less important. In line with the CMIP5 models, HiGEM1.2 shows a band of strong OHU in the mid-latitude Southern Ocean in the 4xCO2 run, which is mostly advective. By contrast, in the high-latitude Southern Ocean regions it is the suppression of convection that leads to OHU. In the enhanced windstress run, convection is strengthened at high Southern latitudes, leading to heat loss, while the magnitude of the OHU in the Southern mid-latitudes is very similar to the 4xCO2 results. Remarkably, there is only very small global OHU in the enhanced windstress run. The wind stress forcing just leads to a redistribution of heat. We relate the ocean changes at high southern latitudes to the effect of climate change on the Antarctic Circumpolar Current (ACC). It weakens in the 4xCO2 run and strengthens in the wind stress run. The weakening is due to a narrowing of the ACC, caused by an expansion of the Weddell Gyre, and a flattening of the isopycnals, which are explained by a combination of the wind stress forcing and increased precipitation.

Probabilistic downscaling of remote sensing data with applications for multi-scale biogeochemical flux modeling

Upscaling ecological information to larger scales in space and downscaling remote sensing observations or model simulations to finer scales remain grand challenges in Earth system science. Downscaling often involves inferring subgrid information from coarse-scale data, and such ill-posed problems are classically addressed using regularization. Here, we apply two-dimensional Tikhonov Regularization (2DTR) to simulate subgrid surface patterns for ecological applications. Specifically, we test the ability of 2DTR to simulate the spatial statistics of high-resolution (4 m) remote sensing observations of the normalized difference vegetation index (NDVI) in a tundra landscape. We find that the 2DTR approach as applied here can capture the major mode of spatial variability of the high-resolution information, but not multiple modes of spatial variability, and that the Lagrange multiplier (γ) used to impose the condition of smoothness across space is related to the range of the experimental semivariogram. We used observed and 2DTR-simulated maps of NDVI to estimate landscape-level leaf area index (LAI) and gross primary productivity (GPP). NDVI maps simulated using a γ value that approximates the range of observed NDVI result in a landscape-level GPP estimate that differs by ca 2% from those created using observed NDVI. Following findings that GPP per unit LAI is lower near vegetation patch edges, we simulated vegetation patch edges using multiple approaches and found that simulated GPP declined by up to 12% as a result. 2DTR can generate random landscapes rapidly and can be applied to disaggregate ecological information and compare of spatial observations against simulated landscapes.

Quantifying particle size and turbulent scale dependence of dust uplift in the Sahara using aircraft measurements

The first size-resolved airborne measurements of dust fluxes and the first dust flux measurements from the central Sahara are presented and compared with a parameterization by Kok (2011a). High-frequency measurements of dust size distribution were obtained from 0.16 to 300 µm diameter, and eddy covariance fluxes were derived. This is more than an order of magnitude larger size range than previous flux estimates. Links to surface emission are provided by analysis of particle drift velocities. Number flux is described by a −2 power law between 1 and 144 µm diameter, significantly larger than the 12 µm upper limit suggested by Kok (2011a). For small particles, the deviation from a power law varies with terrain type and the large size cutoff is correlated with atmospheric vertical turbulent kinetic energy, suggesting control by vertical transport rather than emission processes. The measured mass flux mode is in the range 30–100 µm. The turbulent scales important for dust flux are from 0.1 km to 1–10 km. The upper scale increases during the morning as boundary layer depth and eddy size increase. All locations where large dust fluxes were measured had large topographical variations. These features are often linked with highly erodible surface features, such as wadis or dunes. We also hypothesize that upslope flow and flow separation over such features enhance the dust flux by transporting large particles out of the saltation layer. The tendency to locate surface flux measurements in open, flat terrain means these favored dust sources have been neglected in previous studies.

Separating the influence of projected changes in air temperature and wind on patterns of sea level change and ocean heat content

We present ocean model sensitivity experiments aimed at separating the influence of the projected changes in the “thermal” (near-surface air temperature) and “wind” (near-surface winds) forcing on the patterns of sea level and ocean heat content. In the North Atlantic, the distribution of sea level change is more due to the “thermal” forcing, whereas it is more due to the “wind” forcing in the North Pacific; in the Southern Ocean, the “thermal” and “wind” forcing have a comparable influence. In the ocean adjacent to Antarctica the “thermal” forcing leads to an inflow of warmer waters on the continental shelves, which is somewhat attenuated by the “wind” forcing. The structure of the vertically integrated heat uptake is set by different processes at low and high latitudes: at low latitudes it is dominated by the heat transport convergence, whereas at high latitudes it represents a small residual of changes in the surface flux and advection of heat. The structure of the horizontally integrated heat content tendency is set by the increase of downward heat flux by the mean circulation and comparable decrease of upward heat flux by the subgrid-scale processes; the upward eddy heat flux decreases and increases by almost the same magnitude in response to, respectively, the “thermal” and “wind” forcing. Regionally, the surface heat loss and deep convection weaken in the Labrador Sea, but intensify in the Greenland Sea in the region of sea ice retreat. The enhanced heat flux anomaly in the subpolar Atlantic is mainly caused by the “thermal” forcing.

The relationship between the ITCZ and the Southern Hemispheric eddy-driven jet

We study the effect of a thermal forcing confined to the midlatitudes of one hemisphere on the eddy-driven jet in the opposite hemisphere. We demonstrate the existence of an “interhemispheric teleconnection,” whereby warming (cooling) the Northern Hemisphere causes both the intertropical convergence zone (ITCZ) and the Southern Hemispheric midlatitude jet to shift northward (southward). The interhemispheric teleconnection is effected by a change in the asymmetry of the Hadley cells: as the ITCZ shifts away from the Equator, the cross-equatorial Hadley cell intensifies, fluxing more momentum toward the subtropics and sustaining a stronger subtropical jet. Changes in subtropical jet strength, in turn, alter the propagation of extratropical waves into the tropics, affecting eddy momentum fluxes and the eddy-driven westerlies. The relevance of this mechanism is demonstrated in the context of future climate change simulations, where shifts of the ITCZ are significantly related to shifts of the Southern Hemispheric eddy-driven jet in austral winter. The possible relevance of the proposed mechanism to paleoclimates is discussed, particularly with regard to theories of ice age terminations.

The response of the Southern Hemispheric eddy-driven jet to future changes in shortwave radiation in CMIP5

A strong relationship is found between changes in the meridional gradient of absorbed shortwave radiation (ASR) and Southern Hemispheric jet shifts in 21st century climate simulations of CMIP5 (Coupled Model Intercomparison Project phase 5) coupled models. The relationship is such that models with increases in the meridional ASR gradient around the southern midlatitudes, and therefore increases in midlatitude baroclinicity, tend to produce a larger poleward jet shift. The ASR changes are shown to be dominated by changes in cloud properties, with sea ice declines playing a secondary role. We demonstrate that the ASR changes are the cause, and not the result, of the intermodel differences in jet response by comparing coupled simulations with experiments in which sea surface temperature increases are prescribed. Our results highlight the importance of reducing the uncertainty in cloud feedbacks in order to constrain future circulation changes.

CO2 emission estimation in the urban environment: measurement of the CO2 storage term

Eddy covariance has been used in urban areas to evaluate the net exchange of CO2 between the surface and the atmosphere. Typically, only the vertical flux is measured at a height 2–3 times that of the local roughness elements; however, under conditions of relatively low instability, CO2 may accumulate in the airspace below the measurement height. This can result in inaccurate emissions estimates if the accumulated CO2 drains away or is flushed upwards during thermal expansion of the boundary layer. Some studies apply a single height storage correction; however, this requires the assumption that the response of the CO2 concentration profile to forcing is constant with height. Here a full seasonal cycle (7th June 2012 to 3rd June 2013) of single height CO2 storage data calculated from concentrations measured at 10 Hz by open path gas analyser are compared to a data set calculated from a concurrent switched vertical profile measured (2 Hz, closed path gas analyser) at 10 heights within and above a street canyon in central London. The assumption required for the former storage determination is shown to be invalid. For approximately regular street canyons at least one other measurement is required. Continuous measurements at fewer locations are shown to be preferable to a spatially dense, switched profile, as temporal interpolation is ineffective. The majority of the spectral energy of the CO2 storage time series was found to be between 0.001 and 0.2 Hz (500 and 5 s respectively); however, sampling frequencies of 2 Hz and below still result in significantly lower CO2 storage values. An empirical method of correcting CO2 storage values from under-sampled time series is proposed.

Investigation of the impact of anthropogenic heat flux within an urban land surface model and PILPS-urban

Results from the first international urban model comparison experiment (PILPS-Urban) suggested that models which neglected the anthropogenic heat flux within the surface energy balance performed at least as well as models that include the source term, but this could not be explained. The analyses undertaken show that the results from PILPS-Urban were masked by the signal from including vegetation, which was identified in PILPS-Urban as being important. Including the anthropogenic heat flux does give improved performance, but the benefit is small for the site studied given the relatively small magnitude of this flux relative to other terms in the surface energy balance. However, there is no further benefit from including temporal variations in the flux at this site. The importance is expected to increase at sites with a larger anthropogenic heat flux and greater temporal variations.