Optical characterization of the growth and spatial structure of a substorm onset arc

We present a detailed case study of the characteristics of auroral forms that constitute the first ionospheric signatures of substorm expansion phase onset. Analysis of the optical frequency and along-arc (azimuthal) wave number spectra provides the strongest constraint to date on the potential mechanisms and instabilities in the near-Earth magnetosphere that accompany auroral onset and which precede poleward arc expansion and auroral breakup. We evaluate the frequency and growth rates of the auroral forms as a function of azimuthal wave number to determine whether these wave characteristics are consistent with current models of the substorm onset mechanism. We find that the frequency, spatial scales, and growth rates of the auroral forms are most consistent with the cross-field current instability or a ballooning instability, most likely triggered close to the inner edge of the ion plasma sheet. This result is supportive of a near-Earth plasma sheet initiation of the substorm expansion phase. We also present evidence that the frequency and phase characteristics of the auroral undulations may be generated via resonant processes operating along the geomagnetic field. Our observations provide the most powerful constraint to date on the ionospheric manifestation of the physical processes operating during the first few minutes around auroral substorm onset.

Monte Carlo field-theoretic simulations for melts of symmetric diblock copolymer

Monte Carlo field-theoretic simulations (MCFTS) are performed on melts of symmetric diblock copolymer for invariant polymerization indexes extending down to experimentally relevant values of N̅ ∼ 10^4. The simulations are performed with a fluctuating composition field, W_−(r), and a pressure field, W_+(r), that follows the saddle-point approximation. Our study focuses on the disordered-state structure function, S(k), and the order−disorder transition (ODT). Although shortwavelength fluctuations cause an ultraviolet (UV) divergence in three dimensions, this is readily compensated for with the use of an effective Flory−Huggins interaction parameter, χ_e. The resulting S(k) matches the predictions of renormalized one-loop (ROL) calculations over the full range of χ_eN and N̅ examined in our study, and agrees well with Fredrickson−Helfand (F−H) theory near the ODT. Consistent with the F−H theory, the ODT is discontinuous for finite N̅ and the shift in (χ_eN)_ODT follows the predicted N̅^−1/3 scaling over our range of N̅.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr – Part 2: A new reconstruction of the interplanetary magnetic field

We present a new reconstruction of the interplanetary magnetic field (IMF, B) for 1846–2012 with a full analysis of errors, based on the homogeneously constructed IDV(1d)composite of geomagnetic activity presented in Part 1 (Lockwood et al., 2013a). Analysis of the dependence of the commonly used geomagnetic indices on solar wind parameters is presented which helps explain why annual means of interdiurnal range data, such as the new composite, depend only on the IMF with only a very weak influence of the solar wind flow speed. The best results are obtained using a polynomial (rather than a linear) fit of the form B = χ · (IDV(1d) − β)α with best-fit coefficients χ = 3.469, β = 1.393 nT, and α = 0.420. The results are contrasted with the reconstruction of the IMF since 1835 by Svalgaard and Cliver (2010).

Landscape structure mediates the effects of a stressor on field vole populations

Spatio-temporal landscape heterogeneity has rarely been considered in population-level impact assessments. Here we test whether landscape heterogeneity is important by examining the case of a pesticide applied seasonally to orchards which may affect non-target vole populations, using a validated ecologically realistic and spatially explicit agent-based model. Voles thrive in unmanaged grasslands and untreated orchards but are particularly exposed to applied pesticide treatments during dispersal between optimal habitats. We therefore hypothesised that vole populations do better (1) in landscapes containing more grassland and (2) where areas of grassland are closer to orchards, but (3) do worse if larger areas of orchards are treated with pesticide. To test these hyposeses we made appropriate manipulations to a model landscape occupied by field voles. Pesticide application reduced model population sizes in all three experiments, but populations subsequently wholly or partly recovered. Population depressions were, as predicted, lower in landscapes containing more unmanaged grassland, in landscapes with reduced distance between grassland and orchards, and in landscapes with fewer treated orchards. Population recovery followed a similar pattern except for an unexpected improvement in recovery when the area of treated orchards was increased. Outside the period of pesticide application, orchards increase landscape connectivity and facilitate vole dispersal and so speed population recovery. Overall our results show that accurate prediction of population impact cannot be achieved without taking account of landscape structure. The specifics of landscape structure and habitat connectivity are likely always important in mediating the effects of stressors.

Open solar flux estimates from near-Earth measurements of the interplanetary magnetic field: comparison of the first two perihelion passes of the Ulysses spacecraft

Results from all phases of the orbits of the Ulysses spacecraft have shown that the magnitude of the radial component of the heliospheric field is approximately independent of heliographic latitude. This result allows the use of near- Earth observations to compute the total open flux of the Sun. For example, using satellite observations of the interplanetary magnetic field, the average open solar flux was shown to have risen by 29% between 1963 and 1987 and using the aa geomagnetic index it was found to have doubled during the 20th century. It is therefore important to assess fully the accuracy of the result and to check that it applies to all phases of the solar cycle. The first perihelion pass of the Ulysses spacecraft was close to sunspot minimum, and recent data from the second perihelion pass show that the result also holds at solar maximum. The high level of correlation between the open flux derived from the various methods strongly supports the Ulysses discovery that the radial field component is independent of latitude. We show here that the errors introduced into open solar flux estimates by assuming that the heliospheric field’s radial component is independent of latitude are similar for the two passes and are of order 25% for daily values, falling to 5% for averaging timescales of 27 days or greater. We compare here the results of four methods for estimating the open solar flux with results from the first and second perehelion passes by Ulysses. We find that the errors are lowest (1–5% for averages over the entire perehelion passes lasting near 320 days), for near-Earth methods, based on either interplanetary magnetic field observations or the aa geomagnetic activity index. The corresponding errors for the Solanki et al. (2000) model are of the order of 9–15% and for the PFSS method, based on solar magnetograms, are of the order of 13–47%. The model of Solanki et al. is based on the continuity equation of open flux, and uses the sunspot number to quantify the rate of open flux emergence. It predicts that the average open solar flux has been decreasing since 1987, as Correspondence to: M. Lockwood (m.lockwood@rl.ac.uk) is observed in the variation of all the estimates of the open flux. This decline combines with the solar cycle variation to produce an open flux during the second (sunspot maximum) perihelion pass of Ulysses which is only slightly larger than that during the first (sunspot minimum) perihelion pass.

Relationship between the near-Earth interplanetary field and the coronal source flux: dependence on timescale

The Ulysses spacecraft has shown that the radial component of the heliospheric magnetic field is approximately independent of latitude. This has allowed quantification of the total open solar flux from near-Earth observations of the interplanetary magnetic field. The open flux can also be estimated from photospheric magnetograms by mapping the fields up to the ‘‘coronal source surface’’ where the field is assumed to be radial and which is usually assumed to be at a heliocentric distance r = 2.5R_{S} (a mean solar radius, 1R_{S} = 6.96x10^{8} m). These two classes of open flux estimate will differ by the open flux that threads the heliospheric current sheet(s) inside Earth’s orbit at 2.5R_{S} < r < 1R{1} (where the mean Earth-Sun distance, 1R_{1} = 1 AU = 1.5 x 10^{11} m). We here use near-Earth measurements to estimate this flux and show that at sunspot minimum it causes only a very small (approximately 0.5%) systematic difference between the two types of open flux estimate, with an uncertainty that is of order ±24% in hourly values, ±16% in monthly averages, and between -6% and +2% in annual values. These fractions may be somewhat larger for sunspot maximum because of flux emerging at higher heliographic latitudes.

Dynamic auroral structure in the vicinity of the polar cusp - multipoint observations during southward and northward IMF

The low- and high-latitude boundary layers of the earth's magnetosphere [low-latitude boundary layer (LLBL) and mantle] play important roles in transferring momentum and energy from the solar wind to the magnetosphere-ionosphere system. Particle precipitation, field-aligned current, auroral emission, ionospheric ion drift and ground magnetic perturbations are among the low-altitude parameters that show signatures of various plasma processes in the LLBL and the magnetopause current layer. Magnetic merging events, Kelvin-Helmholtz waves, and pressure pulses excited by the variable solar wind/magnetosheath plasma are examples of boundary phenomena that may be coupled to the ionosphere via field-aligned currents. Optical auroral observation, by photometry and all-sky TV cameras, is a unique technique for investigating the spatial and temporal structure of the electron precipitation associated with such phenomena. However, the distinction between the different boundary layer plasma populations cannot in general be unambiguously determined by optics alone. Additional information, such as satellite observations of particle boundaries and field-aligned currents, is needed in order to identify the plasma source(s) and the magnetosphere-ionosphere coupling mode(s). Two categories of auroral activity/structure in the vicinity of the polar cusp are discussed in this paper, based on combined ground and satellite data. In one case, the quasi-periodic sequence of auroral events at the polar cap boundary involves accelerated electrons (< 1 keV) moving poleward (< 1 km s-1) and azimuthally along the persistent cusp/cleft arc poleward boundary with velocities (< 4 km s-1), comparable to the local ionospheric ion drift during periods of southward IMF. A critical question is whether or not the optical events signify a corresponding plasma flow across the open/closed field line boundary in such cases. Near-simultaneous observations of magnetopause flux transfer events (FTEs) and such optical/ion drift events are reported. The reverse pattern of motion of discrete auroral forms is observed during positive interplanetary magnetic field (IMF) B(Z), i.e. equatorward motion into the cusp/cleft background arc from the poleward edge. Combined satellite and ground-based information for the latter cases indicate a source mechanism, poleward of the cusp at the high-latitude magnetopause or plasma mantle, giving rise to strong momentum transfer and electron precipitation structures within a approximately 200 km-wide latitudinal zone at the cusp/cleft poleward boundary. The striking similarities of auroral electrodynamics in the cleft/mantle region during northward and southward IMF indicate that a qualitatively similar solar wind-magnetosphere coupling mode is operating. It is suggested that, in both cases, the discrete auroral forms represent temporal/spatial structure of larger-scale convection over the polar magnetosphere.

A diabatically-generated potential vorticity structure near the extratropical tropopause in three simulated extratropical cyclones

The structure of near-tropopause potential vorticity (PV) acts as a primary control on the evolution of extratropical cyclones. Diabatic processes such as the latent heating found in ascending moist warm conveyor belts modify PV. A dipole in diabatically-generated PV (hereafter diabatic PV) straddling the extratropical tropopause, with the positive pole above the negative pole, was diagnosed in a recently published analysis of a simulated extratropical cyclone. This PV dipole has the potential to significantly modify the propagation of Rossby waves and the growth of baroclinically-unstable waves. This previous analysis was based on a single case study simulated with 12-km horizontal grid spacing and parameterized convection. Here, the dipole is investigated in three additional cold-season extratropical cyclones simulated in both convection-parameterizing and convection-permitting model configurations. A diabatic PV dipole across the extratropical tropopause is diagnosed in all three cases. The amplitude of the dipole saturates approximately 36 hours from the time diabatic PV is accumulated. The node elevation of the dipole varies between 2-4 PVU in the three cases, and the amplitude of the system-averaged dipole varies between 0.2-0.4 PVU. The amplitude of the negative pole is similar in the convection-parameterizing and convection-permitting simulations. The positive pole, which is generated by long-wave radiative cooling, is weak in the convection-permitting simulations due to the small domain size which limits the accumulation of diabatic tendencies within the interior of the domain. The possible correspondence between the diabatic PV dipole and the extratropical tropopause inversion layer is discussed.

Modelling spectral reflectance of open cork oak woodland: a simulation analysis of the effects of vegetation structure and background

This study analyses the influence of vegetation structure (i.e. leaf area index and canopy cover) and seasonal background changes on moderate-resolution imaging spectrometer (MODIS)-simulated reflectance data in open woodland. Approximately monthly spectral reflectance and transmittance field measurements (May 2011 to October 2013) of cork oak tree leaves (Quercus suber) and of the herbaceous understorey were recorded in the region of Ribatejo, Portugal. The geometric-optical and radiative transfer (GORT) model was used to simulate MODIS response (red, near-infrared) and to calculate vegetation indices, investigating their response to changes in the structure of the overstorey vegetation and to seasonal changes in the understorey using scenarios corresponding to contrasting phenological status (dry season vs. wet season). The performance of normalized difference vegetation index (NDVI), soil-adjusted vegetation index (SAVI), and enhanced vegetation index (EVI) is discussed. Results showed that SAVI and EVI were very sensitive to the emergence of background vegetation in the wet season compared to NDVI and that shading effects lead to an opposing trend in the vegetation indices. The information provided by this research can be useful to improve our understanding of the temporal dynamic of vegetation, monitored by vegetation indices.

Using Distributed Temperature Sensing to monitor field scale dynamics of ground surface temperature and related substrate heat flux

We present one of the first studies of the use of Distributed Temperature Sensing (DTS) along fibre-optic cables to purposely monitor spatial and temporal variations in ground surface temperature (GST) and soil temperature, and provide an estimate of the heat flux at the base of the canopy layer and in the soil. Our field site was at a groundwater-fed wet meadow in the Netherlands covered by a canopy layer (between 0-0.5 m thickness) consisting of grass and sedges. At this site, we ran a single cable across the surface in parallel 40 m sections spaced by 2 m, to create a 40×40 m monitoring field for GST. We also buried a short length (≈10 m) of cable to depth of 0.1±0.02 m to measure soil temperature. We monitored the temperature along the entire cable continuously over a two-day period and captured the diurnal course of GST, and how it was affected by rainfall and canopy structure. The diurnal GST range, as observed by the DTS system, varied between 20.94 and 35.08◦C; precipitation events acted to suppress the range of GST. The spatial distribution of GST correlated with canopy vegetation height during both day and night. Using estimates of thermal inertia, combined with a harmonic analysis of GST and soil temperature, substrate and soil-heat fluxes were determined. Our observations demonstrate how the use of DTS shows great promise in better characterising area-average substrate/soil heat flux, their spatiotemporal variability, and how this variability is affected by canopy structure. The DTS system is able to provide a much richer data set than could be obtained from point temperature sensors. Furthermore, substrate heat fluxes derived from GST measurements may be able to provide improved closure of the land surface energy balance in micrometeorological field studies. This will enhance our understanding of how hydrometeorological processes interact with near-surface heat fluxes.

Spatial and temporal variability of the concentration field from localized releases in a regular building array

Spatial and temporal fluctuations in the concentration field from an ensemble of continuous point-source releases in a regular building array are analyzed from data generated by direct numerical simulations. The release is of a passive scalar under conditions of neutral stability. Results are related to the underlying flow structure by contrasting data for an imposed wind direction of 0 deg and 45 deg relative to the buildings. Furthermore, the effects of distance from the source and vicinity to the plume centreline on the spatial and temporal variability are documented. The general picture that emerges is that this particular geometry splits the flow domain into segments (e.g. “streets” and “intersections”) in each of which the air is, to a first approximation, well mixed. Notable exceptions to this general rule include regions close to the source, near the plume edge, and in unobstructed channels when the flow is aligned. In the oblique (45 deg) case the strongly three-dimensional nature of the flow enhances mixing of a scalar within the canopy leading to reduced temporal and spatial concentration fluctuations within the plume core. These fluctuations are in general larger for the parallel flow (0 deg) case, especially so in the long unobstructed channels. Due to the more complex flow structure in the canyon-type streets behind buildings, fluctuations are lower than in the open channels, though still substantially larger than for oblique flow. These results are relevant to the formulation of simple models for dispersion in urban areas and to the quantification of the uncertainties in their predictions.