On reliability analysis of multi-categorical forecasts

Reliability analysis of probabilistic forecasts, in particular through the rank histogram or Talagrand diagram, is revisited. Two shortcomings are pointed out: Firstly, a uniform rank histogram is but a necessary condition for reliability. Secondly, if the forecast is assumed to be reliable, an indication is needed how far a histogram is expected to deviate from uniformity merely due to randomness. Concerning the first shortcoming, it is suggested that forecasts be grouped or stratified along suitable criteria, and that reliability is analyzed individually for each forecast stratum. A reliable forecast should have uniform histograms for all individual forecast strata, not only for all forecasts as a whole. As to the second shortcoming, instead of the observed frequencies, the probability of the observed frequency is plotted, providing and indication of the likelihood of the result under the hypothesis that the forecast is reliable. Furthermore, a Goodness-Of-Fit statistic is discussed which is essentially the reliability term of the Ignorance score. The discussed tools are applied to medium range forecasts for 2 m-temperature anomalies at several locations and lead times. The forecasts are stratified along the expected ranked probability score. Those forecasts which feature a high expected score turn out to be particularly unreliable.

Toward a new generation of world climate research and computing facilities

The impending threat of global climate change and its regional manifestations is among the most important and urgent problems facing humanity. Society needs accurate and reliable estimates of changes in the probability of regional weather variations to develop science-based adaptation and mitigation strategies. Recent advances in weather prediction and in our understanding and ability to model the climate system suggest that it is both necessary and possible to revolutionize climate prediction to meet these societal needs. However, the scientific workforce and the computational capability required to bring about such a revolution is not available in any single nation. Motivated by the success of internationally funded infrastructure in other areas of science, this paper argues that, because of the complexity of the climate system, and because the regional manifestations of climate change are mainly through changes in the statistics of regional weather variations, the scientific and computational requirements to predict its behavior reliably are so enormous that the nations of the world should create a small number of multinational high-performance computing facilities dedicated to the grand challenges of developing the capabilities to predict climate variability and change on both global and regional scales over the coming decades. Such facilities will play a key role in the development of next-generation climate models, build global capacity in climate research, nurture a highly trained workforce, and engage the global user community, policy-makers, and stakeholders. We recommend the creation of a small number of multinational facilities with computer capability at each facility of about 20 peta-flops in the near term, about 200 petaflops within five years, and 1 exaflop by the end of the next decade. Each facility should have sufficient scientific workforce to develop and maintain the software and data analysis infrastructure. Such facilities will enable questions of what resolution, both horizontal and vertical, in atmospheric and ocean models, is necessary for more confident predictions at the regional and local level. Current limitations in computing power have placed severe limitations on such an investigation, which is now badly needed. These facilities will also provide the world's scientists with the computational laboratories for fundamental research on weather–climate interactions using 1-km resolution models and on atmospheric, terrestrial, cryospheric, and oceanic processes at even finer scales. Each facility should have enabling infrastructure including hardware, software, and data analysis support, and scientific capacity to interact with the national centers and other visitors. This will accelerate our understanding of how the climate system works and how to model it. It will ultimately enable the climate community to provide society with climate predictions, which are based on our best knowledge of science and the most advanced technology.

The role of undisclosed reserves in English open outcry auctions

This study considers the role that reserve prices may play in residential property auctions. In comparison to much of the previous empirical work, this study has access to undisclosed reserve prices from English auctions. Consistent with theoretical arguments in the auction literature, the results obtained illustrate that whilst higher reserve prices increase the revenue obtained for the seller, they also reduce the probability of sale. The findings also highlight the importance of auction participation, with the number of individual bidders and the number of bids significant in most specifications.

Rossby wave-breaking analysis of explosive cyclones in the Euro-Atlantic sector

The two-way relationship between Rossby Wave-Breaking (RWB) and intensification of extra tropical cyclones is analysed over the Euro-Atlantic sector. In particular, the timing, intensity and location of cyclone development are related to RWB occurrences. For this purpose, two potential-temperature based indices are used to detect and classify anticyclonic and cyclonic RWB episodes from ERA-40 Re-Analysis data. Results show that explosive cyclogenesis over the North Atlantic (NA) is fostered by enhanced occurrence of RWB on days prior to the cyclone’s maximum intensification. Under such conditions, the eddy-driven jet stream is accelerated over the NA, thus enhancing conditions for cyclogenesis. For explosive cyclogenesis over the eastern NA, enhanced cyclonic RWB over eastern Greenland and anticyclonic RWB over the sub-tropical NA are observed. Typically only one of these is present in any given case, with the RWB over eastern Greenland being more frequent than its southern counterpart. This leads to an intensification of the jet over the eastern NA and enhanced probability of windstorms reaching Western Europe. Explosive cyclones evolving under simultaneous RWB on both sides of the jet feature a higher mean intensity and deepening rates than cyclones preceded by a single RWB event. Explosive developments over the western NA are typically linked to a single area of enhanced cyclonic RWB over western Greenland. Here, the eddy-driven jet is accelerated over the western NA. Enhanced occurrence of cyclonic RWB over southern Greenland and anticyclonic RWB over Europe is also observed after explosive cyclogenesis, potentially leading to the onset of Scandinavian Blocking. However, only very intense developments have a considerable influence on the large-scale atmospheric flow. Non-explosive cyclones depict no sign of enhanced RWB over the whole NA area. We conclude that the links between RWB and cyclogenesis over the Euro-Atlantic sector are sensitive to the cyclone’s maximum intensity, deepening rate and location.

Impact of I/Q imbalance on the performance of two-way CSI-assisted AF relaying

In this paper, we investigate half-duplex two-way dual-hop channel state information (CSI)-assisted amplify-and-forward (AF) relaying in the presence of in-phase and quadrature-phase (I/Q) imbalance. A compensation approach for the I/Q imbalance is proposed, which employs the received signals together with their conjugations to detect the desired signal. We also derive the average symbol error probability of the considered half-duplex two-way dual-hop CSI-assisted AF relaying networks with and without compensation for I/Q imbalance in Rayleigh fading channels. Numerical results are provided and show that the proposed compensation method mitigates the impact of I/Q imbalance to a certain extent.

Probabilistic physically-based cloud screening of satellite infra-red imagery for operational sea surface temperature retrieval

We propose and demonstrate a fully probabilistic (Bayesian) approach to the detection of cloudy pixels in thermal infrared (TIR) imagery observed from satellite over oceans. Using this approach, we show how to exploit the prior information and the fast forward modelling capability that are typically available in the operational context to obtain improved cloud detection. The probability of clear sky for each pixel is estimated by applying Bayes' theorem, and we describe how to apply Bayes' theorem to this problem in general terms. Joint probability density functions (PDFs) of the observations in the TIR channels are needed; the PDFs for clear conditions are calculable from forward modelling and those for cloudy conditions have been obtained empirically. Using analysis fields from numerical weather prediction as prior information, we apply the approach to imagery representative of imagers on polar-orbiting platforms. In comparison with the established cloud-screening scheme, the new technique decreases both the rate of failure to detect cloud contamination and the false-alarm rate by one quarter. The rate of occurrence of cloud-screening-related errors of >1 K in area-averaged SSTs is reduced by 83%. Copyright © 2005 Royal Meteorological Society.

The effects of corporate social performance on the cost of corporate debt and credit ratings

This study investigates the differential impact that various dimensions of corporate social performance have on the pricing of corporate debt as well as the assessment of the credit quality of specific bond issues. The empirical analysis, based on an extensive longitudinal data set, suggests that overall, good performance is rewarded and corporate social transgressions are penalized through lower and higher corporate bond yield spreads, respectively. Similar conclusions can be drawn when focusing on either the bond rating assigned to a specific debt issue or the probability of it being considered to be an asset of speculative grade.

Immediate auditory repetition of words and nonwords: An ERP study of lexical and sublexical processing

ERPs were elicited to (1) words, (2) pseudowords derived from these words, and (3) nonwords with no lexical neighbors, in a task involving listening to immediately repeated auditory stimuli. There was a significant early (P200) effect of phonotactic probability in the first auditory presentation, which discriminated words and pseudowords from nonwords; and a significant somewhat later (N400) effect of lexicality, which discriminated words from pseudowords and nonwords. There was no reliable effect of lexicality in the ERPs to the second auditory presentation. We conclude that early sublexical phonological processing differed according to phonotactic probability of the stimuli, and that lexically-based redintegration occurred for words but did not occur for pseudowords or nonwords. Thus, in online word recognition and immediate retrieval, phonological and/or sublexical processing plays a more important role than lexical level redintegration.

Centennial variations in sunspot number, open solar flux and streamer belt width: 1. correction of the sunspot number record since 1874

We analyse the widely-used international/ Zürich sunspot number record, R, with a view to quantifying a suspected calibration discontinuity around 1945 (which has been termed the “Waldmeier discontinuity” [Svalgaard, 2011]). We compare R against the composite sunspot group data from the Royal Greenwich Observatory (RGO) network and the Solar Optical Observing Network (SOON), using both the number of sunspot groups, N{sub}G{\sub}, and the total area of the sunspots, A{sub}G{\sub}. In addition, we compare R with the recently developed interdiurnal variability geomagnetic indices IDV and IDV(1d). In all four cases, linearity of the relationship with R is not assumed and care is taken to ensure that the relationship of each with R is the same before and after the putative calibration change. It is shown the probability that a correction is not needed is of order 10{sup}−8{\sup} and that R is indeed too low before 1945. The optimum correction to R for values before 1945 is found to be 11.6%, 11.7%, 10.3% and 7.9% using A{sub}G{\sub}, N{sub)G{\sub}, IDV, and IDV(1d), respectively. The optimum value obtained by combining the sunspot group data is 11.6% with an uncertainty range 8.1-14.8% at the 2σ level. The geomagnetic indices provide an independent yet less stringent test but do give values that fall within the 2σ uncertainty band with optimum values are slightly lower than from the sunspot group data. The probability of the correction needed being as large as 20%, as advocated by Svalgaard [2011], is shown to be 1.6 × 10{sup}−5{\sup}.

Long term changes in EUV and X-ray emissions from the solar corona and chromosphere as measured by the response of the Earth’s ionosphere during total solar eclipses from 1932 to 1999

Measurements of the ionospheric E region during total solar eclipses in the period 1932-1999 have been used to investigate the fraction of Extreme Ultra Violet and soft X-ray radiation, phi, that is emitted from the limb corona and chromosphere. The relative apparent sizes of the Moon and the Sun are different for each eclipse, and techniques are presented which correct the measurements and, therefore, allow direct comparisons between different eclipses. The results show that the fraction of ionising radiation emitted by the limb corona has a clear solar cycle variation and that the underlying trend shows this fraction has been increasing since 1932. Data from the SOHO spacecraft are used to study the effects of short-term variability and it is shown that the observed long-term rise in phi has a negligible probability of being a chance occurrence.

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.

A new source of suprathermal O+ions near the dayside polar cap boundary

A new dayside source of O+ ions for the polar magnetosphere is described, and a statistical survey presented of upward flows of O+ ions using 2 years of data from the retarding ion mass spectrometer (RIMS) experiment on board DE 1, at geocentric distances below 3 RE and invariant latitudes above 40°. The flows are classified according to their spin angle distributions. It is believed that the spacecraft potential near perigee is generally less than +2 V, in which case the entire O+ population at energies below about 60 eV is sampled. Examples are given of field-aligned flow and of transversely accelerated “core” O+ ions; in the latter events a large fraction of the total O+ ion population has been transversely accelerated, and in some extreme cases all the observed ions (of all ion species) have been accelerated, and no residual cold population is observed (“toroidal” distributions). However, by far the most common type of O+ upflow seen by DE RIMS lies near the dayside polar cap boundary (particularly in the prenoon sector) and displays an asymmetric spin angle distribution. In such events the ions carry an upward heat flux, and strong upflow of all species is present (H+, He+, O+, O++, and N+ have all been observed with energies up to about 30 eV, but with the majority of ions below about 2 eV); hence, these have been termed upwelling ion events. The upwelling ions are embedded in larger regions of classical light ion polar wind and are persistently found under the following conditions: at geocentric distances greater than 1.4 RE; at all Kp in summer, but only at high Kp in winter. Low-energy conical ions (<30 eV) are only found near the equatorial edge of the events, the latitude of which moves equatorward with increasing Kp and is highly correlated with the location of field-aligned currents. The RIMS data are fully consistent with a “mass spectrometer effect,” whereby light ions and the more energetic O+ ions flow into the lobes and mantle and hence the far-tail plasma sheet, but lower-energy O+ is swept across the polar cap by the convection electric field, potentially acting as a source for the nightside auroral acceleration regions. The occurrence probability of upwelling ion events, as compared to those of low-altitude transversely accelerated core ions and of field-aligned flow, suggests this could be the dominant mechanism for supplying the nightside auroral acceleration region, and subsequently the ring current and near-earth plasma sheet, with ionospheric O+ ions. It is shown that the total rate of O+ outflow in upwelling ion events (greater than 10^25 s^{−1}) is sufficient for the region near the dayside polar cap boundary to be an important ionospheric heavy ion source.

Mode-specific reactivity of CH4 on Pt(110)-(1Ã2): The concerted role of stretch and bend excitation

The state-resolved reaction probability of CH4 on Pt�110-�1�2 was measured as a function of CH4 translational energy for four vibrational eigenstates comprising different amounts of C-H stretch and bend excitation. Mode-specific reactivity is observed both between states from different polyads and between isoenergetic states belonging to the same polyad of CH4. For the stretch/bend combination states, the vibrational efficacy of reaction activation is observed to be higher than for either pure C-H stretching or pure bending states, demonstrating a concerted role of stretch and bend excitation in C-H bond scission. This concerted role, reflected by the nonadditivity of the vibrational efficacies, is consistent with transition state structures found by ab initio calculations and indicates that current dynamical models of CH4 chemisorption neglect an important degree of freedom by including only C-H stretching motion.

A comparison of methods for treatment selection in seamless phase II/III clinical trials incorporating information on short-term endpoints

In an adaptive seamless phase II/III clinical trial interim analysis, data are used for treatment selection, enabling resources to be focused on comparison of more effective treatment(s) with a control. In this paper, we compare two methods recently proposed to enable use of short-term endpoint data for decision-making at the interim analysis. The comparison focuses on the power and the probability of correctly identifying the most promising treatment. We show that the choice of method depends on how well short-term data predict the best treatment, which may be measured by the correlation between treatment effects on short- and long-term endpoints.

Observations on the flight paths of the day-flying moth Virbia lamae during periods of mate location: do males have a strategy for contacting the pheromone plume?

1. To maximize the probability of rapid contact with a female’s pheromone plume, the trajectories of male foraging flights might be expected to be directed with respect to wind flow and also to be energetically efficient. 2. Flights directed either upwind, downwind, or crosswind have been proposed as optimal strategies for rapid and/or energetically efficient plume contact. Other possible strategies are random and Lévy walks, which have trajectories and turn frequencies that are not dictated by the direction of wind flow. 3. The planar flight paths of males of the day-active moth Virbia lamae were recorded during the customary time of its sexual activity. 4. We found no directional preference in these foraging flights with respect to the direction of contemporaneous wind flow, but, because crosswind encompasses twice the possible orientations of either upwind or downwind, a random orientation is in effect a de facto crosswind strategy. 5. A crosswind preference should be favoured when the plume extends farther downwind than crosswind, and this strategy is realized by V. lamae males by a random orientation of their trajectories with respect to current wind direction