Meteorological factors controlling low-level continental pollutant outflow across a coast

Coastal outflow describes the horizontal advection of pollutants from the continental boundary layer across a coastline into a layer above the marine boundary layer. This process can ventilate polluted continental boundary layers and thus regulate air quality in highly populated coastal regions. This paper investigates the factors controlling coastal outflow and quantifies its importance as a ventilation mechanism. Tracers in the Met Office Unified Model (MetUM) are used to examine the magnitude and variability of coastal outflow over the eastern United States for a 4 week period during summer 2004. Over the 4 week period, ventilation of tracer from the continental boundary layer via coastal outflow occurs with the same magnitude as vertical ventilation via convection and advection. The relative importance of tracer decay rate, cross-coastal advection rate, and a parameter based on the relative continental and marine boundary layer heights, on coastal outflow is assessed by reducing the problem to a time-dependent box-model. The ratio of the advection rate and decay rate is a dimensionless parameter which determines whether tracers are long-lived or short-lived. Long- and short-lived tracers exhibit different behaviours with respect to coastal outflow. For short-lived tracers, increasing the advection rate increases the diurnally averaged magnitude of coastal outflow, but has the opposite effect for very long-lived tracers. Short-lived tracers exhibit large diurnal variability in coastal outflow but long-lived tracers do not. By combining the MetUM and box-model simulations a landwidth is determined which represents the distance inland over which emissions contribute significantly to coastal outflow. A landwidth of between 100 and 400 km is found to be representative for a tracer with a lifetime of 24 h.

A frequency-independent boundary element method for scattering by two-dimensional screens and apertures

We propose and analyse a hybrid numerical–asymptotic hp boundary element method (BEM) for time-harmonic scattering of an incident plane wave by an arbitrary collinear array of sound-soft two-dimensional screens. Our method uses an approximation space enriched with oscillatory basis functions, chosen to capture the high-frequency asymptotics of the solution. We provide a rigorous frequency-explicit error analysis which proves that the method converges exponentially as the number of degrees of freedom N increases, and that to achieve any desired accuracy it is sufficient to increase N in proportion to the square of the logarithm of the frequency as the frequency increases (standard BEMs require N to increase at least linearly with frequency to retain accuracy). Our numerical results suggest that fixed accuracy can in fact be achieved at arbitrarily high frequencies with a frequency-independent computational cost, when the oscillatory integrals required for implementation are computed using Filon quadrature. We also show how our method can be applied to the complementary ‘breakwater’ problem of propagation through an aperture in an infinite sound-hard screen.

Investigation of discrepancies in satellite rainfall estimates over Ethiopia

Tropical Applications of Meteorology Using Satellite and Ground-Based Observations (TAMSAT) rainfall estimates are used extensively across Africa for operational rainfall monitoring and food security applications; thus, regional evaluations of TAMSAT are essential to ensure its reliability. This study assesses the performance of TAMSAT rainfall estimates, along with the African Rainfall Climatology (ARC), version 2; the Tropical Rainfall Measuring Mission (TRMM) 3B42 product; and the Climate Prediction Center morphing technique (CMORPH), against a dense rain gauge network over a mountainous region of Ethiopia. Overall, TAMSAT exhibits good skill in detecting rainy events but underestimates rainfall amount, while ARC underestimates both rainfall amount and rainy event frequency. Meanwhile, TRMM consistently performs best in detecting rainy events and capturing the mean rainfall and seasonal variability, while CMORPH tends to overdetect rainy events. Moreover, the mean difference in daily rainfall between the products and rain gauges shows increasing underestimation with increasing elevation. However, the distribution in satellite–gauge differences demon- strates that although 75% of retrievals underestimate rainfall, up to 25% overestimate rainfall over all eleva- tions. Case studies using high-resolution simulations suggest underestimation in the satellite algorithms is likely due to shallow convection with warm cloud-top temperatures in addition to beam-filling effects in microwave- based retrievals from localized convective cells. The overestimation by IR-based algorithms is attributed to nonraining cirrus with cold cloud-top temperatures. These results stress the importance of understanding re- gional precipitation systems causing uncertainties in satellite rainfall estimates with a view toward using this knowledge to improve rainfall algorithms.

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

A numerical model of the ionospheric signatures of time-varying magnetic reconnection: II. Measuring expansions in the ionospheric flow response

A numerical model embodying the concepts of the Cowley-Lockwood (Cowley and Lockwood, 1992, 1997) paradigm has been used to produce a simple Cowley– Lockwood type expanding flow pattern and to calculate the resulting change in ion temperature. Cross-correlation, fixed threshold analysis and threshold relative to peak are used to determine the phase speed of the change in convection pattern, in response to a change in applied reconnection. Each of these methods fails to fully recover the expansion of the onset of the convection response that is inherent in the simulations. The results of this study indicate that any expansion of the convection pattern will be best observed in time-series data using a threshold which is a fixed fraction of the peak response. We show that these methods used to determine the expansion velocity can be used to discriminate between the two main models for the convection response to a change in reconnection.

A comparison between ion characteristics observed by the POLAR and DMSP spacecraft in the high-latitude magnetosphere

We study here the injection and transport of ions in the convection-dominated region of the Earth’s magnetosphere. The total ion counts from the CAMMICE MICS instrument aboard the POLAR spacecraft are used to generate occurrence probability distributions of magnetospheric ion populations. MICS ion spectra are characterised by both the peak in the differential energy flux, and the average energy of ions striking the detector. The former permits a comparison with the Stubbs et al. (2001) survey of He2+ ions of solar wind origin within the magnetosphere. The latter can address the occurrences of various classifications of precipitating particle fluxes observed in the topside ionosphere by DMSP satellites (Newell and Meng, 1992). The peak energy occurrences are consistent with our earlier work, including the dawn-dusk asymmetry with enhanced occurrences on the dawn flank at low energies, switching to the dusk flank at higher energies. The differences in the ion energies observed in these two studies can be explained by drift orbit effects and acceleration processes at the magnetopause, and in the tail current sheet. Near noon at average ion energies of _1 keV, the cusp and open LLBL occur further poleward here than in the Newell and Meng survey, probably due to convectionrelated time-of-flight effects. An important new result is that the pre-noon bias previously observed in the LLBL is most likely due to the component of this population on closed field lines, formed largely by low energy ions drifting earthward from the tail. There is no evidence here of mass and momentum transfer from the solar wind to the LLBL by nonreconnection coupling. At higher energies (_2–20 keV), we observe ions mapping to the auroral oval and can distinguish between the boundary and central plasma sheets. We show that ions at these energies relate to a transition from dawnward to duskward dominated flow, this is evidence of how ion drift orbits in the tail influence the location and behaviour of the plasma populations in the magnetosphere.

Secondary cyclogenesis along an occluded front leading to damaging wind gusts: windstorm Kyrill, January 2007

Windstorm Kyrill affected large parts of Europe in January 2007 and caused widespread havoc and loss of life. In this study the formation of a secondary cyclone, Kyill II, along the occluded front of the mature cyclone Kyrill and the occurrence of severe wind gusts as Kyrill II passed over Germany are investigated with the help of high-resolution regional climate model simulations. Kyrill underwent an explosive cyclogenesis south of Greenland as the storm crossed polewards of an intense upper-level jet stream. Later in its life cycle secondary cyclogenesis occurred just west of the British Isles. The formation of Kyrill II along the occluded front was associated (a) with frontolytic strain and (b) with strong diabatic heating in combination with a developing upper-level shortwave trough. Sensitivity studies with reduced latent heat release feature a similar development but a weaker secondary cyclone, revealing the importance of diabatic processes during the formation of Kyrill II. Kyrill II moved further towards Europe and its development was favored by a split jet structure aloft, which maintained the cyclone’s exceptionally deep core pressure (below 965 hPa) for at least 36 hours. The occurrence of hurricane force winds related to the strong cold front over North and Central Germany is analyzed using convection-permitting simulations. The lower troposphere exhibits conditional instability, a turbulent flow and evaporative cooling. Simulation at high spatio-temporal resolution suggests that the downward mixing of high momentum (the wind speed at 875 hPa widely exceeded 45 m s-1) accounts for widespread severe surface wind gusts, which is in agreement with observed widespread losses.

Ground-based and satellite observations of high-latitude auroral activity in the dusk sector of the auroral oval

On 7 December 2000, during 13:30-15:30 UT the MIRACLE all-sky camera at Ny Alesund observed auroras at high-latitudes (MLAT similar to 76) simultaneously when the Cluster spacecraft were skimming the magnetopause in the same MLT sector (at similar to 16:00-18:00 MLT). The location of the auroras (near the ionospheric convection reversal boundary) and the clear correlation between their dynamics and IMF variations suggests their close relationship with R1 currents. Consequently, we can assume that the Cluster spacecraft were making observations in the magnetospheric region associated with the auroras, although exact magnetic conjugacy between the ground-based and satellite observations did not exist. The solar wind variations appeared to control both the behaviour of the auroras and the magnetopause dynamics. Auroral structures were observed at Ny Alesund especially during periods of negative IMF B-Z. In addition, the Cluster spacecraft experienced periodic (T similar to 4 - 6 min) encounters between magnetospheric and magnetosheath plasmas. These undulations of the boundary can be interpreted as a consequence of tailward propagating magnetopause surface waves. Simultaneous dusk sector ground-based observations show weak, but discernible magnetic pulsations (Pc 5) and occasionally periodic variations (T - 2 - 3 min) in the high-latitude auroras. In the dusk sector, Pc 5 activity was stronger and had characteristics that were consistent with a field line resonance type of activity. When IMF BZ stayed positive for a longer period, the auroras were dimmer and the spacecraft stayed at the outer edge of the magnetopause where they observed electromagnetic pulsations with T similar to 1 min. We find these observations interesting especially from the viewpoint of previously presented studies relating poleward-moving high-latitude auroras with pulsation activity and MHD waves propagating at the magnetospheric boundary layers.

Coordinated ground-based, low altitude satellite and Cluster observations on global and local scales during a transient post-noon sector excursion of the magnetospheric cusp

On 14 January 2001, the four Cluster spacecraft passed through the northern magnetospheric mantle in close conjunction to the EISCAT Svalbard Radar (ESR) and approached the post-noon dayside magnetopause over Greenland between 13:00 and 14:00 UT During that interval, a sudden reorganisation of the high-latitude dayside convection pattern accurred after 13:20 UT most likely caused by a direction change of the Solar wind magnetic field. The result was an eastward and poleward directed flow-channel, as monitored by the SuperDARN radar network and also by arrays of ground-based magnetometers in Canada, Greenland and Scandinavia. After an initial eastward and later poleward expansion of the flow-channel between 13:20 and 13:40 UT, the four Cluster spacecraft, and the field line footprints covered by the eastward looking scan cycle of the Sondre Stromfjord incoherent scatter radar were engulfed by cusp-like precipitation with transient magnetic and electric field signatures. In addition, the EISCAT Svalbard Radar detected strong transient effects of the convection reorganisation, a poleward moving precipitation, and a fast ion flow-channel in association with the auroral structures that suddenly formed to the west and north of the radar. From a detailed analysis of the coordinated Cluster and ground-based data, it was found that this extraordinary transient convection pattern, indeed, had moved the cusp precipitation from its former pre-noon position into the late post-noon sector, allowing for the first and quite unexpected encounter of the cusp by the Cluster spacecraft. Our findings illustrate the large amplitude of cusp dynamics even in response to moderate solar wind forcing. The global ground-based data proves to be an invaluable tool to monitor the dynamics and width of the affected magnetospheric regions.

Coordinated Cluster, ground-based instrumentation and low-altitude satellite observations of transient poleward-moving events in the ionosphere and in the tail lobe

During the interval between 8:00-9:30 on 14 January 2001, the four Cluster spacecraft were moving from the central magnetospheric lobe, through the dusk sector mantle, on their way towards intersecting the magnetopause near 15:00 MLT and 15:00 UT. Throughout this interval, the EIS-CAT Svalbard Radar (ESR) at Longyearbyen observed a series of poleward-moving transient events of enhanced F-region plasma concentration ("polar cap patches"), with a repetition period of the order of 10 min. Allowing for the estimated solar wind propagation delay of 75 ( 5) min, the interplanetary magnetic field (IMF) had a southward component during most of the interval. The magnetic footprint of the Cluster spacecraft, mapped to the ionosphere using the Tsyganenko T96 model (with input conditions prevailing during this event), was to the east of the ESR beams. Around 09:05 UT, the DMSP-F12 satellite flew over the ESR and showed a sawtooth cusp ion dispersion signature that also extended into the electrons on the equatorward edge of the cusp, revealing a pulsed magnetopause reconnection. The consequent enhanced ionospheric flow events were imaged by the SuperDARN HF backscatter radars. The average convection patterns (derived using the AMIE technique on data from the magnetometers, the EISCAT and SuperDARN radars, and the DMSP satellites) show that the associated poleward-moving events also convected over the predicted footprint of the Cluster spacecraft. Cluster observed enhancements in the fluxes of both electrons and ions. These events were found to be essentially identical at all four spacecraft, indicating that they had a much larger spatial scale than the satellite separation of the order of 600 km. Some of the events show a correspondence between the lowest energy magnetosheath electrons detected by the PEACE instrument on Cluster (10-20 eV) and the topside ionospheric enhancements seen by the ESR (at 400-700 km). We suggest that a potential barrier at the magnetopause, which prevents the lowest energy electrons from entering the magnetosphere, is reduced when and where the boundary-normal magnetic field is enhanced and that the observed polar cap patches are produced by the consequent enhanced precipitation of the lowest energy electrons, making them and the low energy electron precipitation fossil remnants of the magnetopause reconnection rate pulses.

Long-term variations in the magnetic fields of the Sun and the heliosphere: their origin, effects, and implications

Recent studies of the variation of geomagnetic activity over the past 140 years have quantified the "coronal source" magnetic flux F-s that leaves the solar atmosphere and enters the heliosphere and have shown that it has risen, on average, by an estimated 34% since 1963 and by 140% since 1900. This variation of open solar flux has been reproduced by Solanki et al. [2000] using a model which demonstrates how the open flux accumulates and decays, depending on the rate of flux emergence in active regions and on the length of the solar cycle. We here use a new technique to evaluate solar cycle length and find that it does vary in association with the rate of change of F-s in the way predicted. The long-term variation of the rate of flux emergence is found to be very similar in form to that in F-s, which may offer a potential explanation of why F-s appears to be a useful proxy for extrapolating solar total irradiance back in time. We also find that most of the variation of cosmic ray fluxes incident on Earth is explained by the strength of the heliospheric field (quantified by F-s) and use observations of the abundance of the isotope Be-10 (produced by cosmic rays and deposited in ice sheets) to study the decrease in F-s during the Maunder minimum. The interior motions at the base of the convection zone, where the solar dynamo is probably located, have recently been revealed using the helioseismology technique and found to exhibit a 1.3-year oscillation. This periodicity is here reported in observations of the interplanetary magnetic field and geomagnetic activity but is only present after 1940, When present, it shows a strong 22-year variation, peaking near the maximum of even-numbered sunspot cycles and showing minima at the peaks of odd-numbered cycles. We discuss the implications of these long-term solar and heliospheric variations for Earth's environment.

Dawn-dusk asymmetry in particles of solar wind origin within the magnetosphere

Solar wind/magnetosheath plasma in the magnetosphere can be identified using a component that has a higher charge state, lower density and, at least soon after their entry into the magnetosphere, lower energy than plasma from a terrestrial source. We survey here observations taken over 3 years of He2+ ions made by the Magnetospheric Ion Composition Sensor (MICS) of the Charge and Mass Mgnetospheric Ion Composition Experiment (CAMMICE) instrument aboard POLAR. The occurrence probability of these solar wind ions is then plotted as a function of Magnetic Local Time (MLT) and invariant latitude (3) for various energy ranges. For all energies observed by MICS (1.8–21.4 keV) and all solar wind conditions, the occurrence probabilities peaked around the cusp region and along the dawn flank. The solar wind conditions were filtered to see if this dawnward asymmetry is controlled by the Svalgaard-Mansurov effect (and so depends on the BY component of the interplanetary magnetic field, IMF) or by Fermi acceleration of He2+ at the bow shock (and so depends on the IMF ratio BX/BY ). It is shown that the asymmetry remained persistently on the dawn flank, suggesting it was not due to effects associated with direct entry into the magnetosphere. This asymmetry, with enhanced fluxes on the dawn flank, persisted for lower energy ions (below a “cross-over” energy of about 23 keV) but reversed sense to give higher fluxes on the dusk flank at higher energies. This can be explained by the competing effects of gradient/curvature drifts and the convection electric field on ions that are convecting sunward on re-closed field lines. The lower-energy He2+ ions E × B drift dawnwards as they move earthward, whereas the higher energy ions curvature/gradient drift towards dusk. The convection electric field in the tail is weaker for northward IMF. Ions then need less energy to drift to the dusk flank, so that the cross-over energy, at which the asymmetry changes sense, is reduced.

ESR and EISCAT observations of the response of the cusp and cleft to IMF orientation changes

We report observations of the cusp/cleft ionosphere made on December 16th 1998 by the EISCAT (European incoherent scatter) VHF radar at Tromso and the EISCAT Svalbard radar (ESR). We compare them with observations of the dayside auroral luminosity, as seen by meridian scanning photometers at Ny Alesund and of HF radar backscatter, as observed by the CUTLASS radar. We study the response to an interval of about one hour when the interplanetary magnetic field (IMF), monitored by the WIND and ACE spacecraft, was southward. The cusp/cleft aurora is shown to correspond to a spatially extended region of elevated electron temperatures in the VHF radar data. Initial conditions were characterised by a northward-directed IMF and cusp/cleft aurora poleward of the ESR. A strong southward turning then occurred, causing an equatorward motion of the cusp/cleft aurora. Within the equatorward expanding, southward-IMF cusp/cleft, the ESR observed structured and elevated plasma densities and ion and electron temperatures. Cleft ion fountain upflows were seen in association with elevated ion temperatures and rapid eastward convection, consistent with the magnetic curvature force on newly opened held lines for the observed negative IMF B-y. Subsequently, the ESR beam remained immediately poleward of the main cusp/cleft and a sequence of poleward-moving auroral transients passed over it. After the last of these, the ESR was in the polar cap and the radar observations were characterised by extremely low ionospheric densities and downward field-aligned flows. The IMF then turned northward again and the auroral oval contracted such that the ESR moved back into the cusp/cleft region. For the poleward-retreating northward-IMF cusp/cleft, the convection flows were slower, upflows were weaker and the electron density and temperature enhancements were less structured. Following the northward turning, the bands of high electron temperature and cusp/cleft aurora bifurcated, consistent with both subsolar and lobe reconnection taking place simultaneously. The present paper describes the large-scale behaviour of the ionosphere during this interval, as observed by a powerful combination of instruments. Two companion papers, by Lockwood et al. (2000) and Thorolfsson et al. (2000), both in this issue, describe the detailed behaviour of the poleward-moving transients observed during the interval of southward B-z, and explain their morphology in the context of previous theoretical work.

Plasma structure within poleward-moving cusp/cleft auroral transients: EISCAT Svalbard radar observations and an explanation in terms of large local time extent of events

We report high-resolution observations of the southward-IMF cusp/cleft ionosphere made on December 16th 1998 by the EISCAT (European incoherent scatter) Svalbard radar (ESR), and compare them with observations of dayside auroral luminosity, as seen at a wavelength of 630 nm by a meridian scanning photometer at Ny Alesund, and of plasma flows, as seen by the CUTLASS (co-operative UK twin location auroral sounding system) Finland HF radar. The optical data reveal a series of poleward-moving transient red-line (630 nm) enhancements, events that have been associated with bursts in the rate of magnetopause reconnection generating new open flux. The combined observations at this time have strong similarities to predictions of the effects of soft electron precipitation modulated by pulsed reconnection, as made by Davis and Lockwood (1996); however, the effects of rapid zonal flow in the ionosphere, caused by the magnetic curvature force on the newly opened field lines, are found to be a significant additional factor. In particular, it is shown how enhanced plasma loss rates induced by the rapid convection can explain two outstanding anomalies of the 630 nm transients, namely how minima in luminosity form between the poleward-moving events and how events can re-brighten as they move poleward. The observations show how cusp/cleft aurora and transient poleward-moving auroral forms appear in the ESR data and the conditions which cause enhanced 630 nm emission in the transients: they are an important first step in enabling the ESR to identify these features away from the winter solstice when supporting auroral observations are not available.

A multipoint study of a substorm occurring on 7 December, 1992, and its theoretical implications

On 7 December 1992, a moderate substorm was observed by a variety of satellites and ground-based instruments. Ionospheric flows were monitored near dusk by the Goose Bay HF radar and near midnight by the EISCAT radar. The observed flows are compared here with magnetometer observations by the IMAGE array in Scandinavia and the two Greenland chains, the auroral distribution observed by Freja and the substorm cycle observations by the SABRE radar, the SAMNET magnetometer array and LANL geosynchronous satellites. Data from Galileo Earth-encounter II are used to estimate the IMF B-z component. The data presented show that the substorm onset electrojet at midnight was confined to closed field lines equatorward of the preexisting convection reversal boundaries observed in the dusk and midnight regions. No evidence of substantial closure of open flux was detected following this substorm onset. Indeed the convection reversal boundary on the duskside continued to expand equatorward after onset due to the continued presence of strong southward IMF, such that growth and expansion phase features were simultaneously present. Clear indications of closure of open flux were not observed until a subsequent substorm intensification 25 min after the initial onset. After this time, the substorm auroral bulge in the nightside hours propagated well poleward of the pre-existing convection reversal boundary, and strong flow perturbations were observed by the Goose Bay radar, indicative of flows driven by reconnection in the tail.