1 9 th international conference zaragoza-pau on applied mathematics and statistics on heat wave...
TRANSCRIPT
![Page 1: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/1.jpg)
1
9th International Conference Zaragoza-Pau
on Applied Mathematics and Statistics
On heat wave definitionAbaurrea J., Cebrián A.C., Asín J., Centelles A.
19-21 September 2005
![Page 2: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/2.jpg)
2
Introduction (1)Heat waves have not a standard operational definition
Usual approach to define them:
• An excess of daily maximum temperature, Tx, over a fixed threshold (POT)
Other conditions required:
• A minimum duration of the event
• The daily minimum temperature exceeds another threshold
Problems and evidences not considered:
• Greater effects on morbidity and mortality of heat waves occurring in the early summer
• Possibility of longer heat waves including intermediate “cool” days
![Page 3: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/3.jpg)
3
• Kysely (2000):
A group of consecutive days is considered a heat wave if:
a) Tx ≥ T1 for, at least, three days
b) Tx ≥ T2 for every day
c) Mean (Tx) ≥ T1
Tx: Daily maximum temperature
T1: Threshold for hot days
T2: Threshold for warm days
For Central Europe T1 = 30ºC and T2 = 25ºC
Introduction (2)
![Page 4: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/4.jpg)
4
a) and b) as Kysely
c) Mean(Tx) T1 for the whole period and for each partial sequence, HC, HCHC, etc., where H stands for a hot spell and C a cool spell
d) The length and the area (the accumulated sum of differences to T1) for each cool spell included in the wave, must be lower than the corresponding 90th percentile in the reference period
T1 and T2 are, respectively, the 95th and 50th percentiles of Tx values observed in June, July and August, in the reference period 1961-1990
A time period is considered a heat wave if:
Introduction (3)• Abaurrea et al. (2004):
![Page 5: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/5.jpg)
5
Heat wave: A period of arbitrary length where Tx exceeds a “shot temperature”Shot temperature = extraordinary increase of mortality
– Madrid (36.5ºC)– Barcelona (30.3ºC)
Introduction (4)• Díaz et al. (2003):
–Sevilla (41ºC)
–Lisboa (33.5ºC)
These thresholds are the 95th percentiles of the corresponding Tx value distributions in JJAS, 1991-2002
In this way, they obtain the heat wave thresholds for main Spanish towns: Zaragoza (37.3ºC), Huesca (36.1ºC)
![Page 6: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/6.jpg)
6
Data and preliminary analysis
•Zaragoza and Huesca
•Daily Tx and Tn data for 1951-2004
•Daily mortality data for 1975-2002 (people aged 65 or more years)
![Page 7: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/7.jpg)
7
Temperature
• Tx and Tn evolution during the studied period
1951-75 stability 1976-90 increase
1991-96 stability 1997-2004 increase
Data and preliminary analysis
![Page 8: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/8.jpg)
8
• Tx evolution in different summer periods
TemperatureData and preliminary analysis
Lowess (40) of Tx data in 7 overlapping 5-week intervals
![Page 9: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/9.jpg)
9 Lowess (30) of Mr corresponding to seven overlapping 10
year long intervals
Mr decreases between 1975 and 2002
Change in the seasonal profile
MortalityData and preliminary analysis
• Mr: Daily mortality rate per 1000 inhabitants
![Page 10: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/10.jpg)
10
•The effect of hot temperatures on mortality occurs in the short term (1-3 days) (Díaz et al. 2005)
•For daily variables, the correlation between Tx and Mr is maximum with a 24 hour delay
•The greatest correlation between 3 days averaged values is also obtained for a 24 hour delay
Temperature-Mortality relationshipData and preliminary analysis
![Page 11: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/11.jpg)
11
Data and preliminary analysis
The impact on Mr of a fixed high temperature changes in time and along the summer
To show this property we select a temperature value, 33.3ºC, the 90th percentile of daily Tx values in June, 1975-81
•1975-02 is divided into four 7-year periods and we consider observations from June, July and August
Temperature-Mortality relationship
![Page 12: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/12.jpg)
12
Data and preliminary analysis
•Decrease of Mr 90th percentile and mean values in time and along the summer
•33.3ºC is a critical value, regarding the Mr response, for the first period and it is not for the last one
Temperature-Mortality relationship
![Page 13: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/13.jpg)
13
Data and preliminary analysis
•Decrease of Mr 90th percentile and mean values in time and along the summer
•33.3ºC is a critical value, regarding the Mr response, for the first period and it is not for the last one
Temperature-Mortality relationship
![Page 14: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/14.jpg)
14
Mortality ExcessWe define the mortality excess, Mex(t), in day t as the difference between the number of deaths, Mf(t), and its expected value, Me(t)
Mex (t) = Mf(t) –Me(t)
The expected mortality is obtained by fitting a regression model including:
a) Time terms until the second order (long term evolution)
b) Harmonic terms until the fourth order (seasonality)
c) Dummy variables for indicating the periods 75-86, 87-96 and 97-02, in order to fit different seasonal patterns
![Page 15: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/15.jpg)
15
Expected and observed mortality lowess
Mortality Excesss
![Page 16: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/16.jpg)
16
•We try to identify the ‘shot temperature’ for each 7-year period and summer interval, looking for the change point in the lowess smoother of Mex vs. Tx
PROBLEMS
Smoothed curves are frequently erratic because of small sample sizes
A proper shot temperature doesn’t appear in many graphs
Threshold selection
![Page 17: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/17.jpg)
17
Threshold selection
![Page 18: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/18.jpg)
18
Threshold selection
![Page 19: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/19.jpg)
19
High excess crossing temperatures increase in time but remain constant when they are transformed to a percentile scale
Threshold selection
![Page 20: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/20.jpg)
20
Threshold selection processThe process to define T1 threshold needs several steps
a) Analysed interval: 14-May to 16-September in 1975-2002
b) Four 7-year periods (1975-81, 1982-88, 1989-95, 1996-2002)
c) Several divisions of the 14/5-16/9 interval using different
length cells: 3-weeks, 4-weeks, month,...
![Page 21: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/21.jpg)
21
d) Identifying 1.25-excess crossing temperature percentiles
Threshold selection process
![Page 22: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/22.jpg)
22
Threshold selection process
e) Percentile-Threshold allocation to 11 selected dates along the summer
1.25-excess crossing temperature percentile values
![Page 23: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/23.jpg)
23
f) Transformation of the percentile-threshold into its equivalent temperature-threshold (T1) using an adequate probabilistic distribution
Threshold selection process
Probabilistic distributions:
N: Normal
LN: Lognormal
W: Weibull
EV: Extreme Value
L: Logistic
LL: Log-Logistic
![Page 24: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/24.jpg)
24
g) Estimation of the daily T1 threshold for each 7-year period
Threshold selection process
![Page 25: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/25.jpg)
25
•4th June: 32.3-34ºC•27th August: 37-39.8ºC
•The increase of T1 along the period 1975-2002 is about 2ºC
•The bigger slope of the 3rd period is due to temperature warming in August and July, whereas the smaller slope of the 4th period is due to strong temperature increase in June
Threshold selection process
T1 thresholds for the period 1975-2002
![Page 26: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/26.jpg)
26
in comparison with the T1-Díaz performance
Results and conclusions
Evaluation of T1- threshold results
![Page 27: 1 9 th International Conference Zaragoza-Pau on Applied Mathematics and Statistics On heat wave definition Abaurrea J., Cebrián A.C., Asín J., Centelles](https://reader036.vdocument.in/reader036/viewer/2022062712/56649c7b5503460f9492f52d/html5/thumbnails/27.jpg)
27
References Abaurrea, J., et al. (2004). Modelling hot extreme temperature events using a non homogeneous Poisson model. 6th World Congress of Bernoulli Society for Mathematical Statistics and Probability, Barcelona.
Díaz, J., et al. (2002). Effects of extremely hot days on people older than 65 years in Seville (Spain) from 1986 to 1997. Int. J. Biometeorology, 46, 145-149.
Díaz, J., Linares, C., García-Herrera, R. (2005). Impacto de las temperaturas extremas en la salud pública. Futuras actuaciones. Rev. Esp. Salud Pública, 79, 145-157.
Kysely, J. (2002). Temporal fluctuations in heat waves at Prague, the Czech republic, from 1901-97 and their relationship to atmospheric circulation. Int. J. Climatol., 22, 33-50.
Robinson, P. J. (2001). On the definition of a heat wave. J. of Applied Meteorology, 40, 762-75.