Variabilité climatique et statistiques. Etude par simulation de la puissance et de la robustesse de quelques tests utilisés pour vérifier l'homogénéité de chroniques

L'analyse statistique de séries chronologiques de données hydrométéorologiques est un des outils d'identification de variations climatiques. Cette analyse consiste le plus souvent à la mise en œuvre et à l'interprétation de tests statistiques d'homogénéité des séries. Les séries hydrologiques (données de pluie ou de débit) se caractérisent fréquemment par des effectifs faibles, et ne répondent que rarement aux conditions requises par l'application des tests statistiques dont certains sont paramétriques.Nous avons cherché à évaluer, en terme de puissance et de robustesse, le comportement de quelques méthodes statistiques largement employées dans les études de variabilité climatique. Ce travail a été mené dans chaque cas étudié au moyen de procédures de simulations type Monte-Carlo de 100 échantillons de 50 valeurs conformes aux caractéristiques souvent rencontrées dans les séries naturelles. La variabilité simulée est celle d'un changement brutal de la moyenne. Les procédures concernées sont le test de corrélation sur le rang, le test de Pettitt, le test de Buishand, la procédure bayésienne de Lee et Heghinian, et la procédure de segmentation des séries hydrométéorologiques de Hubert et Carbonnel. Des séries artificielles soit stationnaires, soit affectées par une rupture de la moyenne, normales, non-normales, autocorrélées, présentant une tendance linéaire ou un changement brutal de la variance ont été générées.Les conclusions de ce travail doivent être nuancées selon la méthode considérée. D'une manière générale la puissance maximale estimée se situe autour de 50% pour des taux de rupture de la moyenne de l'ordre de 75% de la valeur de l'écart-type. Par ailleurs il apparaît que l'autocorrélation et la présence d'une tendance dans les séries sont les deux caractéristiques qui pénalisent le plus les performances des procédures.Statistical analysis of hydrometeorological time series is often used to identify climatic variations. Most often this analysis consists of applying and interpreting statistical tests of time series homogeneity. Hydrological time series (rainfall and runoff data) are often short and do not always comply with the hypotheses of the statistical methods. Through simulation we have investigated the power and the robustness of some tests which are widely used in the studies dealing with climatic variability. In each case studied, one hundred samples of fifty elements have been generated based on the main characteristics of natural rainfall series. A shift in the mean has been used to represent a possible climatic variation. The procedures used are the rank correlation test, Pettitt's test, Buishand's test, Lee and Heghinian's bayesian procedure, and Hubert and Carbonnel's segmentation procedure for hydrometeorological series.Each simulation of one hundred samples is used to assess the performances of the methods considering a specific characteristic of the series: normality or non-normality, autocorrelation, trend, shift in the variance. First of all, stationary series have been simulated to evaluate the type I error of the tests. Then the series have been simulated with a break in the mean with different levels of amplitude, from 25% to 100% of the standard deviation value. The rank correlation test, Pettitt's test, Buishand's test and the segmentation procedure with a significance level of 1% (significance level of Scheffé's test) reject as heterogeneous less than ten series over one hundred homogeneous simulated series. This result is consistent with the type I error of a statistical test. On the other hand, Lee and Heghinian's bayesian method rejects about 40% of the series. This result means that this latter procedure must only be applied under the hypothesis of heterogeneity. The estimated power of the methods exceeds 40% to 50% when the break in the mean is more than 75% of the standard deviation value.Independent series have been simulated from normal, log-normal and Pearson distributions to compare the performances of the methods requiring normality. The results show that normality has no significant impact on the performances of these methods. However, the simulations do show that the condition of independence of the successive elements of the series is essential to keep performances constant. Otherwise a trend in the series makes the tests inefficient, except for the rank correlation test for which the alternative is a trend. No method seems to be robust against both negative and positive autoregressive dependencies. The procedures requiring a constant variance are robust when the series keep a constant mean, but seem more or less slightly influenced by a break both in the mean and in the standard deviation

Thermodynamic Density Matrix renormalization Group Study of the Magnetic Susceptibility of Half-integer Quantum Spin Chains

It is shown that White's density matrix renormalization group technique can be adapted to obtain thermodynamic quantities. As an illustration, the magnetic susceptibility of Heisenberg S=1/2 and S=3/2 spin chains are computed. A careful finite size analysis is made to determine the range of temperatures where the results are reliable. For the S=1/2 chain, the comparison with the exact Bethe ansatz curve shows an agreement within 1% down to T=0.05J.Comment: 9 pages, 4 figures. To be published in PR

A second catalog of gamma ray bursts: 1978 - 1980 localizations from the interplanetary network

Eighty-two gamma ray bursts were detected between 1978 September 14 and 1980 February 13 by the experiments of the interplanetary network (Prognoz 7, Venera 11 and 12 SIGNE experiments, Pioneer Venus Orbiter, International Sun-Earth Explorer 3, Helios 2, and Vela). Sixty-five of these events have been localized to annuli or error boxes by the method of arrival time analysis. The distribution of sources is consistent with isotropy, and there is no statistically convincing evidence for the detection of more than one burst from any source position. The localizations are compared with those of two previous catalogs

Limits to the burster repetition rate as deduced from the 2nd catalog of the interplanetary network

The burster repetition rate is an important parameter in many gamma ray burst models. The localizations of the interplanetary network, which have a relatively small combined surface area, may be used to estimate the average repetition rate. The method consists of (1) estimating the number of random overlaps between error boxes expected in the catalog and comparing this number to that actually observed; (2) modeling the response of the detectors in the network, so that the probability of detecting a burst can be estimated; and (3) simulating the arrival of bursts at the network assuming that burster repetition is governed by a Poisson process. The application of this method for many different burster luminosity functions shows that (1) the lower limit to the burster repetition rate depends strongly upon the assumed luminosity function; (2) the best lower limit to the repetition period obtainable from the data of the network is about 100 months; and (3) that a luminosity function for all bursters similar to that of the 1979 Mar 5 burster is inconsistent with the data

De différents aspects de la variabilité de la pluviométrie en Afrique de l'Ouest et Centrale non sahélienne

La sécheresse observée depuis plus d'une vingtaine d'années dans les pays sahéliens se fait également ressentir plus au sud dans des régions d'Afrique aux climats plus humides. Cette baisse de la pluviométrie et la diminution des apports en eau de surface qu'elle entraîne y sont de nature à pénaliser les différents projets de développements liés à l'eau. Le programme ICCARE mené par l'ORSTOM a pour objet l'identification et les conséquences de cette variabilité climatique dans l'ensemble de la zone non sahélienne d'Afrique de l'Ouest et Centrale, en s'appuyant sur les données de deux cents postes pluviométriques et sur un ensemble de méthodes alliant représentations cartographiques et procédures statistiques de détection de ruptures dans les séries chronologiques, univariées et multivariées. La simple étude des séries chronologiques de hauteurs précipitées annuelles fait apparaître une nette et brutale fluctuation du régime pluviométrique dans toute la région considérée, à la fin des années 1960 et au début des années 1970. D'une manière générale, il apparaît que ce sont les zones à régime pluviométrique extrême qui ont subi les modifications les plus importantes : les plus arrosées (de la Guinée à la Côte d'Ivoire) et les plus arides (la bordure sahélienne au nord de la zone étudiée). Entre les deux, le phénomène est d'intensité plus nuancée. Les différentes procédures statistiques appliquées aux séries de hauteurs annuelles précipitées soulignent l'existence d'une rupture survenue à la fin des années 1960 ou au début des années 1970, et donc en phase avec ce qui a été observé et étudié au Sahel. D'autres variables permettant une caractérisation plus "qualitative" du phénomène ont également été étudiées. Elles apportent un complément d'information quant aux manifestations de cette variabilité pluviométrique et montrent que la variabilité climatique se traduit à différents niveaux (durée des saisons des pluies, quantités précipitées hors saisons des pluies, etc.). L'examen des séries chronologiques depuis l'origine des stations a permis de resituer l'événement observé dans une perspective historique faite d'alternances de périodes sèches et de périodes humides. Le phénomène observé à la fin des années 1960 et au début des années 1970 apparaît, cependant comme le plus significatif du point de vue statistique. Si les causes premières d'apparition du phénomène sont, à l'heure actuelle, encore insuffisamment expliquées, et ce même si certaines activités humaines y ont, sans aucun doute, contribué, cette baisse de la pluviométrie a, bien entendu, des conséquences importantes sur la disponibilité des ressources en eau dans ces régions. Si la carence pure et simple n'est pas à craindre dans ces régions où les quantités précipitées restent importantes dans l'absolu, les effets de cette variabilité climatique peuvent, malgré tout, se révéler désastreux, en ce sens qu'ils modifient les données d'un équilibre déjà souvent mis à mal par ailleurs (pression anthropique et déforestation par exemple).The drought observed for more than twenty years in the sahelian countries has also affected those located more to the South with more humid climates (SUTCLIFFE and KNOTT, 1987; NICHOLSON et al, 1988; MAHE and OLIVRY, 1991; OLIVRY et al, 1993 PATUREL et al, 1995). The decrease in rainfall and consequently that in runoff might penalise development projects linked with water supply. The ICCARE programme led by ORSTOM aimed at identifying and measuring the consequences of this climatic variability in the non-sahelian parts of the West and Central African region as a whole. The study was based on rainfall data from more than two hundred stations, break detection methods in the time series as well as cartographic tools were used. This study allowed to highlight the manifestations of the climatic variability observed for nearly 25 years in West and Central Africa. Whereas it had been thought for a long time that the rainfall deficit was restricted to the sahelian region, this study showed that the forest covered regions were also affected and generally speaking the so-called 'humid Africa'. The decrease in rainfall entails that in runoff and thus a change in water resources availability which is the cornerstone of a fairly great number of development projects. Hydrological regimes variability and possible modifications of rainfall-runoff relationship are to be the next stages of the ICCARE programme, partial results were already published (SERVAT et al, 1997). A simple study of the time series showed straightforward fluctuations of the rainfall patterns (figures 2 and 3), which happened during the late sixties or the early seventies over the whole region. The mapping of the time series analysis depicts a clear general trend towards a shift of the isohyets to the south-south-west from the 1950's to the 1980's. This shift reflects a sharp drop in the annual rainfall over the whole of the non-sahelian west and central Africa. Generally speaking it appears that the zones with extreme rainfall patterns underwent the most important modifications, namely : the wettest ones from Guinea to Ivory Coast and the driest ones, bordering the sahelian region to the North of the studied area. Elsewhere changes are less drastic (SERVAT et al, 1996). The different statistical procedures applied to the series of annual rainfall showed breaks taking place during the late 1960's or the early 1970's (figure 4), which was in keeping with what had been observed in the sahelian region. Rainfall deficits were in the order of 20% and they could reach values higher than 25% (table 1), in particular along the Atlantic Coast or in the North, which upholds the fact that 'humid Africa' was also severely affected by the rainfall variability. Other variables which allow a more qualitative characterisation of the phenomenon were also studied (table 2) (figures 5 and 6) ( PATUREL et al, 1997, SERVAT et al, 1997). They brought complementary information about the ways the rainfall variability expresses itself. The pattern of the rainy season was slightly different from what it had been before the 1970's, its length was generally shorter either because it started later than before or because it ended earlier. Likewise, the rainfall distribution was modified, which resulted in a more 'homogeneous' pattern for the zone with only one rainy season and in a sensible change in the ratio of the rainfall heights of the two rainy seasons. Some of the regions of the so-called 'wooded' savannah saw a modification of their climate with a shift from a 'guinean' climate toward a 'sudanese' one. In west Africa some regions saw also a decrease of the amount of precipitation occurring outside the rainy season, which led to a strengthening of the dry season and contributed, if need be, to the perception of the phenomenon by local populations. The decrease in the number of rainy days, where it was possible to study it, was in line with the rainfall deficit. A complementary statistical approach was carried out, it consists of a spatio-temporal study using a multidimensional exploratory analysis (KHODJA et al, 1998). This led to using a multivariate test for detecting a shift in the mean value. This approach confirmed the results obtained with univariate analyses whether it is for the time location of the break (late 1960's, early 1970's) or whether for the heterogeneous character of the phenomenon both from a spatial or temporal point of view. A major characteristic of this persisting rainfall deficit seems to be the existence of two axes of heterogeneity along the north-south and east-west directions (table 3) (figure 7). The survey of the rainfall time series from the origin allowed to place the drought in a historical perspective. So, it appears that, since the beginning of the century, the region underwent a succession of dry and wet periods, although it is difficult to speak of cycles. The phenomenon observed during the late 60's and the early 70's appears, however, as the most significant from a statistical point of view. Besides, the still lasting period of deficit has displayed a length and an intensity quite remarkable, in particular in the north and west sectors of the studied zone where the phenomenon presents an even more exceptional character (figures 8,9,10 and 11). Even though what brought about this diminution in rainfall remains, until now, unexplained, certain human activities undoubtedly contributed to the aggravation of the phenomenon. Although deforestation cannot be held entirely responsible for the drought, the fact remains that overlogging helped to increase the rainfall deficit in numerous regions along the Atlantic coast and the Guinean gulf. Of course, this rainfall deficit has important consequences on the availability of water resources in those regions. Agriculture, the filling of dams and therefore the hydroelectric production, to mention only a few domains, are strongly penalised by this decrease in resources. If a real shortage is not to be feared in those regions where the quantities of precipitation remain high in absolute terms, the effects of that variability can still prove to be disastrous, as they modify the elements of a balance that is already threatened by other factors (anthropic pressure and deforestation, for example). The ICCARE programme, which goes on with the study of the modifications of the river hydrological regimes, will give answers as to the effect of the rainfall deficit upon water resources availability

Critical dynamics in thin films

Critical dynamics in film geometry is analyzed within the field-theoretical approach. In particular we consider the case of purely relaxational dynamics (Model A) and Dirichlet boundary conditions, corresponding to the so-called ordinary surface universality class on both confining boundaries. The general scaling properties for the linear response and correlation functions and for dynamic Casimir forces are discussed. Within the Gaussian approximation we determine the analytic expressions for the associated universal scaling functions and study quantitatively in detail their qualitative features as well as their various limiting behaviors close to the bulk critical point. In addition we consider the effects of time-dependent fields on the fluctuation-induced dynamic Casimir force and determine analytically the corresponding universal scaling functions and their asymptotic behaviors for two specific instances of instantaneous perturbations. The universal aspects of nonlinear relaxation from an initially ordered state are also discussed emphasizing the different crossovers that occur during this evolution. The model considered is relevant to the critical dynamics of actual uniaxial ferromagnetic films with symmetry-preserving conditions at the confining surfaces and for Monte Carlo simulations of spin system with Glauber dynamics and free boundary conditions.Comment: 64 pages, 21 figure

Fishmeal supplementation during ovine pregnancy and lactation protects against maternal stress-induced programming of the offspring immune system

Background: Prenatally stressed offspring exhibit increased susceptibility to inflammatory disorders due to in utero programming. Research into the effects of n-3 PUFAs shows promising results for the treatment and prevention of these disorders. The purpose of this study was to investigate whether maternal fishmeal supplementation during pregnancy and lactation protects against programming of the offspring\u27s immune response following simulated maternal infection. Methods: In order to accomplish this, 53 ewes were fed a diet supplemented with fishmeal (FM; rich in n-3 PUFA) or soybean meal (SM; rich in n-6 PUFAs) from day 100 of gestation (gd 100) through lactation. On gd135, half the ewes from each dietary group were challenged with either 1.2 μg/kg Escherichia coli lipopolysaccharide (LPS) endotoxin to simulate a bacterial infection, or saline as the control. At 4.5 months of age the offspring\u27s dermal immune response was assessed by cutaneous hypersensitivity testing with ovalbumin (OVA) and candida albicans (CAA) 21 days after sensitization. Skinfold measurements were taken and serum blood samples were also collected to assess the primary and secondary antibody immune response. Results: Offspring born to SM + LPS mothers had a significantly greater change in skinfold thickness in response to both antigens as well as a greater secondary antibody response to OVA compared to all treatments. Conclusions: Supplementation during pregnancy with FM appears to protect against adverse fetal programming that may occur during maternal infection and this may reduce the risk of atopic disease later in life