The methodology outlined in Chapter 3 has been evaluated in two case studies within the SUDPLAN project; one in Linz (Austria) focusing on continuous time series and one in Wuppertal (Germany) focusing on extreme events. In both cases, short-term precipitation data are extracted from climate projections and analyzed with respect to key statistics and DCF relationships.
4.2. Linz: Continuous Time Series
Linz is the capital of the Austrian federal state Upper Austria and is situated in the northern part of Austria directly at the Danube River. The city of Linz and its agglomeration have approximately 270,000 inhabitants. The Danube is the biggest Austrian river with a mean runoff of about 2000 m³/s at Linz. The average annual rainfall is 830 mm. For the urban drainage pilot in the SUDPLAN project, the whole of Linz and its 39 neighboring communes are of interest as they are drained together to one central wastewater treatment plant (WWTP) by a combined sewer system. In total this area covers approximately 900 km².
In combined sewer systems domestic and industrial sewage as well as storm water is transported together in the same pipe. In case of heavy rainfall, the combined sewage is in general either stored temporarily in the system or spilled out by combined sewer overflows in order to reduce the hydraulic pressure on the system and to limit the discharge to the WWTP. These overflows can have an important ecological impact on the receiving waters and are hygienically problematic. The sewer operator at Linz has put considerable effort and investments in storage facilities and control strategies in order to reduce the overflows and currently meets the Austrian requirements.
As rainfall is the most important driver of the runoff and overflow processes in a combined sewer system, information on possible future rainfall is crucial in order to develop proper mitigation strategies. This is also the major aim for the Linz Pilot in the SUDPLAN project: the SUDPLAN SMS will allow the operator to easily compare the effect of different future scenarios e.g., due to predicted rainfall by means of a simulation model. Details on the model are given in [
24]. One local high resolution rainfall time series for Downtown Linz was available from the Austrian NIEDA tool [
25]. Data is recorded by a tipping bucket rain gauge and is prepared especially for use in urban drainage modeling. Data are available fromthe period 1993–2006 and the 30-year reference period used in this case is 1981–2010.
In the observed time series, both monthly total precipitation and maximum 30-min intensity, here defined as the 99th percentile of the frequency distribution of non-zero values, is highest in summer (
Table 1). In winter, the total precipitation is about half of that in summer and the maximum intensity 6–7 times lower. The characteristics for spring and autumn are rather similar, both being in between the values found for summer and winter. The wet fraction, a measure of occurrence defined as the percentage of 30-min periods with a non-zero precipitation, is consistently close to 10%, slightly lower in summer and slightly higher in winter.
The ECH projection indicates a pronounced relative increase (~25%) of the total autumn and winter precipitation, a smaller increase in spring (14%) and a clear decrease in summer (−14%) (
Table 1). The annual cycle becomes flatter with all months approaching totals of 70–80 mm. The maximum intensity changes in all seasons, somewhat more in autumn and winter (~25–35%) than in spring and summer (~15–20%) but the shape of the annual cycle is largely preserved. The wet fraction shows a minor increase (~5–10%) in all seasons expect summer where it decreases by 16%. In summary, all considered aspects of precipitation (total, maximum, occurrence) increase in autumn, winter and spring. In summer only the maximum intensity increases; the total and the occurrence decreases. This illustrates the need for a flexible Delta Change approach that can manage e.g., contrasting trends for precipitation totals and maxima, respectively, as well as both positive and negative changes in the frequency of occurrence.
Table 1.
Total precipitation (mm/month), maximum intensity (mm/30 min) and wet fraction (%) for all seasons in Linz in observations (OBS) and future projections ECH and HAD. For the projections, the percentage change is added in parentheses.
Table 1.
Total precipitation (mm/month), maximum intensity (mm/30 min) and wet fraction (%) for all seasons in Linz in observations (OBS) and future projections ECH and HAD. For the projections, the percentage change is added in parentheses.
Season | Variable | OBS | ECH | HAD |
---|
Spring | Total | 66.0 | 75.1 (+ 14) | 72.4 (+ 10) |
Maximum | 15.3 | 17.8 (+16) | 17.3 (+ 13) |
Wet fraction | 9.7 | 10.4 (+ 7) | 10.3 (+ 6) |
Summer | Total | 99.2 | 85.3 (− 14) | 95.3 (− 4) |
Maximum | 26.8 | 32.2 (+ 20) | 31.7 (+ 18) |
Wet fraction | 8.9 | 7.5 (− 16) | 8.4 (− 5) |
Autumn | Total | 58.3 | 71.2 (+ 22) | 59.9 (+ 3) |
Maximum | 9.4 | 12.8 (+ 36) | 11.0 (+ 17) |
Wet fraction | 8.6 | 8.9 (+ 4) | 7.9 (− 8) |
Winter | Total | 55.2 | 69.8 (+ 26) | 66.2 (+ 20) |
Maximum | 4.2 | 5.3 (+ 27) | 5.6 (+ 33) |
Wet fraction | 10.9 | 11.8 (+ 9) | 11.3 (+ 3) |
The future changes indicated in the HAD projection are similar overall to the changes in the ECH projection (
Table 1). The sign of the change differs in only one case; the HAD projection indicates a slight decrease of the wet fraction in autumn, in contrast with the minor increase found in the ECH projection. Other notable differences concern the autumn total and maximum, which both increase much less in HAD than in ECH. Further, the decrease of summer total and occurrence are less pronounced in HAD than in ECH. In winter and spring all aspects change similarly in the two projections.
Concerning the DCF functions, these may be divided into two parts with different characteristics (
Figure 2 and
Figure 3). Below the 70th–80th percentile, the DCF functions fluctuate around a relatively constant average value that changes with the season. The average value is closely related to the change in precipitation totals (
Table 1). In winter, with the maximum increase of total precipitation, the average DCF for low and intermediate percentiles is somewhere between 1.1 and 1.2. In summer, when total precipitation decreases in the projections, the average DCF is below 1. The constant DCF value approached for low and intermediate percentiles thus reflects the change found for the majority of intensities, which is manifested in the total precipitation. The notable fluctuations in the DCF functions, which are most pronounced for the lowest percentiles, are caused by the intensity resolution of the RCA3 precipitation data used here (0.01 mm/30 min). Low percentiles are close to this value, which creates sudden jumps in the DCF functions when the percentiles change with 0.01 up or down.
The spread between the different five grid boxes are dominated by the fluctuations related to the intensity resolution and therefore difficult to judge visually, but overall the spread around the mean DCF function (thick black line in
Figure 2 and
Figure 3) is limited and similar for all seasons.
Above the 70th–80th percentile, the DCF functions are smoother and generally increasing towards their highest value at the 100th percentile. This implies that the highest future relative increase is expected for the highest short-term intensities. The general changes in maximum precipitation (
Table 1) are reflected in the DCF value at the 100th percentile, being generally 1.3–1.4 in autumn and winter (
Figure 2c,d and
Figure 3c,d) and 1.1–1.2 in spring and summer (
Figure 2a,b and
Figure 3a,b). It should be emphasised that local precipitation extremes, relevant for, e.g., sewer system design, are generally generated by convective storms in summer. Thus the summer DCF values reflect the expected impact on design precipitation, even if these are lower than DCF values for other seasons. Some sudden bends are sometimes found in the DCF functions at the highest percentiles and these are caused by statistical variability of extreme intensities related to the most intense events. The smoothness of the functions reflects the fact that the intensity resolution has no impact on the DCF values for higher intensities. The spread between grid boxes for large percentiles varies from quite pronounced (e.g.,
Figure 3c) to very limited (e.g.,
Figure 3a) but overall the mean function is considered to represent the general patterns sufficiently well.
Comparing the two projections ECH (
Figure 2) and HAD (
Figure 3) indicates overall similar future changes of not only general aspects of precipitation (
Table 1) but also in terms of the of the DCF functions. An obvious possibility is to parameterize the functions. This has been done in local applications of the approach, using simple expressions. However, it has turned out to be difficult to find a general parameterisation that is guaranteed to work satisfactorily for all attainable shapes of the DCF function and which generally, but not always, looks like the ones in
Figure 2 and
Figure 3. Further, careful weighing is required as the highest percentiles are generally the most important and must be well matched by the fit, but also the remaining part of the function must be satisfactorily reproduced. For these reasons, at least initially in the SUDPLAN tool, no parameterization is done but the empirical DCFs for each percentile are used in the re-scaling.
An event-based analysis of the two projections was made in order to evaluate the assumption behind the proposed approach to frequency adjustment (
Section 3.1). A limit of 2 h was used to separate independent events and events with a total volume below 1 mm were not included. Four types of changes from the reference period to the future period were considered.
ΔWF: relative change in the wet fraction, defined as above.
ΔNE: relative change in the number of events.
ΔEV: relative change in the average event volume.
ΔED: relative change in the average event duration.
Figure 2.
Delta Change (DC) factor as a function of percentile for all seasons in the ECH projection. Thin grey lines represent individual RCA3 grid boxes; thick black line is the average function.
Figure 2.
Delta Change (DC) factor as a function of percentile for all seasons in the ECH projection. Thin grey lines represent individual RCA3 grid boxes; thick black line is the average function.
Figure 3.
Delta Change (DC) factor as a function of percentile for all seasons in the HAD projection. Thin grey lines represent individual RCA3 grid boxes; thick black line is the average function.
Figure 3.
Delta Change (DC) factor as a function of percentile for all seasons in the HAD projection. Thin grey lines represent individual RCA3 grid boxes; thick black line is the average function.
The changes of wet fraction ΔWF are overall similar to the ones reported in
Table 1 but not identical as small events were discarded in the event analysis (
Table 2). The main changes were a decrease in summer and an increase in winter. Concerning summer, in projection ECH the changes in average volume ΔEV and duration ΔED are very small but the projected decrease in wet fraction (−1%) is essentially only caused by a reduced number of events (−22%). The situation is less clear in projection HAD, which has notably fewer future events (−8%) but also a clear decrease of the mean duration (−12%). In both projections, the clear increase of the wet fraction in winter is mainly related to an increased number of events. In spring and autumn only smaller changes in the wet fraction are indicated, which are clearly related to changes in the number of events but also some changes in the event volume are indicated. In total, frequency adjustment by a changed number of events rather than changed event characteristics appears justified.
Table 2.
Seasonal relative changes (%) in the wet fraction (ΔWF), number of events (ΔNE), average volume (ΔEV) and average duration (ΔED) between the reference and the future period in projections ECH and HAD.
Table 2.
Seasonal relative changes (%) in the wet fraction (ΔWF), number of events (ΔNE), average volume (ΔEV) and average duration (ΔED) between the reference and the future period in projections ECH and HAD.
Season | ΔWF | ΔNE | ΔEV | ΔED |
---|
ECH | HAD | ECH | HAD | ECH | HAD | ECH | HAD |
---|
Spring | +9 | +4 | +10 | +4 | +3 | +7 | −2 | +1 |
Summer | −21 | −15 | −22 | −8 | +3 | 0 | −6 | −12 |
Autumn | +8 | −3 | +7 | −1 | +21 | +9 | +2 | −1 |
Winter | +22 | +19 | +23 | +29 | +11 | +4 | +2 | −5 |
4.3. Wuppertal: Extreme Events
The city of Wuppertal, a town with approximately 350,000 residents, is situated in a hill country (~100–350 m.a.s.l.) in the state of North-Rhine-Westphalia in Germany. It is located in the valley of the river Wupper (
Figure 4). There are several creeks on both sides of this valley that open into the storm water system before they finally end in the main receiving river Wupper. During heavy rainfall events the city’s storm water system is quickly blocked by those swollen creeks causing the precipitation to run off on the surface. The storm water runoff may thereby affect valuable public infrastructure and private property, as observed in the recent years in some of the creek catchment areas. Due to the complex geography it is completely unpredictable where a heavy rainfall event might occur and unknown whether there will be flooding.
Up to now the mid- and long-term planning of the storm water sewage system in Wuppertal has been accomplished with iterative model runs of a hydrological model (for the natural creeks) and a hydrodynamic model (for the canalized creeks and sewage system). Surface runoff was not taken into account in detailed simulations but by expert judgement. In the framework of the SUDPLAN-project the newly combined 1D-2D-hydrodynamical model DYNA-GeoCPM (
www.tandler.com and Pecher Software) is used to detect critical spots with a high risk of flooding, e.g., regarding valuable and vulnerable facilities [
26].
The aim is to mitigate the risk of flooding during serious storm events (return period T = 30 years or more) for the detected critical spots. The traditional strategies and measures to achieve this are either the enlargement of the storm water system profiles or the construction of retention basins. Given these two options the potential needs for investments would be immense, considering that the city faces inflow from 350 km of creeks and 650 km of storm water sewer systems. An alternative and much more cost-efficient strategy is to look for localized planning options which are likely to prevent damage. Examples of such structural measures are the alteration of street profiles by means of higher road curbs or the installation of mobile or stationary walls.
Figure 4.
The city limit of Wuppertal with the test catchment Luentenbeck (red area) and the KOSTRA grid (yellow grid).
Figure 4.
The city limit of Wuppertal with the test catchment Luentenbeck (red area) and the KOSTRA grid (yellow grid).
The SUDPLAN-platform will provide the responsible planners and hydraulic modelers in Wuppertal with a tool that enables them to simulate a multitude of scenarios with the model components for the sewer system and surface runoff, both to detect the critical spots and to simulate the effects of different structural measures at the critical spots. The tool should be able to store the parameters and results of such a model run and to visualize the results regarding climate change effects.
EULER II design storms for the Wuppertal pilot were calculated according to DWA-A 118 [
2] based on the IDF curves described in [
27]. The KOSTRA base data for calculating and regionalizing the precipitation amounts for duration periods between 24 h and 72 h consists of the daily 1 km gridded precipitation totals (1951–2000) using the REGNIE method (regionalization of precipitation totals in ~350,000 grid boxes). Short duration periods of the KOSTRA IDF curves are based on a collective of about 200 precipitation stations (1951–2000). The results from the actual analyses are being converted into the habitual KOSTRA 8.45 × 8.45 km grid with ~5350 grid boxes. It may be remarked that the area covered by the six grid boxes in
Figure 4, ~430 km², corresponds to roughly 1/6 of an RCM grid box with the 50 × 50 km resolution used here. The test catchment area “Luentenbeck” for the Wuppertal pilot is located in the KOSTRA grid box 11/51 (
Figure 4) and the corresponding 30-year intensities for a range of durations are given in
Table 3. The reference period used in this case study is 1961–1990.
The Gumbel distribution fits to annual maxima for different durations calculated from time series from individual RCA3 grid boxes are generally very accurate (
Figure 5). Some weak oscillations can be noted for the extracted maxima but in total the log-linear Gumbel fits are very good approximations.
Table 3.
The 30-year precipitation depth (mm) in Wuppertal in KOSTRA data and future projections ECH and HAD. For the projections, the percentage change is added in parentheses.
Table 3.
The 30-year precipitation depth (mm) in Wuppertal in KOSTRA data and future projections ECH and HAD. For the projections, the percentage change is added in parentheses.
Duration | KOSTRA | ECH | HAD |
---|
10 min | 20.1 | 22.4 (+ 11) * | 22.9 (+ 14) * |
30 min | 31.7 | 35.3 (+ 11) | 36.1 (+ 14) |
1 h | 39.8 | 44.5 (+ 12) | 45.4 (+ 14) |
2 h | 46.8 | 52.8 (+ 13) | 53.6 (+ 15) |
6 h | 61.0 | 71.4 (+ 17) | 71.2 (+ 17) |
12 h | 72.2 | 89.1 (+ 23) | 86.5 (+ 20) |
Figure 5.
Example of fitted Gumbel distributions (lines) to annual maxima (points) for different durations calculated from 30-min precipitation in one RCA3 grid box.
Figure 5.
Example of fitted Gumbel distributions (lines) to annual maxima (points) for different durations calculated from 30-min precipitation in one RCA3 grid box.
Concerning the future change of the 30-year intensity in Wuppertal, the projection ECH indicates a clear dependence of the DCF function on the duration and towards larger changes for longer durations (
Figure 6a). In fact, no change at all is projected for the shortest duration (30 min), or even a decrease (although this is conceivably caused by statistical variability in the extreme value analysis rather than reflecting a true signal), whereas an increase exceeding 30% (
i.e., DCF > 1.3) is indicated for the longest durations considered (12–24 h). This type of duration dependence is in some contradiction with, e.g. a study of extreme short-term precipitation in Denmark from another RCM projection [
4], in which an increase of DCF with decreasing duration is indicated. We can only speculate that the qualitative difference found is due to either a different response to global warming in the Wuppertal region or to differences in the descriptions of the precipitation process in the RCMs used. The spread between grid boxes is in the range 0.2–0.3 units; in particular one grid box stands out with consistently ~0.1 higher DCF than the average line. The projection HAD (
Figure 6b) exhibits a principally similar dependence of the future change on duration, but the difference is smaller with DCFs ranging between 1.1 for 30 min to 1.25 for 24 h. The spread between grid boxes is smaller than for the ECH projection, being consistently in the range 0.1–0.2 units. In total, the spread is considered limited enough to meaningfully use the average change in subsequent downscaling. In the remaining part of this section, only changes averaged over all grid boxes are considered.
Figure 6.
Delta Change (DC) factor as a function of duration for the 30-year intensity in Wuppertal in projection ECH (a) and HAD (b). Thin grey lines represent individual RCA3 grid boxes; the thick black line is the average function.
Figure 6.
Delta Change (DC) factor as a function of duration for the 30-year intensity in Wuppertal in projection ECH (a) and HAD (b). Thin grey lines represent individual RCA3 grid boxes; the thick black line is the average function.
Concerning the dependence of the future change on frequency (
i.e., return period), for projection ECH there is a tendency towards a larger increase for long durations as the return period increases (
Figure 7a; point values). For the shortest duration the tendency is inversed towards a smaller future increase as the return period increases. For intermediate durations the dependence on return period is small, which makes the impact of return period small,also in the multiple linear regression used to obtain the final DCF functions (
Figure 7a; lines). The DCF functions corresponding to return periods 1 and 10 years, respectively, are virtually identical whereas the 100-year function is located slightly higher. From visual inspection only, the fit of the final DCF functions to the individual curves may not seem impressive, but this is mainly due to the pronounced variability for durations of 30 min and 12 h. Overall, the magnitude and dependencies of the DCF values appear reasonably well captured. More advanced regression may be considered to obtain improved fits, but likely at the expense of a reduced robustness.
Figure 7.
Delta Change (DC) factor as a function of duration for the T-year intensity in Wuppertal in projection ECH (a) and HAD (b). Point values represent original values; lines represent fits of Equation (1). The colors represent different return periods T = 1, 10 and 100 years.
Figure 7.
Delta Change (DC) factor as a function of duration for the T-year intensity in Wuppertal in projection ECH (a) and HAD (b). Point values represent original values; lines represent fits of Equation (1). The colors represent different return periods T = 1, 10 and 100 years.
For projection HAD, the dependence on return period is less clear, although the same tendency as in the ECH projection is found for short durations (
Figure 7b; point values). However, for the fitted DCF functions, contrary to the ECH projection, the 100-year function is located slightly below the other functions. Potentially, the dependence on the return period is unnecessary in the IDF curve downscaling, but this needs to be verified by further analyses. In the case of HAD, the fits are rather close to the individual curves for the entire range of durations. The weak dependence of DCF on the return period is in some contradiction with [
4], in which a clear tendency towards higher DCF with longer return period was indicated for Denmark.
The future estimated changes of the 30-year intensities obtained from the fitted DCF functions are listed in
Table 3. The value for a duration of 10 min has been extrapolated using the final functions. Overall, the two projections correspond well with respect to both the magnitude of change and the dependency on duration, as also seen in
Figure 7. The 1-h values may be compared with the European analysis performed in [
4], where an increase of ~20% was found for the same region between the same 30-year periods in another climate projection.
Finally, the 1-h EULER II design storm with a 30-year return period was rescaled in line with the projected future changes (
Figure 8). All the individual 5-min values have been multiplied by the 1-hour 30-year DCF obtained from ECH (1.12) and HAD (1.14). As discussed in
Section 3.2., in light of the weak duration dependence for short durations, this DCF value may also be assumed valid for shorter durations, down to the scale corresponding to the temporal resolution of the design storm (5 min).
Figure 8.
Historical 30-year 1-h EULER II design storm for Wuppertal (OBS) and downscaled versions based on future projection ECH and HAD.
Figure 8.
Historical 30-year 1-h EULER II design storm for Wuppertal (OBS) and downscaled versions based on future projection ECH and HAD.