STATISTICAL DOWNSCALING OF DAILY PRECIPITATION IN THE ARGENTINE PAMPAS REGION Bettolli ML- Penalba OC Department of Atmospheric and Ocean Sciences, University of Buenos Aires, Argentina National Council of Research and Development (CONICET), Argentina bettolli@at.fcen.uba.ar ICRC-CORDEX 2016
MOTIVATION The Pampas covers the most productive agricultural land in Argentina The main crops are: Soybean, corn, wheat and sunflower 10 0-10 -20-50 Pampas Region -60-90 -80-70 -60-50
MOTIVATION Argentina is the largest exporter of soybean products and in terms of production size comes third in the world after the United States and Brazil 10 0-10 -20-50 Pampas Region -60-90 -80-70 -60-50
MOTIVATION Crops grow under rain-fed conditions. Precipitation is one of the principal climatic variables of influence for the production and management of crops. 10 0-10 -20-50 Pampas Region -60-90 -80-70 -60-50
MOTIVATION CORDEX-ESD-Experiment 1 was set and performed in the La Plata Basin region. Prior to this experiment, empirical statistical downscaling on a daily basis had not been widely developed or applied over the region. -15-20 -45-50 -55 La Plata Basin Pampas Region South America -75-70 -65-60 -55-50 -45
MOTIVATION There is a need to promote ESD applications in southern South American regions. -15-20 La Plata Basin -45-50 Pampas Region South America -55-75 -70-65 -60-55 -50-45
OBJECTIVE To calibrate and validate a statistical method to downscale daily precipitation in the Argentine Pampas region.
DATA AND METHODOLOGY 0 Predictand: Daily precipitation. 28 weather stations. Predictors: NCEP R2 daily fields. Three different domains (orange rectangles) -10-20 -50 MSLP hgt 500 hpa hgt 850 hpa u 850 hpa v 850 hpa v 200 hpa T 850 hpa rhum 850 hpa -60
DATA AND METHODOLOGY 0-10 -20-50 -60-90 -80-70 -60-50 Analogue Method: Unrestricted analogues were chosen in a different calendar + year - 1 year Euclidean distance Different domains and combinations of variables were evaluated Calibration: 1979-2000 Validation: 2001-2010
DATA AND METHODOLOGY 0-10 -20 Different aspects of daily precipitation were evaluated through different metrics: Precipitation ocurrence Precipitation amount -50-60 -90-80 -70-60 -50
Occurrence Index I: (Timbal et al. 2003) Assesses false alarms, missed forecasts and hits of rainy days The higher the index, the better the method performance
45 Occurrence 40 Index I 35 30 25 20 15 10 5 0 MSLPLD hgt500ld u850ld v850ld u850ld_v850ld T850LD rhum850ld hgt500ld_t850ld MSLPLD_T850LD hgt500ld_rhum850ld u850ld_v850ld_t850ld hgt500md MSLPMD u850md v850md T850MD rhum850md hgt500md_t850md Combinations of Predictor Variables and Domains Median 25%-75% Min-Max hgt500md_rhum850md MSLPMD_T850MD MSLPMD_rhum850MD u850md_v850md_t850md u850md_v850md_rhum850md hgt500ld_t850md hgt500ld_rhum850md MSLPLD_T850MD MSLPSD_T850SD hgt500sd_t850sd u850sd_v850sd_t850sd u850sd_v850sd_rhum850sd u850sd_v850sd_rhum850sd_t850sd MSLPSD_T850SD_rhum850SD MSLPMD_T850SD_rhum850SD u850md_v850md_t850sd_rhum850sd
Each group of 5 boxplots correspond to Summer, Fall, Winter, Spring and Annual Index I. 45 Occurrence 40 Index I Boxplots of the index I values between daily estimated and observed values over the 28 meteorological stations. 35 30 25 20 15 10 5 0 MSLPLD hgt500ld u850ld v850ld u850ld_v850ld Median 25%-75% Min-Max T850LD rhum850ld hgt500ld_t850ld MSLPLD_T850LD hgt500ld_rhum850ld u850ld_v850ld_t850ld hgt500md MSLPMD u850md v850md T850MD rhum850md hgt500md_t850md Combinations of Predictor Variables and Domains hgt500md_rhum850md MSLPMD_T850MD MSLPMD_rhum850MD u850md_v850md_t850md u850md_v850md_rhum850md hgt500ld_t850md hgt500ld_rhum850md MSLPLD_T850MD MSLPSD_T850SD hgt500sd_t850sd u850sd_v850sd_t850sd u850sd_v850sd_rhum850sd u850sd_v850sd_rhum850sd_t850sd MSLPSD_T850SD_rhum850SD MSLPMD_T850SD_rhum850SD u850md_v850md_t850sd_rhum850sd
Each group of 5 boxplots correspond to Summer, Fall, Winter, Spring and Annual Index I. 45 Occurrence 40 Index I Boxplots of the index I values between daily estimated and observed values over the 28 meteorological stations. 35 30 25 20 15 10 5 0 MSLPLD hgt500ld u850ld v850ld u850ld_v850ld Median 25%-75% Min-Max T850LD rhum850ld hgt500ld_t850ld MSLPLD_T850LD hgt500ld_rhum850ld u850ld_v850ld_t850ld hgt500md MSLPMD u850md v850md T850MD rhum850md hgt500md_t850md Combinations of Predictor Variables and Domains hgt500md_rhum850md MSLPMD_T850MD MSLPMD_rhum850MD u850md_v850md_t850md u850md_v850md_rhum850md hgt500ld_t850md hgt500ld_rhum850md MSLPLD_T850MD MSLPSD_T850SD hgt500sd_t850sd u850sd_v850sd_t850sd u850sd_v850sd_rhum850sd u850sd_v850sd_rhum850sd_t850sd MSLPSD_T850SD_rhum850SD MSLPMD_T850SD_rhum850SD u850md_v850md_t850sd_rhum850sd
Each group of 5 boxplots correspond to Summer, Fall, Winter, Spring and Annual Index I. 45 Occurrence 40 Index I Boxplots of the index I values between daily estimated and observed values over the 28 meteorological stations. 35 30 25 20 15 10 5 0 The skill is improved when zonal and meridional wind components are considered east of the Andes range: The Andes range canalizes the moisture advection to the region of interest. Median The reanalyisis performance 25%-75% near the Andes range is affected particularly at low levels. Min-Max (Timbal 2004) MSLPLD hgt500ld u850ld v850ld u850ld_v850ld T850LD rhum850ld hgt500ld_t850ld MSLPLD_T850LD hgt500ld_rhum850ld u850ld_v850ld_t850ld hgt500md MSLPMD u850md v850md T850MD rhum850md hgt500md_t850md hgt500md_rhum850md MSLPMD_T850MD MSLPMD_rhum850MD u850md_v850md_t850md u850md_v850md_rhum850md hgt500ld_t850md hgt500ld_rhum850md MSLPLD_T850MD MSLPSD_T850SD hgt500sd_t850sd u850sd_v850sd_t850sd u850sd_v850sd_rhum850sd u850sd_v850sd_rhum850sd_t850sd MSLPSD_T850SD_rhum850SD MSLPMD_T850SD_rhum850SD u850md_v850md_t850sd_rhum850sd
Each group of 5 boxplots correspond to Summer, Fall, Winter, Spring and Annual Index I. 45 Occurrence 40 Index I Boxplots of the index I values between daily estimated and observed values over the 28 meteorological stations. 35 30 25 20 15 10 5 0 In all cases: the highest Index I values are found in Winter and the lowest in Summer (Timbal 2004) MSLPLD hgt500ld u850ld v850ld u850ld_v850ld Median 25%-75% Min-Max T850LD rhum850ld hgt500ld_t850ld MSLPLD_T850LD hgt500ld_rhum850ld u850ld_v850ld_t850ld hgt500md MSLPMD u850md v850md T850MD rhum850md hgt500md_t850md hgt500md_rhum850md MSLPMD_T850MD MSLPMD_rhum850MD u850md_v850md_t850md u850md_v850md_rhum850md hgt500ld_t850md hgt500ld_rhum850md MSLPLD_T850MD MSLPSD_T850SD hgt500sd_t850sd u850sd_v850sd_t850sd u850sd_v850sd_rhum850sd u850sd_v850sd_rhum850sd_t850sd MSLPSD_T850SD_rhum850SD MSLPMD_T850SD_rhum850SD u850md_v850md_t850sd_rhum850sd
Each group of 5 boxplots correspond to Summer, Fall, Winter, Spring and Annual RMSE. Amount RMSE Boxplots of RMSE values between daily estimated and observed values over the 28 meteorological stations. 22 20 18 16 14 12 10 8 6 4 2 MSLPLD hgt500ld u850ld v850ld u850ld_v850ld T850LD rhum850ld hgt500ld_t850ld MSLPLD_T850LD hgt500ld_rhum850ld u850ld_v850ld_t850ld hgt500md MSLPMD u850md v850md T850MD rhum850md hgt500md_t850md Combinations of Predictor Variables and Domains Median 25%-75% Min-Max hgt500md_rhum850md MSLPMD_T850MD MSLPMD_rhum850MD u850md_v850md_t850md u850md_v850md_rhum850md hgt500ld_t850md hgt500ld_rhum850md MSLPLD_T850MD MSLPSD_T850SD hgt500sd_t850sd u850sd_v850sd_t850sd u850sd_v850sd_rhum850sd u850sd_v850sd_rhum850sd_t850sd MSLPSD_T850SD_rhum850SD MSLPMD_T850SD_rhum850SD u850md_v850md_t850sd_rhum850sd
Each group of 5 boxplots correspond to Summer, Fall, Winter, Spring and Annual RMSE. Amount RMSE Boxplots of RMSE values between daily estimated and observed values over the 28 meteorological stations. 22 20 18 16 14 12 10 8 6 4 2 RMSE decrease when zonal and meridional wind components are considered Combinations of Predictor Variables and Domains east of the Andes range. The lowest RMSE values are found in Winter MSLPLD hgt500ld u850ld v850ld u850ld_v850ld Median 25%-75% Min-Max T850LD rhum850ld hgt500ld_t850ld MSLPLD_T850LD hgt500ld_rhum850ld u850ld_v850ld_t850ld hgt500md MSLPMD u850md v850md T850MD rhum850md hgt500md_t850md hgt500md_rhum850md MSLPMD_T850MD MSLPMD_rhum850MD u850md_v850md_t850md u850md_v850md_rhum850md hgt500ld_t850md hgt500ld_rhum850md MSLPLD_T850MD MSLPSD_T850SD hgt500sd_t850sd u850sd_v850sd_t850sd u850sd_v850sd_rhum850sd u850sd_v850sd_rhum850sd_t850sd MSLPSD_T850SD_rhum850SD MSLPMD_T850SD_rhum850SD u850md_v850md_t850sd_rhum850sd
Combination: u + v + rhum + T 850 hpa Small Domain Validation Period: 2001-2010
Interannual variability Summer precipitation frequency Correlation values 0.49 0.44 0.63 0.75 0.52 0.6 0.63 0.52 0.5 0.62 0.67 0.54 0.6 0.29 0.78 0.58 0.54 0.61 0.59 0.53 0.67 0.56 0.52 0.480.66 0.68 0.47-70 -65-60 -55
Interannual variability Summer precipitation frequency Correlation values 0.49 0.44 0.63 0.75 0.52 0.6 0.63 0.52 0.5 0.62 0.67 0.54 0.6 0.29 0.78 0.58 0.54 0.61 0.59 0.53 0.67 0.56 0.52 0.480.66 0.68 0.47 Good general agreement in the validation period -70-65 -60-55
Interannual variability Summer precipitation amount Correlation values 0.43 0.61 0.4 0.26 0.42 0.57 0.48 0.57 0.16 0.51 0.56 0.42 0.55 0.4 0.62 0.53 0.47 0.54 0.54 0.51 0.43 0.55 0.27 0.24 0.220.4 0.28 0.37-70 -65-60 -55
Interannual variability Summer precipitation amount Correlation values 0.43 0.61 0.4 0.26 0.42 0.57 0.48 0.57 0.16 0.51 0.56 0.42 0.55 0.4 0.62 0.53 0.47 0.54 0.54 0.51 0.43 0.55 0.27 0.24 0.220.4 0.28 0.37 The agreement depends on the station -70-65 -60-55
Wet spells: Winter Summer -70-65 -60-55
Wet spells: Winter Summer The probability of wet spells is well captured except for Santa Rosa St in Winter (low precipitation in Winter) -70-65 -60-55
Distributions: daily precipiation amount Winter 600 500 K-S Test, no difference 95% confidence Concordia St Observed Concordia St Estimated 600 500 Summer Concodia St Observed Concordia St Estimated 400 400 No of obs 300 No of obs 300 200 200 100 100 600 500 0-20 0 20 40 60 80 100 120 140 160 180 200 220 Daily Precipitation (mm) Santa Rosa St Observed Santa Rosa St Estimated 0 600 500-20 0 20 40 60 80 100 120 140 160 180 200 220 Daily Precipitation (mm) Santa Rosa St Observed Santa Rosa St Estimated No of obs 400 300-70 -65-60 -55 No of obs 400 300 200 200 100 100 0-10 0 10 20 30 40 50 60 70 80 90 0-10 0 10 20 30 40 50 60 70 80 90 100
Cold Season Warm Season RESULTS 75th Percentile Observed Estimated AN NCEP R2-70 -65-60 -55-70 -65-60 -55-70 -65-60 -55 mm/day 24 22 20 18 16 14 12 10 8 18 16 14 12 10 8-70 -65-60 -55-70 -65-60 -55-70 -65-60 -55 6
Cold Season Warm Season RESULTS 75th Percentile -70-65 -60-55 Observed Estimated AN NCEP R2 Good agreement in the spatial -70-65 -60-55 -70-65 -60-55 -70-65 -60-55 distribution and values of the 75th underestimates the percentile in both seasons. 75th Percentile -70-65 -60-55 NCEP R2-70 -65-60 -55 mm/day 24 22 20 18 16 14 12 10 8 18 16 14 12 10 8 6
Concluding Remarks The ESD performance depends on the season and station: the lowest skill was found in Summer related to small scale processes that lead to precipitation. However, the results show the great potential of the method which is able to reproduce different aspects of daily precipitation with a very good level of accuracy.
Concluding Remarks These results encurage us to: continue exploring the advantages of ESD to produce actionable regional climate information evaluate different ESD methods apart from the analogue method
Thanks!