Influence of the seasonal cycle on the termination of El Nino events in a coupled general circulation model

In this study, the mechanisms leading to the El Nino peak and demise are explored through a coupled general circulation model ensemble approach evaluated against observations. The results here suggest that the timing of the peak and demise for intense El Nino events is highly predictable as the evolution of the coupled system is strongly driven by a southward shift of the intense equatorial Pacific westerly anomalies during boreal winter. In fact, this systematic late-year shift drives an intense eastern Pacific thermocline shallowing, constraining a rapid El Nino demise in the following months. This wind shift results from a southward displacement in winter of the central Pacific warmest SSTs in response to the seasonal evolution of solar insolation. In contrast, the intensity of this seasonal feedback mechanism and its impact on the coupled system are significantly weaker in moderate El Nino events, resulting in a less pronounced thermocline shallowing. This shallowing transfers the coupled system into an unstable state in spring but is not sufficient to systematically constrain the equatorial Pacific evolution toward a rapid El Nino termination. However, for some moderate events, the occurrence of intense easterly wind anomalies in the eastern Pacific during that period initiate a rapid surge of cold SSTs leading to La Nina conditions. In other cases, weaker trade winds combined with a slightly deeper thermocline allow the coupled system to maintain a broad warm phase evolving through the entire spring and summer and a delayed El Nino demise, an evolution that is similar to the prolonged 1986/87 El Nino event. La Nina events also show a similar tendency to peak in boreal winter, with characteristics and mechanisms mainly symmetric to those described for moderate El Nino cases.

Cone model-based SEP event calculations for applications to multipoint observations

The problem of modeling solar energetic particle (SEP) events is important to both space weather research and forecasting, and yet it has seen relatively little progress. Most important SEP events are associated with coronal mass ejections (CMEs) that drive coronal and interplanetary shocks. These shocks can continuously produce accelerated particles from the ambient medium to well beyond 1 AU. This paper describes an effort to model real SEP events using a Center for Integrated Space weather Modeling (CISM) MHD solar wind simulation including a cone model of CMEs to initiate the related shocks. In addition to providing observation-inspired shock geometry and characteristics, this MHD simulation describes the time-dependent observer field line connections to the shock source. As a first approximation, we assume a shock jump-parameterized source strength and spectrum, and that scatter-free transport occurs outside of the shock source, thus emphasizing the role the shock evolution plays in determining the modeled SEP event profile. Three halo CME events on May 12, 1997, November 4, 1997 and December 13, 2006 are used to test the modeling approach. While challenges arise in the identification and characterization of the shocks in the MHD model results, this approach illustrates the importance to SEP event modeling of globally simulating the underlying heliospheric event. The results also suggest the potential utility of such a model for forcasting and for interpretation of separated multipoint measurements such as those expected from the STEREO mission.

An event-based approach to validating solar wind speed predictions: high-speed enhancements in the Wang-Sheeley-Arge model

One of the primary goals of the Center for Integrated Space Weather Modeling (CISM) effort is to assess and improve prediction of the solar wind conditions in near‐Earth space, arising from both quasi‐steady and transient structures. We compare 8 years of L1 in situ observations to predictions of the solar wind speed made by the Wang‐Sheeley‐Arge (WSA) empirical model. The mean‐square error (MSE) between the observed and model predictions is used to reach a number of useful conclusions: there is no systematic lag in the WSA predictions, the MSE is found to be highest at solar minimum and lowest during the rise to solar maximum, and the optimal lead time for 1 AU solar wind speed predictions is found to be 3 days. However, MSE is shown to frequently be an inadequate “figure of merit” for assessing solar wind speed predictions. A complementary, event‐based analysis technique is developed in which high‐speed enhancements (HSEs) are systematically selected and associated from observed and model time series. WSA model is validated using comparisons of the number of hit, missed, and false HSEs, along with the timing and speed magnitude errors between the forecasted and observed events. Morphological differences between the different HSE populations are investigated to aid interpretation of the results and improvements to the model. Finally, by defining discrete events in the time series, model predictions from above and below the ecliptic plane can be used to estimate an uncertainty in the predicted HSE arrival times.

Transient teleconnection event at the onset of a planet-encircling dust storm on Mars

We use proper orthogonal decomposition (POD) to study a transient teleconnection event at the onset of the 2001 planet-encircling dust storm on Mars, in terms of empirical orthogonal functions (EOFs). There are several differences between this and previous studies of atmospheric events using EOFs. First, instead of using a single variable such as surface pressure or geopotential height on a given pressure surface, we use a dataset describing the evolution in time of global and fully three-dimensional atmospheric fields such as horizontal velocity and temperature. These fields are produced by assimilating Thermal Emission Spectrometer observations from NASA's Mars Global Surveyor spacecraft into a Mars general circulation model. We use total atmospheric energy (TE) as a physically meaningful quantity which weights the state variables. Second, instead of adopting the EOFs to define teleconnection patterns as planetary-scale correlations that explain a large portion of long time-scale variability, we use EOFs to understand transient processes due to localised heating perturbations that have implications for the atmospheric circulation over distant regions. The localised perturbation is given by anomalous heating due to the enhanced presence of dust around the northern edge of the Hellas Planitia basin on Mars. We show that the localised disturbance is seemingly restricted to a small number (a few tens) of EOFs. These can be classified as low-order, transitional, or high-order EOFs according to the TE amount they explain throughout the event. Despite the global character of the EOFs, they show the capability of accounting for the localised effects of the perturbation via the presence of specific centres of action. We finally discuss possible applications for the study of terrestrial phenomena with similar characteristics.

Numerical simulation of the 12 May 1997 CME Event: The role of magnetic reconnection

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the 12 May 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance (≈1R ) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

Action following the discovery of a global association between the whole genome and adverse event risk in a clinical drug-development programme

Observation of adverse drug reactions during drug development can cause closure of the whole programme. However, if association between the genotype and the risk of an adverse event is discovered, then it might suffice to exclude patients of certain genotypes from future recruitment. Various sequential and non-sequential procedures are available to identify an association between the whole genome, or at least a portion of it, and the incidence of adverse events. In this paper we start with a suspected association between the genotype and the risk of an adverse event and suppose that the genetic subgroups with elevated risk can be identified. Our focus is determination of whether the patients identified as being at risk should be excluded from further studies of the drug. We propose using a utility function to? determine the appropriate action, taking into account the relative costs of suffering an adverse reaction and of failing to alleviate the patient's disease. Two illustrative examples are presented, one comparing patients who suffer from an adverse event with contemporary patients who do not, and the other making use of a reference control group. We also illustrate two classification methods, LASSO and CART, for identifying patients at risk, but we stress that any appropriate classification method could be used in conjunction with the proposed utility function. Our emphasis is on determining the action to take rather than on providing definitive evidence of an association. Copyright (C) 2008 John Wiley & Sons, Ltd.

Changes in milk composition and processing properties during the spring flush period

The period, known to UK farmers and processors as the "spring flush", when the cows' diet changes from dry feed to spring pasture, has long been established as a time of change in milk properties and processing characteristics. Although it is believed to be a time when problems in processing are most likely to occur (e.g. milk that does not form clots or forms weak gels during cheesemaking), there is little evidence in the literature of detailed changes in milk composition and their impact on product manufacture. In this study, a range of physicochemical properties were analysed in milk collected from five commercial dairy herds before, during and after the spring flush period of 2006. In particular, total and ionic calcium contents of milk were studied in relation to other parameters including rennet clotting, acid gel properties, heat coagulation, alcohol stability, micelle size and zeta potential. Total divalent cations were significantly reduced from 35.4 to 33.4 mmol.L-1 during the study, while ionic calcium was reduced from 1.48 to 1.40 mmol.L-1 over the same period. Many parameters varied significantly between the sample dates. However, there was no evidence to suggest that any of the milk samples would have been unsuitable for processing - e.g. there were no samples that did not form clots with chymosin within a reasonable time or formed especially weak rennet or acid gels. A number of statistically significant correlations were found within the data, including ionic calcium concentration and pH; rennet clotting time (RCT) and micelle diameter; and RCT and ethanol stability. Overall, while there were clear variations in milk composition and properties over this period, there was no evidence to support the view that serious processing problems are likely during the change from dry feed to spring pasture.