An infrared desert dust index for the Along-Track Scanning Radiometers

A new aerosol index for the Along-Track Scanning Radiometers (ATSRs) is presented that provides a means to detect desert dust contamination in infrared SST retrievals. The ATSR Saharan dust index (ASDI) utilises only the thermal infrared channels and may therefore be applied consistently to the entire ATSR data record (1991 to present), for both day time and night time observations. The derivation of the ASDI is based on a principal component (PC) analysis (PCA) of two unique pairs of channel brightness temperature differences (BTDs). In 2-D space (i.e. BTD vs BTD), it is found that the loci of data unaffected by aerosol are confined to a single axis of variability. In contrast, the loci of aerosol-contaminated data fall off-axis, shifting in a direction that is approximately orthogonal to the clear-sky axis. The ASDI is therefore defined to be the second PC, where the first PC accounts for the clear-sky variability. The primary ASDI utilises the ATSR nadir and forward-view observations at 11 and 12 μm (ASDI2). A secondary, three-channel nadir-only ASDI (ASDI3) is also defined for situations where data from the forward view are not available. Empirical and theoretical analyses suggest that ASDI is well correlated with aerosol optical depth (AOD: correlation r is typically > 0.7) and provides an effective tool for detecting desert mineral dust. Overall, ASDI2 is found to be more effective than ASDI3, with the latter being sensitive only to very high dust loading. In addition, use of ASDI3 is confined to night time observations as it relies on data from the 3.7 μm channel, which is sensitive to reflected solar radiation. This highlights the benefits of having data from both a nadir- and a forward-view for this particular approach to aerosol detection.

Modeling ocean primary production: sensitivity to spectral resolution of attenuation and absorption of light

Modeling the vertical penetration of photosynthetically active radiation (PAR) through the ocean, and its utilization by phytoplankton, is fundamental to simulating marine primary production. The variation of attenuation and absorption of light with wavelength suggests that photosynthesis should be modeled at high spectral resolution, but this is computationally expensive. To model primary production in global 3d models, a balance between computer time and accuracy is necessary. We investigate the effects of varying the spectral resolution of the underwater light field and the photosynthetic efficiency of phytoplankton (α∗), on primary production using a 1d coupled ecosystem ocean turbulence model. The model is applied at three sites in the Atlantic Ocean (CIS (∼60°N), PAP (∼50°N) and ESTOC (∼30°N)) to include the effect of different meteorological forcing and parameter sets. We also investigate three different methods for modeling α∗ – as a fixed constant, varying with both wavelength and chlorophyll concentration [Bricaud, A., Morel, A., Babin, M., Allali, K., Claustre, H., 1998. Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic (case 1) waters. Analysis and implications for bio-optical models. J. Geophys. Res. 103, 31033–31044], and using a non-spectral parameterization [Anderson, T.R., 1993. A spectrally averaged model of light penetration and photosynthesis. Limnol. Oceanogr. 38, 1403–1419]. After selecting the appropriate ecosystem parameters for each of the three sites we vary the spectral resolution of light and α∗ from 1 to 61 wavebands and study the results in conjunction with the three different α∗ estimation methods. The results show modeled estimates of ocean primary productivity are highly sensitive to the degree of spectral resolution and α∗. For accurate simulations of primary production and chlorophyll distribution we recommend a spectral resolution of at least six wavebands if α∗ is a function of wavelength and chlorophyll, and three wavebands if α∗ is a fixed value.

Deriving a sea surface temperature record suitable for climate change research from the along-track scanning radiometers

We describe the approach to be adopted for a major new initiative to derive a homogeneous record of sea surface temperature for 1991–2007 from the observations of the series of three along-track scanning radiometers (ATSRs). This initiative is called (A)RC: (Advanced) ATSR Re-analysis for Climate. The main objectives are to reduce regional biases in retrieved sea surface temperature (SST) to less than 0.1 K for all global oceans, while creating a very homogenous record that is stable in time to within 0.05 K decade−1, with maximum independence of the record from existing analyses of SST used in climate change research. If these stringent targets are achieved, this record will enable significantly improved estimates of surface temperature trends and variability of sufficient quality to advance questions of climate change attribution, climate sensitivity and historical reconstruction of surface temperature changes. The approach includes development of new, consistent estimators for SST for each of the ATSRs, and detailed analysis of overlap periods. Novel aspects of the approach include generation of multiple versions of the record using alternative channel sets and cloud detection techniques, to assess for the first time the effect of such choices. There will be extensive effort in quality control, validation and analysis of the impact on climate SST data sets. Evidence for the plausibility of the 0.1 K target for systematic error is reviewed, as is the need for alternative cloud screening methods in this context.

A cross-calibration of GMS-5 thermal channels against ATSR-2

Motivated by the importance to weather and climate of the Indo-Pacific seas, we present a new calibration of the Visible Infrared Spin-Scan Radiometer (VISSR) on the geostationary meteorological satellite, GMS-5. VISSR imagery has significant potential for exploring the dynamics of the ocean and air–sea interactions in this poorly characterized region, by virtue of the VISSR's surface temperature retrieval capability and hourly sampling. However, the calibration of the thermal imagery supplied by the Japanese Meteorological Agency (JMA) is inconsistent with the spectral characteristics of the channels, and published details of the JMA calibration procedure are scant. We use the well-characterized Along-Track Scanning Radiometer 2 (ATSR-2) as a reference, and determine calibration corrections for GMS-5 VISSR. We obtain more credible VISSR brightness temperatures and demonstrate sea surface temperature (SST) retrieval that validates well against in situ measurements (bias ∼0.3 and scatter ∼0.4 K). Comparison with a widely used sea surface temperature analysis shows that the GMS-5 VISSR SST fields capture important spatial structure, absent in the analysis.

Equator-to-pole temperature differences and the extra-tropical storm track responses of the CMIP5 climate models

This paper aims to understand the physical processes causing the large spread in the storm track projections of the CMIP5 climate models. In particular, the relationship between the climate change responses of the storm tracks, as measured by the 2–6 day mean sea level pressure variance, and the equator-to-pole temperature differences at upper- and lower-tropospheric levels is investigated. In the southern hemisphere the responses of the upper- and lower-tropospheric temperature differences are correlated across the models and as a result they share similar associations with the storm track responses. There are large regions in which the storm track responses are correlated with the temperature difference responses, and a simple linear regression model based on the temperature differences at either level captures the spatial pattern of the mean storm track response as well explaining between 30 and 60 % of the inter-model variance of the storm track responses. In the northern hemisphere the responses of the two temperature differences are not significantly correlated and their associations with the storm track responses are more complicated. In summer, the responses of the lower-tropospheric temperature differences dominate the inter-model spread of the storm track responses. In winter, the responses of the upper- and lower-temperature differences both play a role. The results suggest that there is potential to reduce the spread in storm track responses by constraining the relative magnitudes of the warming in the tropical and polar regions.

CP-OQAM-OFDM based SC-FDMA: adjustable user bandwidth and space-time coding

The discrete Fourier transmission spread OFDM DFTS-OFDM) based single-carrier frequency division multiple access (SC-FDMA) has been widely adopted due to its lower peak-to-average power ratio (PAPR) of transmit signals compared with OFDM. However, the offset modulation, which has lower PAPR than general modulation, cannot be directly applied into the existing SC-FDMA. When pulse-shaping filters are employed to further reduce the envelope fluctuation of transmit signals of SC-FDMA, the spectral efficiency degrades as well. In order to overcome such limitations of conventional SC-FDMA, this paper for the first time investigated cyclic prefixed OQAMOFDM (CP-OQAM-OFDM) based SC-FDMA transmission with adjustable user bandwidth and space-time coding. Firstly, we propose CP-OQAM-OFDM transmission with unequally-spaced subbands. We then apply it to SC-FDMA transmission and propose a SC-FDMA scheme with the following features: a) the transmit signal of each user is offset modulated single-carrier with frequency-domain pulse-shaping; b) the bandwidth of each user is adjustable; c) the spectral efficiency does not decrease with increasing roll-off factors. To combat both inter-symbolinterference and multiple access interference in frequencyselective fading channels, a joint linear minimum mean square error frequency domain equalization using a prior information with low complexity is developed. Subsequently, we construct space-time codes for the proposed SC-FDMA. Simulation results confirm the powerfulness of the proposed CP-OQAM-OFDM scheme (i.e., effective yet with low complexity).

Variability of the North Atlantic summer storm track: mechanisms and impacts on European climate

The summertime variability of the extratropical storm track over the Atlantic sector and its links to European climate have been analysed for the period 1948–2011 using observations and reanalyses. The main results are as follows. (1) The dominant mode of the summer storm track density variability is characterized by a meridional shift of the storm track between two distinct paths and is related to a bimodal distribution in the climatology for this region. It is also closely related to the Summer North Atlantic Oscillation (SNAO). (2) A southward shift is associated with a downstream extension of the storm track and a decrease in blocking frequency over the UK and northwestern Europe. (3) The southward shift is associated with enhanced precipitation over the UK and northwestern Europe and decreased precipitation over southern Europe (contrary to the behaviour in winter). (4) There are strong ocean–atmosphere interactions related to the dominant mode of storm track variability. The atmosphere forces the ocean through anomalous surface fluxes and Ekman currents, but there is also some evidence consistent with an ocean influence on the atmosphere, and that coupled ocean–atmosphere feedbacks might play a role. The ocean influence on the atmosphere may be particularly important on decadal timescales, related to the Atlantic Multidecadal Oscillation (AMO).

Further evidence for variable synchrotron emission in XTE J1118+480 in outburst

We present simultaneous multicolor infrared and optical photometry of the black hole X-ray transient XTE J1118+480 during its short 2005 January outburst, supported by simultaneous X-ray observations. The variability is dominated by short timescales, ~10 s, although a weak superhump also appears to be present in the optical. The optical rapid variations, at least, are well correlated with those in X-rays. Infrared JHKs photometry, as in the previous outburst, exhibits especially large-amplitude variability. The spectral energy distribution (SED) of the variable infrared component can be fitted with a power law of slope α=-0.78+/-0.07, where F_ν~ν^α. There is no compelling evidence for evolution in the slope over five nights, during which time the source brightness decayed along almost the same track as seen in variations within the nights. We conclude that both short-term variability and longer timescale fading are dominated by a single component of constant spectral shape. We cannot fit the SED of the IR variability with a credible thermal component, either optically thick or thin. This IR SED is, however, approximately consistent with optically thin synchrotron emission from a jet. These observations therefore provide indirect evidence to support jet-dominated models for XTE J1118+480 and also provide a direct measurement of the slope of the optically thin emission, which is impossible, based on the average spectral energy distribution alone.

Effects of non-uniform crosswind fields on scintillometry measurements

The effects of a non-uniform wind field along the path of a scintillometer are investigated. Theoretical spectra are calculated for a range of scenarios where the crosswind varies in space or time and compared to the ‘ideal’ spectrum based on a constant uniform crosswind. It is verified that the refractive-index structure parameter relation with the scintillometer signal remains valid and invariant for both spatially and temporally-varying crosswinds. However, the spectral shape may change significantly preventing accurate estimation of the crosswind speed from the peak of the frequency spectrum and retrieval of the structure parameter from the plateau of the power spectrum. On comparison with experimental data, non-uniform crosswind conditions could be responsible for previously unexplained features sometimes seen in observed spectra. By accounting for the distribution of crosswind, theoretical spectra can be generated that closely replicate the observations, leading to a better understanding of the measurements. Spatial variability of wind speeds should be expected for paths other than those that are parallel to the surface and over flat, homogenous areas, whilst fluctuations in time are important for all sites.

Spectral asymmetry of the massless Dirac operator on a 3-torus

Consider the massless Dirac operator on a 3-torus equipped with Euclidean metric and standard spin structure. It is known that the eigenvalues can be calculated explicitly: the spectrum is symmetric about zero and zero itself is a double eigenvalue. The aim of the paper is to develop a perturbation theory for the eigenvalue with smallest modulus with respect to perturbations of the metric. Here the application of perturbation techniques is hindered by the fact that eigenvalues of the massless Dirac operator have even multiplicity, which is a consequence of this operator commuting with the antilinear operator of charge conjugation (a peculiar feature of dimension 3). We derive an asymptotic formula for the eigenvalue with smallest modulus for arbitrary perturbations of the metric and present two particular families of Riemannian metrics for which the eigenvalue with smallest modulus can be evaluated explicitly. We also establish a relation between our asymptotic formula and the eta invariant.

Linking Northern Hemisphere blocking and storm track biases in the CMIP5 climate models

The relationship between biases in Northern Hemisphere (NH) atmospheric blocking frequency and extratropical cyclone track density is investigated in 12 CMIP5 climate models to identify mechanisms underlying climate model biases and inform future model development. Biases in the Greenland blocking and summer Pacific blocking frequencies are associated with biases in the storm track latitudes while biases in winter European blocking frequency are related to the North Atlantic storm track tilt and Mediterranean cyclone density. However, biases in summer European and winter Pacific blocking appear less related with cyclone track density. Furthermore, the models with smaller biases in winter European blocking frequency have smaller biases in the cyclone density in Europe, which suggests that they are different aspects of the same bias. This is not found elsewhere in the NH. The summer North Atlantic and the North Pacific mean CMIP5 track density and blocking biases might therefore have different origins.

Cloud clearing techniques over land for land surface temperature retrieval from the Advanced Along Track Scanning Radiometer

We present five new cloud detection algorithms over land based on dynamic threshold or Bayesian techniques, applicable to the Advanced Along Track Scanning Radiometer (AATSR) instrument and compare these with the standard threshold based SADIST cloud detection scheme. We use a manually classified dataset as a reference to assess algorithm performance and quantify the impact of each cloud detection scheme on land surface temperature (LST) retrieval. The use of probabilistic Bayesian cloud detection methods improves algorithm true skill scores by 8-9 % over SADIST (maximum score of 77.93 % compared to 69.27 %). We present an assessment of the impact of imperfect cloud masking, in relation to the reference cloud mask, on the retrieved AATSR LST imposing a 2 K tolerance over a 3x3 pixel domain. We find an increase of 5-7 % in the observations falling within this tolerance when using Bayesian methods (maximum of 92.02 % compared to 85.69 %). We also demonstrate that the use of dynamic thresholds in the tests employed by SADIST can significantly improve performance, applicable to cloud-test data to provided by the Sea and Land Surface Temperature Radiometer (SLSTR) due to be launched on the Sentinel 3 mission (estimated 2014).

A nonlinear oscillator describing storm track variability

We construct a two-variable model which describes the interaction between local baroclinicity and eddy heat flux in order to understand aspects of the variance in storm tracks. It is a heuristic model for diabatically forced baroclinic instability close to baroclinic neutrality. The two-variable model has the structure of a nonlinear oscillator. It exhibits some realistic properties of observed storm track variability, most notably the intermittent nature of eddy activity. This suggests that apparent threshold behaviour can be more accurately and succinctly described by a simple nonlinearity. An analogy is drawn with triggering of convective events.

The real-time optimisation of DNO owned storage devices on the LV network for peak reduction

Energy storage is a potential alternative to conventional network reinforcementof the low voltage (LV) distribution network to ensure the grid’s infrastructure remainswithin its operating constraints. This paper presents a study on the control of such storagedevices, owned by distribution network operators. A deterministic model predictive control (MPC) controller and a stochastic receding horizon controller (SRHC) are presented, wherethe objective is to achieve the greatest peak reduction in demand, for a given storagedevice specification, taking into account the high level of uncertainty in the prediction of LV demand. The algorithms presented in this paper are compared to a standard set-pointcontroller and bench marked against a control algorithm with a perfect forecast. A specificcase study, using storage on the LV network, is presented, and the results of each algorithmare compared. A comprehensive analysis is then carried out simulating a large number of LV networks of varying numbers of households. The results show that the performance of each algorithm is dependent on the number of aggregated households. However, on a typical aggregation, the novel SRHC algorithm presented in this paper is shown to outperform each of the comparable storage control techniques.

Synergistic angular and spectral estimation of aerosol properties using CHRIS/PROBA-1 and simulated Sentinel-3~data

A method has been developed to estimate Aerosol Optical Depth (AOD), Fine Mode Fraction (FMF) and Single Scattering Albedo (SSA) over land surfaces using simulated Sentinel-3 data. The method uses inversion of a coupled surface/atmosphere radiative transfer model, and includes a general physical model of angular surface reflectance. An iterative process is used to determine the optimum value of the aerosol properties providing the best fit of the corrected reflectance values for a number of view angles and wavelengths with those provided by the physical model. A method of estimating AOD using only angular retrieval has previously been demonstrated on data from the ENVISAT and PROBA-1 satellite instruments, and is extended here to the synergistic spectral and angular sampling of Sentinel-3 and the additional aerosol properties. The method is tested using hyperspectral, multi-angle Compact High Resolution Imaging Spectrometer (CHRIS) images. The values obtained from these CHRIS observations are validated using ground based sun-photometer measurements. Results from 22 image sets using the synergistic retrieval and improved aerosol models show an RMSE of 0.06 in AOD, reduced to 0.03 over vegetated targets.