Search

We develop high-resolution (1 km) estimates of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) over the Australian continent for the period January 2003 to June 2022 by empirical upscaling of flux tower measurements. We compare our estimates with nine other products that cover the three broad categories that define current methods for estimating the terrestrial carbon cycle and assess if consiliences between datasets can point to the correct dynamics of Australia's carbon cycle. Our results indicate that regional empirical upscaling greatly improves upon the existing global empirical upscaling efforts, outperforms process-based models, and agrees much better with the dynamics of CO 2 flux over Australia as estimated by two regional atmospheric inversions. Our nearly 20-year estimates of terrestrial carbon fluxes revealed that Australia is a strong net carbon sink of -0.44 PgC yr -1 (interquartile range, IQR = 0.42 PgC yr -1) on average, with an inter-annual variability of 0.18 PgC yr -1 and an average seasonal amplitude of 0.85 PgC yr -1. Annual mean carbon uptake estimated from other methods ranged considerably, while carbon flux anomalies showed much better agreement between methods. NEE anomalies were predominately driven by cumulative rainfall deficits and surpluses, resulting in larger anomalous responses from GPP than ER. In contrast, we show that the long-term average seasonal cycle is dictated more by the variability in ER than GPP, resulting in peak carbon uptake typically occurring during the cooler, drier austral autumn and winter months. This new estimate of Australia's terrestrial carbon cycle provides a benchmark for assessment against land surface model simulations and a means for monitoring of Australia's terrestrial carbon cycle at an unprecedented high resolution. We call this new estimate of Australia's terrestrial carbon cycle "AusEFlux" (Australian Empirical Fluxes). [ABSTRACT FROM AUTHOR]

The McArthur grassland and forest fire danger indices, widely used in Australia, predict six fire danger classes from 'Low-Moderate' to 'Catastrophic.' These classes were linked to the rate of fire spread and difficulty of suppression. However, the lack of rate of fire spread data, especially for elevated fire danger classes, has hindered improvement of the McArthur methodology or an alternate approach. We explored the relationship between fire danger classes and burned areas (derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument) within six climate zones during the 2000–2016 Australian fire seasons. A negative binomial linear regression model was used to explore this relationship. The fire danger classes demonstrated a corresponding increase in burned area from 'Low-Moderate' to 'Very High' classes in Australia's inland regions. The elevated fire danger classes did not contribute to this trend. In coastal regions, the satellite-derived burned area showed no relationship between fire danger classes and satellite-derived burned area. We used accumulated burned area from the daily MODIS product, which could be subjected to lagged detection as observed in the Kilmore East fire. Thus, the satellite-derived total burned area may not be a suitable metric for informing the McArthur fire danger classes across Australia. We explored the relationship of the McArthur fire danger classes with the satellite-derived daily burned area The findings suggest that satellite-derived daily burned area may not be a suitable metric for informing McArthur fire danger classes across Australia because it can be subjected to lagged detection. [ABSTRACT FROM AUTHOR]

Studies using remote sensing data for fire danger prediction have primarily relied on fire ignitions data to develop fire danger indices (FDIs). However, these data may only represent conditions suitable for ignition but may not represent fire danger conditions causing escalating fire size. The fire-related response variable's scalability is a key factor that forms a basis for an FDI to include a broader range of fire danger conditions. Remote sensing derived fire size is a scalable fire characteristic encapsulating all possible fire sizes that previously occurred in the landscape, including extreme fire events. Consequently, we propose a new FDI that uses remote sensing derived fire size as a response variable. We computed fire sizes from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument burned area. We applied random forest (RF) and logistic regression (LR) to develop the FDI for Australia. RF models performed better than LR, and the higher predicted probabilities demonstrated higher chances for ignited fires to be escalated to larger fire sizes at a regional scale across Australia. However, the predicted probabilities cannot be related to the specific range of fire sizes due to data limitations. Further research with higher temporal and spatial resolution data of both the response and predictor variables can help establish a better relationship between a specific range of fire sizes and the predicted probabilities. [ABSTRACT FROM AUTHOR]

Spectral indices derived from satellite optical remote sensing data have typically been used for fire severity estimation, although other remote sensing systems such as Light Detection and Ranging (LiDAR) are increasingly applied. Despite a multitude of remotely sensed fire severity estimation methods, comparisons of method performance are few. Insights into the merits and limitations of remotely sensed fire severity methods help develop appropriate spatial tools for the management of fire-affected areas. We evaluated the performance of seven passive (optical) and active (LiDAR) remotely sensed fire severity estimation methods in classifying and explaining variation in a field-estimated modified Composite Burn Index (MCBI) for a recent large wildfire in south-eastern Australia. Our evaluation included three commonly applied indices; the differenced Normalized Burn Ratio (dNBR), Relative dNBR (RdNBR) and Relative Burn Ratio (RBR). We compared these NBR indices against two recently proposed fire severity estimation methods that have not previously been evaluated with CBI field data – the Vegetation Structure Perpendicular Index (VSPI) spectral index and the LiDAR point cloud-derived Profile Area Change (PAC), along with experimental relativized forms of these indices (RVSPI and RPAC, respectively). The RVSPI (κ = 0.47) demonstrated similar overall classification accuracy (N classes = 4) to the PAC (κ = 0.48), however both indices had lower classification accuracy than the dNBR (κ = 0.59), RdNBR (κ = 0.59) and RBR (κ = 0.61). The VSPI and PAC were unable to accurately represent non-structural changes caused by lower severity fire. Application of these optical and LiDAR indices should consider their discussed limitations in relation to the objectives of their application. [ABSTRACT FROM AUTHOR]

Topographic models provide a useful tool for understanding gully occurrence in the landscape but require reliable estimates of gully head drainage areas. Modern high‐resolution topography data (collected using structure from motion photogrammetry or light detection and ranging) is increasingly used for topographic studies of gullies, but little work has been done to assess the variability of gully head drainage area estimates using different methods. This study evaluated alternative approaches to using high‐resolution digital elevation models (DEMs) so that gully topographic models can be more readily applied to any area with suitably high‐resolution data. Specifically, we investigated the impact of single‐ or multiple‐direction flow routing algorithms, DEM hydrologic‐enforcement procedures and spatial resolution on gully head drainage area estimation. We tested these methods on a 40 km2 site centred on Weany Creek, a low‐relief semi‐arid landscape draining towards the Great Barrier Reef, Australia. Using a subroutine to separate gully heads into those with divergent or convergent flow patterns upslope, we found that divergent flow conditions occurred at half of 484 studied gullies. Drainage areas estimated by different flow routing algorithms were more variable in these divergent cases than for convergent cases. This variation caused a significant difference between topographic threshold parameters (slope b and intercept k) derived from single‐ or multiple‐direction flow routing algorithms, respectively. Different methods of hydrologic enforcement (filling or breaching) also affected threshold analysis, resulting in estimates of the exponent b being ~188% higher if the DEM was filled than if breached. The testing of the methods to date indicates that a finer resolution (≤2 m) DEM and a multiple‐direction flow routing algorithm achieve the most realistic drainage area estimates in low‐relief landscapes. For Weany Creek we estimated threshold parameters k = 0.033 and b = 0.189, indicating that it is highly susceptible to gully erosion. [ABSTRACT FROM AUTHOR]

There is a growing interest globally in the spatial distribution and temporal dynamics of intermittently flowing streams and rivers, and how this varies in relation to climatic and other environmental factors. However, biases in the distribution of stream gauges may give a misleading impression of spatial-temporal variations in streamflow intermittency within river networks. Here, we developed an approach to quantify catchment-wide streamflow intermittency over long time frames and in a spatially explicit manner, using readily accessible and spatially contiguous daily runoff data from a national-scale water balance model. We examined the ability of the water balance model to simulate streamflow in two hydro-climatically distinctive (subtropical and temperate) regions in Australia, with a particular focus on low-flow simulations. We also evaluated the effect of model time step (daily vs. monthly) on flow intermittency estimation to inform future model selection. The water balance model showed better performance in the temperate region characterised by steady baseflow than in the subtropical region with flashy hydrographs and frequent cease-to-flow periods. The model tended to overestimate low-flow magnitude mainly due to overestimation of gains (e.g. groundwater release to baseflow) during low-flow periods. Modelled patterns of flow intermittency revealed highly dynamic behaviour in space and time, with cease-to-flow events affecting between 29 and 80 % of the river network over the period of 1911–2016, using a daily streamflow model. The daily flow model did not perform better than the monthly flow model in quantifying flow intermittency at a monthly time step, and model selection should depend on the intended application of the model outputs. Our general approach to quantifying spatio-temporal patterns of flow intermittency is transferable to other parts of the world, and it can inform hydro-ecological understanding and management of intermittent streams where limited gauging data are available. [ABSTRACT FROM AUTHOR]

The low spatial density of streamflow gauging stations limits the accuracy of spatial streamflow estimates in many parts of the world. Strategies to improve runoff estimates in the absence of dense measurements have tended to focus on estimating parameters of runoff models in ungauged regions, through so‐called parameter regionalization methods. However, parameter regionalization can be affected by overdependence on calibration at gauged sites, model parameter equifinality, and ensuing estimation errors. As a result, spatial model runoff estimates typically exhibit spatially correlated biases. This analysis attempts to enhance the use of observations in spatial runoff estimation. Specifically, we assessed the potential to reduce systematic errors by spatially interpolating residuals (i.e., errors) between prior grid‐based streamflow estimates for Australia at 0.05° × 0.05° grid from the Australian Bureau of Meteorology's calibrated, operational Australian Water Resources Assessment Landscape model (AWRA‐L) and streamflow gauging records from 780 unimpeded, relatively small catchments. We analyzed spatial autocorrelation in residuals and tested an efficient two‐step correction approach involving a uniform correction and subsequent kriging of residuals. The approach removed an average of 41% of systematic bias in the model estimates and also improved other model performance measures. Further reduction in errors at shorter timescales may be achievable through a temporally hierarchical correction scheme. Key Points: We developed a method to reduce spatially variable bias in modeled runoff at the continental scaleThe two‐step regression and kriging scheme removed an average 41% of bias in cross validationCorrections were applied to continuous daily fields of runoff at a 0.05 × 0.05° grid across Australia [ABSTRACT FROM AUTHOR]

Phenology is the study of periodic biological occurrences and can provide important insights into the influence of climatic variability and change on ecosystems. Understanding Australia's vegetation phenology is a challenge due to its diverse range of ecosystems, from savannas and tropical rainforests to temperate eucalypt woodlands, semiarid scrublands, and alpine grasslands. These ecosystems exhibit marked differences in seasonal patterns of canopy development and plant life-cycle events, much of which deviates from the predictable seasonal phenological pulse of temperate deciduous and boreal biomes. Many Australian ecosystems are subject to irregular events (i.e. drought, flooding, cyclones, and fire) that can alter ecosystem composition, structure, and functioning just as much as seasonal change. We show how satellite remote sensing and ground-based digital repeat photography (i.e. phenocams) can be used to improve understanding of phenology in Australian ecosystems. First, we examine temporal variation in phenology on the continental scale using the enhanced vegetation index (EVI), calculated from MODerate resolution Imaging Spectroradiometer (MODIS) data. Spatial gradients are revealed, ranging from regions with pronounced seasonality in canopy development (i.e. tropical savannas) to regions where seasonal variation is minimal (i.e. tropical rainforests) or high but irregular (i.e. arid ecosystems). Next, we use time series colour information extracted from phenocam imagery to illustrate a range of phenological signals in four contrasting Australian ecosystems. These include greening and senescing events in tropical savannas and temperate eucalypt understorey, as well as strong seasonal dynamics of individual trees in a seemingly static evergreen rainforest. We also demonstrate how phenology links with ecosystem gross primary productivity (from eddy covariance) and discuss why these processes are linked in some ecosystems but not others. We conclude that phenocams have the potential to greatly improve the current understanding of Australian ecosystems. To facilitate the sharing of this information, we have formed the Australian Phenocam Network (http://phenocam.org.au/). [ABSTRACT FROM AUTHOR]

Australia is frequently subject to droughts and floods. Its hydrology is strongly connected to oceanic and atmospheric oscillations (climate modes) such as the El Niño-Southern Oscillation (ENSO), Indian Ocean dipole (IOD), and southern annular mode (SAM). A global 32-yr dataset of remotely sensed surface soil moisture (SSM) was used to examine hydrological variations in mainland Australia for the period 1978-2010. Complex empirical orthogonal function (CEOF) analysis was applied to extract independent signals and to investigate their relationships to climate modes. The annual cycle signal represented 46.3% of the total variance and a low but highly significant connection with SAM was found. Two multiannual signals with a lesser share in total variance (6.3% and 4.2%) were identified. The first one had an unstable period of 2-5 yr and reflected an east-west pattern that can be associated with ENSO and SAM but not with IOD. The second one, a 1- to 5-yr oscillation, formed a dipole pattern between the west and north and can be linked to ENSO and IOD. As expected, relationships with ENSO were found throughout the year and are especially strong during southern spring and summer in the east and north. Somewhat unexpectedly, SAM impacts strongest in the north and east during summer and is proposed as the key driver of the annual SSM signal. The IOD explains SSM variations in the north, east, and southeast during spring and also in the west during winter. [ABSTRACT FROM AUTHOR]

Precipitation estimates from reanalyses and satellite observations are routinely used in hydrologie applications, but their accuracy is seldom systematically evaluated. This study used high-resolution gauge-only daily precipitation analyses for Australia (SILO) and South and East Asia [Asian Precipitation--Highly-Resolved Observational Data Integration Towards Evaluation (APHRODITE)] to calculate the daily detection and accuracy metrics for three reanalyses [ECMWF Re-Analysis Interim (ERA-Interim), Japanese 25-yr Reanalysis (JRA-25), and NCEP-Department of Energy (DOE) Global Reanalysis 2] and three satellite-based precipitation products [Tropical Rainfall Measuring Mission (TRMM) 3B42V6, Climate Prediction Center morphing technique (CMORPH), and Precipitation Estimation from Remotely Sensed Imagery Using Artificial Neural Networks (PERSIANN)]. A depth-frequency-adjusted ensemble mean of the reanalyses and satellite products was also evaluated. Re-analyses precipitation from ERA-Interim in southern Australia (SAu) and northern Australasia (NAu) showed higher detection performance. JRA-25 had a better performance in South and East Asia (SEA) except for the monsoon period, in which satellite estimates from TRMM and CMORPH outperformed the reanalyses. In terms of accuracy metrics (correlation coefficient, root-mean-square difference, and a precipitation intensity proxy, which is the ratio of monthly precipitation amount to total days with precipitation) and over the three subdomains, the depth-frequency-adjusted ensemble mean generally outperformed or was nearly as good as any of the single members. The results of the ensemble show that additional information is captured from the different precipitation products. This finding suggests that, depending on precipitation regime and location, combining (re)analysis and satellite products can lead to better precipitation estimates and, thus, more accurate hydrological applications than selecting any single product. [ABSTRACT FROM AUTHOR]

The "Millennium Drought" (2001-2009) can be described as the worst drought on record for southeast Australia. Adaptation to future severe droughts requires insight into the drivers of the drought and its impacts. These were analyzed using climate, water, economic, and remote sensing data combined with biophysical modeling. Prevailing El Niño conditions explained about two thirds of rainfall deficit in east Australia. Results for south Australia were inconclusive; a contribution from global climate change remains plausible but unproven. Natural processes changed the timing and magnitude of soil moisture, streamflow, and groundwater deficits by up to several years, and caused the amplification of rainfall declines in streamflow to be greater than in normal dry years. By design, river management avoided impacts on some categories of water users, but did so by exacerbating the impacts on annual irrigation agriculture and, in particular, river ecosystems. Relative rainfall reductions were amplified 1.5-1.7 times in dryland wheat yields, but the impact was offset by steady increases in cropping area and crop water use efficiency (perhaps partly due to CO2 fertilization). Impacts beyond the agricultural sector occurred (e.g., forestry, tourism, utilities) but were often diffuse and not well quantified. Key causative pathways from physical drought to the degradation of ecological, economic, and social health remain poorly understood and quantified. Combined with the multiple dimensions of multiyear droughts and the specter of climate change, this means future droughts may well break records in ever new ways and not necessarily be managed better than past ones. [ABSTRACT FROM AUTHOR]

On average, Australia is a dry continent with many competing uses for water. Hence, there is an urgent need to know actual evapotranspiration (ETa) patterns across wide areas of agricultural and natural ecosystems, as opposed to just point measurements of ETa. The Australian Government has tasked the science agencies with operationally developing monthly and annual estimates of ETa and other hydrological variables, and with forecasting water availability over periods of days to decades, as part of its national water assessment programme. To meet these needs, Australian researchers have become leaders in developing large-area methods for estimating ETa at regional and continental scales. Ground methods include meteorological models, eddy covariance towers, sap flow sensors and catchment water balance models. Remote sensing methods use thermal infrared, mid infrared and/or vegetation indices usually combined with meteorological data to estimate ETa. Ground and remote sensing ETa estimates are assimilated into the Australian Water Resource Assessment, which issues annual estimates of the state of the continental water balance for policy and planning purposes. The best ETa models are estimated to have an error or uncertainty of 10% to 20% in Australia. Developments in Australian ETa research over the past 20 years are reviewed, and sources of error and uncertainty in current methods and models are discussed. Copyright © 2011 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

Evaluation of fire severity reduction strategies requires the quantification of intervention outcomes and, more broadly, the extent to which fuel characteristics affect fire severity. However, investigations are currently limited by the availability of accurate data on fire severity predictors, particularly relating to fuel. Here, we used airborne LiDAR data collected before the 2019-20 Australian Black Summer fires to investigate the contribution of fuel structure to fire severity under a range of weather conditions. Fire severity was estimated using the Relative Burn Ratio calculated from Sentinel-2 optical remote sensing imagery. We modelled the effects of various fuel structure estimates and other environmental predictors using Random Forest models. In addition to variables estimated at each observation point, we investigated the influence of surrounding landscape characteristics using an innovative method to estimate fireline progression direction. Our models explained 63-76% of fire severity variance using parsimonious predictor sets. Fuel cover in the understorey and canopy, and vertical vegetation heterogeneity, were positively associated with fire severity. Up-fire burnt area and recent planned and unplanned fire reduced fire severity, whereby unplanned fire provided a longer-lasting reduction of fire severity (up to 15 years) than planned fire (up to 10 years). Although fuel structure and land management effects were important predictors, weather and canopy height effects were dominant. By mapping continuous interactions between weather and fuel-related variables, we found strong evidence of diminishing fuel effects below 20-40% relative air humidity. While our findings suggest that land management interventions can provide meaningful fire severity reduction, they also highlight the risk of warmer and drier future climates constraining these advantages., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Access to water is a critical aspect of livestock production, although the relationship between livestock weight gain and water quality remains poorly understood. Previous work has shown that water quality of poorly managed farm dams can be improved by fencing and constructing hardened watering points to limit stock access to the dam, and revegetation to filter contaminant inflow. Here we use cattle weight gain data from three North American studies to develop a cost-benefit analysis for the renovation of farm dams to improve water quality and, in turn, promote cattle weight gain on farms in south-eastern Australia. Our analysis indicated a strong likelihood of positive results and suggested there may be substantial net economic benefit from renovating dams in poor condition to improve water quality. The average per-farm Benefit-Cost Ratios based on deterministic assumptions was 1.5 for New South Wales (NSW) and 3.0 for Victoria in areas where rainfall exceeds 600mm annually. Our analyses suggested that cattle on farms in NSW and Victoria would need to experience additional weight gain from switching to clean water of at least 6.5% and 1.8% per annum respectively, to break even in present value terms. Monte Carlo simulation based on conservative assumptions indicated that the probability of per-farm benefits exceeding costs was greater than 70%. We recommend localised experiments to assess the impact of improved water quality on livestock weight gain in Australian conditions to confirm these expectations empirically., Competing Interests: The authors have declared that no competing interests exist.