Search

Climate change is typically modeled using sophisticated mathematical models (climate models) of physical processes that range in temporal and spatial scales. Multi-model ensemble means of climate models show better correlation with the observations than any of the models separately. Currently, an open research question is how climate models can be combined to create an ensemble mean in an optimal way. We present a novel stochastic approach based on Markov chains to estimate model weights in order to obtain ensemble means. The method was compared to existing alternatives by measuring its performance on training and validation data, as well as model-as-truth experiments. The Markov chain method showed improved performance over those methods when measured by the root mean squared error in validation and comparable performance in model-as-truth experiments. The results of this comparative analysis should serve to motivate further studies in applications of Markov chain and other nonlinear methods that address the issues of finding optimal model weight for constructing ensemble means.

The World Climate Research Programme (WCRP) has an international initiative called the COordinated Regional climate Downscaling EXperiment (CORDEX). The goal of the initiative is to provide regionally downscaled climate projections for most land regions of the globe, as a compliment to the global climate model projections performed within the Coupled Model Intercomparison Projects (CMIP). CORDEX includes data from both dynamical and statistical downscaling. It is anticipated that the CORDEX dataset will provide a link to the impacts and adaptation community through its better resolution and regional focus. Participation in CORDEX is open and any researchers performing climate downscaling are encourage to engage with the initiative. Here I present the current status, evaluation and future projections for the CORDEX-AustralAsia ensemble.The CORDEX-Australasia ensemble is the largest regional climate projection ensemble ever created for the region. It is a 20-member ensemble made by 6 regional climate models downscaling 11 global climate models. Overall the ensemble produces a good representation of recent climate. Consistent biases within the ensemble include an underestimation of the diurnal temperature range and an underestimation of precipitation across much of southern Australia. Under a high emissions scenario projected temperature changes by the end of the twenty-first century reach ~ 5 K in the interior of Australia with smaller increases found toward the coast. Projected precipitation changes are towards drying, particularly in the most populated areas of the southwest and southeast of the continent. The projected precipitation change is very seasonal with summer projected to see little change leaning toward an increase. These results provide a foundation enabling future studies of regional climate changes, climate change impacts, and adaptation options for Australia.

Coarse resolution global climate models (GCM) cannot resolve fine-scale drivers of regional climate, which is the scale where climate adaptation decisions are made. Regional climate models (RCMs) generate high-resolution projections by dynamically downscaling GCM outputs. However, evidence of where and when downscaling provides new information about both the current climate (added value, AV) and projected climate change signals, relative to driving data, is lacking. Seasons and locations where CORDEX-Australasia ERA-Interim and GCM-driven RCMs show AV for mean and extreme precipitation and temperature are identified. A new concept is introduced, ‘realised added value’, that identifies where and when RCMs simultaneously add value in the present climate and project a different climate change signal, thus suggesting plausible improvements in future climate projections by RCMs. ERA-Interim-driven RCMs add value to the simulation of summer-time mean precipitation, especially over northern and eastern Australia. GCM-driven RCMs show AV for precipitation over complex orography in south-eastern Australia during winter and widespread AV for mean and extreme minimum temperature during both seasons, especially over coastal and high-altitude areas. RCM projections of decreased winter rainfall over the Australian Alps and decreased summer rainfall over northern Australia are collocated with notable realised added value. Realised added value averaged across models, variables, seasons and statistics is evident across the majority of Australia and shows where plausible improvements in future climate projections are conferred by RCMs. This assessment of varying RCM capabilities to provide realised added value to GCM projections can be applied globally to inform climate adaptation and model development.

Climate change is typically modelled using sophisticated mathematical models (Climate Models) of physical processes taking place over long periods of time. Multi-model ensembles of climate models show better correlation with the observations than any of the models separately. Currently, an open research question is how climate models can be combined to create an ensemble in an optimal way. We present a novel approach based on Markov chains to estimate model weights in order to obtain ensemble means. The method was compared to existing alternatives by measuring its performance on training and validation data. The Markov chain method showed improved performance over those methods when measured by the root mean squared error and the R-squared metrics. The results of this comparative analysis should serve to motivate further studies in Markov chain and other nonlinear methods application, that address the issues of finding optimal model weight for constructing ensemble means.

South East Australia is characterised by a diverse climate ranging from lush, temperate mountain ranges to hot and arid grasslands. The region is home to Australia's largest river system, the Murray-Darling. The Murray-Darling Basin is an important agricultural region, generating almost 50% of Australia's total irrigated agricultural production in 2018. Rainfall in this region is typically highly variable and subject to severe drought. The Millennium Drought (2001-2009), widely known as the worst drought on record and one of the most severe in the world, has now been superseded by a worse drought (2017-present), setting a new extreme in the drought record. During the current drought, rainfall, root zone soil moisture and water storages have reached record-breaking low levels. High temperatures have also broken historical records on multiple occasions since the drought began. Drought conditions and exceptionally high temperatures have dried the landscape, which has led to intense bushfires that have so far ravaged over 5 million hectares.Yet the degree to which the land surface exacerbates drought in the Murray-Darling Basin remains unknown. In other words, the relative importance of local versus remote processes affecting rainfall, particularly during drought, is uncertain. Where does the moisture come from, and how strongly do local land surface processes attenuate or amplify this atmospheric moisture to affect local rainfall? Establishing the evaporative source regions that supply moisture for rainfall can help reveal the mechanisms driving anomalously low rainfall. In the case of drought, it can help reveal whether anomalous rainfall was due to a reduction in source evaporation, anomalous atmospheric circulation (i.e., the moisture was generated but transported somewhere else), land surface control on the atmosphere through feedbacks, or a combination of factors.We used a Lagrangian back-trajectory approach to determine the long-term average evaporative source regions that supply Australia's rainfall, and the level of recycling that rainfall undergoes. The back-trajectory model tracked water vapour from the location of rainfall events backward in time and space and identified the evaporative origin. From this, we calculated the proportion of rainfall falling across the Murray-Darling Basin that originated as evapotranspiration from the Basin itself; that is, the rainfall recycling ratio.By combining this long-term baseline of source region and rainfall recycling with anomalies of source region evaporation and local atmospheric boundary layer properties, we found that the drivers of low rainfall changed through time during the Millennium Drought. At the peak of the Drought the anomalously low rainfall was driven by a lack of atmospheric moisture advected from the identified typical source region; at other times the low rainfall was due to local conditions unfavorable for the precipitation of available moisture. Overall we found that land surface control on the atmosphere exacerbated the Millennium Drought by approximately 10%.

Wildfire can become a catastrophic natural hazard, especially during dry summer seasons in Australia. Severity is influenced by various meteorological, geographical, and fuel characteristics. Modified Mark 4 McArthur's Grassland Fire Danger Index (GFDI) is a commonly used approach to determine the fire danger level in grassland ecosystems. The degree of curing (DOC, i.e. proportion of dead material) of the grass is one key ingredient in determining the fire danger. It is difficult to collect accurate DOC information in the field, and therefore ground-observed measurements are rather limited. In this study, we explore the possibility of whether adding satellite-observed data responding to vegetation water content (vegetation optical depth, VOD) will improve DOC prediction when compared with the existing satellite-observed data responding to DOC prediction models based on vegetation greenness (normalised difference vegetation index, NDVI). First, statistically significant relationships are established between selected ground-observed DOC and satellite-observed vegetation datasets (NDVI and VOD) with an r2 up to 0.67. DOC levels estimated using satellite observations were then evaluated using field measurements with an r2 of 0.44 to 0.55. Results suggest that VOD-based DOC estimation can reasonably reproduce ground-based observations in space and time and is comparable to the existing NDVI-based DOC estimation models.

Wildfire can become a catastrophic natural hazard, especially during dry summer seasons in Australia. Severity is influenced by various meteorological, geographical, and fuel characteristics. Modified Mark 4 McArthur's Grassland Fire 10 Danger Index (GFDI) is a commonly used approach to determine the fire danger level in grassland ecosystems. The degree of curing (DOC, i.e. proportion of dead material) of the grass is one key ingredient in determining the fire danger. It is difficult to collect accurate DOC information in the field, therefore, ground observed measurements are rather limited. In this study, we used satellite observed vegetation greenness (Normalised Difference Vegetation Index, NDVI) and vegetation water content (Vegetation Optical Depth, VOD) information to improve the accuracy of the DOC estimation. First, a statistically 15 significant relationship is established between selected ground observed DOC and satellite observed vegetation datasets (NDVI and VOD) with an r2 of 0.67. DOC levels estimated using satellite observations were then evaluated using field measurements with an r2 of 0.55. Results suggest that satellite based DOC estimation can reasonably reproduce ground based observations in space and time. Comparison with currently available satellite based DOC products shows that our model has a comparable and arguably more balanced performance.

Land surface albedo, the fraction of incoming solar radiation reflected by the land surface, is a key component of the Earth system. This study evaluates snow-free surface albedo simulations by the Community Atmosphere Biosphere Land Exchange (CABLEv1.4b) model with the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Satellite Pour L'Observation de la Terre (SPOT) albedo. We compare results from offline simulations over the Australian continent. The control simulation has prescribed background snow-free and vegetation-free soil albedo derived from MODIS whilst the experiments use a simple parameterisation based on soil moisture and colour, originally from the Biosphere Atmosphere Transfer Scheme (BATS), and adopted in the Common Land Model (CLM). The control simulation, with prescribed soil albedo, shows that CABLE simulates overall albedo over Australia reasonably well, with differences compared to MODIS and SPOT albedos within ±0.1. Application of the original BATS scheme, which uses an eight-class soil classification, resulted in large differences of up to −0.25 for the near-infrared (NIR) albedo over large parts of the desert regions of central Australia. The use of a recalibrated 20-class soil colour classification from the CLM, which includes a higher range for saturated and VIS (visible) and NIR soil albedos, reduced the underestimation of the NIR albedo. However, this soil colour mapping is tuned to CLM soil moisture, a quantity which is not necessarily transferrable between land surface models. We therefore recalibrated the soil color map using CABLE's climatological soil moisture, which further reduced the underestimation of the NIR albedo to within ±0.15 over most of the continent as compared to MODIS and SPOT albedos. Small areas of larger differences of up to −0.25 remained within the central arid parts of the continent during summer; however, the spatial extent of these large differences is substantially reduced as compared to the simulation using the default eight-class uncalibrated soil colour map. It is now possible to use CABLE coupled to atmospheric models to investigate soil-moisture–albedo feedbacks, an important enhancement of the model.

Vorticity-driven lateral fire spread (VLS) is a form of dynamic fire behaviour, during which a wildland fire spreads rapidly across a steep leeward slope in a direction approximately transverse to the background winds. VLS is often accompanied by a downwind extension of the active flaming region and intense pyro-convection. In this study, the WRF-Fire (WRF stands for Weather Research and Forecasting) coupled atmosphere–fire model is used to examine the sensitivity of resolving VLS to both the horizontal and vertical grid spacing, and the fire-to-atmosphere coupling from within the model framework. The atmospheric horizontal and vertical grid spacing are varied between 25 and 90 m, and the fire-to-atmosphere coupling is either enabled or disabled. At high spatial resolutions, the inclusion of fire-to-atmosphere coupling increases the upslope and lateral rate of spread by factors of up to 2.7 and 9.5, respectively. This increase in the upslope and lateral rate of spread diminishes at coarser spatial resolutions, and VLS is not modelled for a horizontal and vertical grid spacing of 90 m. The lateral fire spread is driven by fire whirls formed due to an interaction between the background winds and the vertical circulation generated at the flank of the fire front as part of the pyro-convective updraft. The laterally advancing fire fronts become the dominant contributors to the extreme pyro-convection. The results presented in this study demonstrate that both high spatial resolution and two-way atmosphere–fire coupling are required to model VLS with WRF-Fire.

We present a novel Bayesian statistical approach to computing model weights in climate change We present a novel Bayesian statistical approach to computing model weights in climate change projection ensembles. The weight of each climate model is obtained by weighting the current day observed data under the posterior distribution admitted under competing climate models. We use a linear model to describe the model output and observations. The approach accounts for uncertainty in model bias, trend and internal variability, as well as including error in the observations used. Our framework is general, requires very little problem specific input, and works well with default priors. We carry out cross-validation checks that confirm that the method produces the correct coverage.

In this study, we have examined the ability of a regional climate model (RCM) to simulate the extended drought that occurred throughout the period of 2002 through 2007 in south-east Australia. In particular, the ability to reproduce the two drought peaks in 2002 and 2006 was investigated. Overall, the RCM was found to reproduce both the temporal and the spatial structure of the drought-related precipitation anomalies quite well, despite using climatological seasonal surface characteristics such as vegetation fraction and albedo. This result concurs with previous studies that found that about two-thirds of the precipitation decline can be attributed to the El Niño–Southern Oscillation (ENSO). Simulation experiments that allowed the vegetation fraction and albedo to vary as observed illustrated that the intensity of the drought was underestimated by about 10 % when using climatological surface characteristics. These results suggest that in terms of drought development, capturing the feedbacks related to vegetation and albedo changes may be as important as capturing the soil moisture–precipitation feedback. In order to improve our modelling of multi-year droughts, the challenge is to capture all these related surface changes simultaneously, and provide a comprehensive description of land surface–precipitation feedback during the droughts development.

Increases in greenhouse gas concentrations are expected to impact the terrestrial hydrologic cycle through changes in radiative forcings and plant physiological and structural responses. Here, we investigate the nature and frequency of non-stationary hydrological response as evidenced through water balance studies over 166 anthropogenically unaffected catchments in Australia. Non-stationarity of hydrologic response is investigated through analysis of long-term trend in annual runoff ratio (1984–2005). Results indicate that a significant trend (p

Combining information derived from satellite-based passive and active microwave sensors has the potential to offer improved retrievals of surface soil moisture variations at global scales. Here we propose a technique to take advantage of retrieval characteristics of passive (AMSR-E) and active (ASCAT) microwave satellite estimates over sparse-to-moderately vegetated areas to obtain an improved soil moisture product. To do this, absolute soil moisture values from AMSR-E and relative soil moisture derived from ASCAT are rescaled against a reference land surface model date set using a cumulative distribution function (CDF) matching approach. While this technique imposes the bias of the reference to the rescaled satellite products, it adjusts both satellite products to the same range and almost preserves the correlation between satellite products and in situ measurements. Comparisons with in situ data demonstrated that over the regions where the correlation coefficient between rescaled AMSR-E and ASCAT is above 0.65 (hereafter referred to as transitional regions), merging the different satellite products together increases the number of observations while minimally changing the accuracy of soil moisture retrievals. These transitional regions also delineate the boundary between sparsely and moderately vegetated regions where rescaled AMSR-E and ASCAT are respectively used in the merged product. Thus the merged product carries the advantages of better spatial coverage overall and increased number of observations particularly for the transitional regions. The combination approach developed in this study has the potential to be applied to existing microwave satellites as well as to new microwave missions. Accordingly, a long-term global soil moisture dataset can be developed and extended, enhancing basic understanding of the role of soil moisture in the water, energy and carbon cycles.

Many studies bias correct daily precipitation from climate models to match the observed precipitation statistics, and the bias corrected data are then used for various modelling applications. This paper presents a review of recent methods used to bias correct precipitation from regional climate models (RCMs). The paper then assesses four bias correction methods applied to the weather research and forecasting (WRF) model simulated precipitation, and the follow-on impact on modelled runoff for eight catchments in southeast Australia. Overall, the best results are produced by either quantile mapping or a newly proposed two-state gamma distribution mapping method. However, the difference between the tested methods is small in the modelling experiments here (and as reported in the literature), mainly because of the substantial corrections required and inconsistent errors over time (non-stationarity). The errors remaining in bias corrected precipitation are typically amplified in modelled runoff. The tested methods cannot overcome limitation of RCM in simulating precipitation sequence, which affects runoff generation. Results further show that whereas bias correction does not seem to alter change signals in precipitation means, it can introduce additional uncertainty to change signals in high precipitation amounts and, consequently, in runoff. Future climate change impact studies need to take this into account when deciding whether to use raw or bias corrected RCM results. Nevertheless, RCMs will continue to improve and will become increasingly useful for hydrological applications as the bias in RCM simulations reduces.

The strength of land–atmosphere coupling during the onset (September) through to the peak (February) of the wet season over Northern Australia is statistically diagnosed using ensembles of land surface model simulations that produce a range of different background soil moisture states. We derive coupling strength between the soil moisture and the planetary boundary layer via a statistical measure of association. The simulated evaporative fraction and the boundary layer are shown to be strongly coupled during both SON and DJF despite the differing background soil moisture states between the two seasons as among the ensemble members. The sign and magnitude of the surface layer soil moisture based coupling strength during the onset of the wet season (SON) differs from the coupling between the evaporative fraction and boundary layer from the same season, and the coupling between the surface soil moisture and boundary layer coupling during DJF. The patterns and magnitude of the surface flux-boundary layer coupling are not captured when coupling is diagnosed using the surface layer soil moisture alone. The conflicting results arise because the surface layer soil moisture lacks strong association with the atmosphere during the monsoon onset because the evapotranspiration is dominated by transpiration. Our results indicate that accurately diagnosing coupling strength in seasonally dry regions, such as Northern Australia, requires root zone soil moisture to be included.