Marie Ekström

Evaluating the performance of a WRF physics ensemble over South-East Australia

When using the Weather Research and Forecasting (WRF) modelling system it is necessary to choose between many parametrisations for each physics option. This study examines the performance of various physics scheme combinations on the simulation of a series of rainfall events near the south-east coast of Australia known as East Coast Lows. A thirty-six member multi-physics ensemble was created such that each member had a unique set of physics parametrisations. No single ensemble member was found to perform best for all events, variables and metrics. This is reflected in the fact that different climate variables are found to be sensitive to different physical parametrisations. While a standardised super-metric can be used to identify best performers, a step-wise decision approach described here, allows explicit recognition of the “robustness” of choosing one parameterisation over another, allowing the identification of a group of “equally robustly” performing physics combinations. These results suggest that the Mellor-Yamada-Janjic planetary boundary layer scheme and the Betts-Miller-Janjic cumulus scheme can be chosen with some robustness. Possibly with greater confidence, the results also suggest that the Yonsei University planetary boundary layer scheme, Kain-Fritsch cumulus scheme and RRTMG radiation scheme should not be used in combination in this region. Results further indicate that the selection of physics scheme options has larger impact on model performance during the more intensive rainfall events.

Optimally choosing small ensemble members to produce robust climate simulations

This study examines the subset climate model ensemble size required to reproduce certain statistical characteristics from a full ensemble. The ensemble characteristics examined are the root mean square error, the ensemble mean and standard deviation. Subset ensembles are created using measures that consider the simulation performance alone or include a measure of simulation independence relative to other ensemble members. It is found that the independence measure is able to identify smaller subset ensembles that retain the desired full ensemble characteristics than either of the performance based measures. It is suggested that model independence be considered when choosing ensemble subsets or creating new ensembles.

Evaluating rainfall patterns using physics scheme ensembles from a regional atmospheric model

This study evaluated the ability of Weather Research and Forecasting (WRF) multi-physics ensembles to simulate storm systems known as East Coast Lows (ECLs). ECLs are intense low-pressure systems that develop off the eastern coast of Australia. These systems can cause significant damage to the region. On the other hand, the systems are also beneficial as they generate the majority of high inflow to coastal reservoirs. It is the common interest of both hazard control and water management to correctly capture the ECL features in modeling, in particular, to reproduce the observed spatial rainfall patterns. We simulated eight ECL events using WRF with 36 model configurations, each comprising physics scheme combinations of two planetary boundary layer (pbl), two cumulus (cu), three microphysics (mp), and three radiation (ra) schemes. The performance of each physics scheme combination and the ensembles of multiple physics scheme combinations were evaluated separately. Results show that using the ensemble average gives higher skill than the median performer within the ensemble. More importantly, choosing a composite average of the better performing pbl and cu schemes can substantially improve the representation of high rainfall both spatially and quantitatively.

Linking climate projections to performance: A yield‐based decision scaling assessment of a large urban water resources system

Despite a decade of research into climate change impacts on water resources, the scientific community has delivered relatively few practical methodological developments for integrating uncertainty into water resources system design. This paper presents an application of the “decision scaling” methodology for assessing climate change impacts on water resources system performance and asks how such an approach might inform planning decisions. The decision scaling method reverses the conventional ethos of climate impact assessment by first establishing the climate conditions that would compel planners to intervene. Climate model projections are introduced at the end of the process to characterize climate risk in such a way that avoids the process of propagating those projections through hydrological models. Here we simulated 1000 multisite synthetic monthly streamflow traces in a model of the Melbourne bulk supply system to test the sensitivity of system performance to variations in streamflow statistics. An empirical relation was derived to convert decision-critical flow statistics to climatic units, against which 138 alternative climate projections were plotted and compared. We defined the decision threshold in terms of a system yield metric constrained by multiple performance criteria. Our approach allows for fast and simple incorporation of demand forecast uncertainty and demonstrates the reach of the decision scaling method through successful execution in a large and complex water resources system. Scope for wider application in urban water resources planning is discussed.

Natural hazards in Australia: floods

Floods are caused by a number of interacting factors, making it remarkably difficult to explain changes in flood hazard. This paper reviews the current understanding of historical trends and variability in flood hazard across Australia. Links between flood and rainfall trends cannot be made due to the influence of climate processes over a number of spatial and temporal scales as well as landscape changes that affect the catchment response. There are also still considerable uncertainties in future rainfall projections, particularly for sub-daily extreme rainfall events. This is in addition to the inherent uncertainty in hydrological modelling such as antecedent conditions and feedback mechanisms.Research questions are posed based on the current state of knowledge. These include a need for high-resolution climate modelling studies and efforts in compiling and analysing databases of sub-daily rainfall and flood records. Finally there is a need to develop modelling frameworks that can deal with the interaction between climate processes at different spatio-temporal scales, so that historical flood trends can be better explained and future flood behaviour understood.

The Lund instrumental record of meteorological observations: reconstruction of monthly sea-level pressure 1780–1997

The reconstructed surface air pressure series from Lund, southern Sweden, covers the period 1780-1997 and comprises mon than 234000 valid observations (three observations per day), i.e. > 98% of all possible observation occasions. For the Early Instrumental Period (EIP; 1780-1860) data were digitised from the original records. For most of the Modern Instrumental Period (MIP; 1861-) a series was compiled from various databases containing instrument corrected data. During EIP, the series of raw monthly means show several substantial inhomogeneities. With the aid of a detailed reconstruction of the station history, it was possible to remove almost all inhomogeneities during EIP by applying the correct instrument corrections (for barometer temperature, to standard gravity and to mean sea-level pressure) to the series of original observations. In particular, corrections for the temperature and altitude of the barometer eliminated several inhomogeneities. A prerequisite for applying these corrections is the availability of high-resolution data (actual raw observations or daily averages). Further homogenisation was attained by intercomparison of the monthly mean pressure with acknowledged homogeneous series (mainly the UKMO monthly grid, station records from Copenhagen and Edinburgh). Statistical tests of homogeneity showed that no substantial inhomogeneities remain in the final version. The modern part of the final monthly pressure series largely follows that of the southern Baltic Sea region. Furthermore, it shows relatively high pressure during spring (MAM) in the period 1780-1820, which was paralleled by severe wind erosion in southern Scandinavia during this time. Relatively high pressure throughout the year is also notable during a period of precipitation deficit in 1970s. Copyright (C) 1999 Royal Meteorological Society.

Natural hazards in Australia: heatwaves

As part of a special issue on natural hazards, this paper reviews the current state of scientific knowledge of Australian heatwaves. Over recent years, progress has been made in understanding both the causes of and changes to heatwaves. Relationships between atmospheric heatwaves and large-scale and synoptic variability have been identified, with increasing trends in heatwave intensity, frequency and duration projected to continue throughout the 21st century. However, more research is required to further our understanding of the dynamical interactions of atmospheric heatwaves, particularly with the land surface. Research into marine heatwaves is still in its infancy, with little known about driving mechanisms, and observed and future changes. In order to address these knowledge gaps, recommendations include: focusing on a comprehensive assessment of atmospheric heatwave dynamics; understanding links with droughts; working towards a unified measurement framework; and investigating observed and future trends in marine heatwaves. Such work requires comprehensive and long-term collaboration activities. However, benefits will extend to the international community, thus addressing global grand challenges surrounding these extreme events.

Rainfall in Australia's eastern seaboard: a review of confidence in projections based on observations and physical processes

There is a high degree of variation in rainfall projections for later this century for Australias eastern seaboard, partly because of how different climate models represent the relevant physical processes. These processes include local environmental conditions, synoptic phenomena and large-scale atmospheric and oceanic modes of variability. We review these processes using a wide range of analyses from observations and modelling. A synthesis of this review is used to produce likelihood and confidence measures for the rainfall projections, intended to inform planning and adaptation.

Comparison of various climate change projections of eastern Australian rainfall

The Australian eastern seaboard is a distinct climate entity from the interior of the continent, with different climatic influences on each side of the Great Dividing Range. Therefore, it is plausible that downscaling of global climate models could reveal meaningful regional detail, or ‘added value’, in the climate change signal of mean rainfall change in eastern Australia un-der future scenarios. However, because downscaling is typically done using a limited set of global climate models and downscaling methods, the results from a downscaling study may not represent the range of uncertainty in plausible projected change for a region suggested by the ensemble of host global climate models. A complete and unbiased representation of the plausible changes in the climate is essential in producing climate projections useful for future planning. As part of this aim it is important to quantify any differences in the change signal between global climate models and downscaling, and understand the cause of these differ-ences in terms of plausible added regional detail in the climate change signal, the impact of sub-sampling global climate models and the effect of the downscaling models themselves. Here we examine rainfall projections in eastern Australia under a high emissions scenario by late in the century from ensembles of global climate models, two dynamical downscaling models and one statistical downscaling model. We find no cases where all three downscaling methods show the same clear regional spatial detail in the change signal that is distinct from the host models. However, some downscaled projections suggest that the eastern seaboard could see little change in spring rainfall, in contrast to the substantial rainfall decrease inland. The change signal in the downscaled outputs is broadly similar at the large scale in the various model outputs, with a few notable exceptions. For example, the model median from dynamical downscaling projects a rainfall increase over the entirety of eastern Australia in autumn that is greater than the global models. Also, there are some instances where a downscaling method produces changes outside the range of host models over eastern Australia as a whole, thus ex-panding the projected range of uncertainty. Results are particularly uncertain for summer, where no two downscaling studies clearly agree. There are also some confounding factors from the model configuration used in downscaling, where the particular zones used for statis-tical models and the model components used in dynamical models have an influence on results and produce additional uncertainty.

Extreme precipitation in WRF during the Newcastle East Coast Low of 2007

In the context of regional downscaling, we study the representation of extreme precipitation in the Weather Research and Forecasting (WRF) model, focusing on a major event that occurred on the 8th of June 2007 along the coast of eastern Australia (abbreviated “Newy”). This was one of the strongest extra-tropical low-pressure systems off eastern Australia in the last 30 years and was one of several storms comprising a test bed for the WRF ensemble that underpins the regional climate change projections for eastern Australia (New South Wales/Australian Capital Territory Regional Climate Modelling Project, NARCliM). Newy provides an informative case study for examining precipitation extremes as simulated by WRF set up for regional downscaling. Here, simulations from the NARCliM physics ensemble of Newy available at ∼10 km grid spacing are used. Extremes and spatio-temporal characteristics are examined using land-based daily and hourly precipitation totals, with a particular focus on hourly accumulations. Of the different physics schemes assessed, the cumulus and the boundary layer schemes cause the largest differences. Although the Betts-Miller-Janjic cumulus scheme produces better rainfall totals over the entire storm, the Kain-Fritsch cumulus scheme promotes higher and more realistic hourly extreme precipitation totals. Analysis indicates the Kain-Fritsch runs are correlated with larger resolved grid-scale vertical moisture fluxes, which are produced through the influence of parameterized convection on the larger-scale circulation and the subsequent convergence and ascent of moisture. Results show that WRF qualitatively reproduces spatial precipitation patterns during the storm, albeit with some errors in timing. This case study indicates that whilst regional climate simulations of an extreme event such as Newy in WRF may be well represented at daily scales irrespective of the physics scheme used, the representation at hourly scales is likely to be physics scheme dependent.