Climate change impacts in the design of drainage systems: case study of Portugal

This paper aims to assess the necessity of updating the intensity-duration-frequency (IDF) curves used in Portugal to design building storm-water drainage systems. A comparative analysis of the design was performed for the three predefined rainfall regions in Portugal using the IDF curves currently in use and estimated for future decades. Data for recent and future climate conditions simulated by a global and regional climate model chain are used to estimate possible changes of rainfall extremes and its implications for the drainage systems. The methodology includes the disaggregation of precipitation up to subhourly scales, the robust development of IDF curves, and the correction of model bias. Obtained results indicate that projected changes are largest for the plains in southern Portugal (5–33%) than for mountainous regions (3–9%) and that these trends are consistent with projected changes in the long-term 95th percentile of the daily precipitation throughout the 21st century. The authors conclude there is a need to review the current precipitation regime classification and change the new drainage systems towards larger dimensions to mitigate the projected changes in extreme precipitation.

The XWS open access catalogue of extreme European windstorms from 1979 to 2012

The XWS (eXtreme WindStorms) catalogue consists of storm tracks and model-generated maximum 3 s wind-gust footprints for 50 of the most extreme winter windstorms to hit Europe in the period 1979–2012. The catalogue is intended to be a valuable resource for both academia and industries such as (re)insurance, for example allowing users to characterise extreme European storms, and validate climate and catastrophe models. Several storm severity indices were investigated to find which could best represent a list of known high-loss (severe) storms. The best-performing index was Sft, which is a combination of storm area calculated from the storm footprint and maximum 925 hPa wind speed from the storm track. All the listed severe storms are included in the catalogue, and the remaining ones were selected using Sft. A comparison of the model footprint to station observations revealed that storms were generally well represented, although for some storms the highest gusts were underestimated. Possible reasons for this underestimation include the model failing to simulate strong enough pressure gradients and not representing convective gusts. A new recalibration method was developed to estimate the true distribution of gusts at each grid point and correct for this underestimation. The recalibration model allows for storm-to-storm variation which is essential given that different storms have different degrees of model bias. The catalogue is available at www.europeanwindstorms.org.

Secondary cyclogenesis along an occluded front leading to damaging wind gusts: windstorm Kyrill, January 2007

Windstorm Kyrill affected large parts of Europe in January 2007 and caused widespread havoc and loss of life. In this study the formation of a secondary cyclone, Kyill II, along the occluded front of the mature cyclone Kyrill and the occurrence of severe wind gusts as Kyrill II passed over Germany are investigated with the help of high-resolution regional climate model simulations. Kyrill underwent an explosive cyclogenesis south of Greenland as the storm crossed polewards of an intense upper-level jet stream. Later in its life cycle secondary cyclogenesis occurred just west of the British Isles. The formation of Kyrill II along the occluded front was associated (a) with frontolytic strain and (b) with strong diabatic heating in combination with a developing upper-level shortwave trough. Sensitivity studies with reduced latent heat release feature a similar development but a weaker secondary cyclone, revealing the importance of diabatic processes during the formation of Kyrill II. Kyrill II moved further towards Europe and its development was favored by a split jet structure aloft, which maintained the cyclone’s exceptionally deep core pressure (below 965 hPa) for at least 36 hours. The occurrence of hurricane force winds related to the strong cold front over North and Central Germany is analyzed using convection-permitting simulations. The lower troposphere exhibits conditional instability, a turbulent flow and evaporative cooling. Simulation at high spatio-temporal resolution suggests that the downward mixing of high momentum (the wind speed at 875 hPa widely exceeded 45 m s-1) accounts for widespread severe surface wind gusts, which is in agreement with observed widespread losses.

Ionospheric and geomagnetic responses to changes in IMF B_Z: a superposed epoch study

Superposed epoch studies have been carried out in order to determine the ionospheric response at mid-latitudes to southward turnings of the interplanetary magnetic field (IMF). This is compared with the geomagnetic response, as seen in the indices K p, AE and Dst. The solar wind, IMF and geomagnetic data used were hourly averages from the years 1967–1989 and thus cover a full 22-year cycle in the solar magnetic field. These data were divided into subsets, determined by the magnitudes of the southward turnings and the concomitant increase in solar wind pressure. The superposed epoch studies were carried out using the time of the southward turning as time zero. The response of the mid-latitude ionosphere is studied by looking at the F-layer critical frequencies, f o F2, from hourly soundings by the Slough ionosonde and their deviation from the monthly median values, δf o F2. For the southward turnings with a change in B z of δB z > 11.5 nT accompanied by a solar wind dynamic pressure P exceeding 5 nPa, the F region critical frequency, f o F2, shows a marked decrease, reaching a minimum value about 20 h after the southward turning. This recovers to pre-event values over the subsequent 24 h, on average. The Dst index shows the classic storm-time decrease to about −60 nT. Four days later, the index has still to fully recover and is at about −25 nT. Both the K p and AE indices show rises before the southward turnings, when the IMF is strongly northward but the solar wind dynamic pressure is enhanced. The average AE index does register a clear isolated pulse (averaging 650 nT for 2 h, compared with a background peak level of near 450 nT at these times) showing enhanced energy deposition at high latitudes in substorms but, like K p, remains somewhat enhanced for several days, even after the average IMF has returned to zero after 1 day. This AE background decays away over several days as the Dst index recovers, indicating that there is some contamination of the currents observed at the AE stations by the continuing enhanced equatorial ring current. For data averaged over all seasons, the critical frequencies are depressed at Slough by 1.3 MHz, which is close to the lower decile of the overall distribution of δf o Fl values. Taking 30-day periods around summer and winter solstice, the largest depression is 1.6 and 1.2 MHz, respectively. This seasonal dependence is confirmed by a similar study for a Southern Hemisphere station, Argentine Island, giving peak depressions of 1.8 MHz and 0.5 MHz for summer and winter. For the subset of turnings where δB z > 11.5 nT and P ≤ 5 nPa, the response of the geomagnetic indices is similar but smaller, while the change in δf o F2 has all but disappeared. This confirms that the energy deposited at high latitudes, which leads to the geomagnetic and ionospheric disturbances following a southward turning of the IMF, increases with the energy density (dynamic pressure) of the solar wind flow. The magnitude of all responses are shown to depend on δB z . At Slough, the peak depression always occurs when Slough rotates into the noon sector. The largest ionospheric response is for southward turnings seen between 15–21 UT.

Field-aligned plasma flow in the quiet, mid-latitude ionosphere deduced from topside soundings

A method for quantifying diffusive flows of O+ ions in the topside ionosphere from satellite soundings is described. A departure from diffusive equilibrium alters the shape of the plasma scale-height profile near the F2-peak where ion-neutral frictional drag is large. The effect enables the evaluation of , the field-aligned flux of O+ ions relative to the neutral oxygen atom gas, using MSIS model values for the neutral thermospheric densities and temperature. Upward flow values are accurate to within about 10%, the largest sources of error being the MSIS prediction for the concentration of oxygen atoms and the plasma temperature gradient deduced from the sounding. Downward flux values are only determined to within 20%. From 60,000 topside soundings, taken at the minimum and rising phase of the solar cycle, a total of 1098 mean scale-height profiles are identified for which no storm sudden commencement had occurred in the previous 12 days and for which Kp was less than 2o, each mean profile being an average of about six soundings. A statistical study ofdeduced from these profiles shows the diurnal cycle of O+ flow in the quiet, topside ionosphere at mid-latitudes and its seasonal variations. The differences betweenand ion flux observations from incoherent scatter radars are considered using the meridional thermospheric winds predicted by a global, three-dimensional model. The mean interhemispheric flow from summer to winter is compared with predictions by a numerical model of the protonospheric coupling of conjugate ionospheres for up to 6 days following a geomagnetic storm. The observed mean (of order 3 × 1016 ions day−1 along a flux tube of area 1 m2 at 1000 km) is larger than predicted for day 6 and the suggested explanation is a decrease in upward flows from the winter, daytime ionosphere between the sixth and twelfth days.

Characterising phosphorus and nitrate inputs to a rural river using high-frequency concentration-flow relationships

The total reactive phosphorus (TRP) and nitrate concentrations of the River Enborne, southern England, were monitored at hourly interval between January 2010 and December 2011. The relationships between these high-frequency nutrient concentration signals and flow were used to infer changes in nutrient source and dynamics through the annual cycle and each individual storm event, by studying hysteresis patterns. TRP concentrations exhibited strong dilution patterns with increasing flow, and predominantly clockwise hysteresis through storm events. Despite the Enborne catchment being relatively rural for southern England, TRP inputs were dominated by constant, non-rain-related inputs from sewage treatment works (STW) for the majority of the year, producing the highest phosphorus concentrations through the spring–summer growing season. At higher river flows, the majority of the TRP load was derived from within-channel remobilisation of phosphorus from the bed sediment, much of which was also derived from STW inputs. Therefore, future phosphorus mitigation measures should focus on STW improvements. Agricultural diffuse TRP inputs were only evident during storms in the May of each year, probably relating to manure application to land. The nitrate concentration–flow relationship produced a series of dilution curves, indicating major inputs from groundwater and to a lesser extent STW. Significant diffuse agricultural inputs with anticlockwise hysteresis trajectories were observed during the first major storms of the winter period. The simultaneous investigation of high-frequency time series data, concentration–flow relationships and hysteresis behaviour through multiple storms for both phosphorus and nitrate offers a simple and innovative approach for providing new insights into nutrient sources and dynamics.

Return periods of losses associated with European windstorm series in a changing climate

Possible future changes of clustering and return periods (RPs) of European storm series with high potential losses are quantified. Historical storm series are identified using 40 winters of reanalysis. Time series of top events (1, 2 or 5 year return levels (RLs)) are used to assess RPs of storm series both empirically and theoretically. Additionally, 800 winters of general circulation model simulations for present (1960–2000) and future (2060–2100) climate conditions are investigated. Clustering is identified for most countries, and estimated RPs are similar for reanalysis and present day simulations. Future changes of RPs are estimated for fixed RLs and fixed loss index thresholds. For the former, shorter RPs are found for Western Europe, but changes are small and spatially heterogeneous. For the latter, which combines the effects of clustering and event ranking shifts, shorter RPs are found everywhere except for Mediterranean countries. These changes are generally not statistically significant between recent and future climate. However, the RPs for the fixed loss index approach are mostly beyond the range of pre-industrial natural climate variability. This is not true for fixed RLs. The quantification of losses associated with storm series permits a more adequate windstorm risk assessment in a changing climate.

Large-scale dynamics associated with clustering of extratropical cyclones affecting Western Europe

Some recent winters in Western Europe have been characterized by the occurrence of multiple extratropical cyclones following a similar path. The occurrence of such cyclone clusters leads to large socio-economic impacts due to damaging winds, storm surges, and floods. Recent studies have statistically characterized the clustering of extratropical cyclones over the North Atlantic and Europe and hypothesized potential physical mechanisms responsible for their formation. Here we analyze 4 months characterized by multiple cyclones over Western Europe (February 1990, January 1993, December 1999, and January 2007). The evolution of the eddy driven jet stream, Rossby wave-breaking, and upstream/downstream cyclone development are investigated to infer the role of the large-scale flow and to determine if clustered cyclones are related to each other. Results suggest that optimal conditions for the occurrence of cyclone clusters are provided by a recurrent extension of an intensified eddy driven jet toward Western Europe lasting at least 1 week. Multiple Rossby wave-breaking occurrences on both the poleward and equatorward flanks of the jet contribute to the development of these anomalous large-scale conditions. The analysis of the daily weather charts reveals that upstream cyclone development (secondary cyclogenesis, where new cyclones are generated on the trailing fronts of mature cyclones) is strongly related to cyclone clustering, with multiple cyclones developing on a single jet streak. The present analysis permits a deeper understanding of the physical reasons leading to the occurrence of cyclone families over the North Atlantic, enabling a better estimation of the associated cumulative risk over Europe.

Climate change projections and selected impacts for Vietnam

In September 2013, the 5th Assessment Report (5AR) of the International Panel on Climate Change (IPCC) has been released. Taking the 5AR cli-mate change scenarios into account, the World Bank published an earli-er report on climate change and its impacts on selected hot spot re-gions, including Southeast Asia. Currently, dynamical and statistical-dynamical downscaling efforts are underway to obtain higher resolution and more robust regional climate change projections for tropical South-east Asia, including Vietnam. Such initiatives are formalized under the World Meteorological Organization (WMO) Coordinated Regional Dynamic Downscaling Experiment (CORDEX) East Asia and Southeast Asia and also take place in climate change impact projects such as the joint Vietnam-ese-German project “Environmental and Water Protection Technologies of Coastal Zones in Vietnam (EWATEC-COAST)”. In this contribution, the lat-est assessments for changes in temperature, precipitation, sea level, and tropical cyclones (TCs) under the 5AR Representative Concentration Pathway (RCP) scenarios 4.5 and 8.5 are reviewed. Special emphasis is put on changes in extreme events like heat waves and/or heavy precipita-tion. A regional focus is Vietnam south of 16°N. A continued increase in mean near surface temperature is projected, reaching up to 5°C at the end of this century in northern Vietnam un-der the high greenhouse-gas forcing scenario RCP8.5. Overall, project-ed changes in annual precipitation are small, but there is a tendency of more rainfall in the boreal winter dry season. Unprecedented heat waves and an increase in extreme precipitation events are projected by both global and regional climate models. Globally, TCs are projected to decrease in number, but an increase in intensity of peak winds and rain-fall in the inner core region is estimated. Though an assessment of changes in land-falling frequency in Vietnam is uncertain due to difficul-ties in assessing changes in TC tracks, some work indicates a reduction in the number of land-falling TCs in Vietnam. Sea level may rise by 75-100 cm until the end of the century with the Vietnamese coastline experienc-ing 10-15% higher rise than on global average. Given the large rice and aquaculture production in the Mekong and Red River Deltas, that are both prone to TC-related storm surges and flooding, this poses a challenge to foodsecurity and protection of coastal population and assets.

The DYMECS project: a statistical approach for the evaluation of convective storms in high-resolution NWP models

A new frontier in weather forecasting is emerging by operational forecast models now being run at convection-permitting resolutions at many national weather services. However, this is not a panacea; significant systematic errors remain in the character of convective storms and rainfall distributions. The DYMECS project (Dynamical and Microphysical Evolution of Convective Storms) is taking a fundamentally new approach to evaluate and improve such models: rather than relying on a limited number of cases, which may not be representative, we have gathered a large database of 3D storm structures on 40 convective days using the Chilbolton radar in southern England. We have related these structures to storm life-cycles derived by tracking features in the rainfall from the UK radar network, and compared them statistically to storm structures in the Met Office model, which we ran at horizontal grid length between 1.5 km and 100 m, including simulations with different subgrid mixing length. We also evaluated the scale and intensity of convective updrafts using a new radar technique. We find that the horizontal size of simulated convective storms and the updrafts within them is much too large at 1.5-km resolution, such that the convective mass flux of individual updrafts can be too large by an order of magnitude. The scale of precipitation cores and updrafts decreases steadily with decreasing grid lengths, as does the typical storm lifetime. The 200-m grid-length simulation with standard mixing length performs best over all diagnostics, although a greater mixing length improves the representation of deep convective storms.

Rainfall estimation by rain gauge-radar combination: a concurrent multiplicative-additive approach

A procedure (concurrent multiplicative-additive objective analysis scheme [CMA-OAS]) is proposed for operational rainfall estimation using rain gauges and radar data. On the basis of a concurrent multiplicative-additive (CMA) decomposition of the spatially nonuniform radar bias, within-storm variability of rainfall and fractional coverage of rainfall are taken into account. Thus both spatially nonuniform radar bias, given that rainfall is detected, and bias in radar detection of rainfall are handled. The interpolation procedure of CMA-OAS is built on Barnes' objective analysis scheme (OAS), whose purpose is to estimate a filtered spatial field of the variable of interest through a successive correction of residuals resulting from a Gaussian kernel smoother applied on spatial samples. The CMA-OAS, first, poses an optimization problem at each gauge-radar support point to obtain both a local multiplicative-additive radar bias decomposition and a regionalization parameter. Second, local biases and regionalization parameters are integrated into an OAS to estimate the multisensor rainfall at the ground level. The procedure is suited to relatively sparse rain gauge networks. To show the procedure, six storms are analyzed at hourly steps over 10,663 km2. Results generally indicated an improved quality with respect to other methods evaluated: a standard mean-field bias adjustment, a spatially variable adjustment with multiplicative factors, and ordinary cokriging.

An assessment of the impact of local processes on dust lifting in martian climate models

Simulation of the lifting of dust from the planetary surface is of substantially greater importance on Mars than on Earth, due to the fundamental role that atmospheric dust plays in the former’s climate, yet the dust emission parameterisations used to date in martian global climate models (MGCMs) lag, understandably, behind their terrestrial counterparts in terms of sophistication. Recent developments in estimating surface roughness length over all martian terrains and in modelling atmospheric circulations at regional to local scales (less than O(100 km)) presents an opportunity to formulate an improved wind stress lifting parameterisation. We have upgraded the conventional scheme by including the spatially varying roughness length in the lifting parameterisation in a fully consistent manner (thereby correcting a possible underestimation of the true threshold level for wind stress lifting), and used a modification to account for deviations from neutral stability in the surface layer. Following these improvements, it is found that wind speeds at typical MGCM resolution never reach the lifting threshold at most gridpoints: winds fall particularly short in the southern midlatitudes, where mean roughness is large. Sub-grid scale variability, manifested in both the near-surface wind field and the surface roughness, is then considered, and is found to be a crucial means of bridging the gap between model winds and thresholds. Both forms of small-scale variability contribute to the formation of dust emission ‘hotspots’: areas within the model gridbox with particularly favourable conditions for lifting, namely a smooth surface combined with strong near-surface gusts. Such small-scale emission could in fact be particularly influential on Mars, due both to the intense positive radiative feedbacks that can drive storm growth and a strong hysteresis effect on saltation. By modelling this variability, dust lifting is predicted at the locations at which dust storms are frequently observed, including the flushing storm sources of Chryse and Utopia, and southern midlatitude areas from which larger storms tend to initiate, such as Hellas and Solis Planum. The seasonal cycle of emission, which includes a double-peaked structure in northern autumn and winter, also appears realistic. Significant increases to lifting rates are produced for any sensible choices of parameters controlling the sub-grid distributions used, but results are sensitive to the smallest scale of variability considered, which high-resolution modelling suggests should be O(1 km) or less. Use of such models in future will permit the use of a diagnosed (rather than prescribed) variable gustiness intensity, which should further enhance dust lifting in the southern hemisphere in particular.

The slowly evolving background state of the atmosphere

The theory of wave–mean flow interaction requires a partition of the atmospheric flow into a notional background state and perturbations to it. Here, a background state, known as the Modified Lagrangian Mean (MLM), is defined as the zonally symmetric state obtained by requiring that every potential vorticity (PV) contour lying within an isentropic layer encloses the same mass and circulation as in the full flow. For adiabatic and frictionless flow, these two integral properties are time-invariant and the MLM state is a steady solution of the primitive equations. The time dependence in the adiabatic flow is put into the perturbations, which can be described by a wave-activity conservation law that is exact even at large amplitude. Furthermore, the effects of non-conservative processes on wave activity can be calculated from the conservation law. A new method to calculate the MLM state is introduced, where the position of the lower boundary is obtained as part of the solution. The results are illustrated using Northern Hemisphere ERA-Interim data. The MLM state evolves slowly, implying that the net non-conservative effects are weak. Although ‘adiabatic eddy fluxes’ cannot affect the MLM state, the effects of Rossby-wave breaking, PV filamentation and subsequent dissipation result in sharpening of the polar vortex edge and meridional shifts in the MLM zonal flow, both at tropopause level and on the winter stratospheric vortex. The rate of downward migration of wave activity during stratospheric sudden warmings is shown to be given by the vertical scale associated with polar vortex tilt divided by the time-scale for wave dissipation estimated from the wave-activity conservation law. Aspects of troposphere–stratosphere interaction are discussed. The new framework is suitable to examine the climate and its interactions with disturbances, such as midlatitude storm tracks, and makes a clean partition between adiabatic and non-conservative processes.

Clouds, circulation and climate sensitivity

Fundamental puzzles of climate science remain unsolved because of our limited understanding of how clouds, circulation and climate interact. One example is our inability to provide robust assessments of future global and regional climate changes. However, ongoing advances in our capacity to observe, simulate and conceptualize the climate system now make it possible to fill gaps in our knowledge. We argue that progress can be accelerated by focusing research on a handful of important scientific questions that have become tractable as a result of recent advances. We propose four such questions below; they involve understanding the role of cloud feedbacks and convective organization in climate, and the factors that control the position, the strength and the variability of the tropical rain belts and the extratropical storm tracks.

Sea level extremes at the coasts of China

Hourly sea level records from 1954 to 2012 at 20 tide gauges at and adjacent to the Chinese coasts are used to analyze extremes in sea level and in tidal residual. Tides and tropical cyclones determine the spatial distribution of sea level maxima. Tidal residual maxima are predominantly determined by tropical cyclones. The 50 year return level is found to be sensitive to the number of extreme events used in the estimation. This is caused by the small number of tropical cyclone events happening each year which lead to other local storm events included thus significantly affecting the estimates. Significant increase in sea level extremes is found with trends in the range between 2.0 and 14.1 mm yr−1. The trends are primarily driven by changes in median sea level but also linked with increases in tidal amplitudes at three stations. Tropical cyclones cause significant interannual variations in the extremes. The interannual variability in the sea level extremes is also influenced by the changes in median sea level at the north and by the 18.6 year nodal cycle at the South China Sea. Neither of PDO and ENSO is found to be an indicator of changes in the size of extremes, but ENSO appears to regulate the number of tropical cyclones that reach the Chinese coasts. Global mean atmospheric temperature appears to be a good descriptor of the interannual variability of tidal residual extremes induced by tropical cyclones but the trend in global temperature is inconsistent with the lack of trend in the residuals.