Articles | Volume 12, issue 5
Model experiment description paper
24 May 2019
Model experiment description paper |  | 24 May 2019

BARRA v1.0: the Bureau of Meteorology Atmospheric high-resolution Regional Reanalysis for Australia

Chun-Hsu Su, Nathan Eizenberg, Peter Steinle, Dörte Jakob, Paul Fox-Hughes, Christopher J. White, Susan Rennie, Charmaine Franklin, Imtiaz Dharssi, and Hongyan Zhu

The Bureau of Meteorology Atmospheric high-resolution Regional Reanalysis for Australia (BARRA) is the first atmospheric regional reanalysis over a large region covering Australia, New Zealand, and Southeast Asia. The production of the reanalysis with approximately 12 km horizontal resolution – BARRA-R – is well underway with completion expected in 2019. This paper describes the numerical weather forecast model, the data assimilation methods, the forcing and observational data used to produce BARRA-R, and analyses results from the 2003–2016 reanalysis. BARRA-R provides a realistic depiction of the meteorology at and near the surface over land as diagnosed by temperature, wind speed, surface pressure, and precipitation. Comparing against the global reanalyses ERA-Interim and MERRA-2, BARRA-R scores lower root mean square errors when evaluated against (point-scale) 2 m temperature, 10 m wind speed, and surface pressure observations. It also shows reduced biases in daily 2 m temperature maximum and minimum at 5 km resolution and a higher frequency of very heavy precipitation days at 5 and 25 km resolution when compared to gridded satellite and gauge analyses. Some issues with BARRA-R are also identified: biases in 10 m wind, lower precipitation than observed over the tropical oceans, and higher precipitation over regions with higher elevations in south Asia and New Zealand. Some of these issues could be improved through dynamical downscaling of BARRA-R fields using convective-scale (<2 km) models.

1 Introduction

Reanalyses are widely used for climate monitoring and studying climate change as they provide long-term spatially complete records of the atmosphere. This is achieved by using data assimilation techniques that produce an observation-constrained model estimate of the atmosphere. They draw short-term model states towards observations from multiple, disparate sources to form an atmospheric analysis. A physically realistic model provides the means to infer atmospheric states at locations without observations from the limited collection of irregularly distributed observations.

Global-scale reanalyses using global atmospheric general circulation models (GCMs) have advanced in quality and quantity during the past 2 decades (Dee et al., 2014; Hartmann et al., 2013). At present, the available global reanalyses established for the satellite era include the NCEP–NCAR reanalysis at 210 km horizontal resolution (Kalnay et al., 1996), the Japanese 55-year Reanalysis (JRA-55) at 60 km (Ebita et al., 2011), the Modern-Era Retrospective analysis for Research and Applications-2 (MERRA-2) at about 50 km (Gelaro et al., 2017), and the European Centre for Medium-Range Weather Forecasts (ECMWF) ReAnalysis Interim (ERA-Interim) at ∼79 km (Dee et al., 2011). The latter is currently being replaced by the new ERA5 ∼31 km reanalysis (Hersbach and Dee, 2016). These global reanalyses have the advantage of providing globally consistent information, but at the expense of spatial resolution. With resolutions typically greater than 50 km, they may be deficient in accounting for important sub-grid variations in meteorology over heterogeneous terrains and islands, across irregular coastlines, and in terms of other small-scale processes (Mesinger et al., 2006; Randall et al., 2007, and references therein).

To address these shortcomings, the development in global reanalysis has also driven concurrent efforts in statistical approaches and dynamical downscaling (e.g. Dickinson et al., 1989; Fowler et al., 2007; Evans and McCabe, 2013). The latter typically embeds a high-resolution meteorological model within a global reanalysis, whereby effects of small-scale forcing and processes such as convection are modelled. Such development is supported by improvements in non-hydrostatic models that run at high resolution in operational numerical weather prediction (NWP) (e.g. Clark et al., 2016). Regional reanalyses are emerging as a step further in this direction. The first regional reanalysis was the North America Regional Reanalysis (NARR; Mesinger et al., 2006). More recent examples include the Arctic System Reanalysis (ASR; Bromwich et al., 2018), Indian Monsoon Data Assimilation and Analysis (IMDAA; Mahmood et al., 2018), and Uncertainties in Ensembles of Regional Reanalyses (UERRA) in Europe (Borsche et al., 2015, and references therein). In contrast to dynamically downscaled global reanalyses, observations are used in regional reanalyses in the same way as in the global ones to reduce model errors in high-resolution simulations (Bollmeyer et al., 2015). The resulting reanalyses are expected to have better representations of frequency distributions, extremes, and actual space- and time-dependent variability (particularly for near-ground variables). UERRA consists of four regional reanalyses developed by the Swedish Meteorological and Hydrological Institute (SMHI), Météo France, the Deutscher Wetterdienst (DWD), and the UK Met Office (UKMO), producing an ensemble of high-resolution (5–25 km) regional reanalyses of essential climate variables. The SMHI's HARMONIE (High-Resolution Limited-Area Model – HIRLAM – Aire Limitée Adaptation Dynamique Développement International – ALADIN – Regional/Mesoscale Operational NWP in Europe) reanalysis has entered production for the Copernicus Climate Change Service (Ridal et al., 2017).

Regional reanalyses provide significant added value to their global counterparts in diverse applications ranging from traditional climate studies to industry applications, including regional climate change assessments that include local impact studies (e.g. Fall et al., 2010) and extreme event reconstructions (e.g. Zick and Matyas, 2015). As regional reanalyses are generally produced with high spatial and temporal resolution, the extremes of variables at local scales may be quantified more accurately. They are also an alternative reference to evaluate climate projections (e.g. Ruiz-Barradas and Nigam, 2006; Radic and Clarke, 2011). At the same time, embedded forecast models can be used within the framework of the Coordinated Regional Climate Downscaling Experiment (CORDEX; Martynov et al., 2013) within a seamless framework for weather and climate prediction, whereby model deficiencies that differ in spatial scales and timescales can be more readily understood (Brown et al., 2012). They also offer useful datasets for designing new infrastructure, particularly if they are sufficiently long and spatially relevant to define the likelihood of extremes. For renewable energy production, they can provide valuable information on the intermittency (e.g. wind lull) and covariability (e.g. correlation spatially or between variables) of phenomena. For instance, the COSMO (Consortium for Small-scale Modelling) 6 km reanalysis has shown the potential to provide realistic sub-daily representations of winds at 10 to 40 m of height (Borsche et al., 2016) and to resolve small-scale cloud structures (Bollmeyer et al., 2015). NARR was used to define a climatology of surface wind extremes (Malloy et al., 2015) and 30-year trends in wind at hub height (Holt and Wang, 2012) over northern America.

To date, while regional reanalyses exist for North America, Europe, and India, no atmospheric regional reanalysis for the Australasian region has been produced. To close this gap, the Bureau of Meteorology Atmospheric high-resolution Regional Reanalysis for Australia (BARRA; Jakob et al., 2017) has been produced. BARRA is the first atmospheric regional reanalysis that covers Australia, New Zealand, Southeast Asia, and south to the Antarctic ice edge (Fig. 1). It is produced by the Australian Bureau of Meteorology (Bureau), with sponsorship from state fire and governmental agencies across Australia, because of the important advantages it provides for planning and management to reduce risks due to extreme weather events, including bushfires. For instance, BARRA will address the lack of accurate climate information on highly variable surface winds over large areas of Australia due to the low density of the surface observation network in remote areas. BARRA covers a 29-year period from 1990 to 2018, with possible further extensions back and forward in time. The BARRA project delivers a whole-domain reanalysis (identified as BARRA-R) with approximately 12 km horizontal resolution and additional convective-scale (1.5 km horizontal grid-length) downscaling (BARRA-x), nested within BARRA-R, centred on major Australian cities to generate additional high-resolution information needed for local-scale applications and studies. These resulting gridded (12 and 1.5 km) products include a variety of 10 min to hourly surface parameters describing weather and land-surface conditions and hourly upper-air parameters covering the troposphere and stratosphere. The fields on standard pressure levels are generated from vertical interpolation of model-level fields. BARRA serves to lay the foundation for future generations of reanalyses at the Bureau and to further develop its capabilities to produce seamless climate information that integrates its observational networks and NWP programme.

Figure 1The BARRA-R domain is enclosed by the dashed box, with the model orography shown.


In this paper, we describe the forecast model, data assimilation methods, and the forcing and observational data used to produce BARRA-R in Sect. 2. Section 3 provides an initial assessment of the reanalysis system over the first 14 years (2003–2016), with a focus on analysing the quality at or near the surface; Sect. 4 concludes with a brief summary of our findings.

2 The BARRA-R reanalysis

The development of BARRA builds on the Bureau's experience in operational (deterministic) NWP forecasting over the Australian region using the Australian Community Climate and Earth-System Simulator (ACCESS-R) system (Bureau of Meteorology, 2010, 2013; Puri et al., 2013), and BARRA-R is produced using the UKMO's system in UERRA (based on Jermey and Renshaw, 2016) but without the ensemble component. An ensemble NWP forecast system is currently under development at the Bureau. BARRA-R is produced by running a limited-area meteorological forecast model forced with global reanalysis boundary conditions, drawn closer to observations via data assimilation. This section provides an overview of these components, while more technical details are included in the references.

2.1 Forecast model

The Unified Model (UM; Davies et al., 2005) is the grid-point atmospheric model used in BARRA-R and ACCESS. It uses a non-hydrostatic, fully compressible, deep-atmosphere formulation and its dynamical core (Even Newer Dynamics for General atmospheric modelling of the environment, ENDGame) solves the equations of motion using mass-conserving, semi-implicit, semi-Lagrangian, time-integration methods (Wood et al., 2014). The model includes a comprehensive set of parameterizations, including a modified boundary layer scheme based on Lock et al. (2000), a variant of Wilson and Ballard (1999) for mixed-phase cloud microphysics, the mass flux convection scheme of Gregory and Rowntree (1990), and the radiation scheme of Edwards and Slingo (1996), which have all since been improved. Other parameterized sub-grid-scale processes include, fractional cloud cover and orographic drag. More details on all of the physics schemes can be found in Walters et al. (2017a).

The prognostic variables are three-dimensional wind components, virtual dry potential temperature and Exner pressure, dry density, and mixing ratios of moist quantities. The model is discretized on a horizontally staggered Arakawa C grid (Arakawa and Lamb, 1977) and a vertically staggered Charney–Phillips grid (Charney and Phillips, 1953). The staggered arrangement of grid points allows for accurate finite differencing but results in different model fields located on staggered grids displaced by half a grid of spacing along both axes. Data have been left on the staggered grids to allow users to apply the most appropriate re-gridding methods suited for given applications. The vertical levels smoothly transition from terrain-following coordinates near the surface to constant height surfaces in the upper atmosphere (Davies et al., 2005).

BARRA-R uses version 10.2 of the UM and is configured with 70 vertical levels extending from near the surface to 80 km above sea level: 50 model levels below 18 km and 20 levels above this. While configured with this height based on ACCESS-R, we have more confidence in the data up to a height of 25–30 km where we have the most information from observations. The horizontal domain of BARRA-R spans from 65.0 to 196.9 E and −65.0 to 19.4 N (Fig. 1), with constant latitude and longitude increments of 0.11× 0.11 and 1200 × 768 grid points in the horizontal. Our choice of the horizontal resolution follows the deterministic component of the UKMO reanalysis and the IMDAA reanalyses. The model was run to produce 12 h (12 h) forecasts in each 6-hourly cycle (see Sect. 2.2) to give extra data for driving dynamical downscaling within the domain.

The model parameterizations in BARRA-R are inherited from the UKMO Global Atmosphere (GA) 6.0 configurations described in Walters et al. (2017a). The GA6 configurations are also suited for limited-area models with resolutions >10 km but with some modifications, as described below.

  • i.

    A variable Charnock coefficient is used in surface heat exchange over the sea to improve the tropical Pacific air–sea exchange (Ma et al., 2015).

  • ii.

    The heat capacity of the “inland water canopy” is set to 2.11×107 J K−1 m−2 for modelling lakes. This improves the diurnal cycle over inland waters. By contrast, grid cells containing salt lakes in Australia are modelled as bare soil surface (for Lake Eyre and Lake Frome) and vegetated surface (e.g. Lake Lefroy, Lake Ballard).

  • iii.

    For its deep convective mass flux scheme, a grid-box-dependent convective available potential energy (CAPE) closure scheme is chosen to limit the role of parameterized convection. When vertical velocity exceeds the given threshold of 1 m s−1, the vertical-velocity-dependent CAPE closure is chosen to release the convective instability efficiently (Zhu and Dietachmayer, 2015). These changes aim to improve the model stability.

  • iv.

    The river-routing scheme has been turned off because it is not designed for a limited-area model. Therefore, there is no routing of runoff from inland grid points out to sea and inland water bodies, and soil moisture is not affected by this hydrological process.

The characteristics of the lower boundary, climatological fields, and natural and anthropogenic emissions are specified using static ancillary fields. These are created as per Walters et al. (2017a; Table 1), with the exceptions of the land–sea mask and canopy tree heights. The land–sea mask is created from the 1 km resolution International Geosphere–Biosphere Programme (IGBP) land cover data (Loveland et al., 2000), and the canopy tree heights are derived from satellite light detection and ranging (lidar; Simard et al., 2011; Dharssi et al., 2015). Climatological aerosol fields (ammonium sulfate, mineral dust, sea salt, biomass burning, fossil-fuel black carbon, fossil-fuel organic carbon, and secondary organic (biogenic) aerosols) are used to derive the cloud droplet number concentration. Absorption and scattering by aerosols are included in both the shortwave and longwave.

2.1.1 Land surface

The UM uses a community land-surface model, the Joint UK Land Environment Simulator (JULES; Best et al., 2011). It models the partitioning of rainfall into canopy interception, surface runoff, and infiltration and uses the Richards' equation and Darcy's law to model soil hydrology. The sub-grid-scale heterogeneity of soil moisture is represented by the Probability Distributed Moisture (PDM) model (Moore, 2007). A nine-tile approach is used to represent sub-grid-scale heterogeneity in land cover, with the surface of each land point subdivided into five vegetation types (broadleaf trees, needle-leaved trees, temperate C3 grass, tropical C4 grass, and shrubs) and four non-vegetated surface types (urban, inland water, bare soil, and land ice). It describes a 3 m soil column with a four-layer soil scheme with soil thicknesses of 0.1, 0.25, 0.65, and 2.0 m and models vertical heat and water transfer within the column with van Genuchten hydraulic parameters. The JULES urban parameters are optimized for Australia as described by Dharssi et al. (2015).

2.1.2 Soil moisture

For the 1990–2014 period, soil moisture fields in BARRA-R are initialized daily at 06:00 UTC using soil moisture analyses from an offline simulation of JULES at 60 km resolution, driven by bias-corrected ERA-Interim atmosphere forcing data, using methods described in Dharssi and Vinodkumar (2017) and Zhao et al. (2017). The simulation used a 10-year-long spin-up period and was then run continuously for the 1990 to 2014 period. The near-surface soil moisture analyses are found to have good skill for the Australian region when validated against ground-based soil moisture observations (Dharssi and Vinodkumar, 2017). As the offline runs were terminated at the end of December 2014, the daily initialization scheme is continued with soil moisture analyses from the Bureau's global NWP system – ACCESS-G (Bureau of Meteorology, 2016). These external soil moisture analyses are downscaled to the BARRA-R grid using a simple method that takes into account differences in soil texture. In each 6-hourly cycle, a land-surface analysis is conducted within BARRA (see Sect. 2.2). The daily initialization was conducted with the purpose of avoiding spurious drift in the BARRA moisture fields and reducing the time needed to spin up from ERA-Interim initial conditions. However, as multiple parallel production streams are needed to produce the reanalysis (see Sect. 2.2), there is a discontinuity in soil moisture in the bottom two layers between successive production streams, although soil moisture in the top two layers becomes stable after 1 month of runs. A discontinuity occurring at the 2014–2015 changeover has recently been reported by BARRA data users. These impacts, particularly on forested regions where trees extract water from the deep soil layers, are under investigation.

2.1.3 Boundary conditions

The BARRA-R sequential data assimilation process is initialized using ERA-Interim analysis fields (see Sect. 2.2), after which the only relationship with ERA-Interim is solely through the lateral boundary conditions. Hourly lateral boundary conditions for BARRA-R are interpolated from ERA-Interim's 6-hourly analysis fields at 0.75× 0.75 resolution. The rim width of the boundary frame is 0.88.

The land boundary is provided by a land-surface analysis (Sect. 2.2). Daily sea surface temperature (SST) and sea ice (SIC) analysis at 0.05× 0.05 resolution from reprocessed (1985–2007; Roberts-Jones et al., 2012) and near-real-time (NRT) Operational Sea Surface Temperature and Ice Analysis (OSTIA; Donlon et al., 2012) are used as lower boundaries over water after being interpolated to the UM grid. The NRT data are used from January 2007. OSTIA is widely used by NWP centres and operational ocean forecasting systems owing to its short real-time latency. Even though the re-processed and NRT data do not constitute a homogeneous time series, OSTIA is favoured over other SST reanalyses owing to its higher spatial resolution. Masunaga et al. (2015, 2018) have shown that steep SST gradients unresolved by coarse SST reanalyses can influence the organization of long-lived rain bands and the enhancement or reduction of surface convergence; this is particularly problematic for atmosphere-only reanalyses as thermal structure and motions in the marine atmospheric boundary layer are not well constrained by data assimilation.

2.2 Data assimilation system

The BARRA-R analysis scheme is based on fixed deterministic atmospheric and land-surface assimilation systems used by the UKMO in UERRA (Jermey and Renshaw, 2016) and IMDAA (Mahmood et al., 2018). BARRA-R uses a sequential data assimilation scheme, advancing forward in time using 6-hourly analysis cycles centred at synoptic hours t0= 00:00, 06:00, 12:00, and 18:00 UTC and 12 h forecast cycles from t0−3 h (Fig. 2).

Figure 2Cycling set-up of BARRA-R at base time t0= 00:00, 06:00, 12:00, and 18:00 UTC. Each UM forecast is initialized at t0−3 h by the previous forecast (grey arrows) with increments from the current analysis (red arrows). The purple bars indicate the time steps of the model states that have been archived.


In each analysis cycle, available observations, distributed across a 6h analysis window t0−3 h t<t0+3 h, are combined with the prior information of the model forecast from the previous cycle (the background state) to provide a more accurate estimate of the atmosphere over this window. This first involves a four-dimensional variational (4D-Var) analysis of the basic upper-air atmospheric fields (wind, temperature, specific humidity, pressure) with conventional and satellite observations (see below). The 4D-Var is favoured over 3D-Var as it takes account of time tendency information in the observations and this has a positive impact on the resulting forecasts (Rawlins et al., 2007). The UKMO's VAR assimilation system (version 2016.03.0) is used. The 4D-Var uses a linear perturbation forecast (PF) model (Lorenc, 2003; Rawlins et al., 2007; Lorenc and Payne, 2007), which uses a simpler model state linearized about a “guess” trajectory (i.e. tangent linear model) with a lower resolution (0.33 cf. 0.11) than the full forecast model. The lower resolution is chosen to limit computational costs. The PF model uses a simplified set of physical parameterizations including a simple boundary layer, cloud latent heat release, large-scale precipitation, and convection. In other words, it is assumed that the lower-resolution corrections to the background state (i.e. increments), interpolated to a higher resolution, are suitable corrections for the full model. The analysis increments from 4D-Var valid at t0−3 h are added to the background state at t0−3 h to produce an improved initial condition for the forecast model to perform the next 12 h forecast from t0−3 h to t0+9 h. A constraint of zero analysis increments is specified at the model boundary such that BARRA-R relies on the driving model ERA-Interim to define large-scale flow and other atmospheric conditions (Sect. 2.1.3). The observation departure statistics of the analysis, which are differences between the analysis and observations, are shown to be less than those of the model background (Supplement, Table S1). The assimilation is therefore behaving as desired by drawing the model towards observations for nearly all observational types.

The variational method minimizes a cost function whose two principal terms penalize distance to the background state and distance to the observations. The two terms are squared differences weighted by the inverse of their corresponding error covariances. In BARRA-R, the background error covariance has been estimated by a smooth parameterized approximation to climatology tuned by forecast differences (Ingleby, 2001). Accordingly, the estimated background error covariance is invariant between successive analysis windows but is time varying within the analysis window. The cost function also includes a pressure-based energy norm that serves as a weak-constraint digital filter to suppress spurious fast oscillations associated with gravity-inertia waves produced in model forecasts when analysis increments are added to the background state (Gauthier and Thépaut, 2001).

The initial land-surface state can have a significant impact on short-term forecasts of screen-level temperature and humidity, and its quality can also be improved through data assimilation. An extended Kalman filter (EKF) using observations of 2 m temperature and humidity is used to analyse the BARRA land state at every 6 h cycle and provide analyses of soil moisture, soil temperature, and skin temperature as described by Dharssi et al. (2012). The assimilation of satellite-retrieved soil moisture is not attempted here as it has not been realized in ACCESS. The UKMO's SURF analysis system (version 2016.07.0) is used to perform EKF. The Jacobian, which relates observed variables to model variables, for the Kalman gain matrix is estimated using finite difference by perturbing each model variable to be analysed in 40 perturbations and performing short 3 h forecasts. Here JULES (version 3.0) is run in the stand-alone mode, decoupled from the UM. The BARRA-R land state is reconfigured with EKF-derived surface analyses at every t0.

Note that the last 6 h forecast of a model run represents the prior state estimates needed for the next analysis cycle. The forecast fields valid at t0−3 h, t0−2 h, and t0−1 h are discarded, as these fields may still be influenced by transient artefacts due to the slight imbalance introduced by the addition of the analysis increments. It is already noted that this effect is also mitigated with the energy norm in the 4D-Var's cost function that penalizes the unbalanced structure in the increments.

The reanalysis is produced with multiple parallel production streams to speed up production. Each stream has a month of spin-up time from the ERA-Interim initial conditions before production data are archived, with most streams producing 1 year of reanalyses. Trials have shown that a 1-month period is a sufficient spin-up for the atmosphere (Renshaw et al., 2013) and top levels of soil moisture but insufficient for soil moisture in the deeper layers.

2.3 Observations

Conventional observations from land-surface stations, ships, drifting buoys, aircrafts, radiosondes, wind profilers, and satellite observations, namely retrieved wind, radiances, and bending angle, are assimilated in BARRA-R. The various observational types are chosen as they have been assimilated in the Bureau's operational NWP systems; other observational types, such as clear-sky radiances, have not been assimilated due to resource constraints. Rain observations from radar and gauges are also not assimilated as their assimilation schemes are still being tested for operational NWP. As listed in Table 1, the datasets are pragmatically taken from multiple sources, as they are being prepared during the production runs. Most of the observations prior to 2009 are supplied by the ECMWF, and the satellite radiance data from 2017 and onwards are extracted from the UKMO operational archive.

Table 1Observations assimilated in BARRA. Only the period concurrent with the reanalysis period is used. The various datasets were retrieved during the production, and thus the exact periods of each set used may differ.

Download Print Version | Download XLSX

The Bureau's archived observational data are also used to support this work, especially for the cycles from 2010 onwards. BARRA-R also assimilates additional high-frequency (10 min) land-surface observations from automatic weather stations in Australia and locally derived satellite atmospheric motion vectors (AMV). Ground positioning system (GPS) radio occultation bending angle data up to 2009 are provided by the Radio Occultation Meteorology Satellite Application Facility (ROM SAF). Additional land-surface observations over New Zealand are extracted from their National Climate Database (CliFlo, 2017). The 4D-Var assimilation of local AMV (Le Marshall et al., 2013) and GPSRO (Le Marshall et al., 2010) has been shown to improve operational forecasts.

Before being assimilated, observations are screened to select the best-quality observations, remove duplicates, and reduce data redundancy via thinning using the UKMO's Observing Processing System (OPS; based on version 2016.03.0) (Rawlins et al., 2007). There are per-cycle quality controls performed based on the method of Lorenc and Hammon (1988). Observations significantly different from the model background are rejected when exceeding a threshold calculated by a Bayesian scheme unless they are consistent with other observations nearby. The observational error variances and thinning distances are established at the UKMO and the Bureau for their NWP systems. For the surface, sonde, and aircraft observations, an observation automatic monitoring system performs monthly blacklisting of sites that show consistently large differences with BARRA-R's forecast over a 1-month period. The system also calculates bias corrections for surface pressure and for aircraft and sonde temperature.

For the satellite data, instruments and their individual channels are rejected when they become unreliable. The blacklisting is informed by the work of the ECMWF and MERRA-2 reanalysis teams. Further, air-mass-dependent variational bias correction is applied to satellite radiances as part of the assimilation process, allowing the time-varying corrections to fit drifts in instrumental bias (Harris and Kelly, 2001; Dee and Uppala, 2009). The bias corrections were calculated monthly, with the satellite radiances during the first month of each production stream not assimilated. There are abrupt changes to the amount of satellite data assimilated at the start and end of satellite missions and the various observational data archives. In some cases, changes occur when corrections were made to the observation screening and thinning rules mid-production for the 2010–2015 reanalyses. The impacts of such changes, known to cause artificial shifts and spurious trends in a reanalysis (e.g. Thorne and Vose, 2010; Dee et al., 2011), are still to be investigated for BARRA-R.

3 Preliminary evaluation

Our evaluation focuses on three areas: surface variables, pressure-level temperature and wind, and precipitation. For the surface variables, we compare BARRA-R against point-scale observations and gridded analyses of observations for 2 m temperature. For the pressure levels, we evaluate BARRA-R against point-scale observations of temperature and wind and examine the time series of the bias between BARRA-R and the global reanalyses. Finally, as rain observations are not assimilated in BARRA-R, gridded analyses of rain observations from gauges and satellites are used to provide the best independent reference in this study.

3.1 Surface

3.1.1 Point-scale evaluation of 2 m temperature, 10 m wind speed, and surface pressure

The t0+6 h model forecasts of 2 m (screen) temperature, 10 m wind speed, and surface pressure are evaluated against land observations. These observations have only an indirect relation to the forecasts as they are not used in the analysis for the associated cycle t0. Since errors tend to grow with the forecast range, the assessment places an upper bound on the true errors of the analysis fields between time t0 and t0+3 h. These fields are interpolated from the model levels using surface similarity theory (Walters et al., 2017a). The ERA-Interim t0+6 h forecasts from 00:00 and 12:00 UTC and the MERRA-2 hourly time-averaged forecast fields (M2T1NXSLV) are also evaluated to serve as benchmarks. It is not ideal to directly compare reanalyses with different resolutions, and interpolating them onto common (observed) locations diminishes some of the improvement achieved by BARRA-R relative to coarser reanalyses. Nonetheless, we undertake the latter to assess whether the models contain finer-scale information captured by point measurements; it therefore does not provide an assessment of the true quality of the reanalyses at their native resolutions.

To correct representativity errors in both reanalyses, their model values at (modelled) land grid cells are interpolated to the observation times and the station locations via bilinear interpolation in time and in the horizontal direction. Height corrections are applied to the interpolated fields to match the station heights: the corrections to the screen temperature are based on dry adiabatic lapse rate (Sheridan et al., 2010), 10 m wind speed is based on Howard and Clark (2007), and the correction to surface pressure is based on the hydrostatic equation under a constant lapse rate. As the observations are irregularly distributed in time, we consider all observations within a t0+5 h to t0+7 h time window, with t0 being 00:00 and 12:00 UTC, and the model grids are linearly interpolated to the observation times. Root mean square difference (RMSD), Pearson's linear correlation, additive bias, and variance bias are calculated at each station, with bias=meandm-mean(do) and the variance bias as Mbias=vardm/var(do)-1, to capture differences in the dispersion, where var(∗) computes the variance in time.

Box plots in Fig. 3 show the distribution of scores across 900–1500 stations in the BARRA-R domain. BARRA-R shows better agreement with the point observations than the global reanalyses for all three surface variables by most of the measures. This result is expected since BARRA-R resolves near-surface features below a 50 km horizontal scale and assimilates more surface observations over Australia and New Zealand. In particular, BARRA-R shows lower RMSD at about 80 % of the stations for screen temperature and 10 m wind speed and at 70 % of stations for surface pressure (see Fig. S1). At closer inspection in Fig. 4a, a percentile comparison plot of screen temperature deviation from monthly mean indicates that the frequency distribution of BARRA-R temperature is closer to that of the observations than ERA-Interim, particularly in regimes below percentile 25 and above percentile 90.

Figure 3Box plots showing the distribution of ERA-Interim, MERRA-2, and BARRA-R evaluation scores for (a) 2 m temperature, (b) 10 m wind speed, and (c) surface pressure over all stations in the BARRA-R domain. The scores are calculated on model forecasts valid between t0+5 h and t0+7 h against observations during 2007–2016. Individual boxes show the interquartile range of the scores, medians are marked in each box, and “whiskers” cover the 5th–95th percentile range.


Figure 4Comparisons of percentile values between observations and reanalyses for (a) 2 m temperature, and (b) 10 m wind speed during 2010–2013. The values for percentiles 0.05 to 99.95 are calculated using values derived from monthly means. The vertical blue dashed lines indicate the corresponding percentiles of the observations.


For 10 m wind speed, there are negative biases for variance in all the reanalyses assessed in this paper. Figure 4b shows that 10 m wind speeds are positively biased during light wind conditions and vice versa during strong wind speeds. There are many possible reasons for underestimating strong winds: inaccurate descriptions of boundary layer mixing, form drag for sub-grid orography, and surface properties such as land cover and vegetation types. Changing the fractional area of the vegetation canopy modifies the scalar roughness of the vegetated tiles, affecting the wind speed. The seemingly linear variation in wind speed is known in the global reanalyses (e.g. Carvalho et al., 2014), and Rose and Apt (2016) attributed the problem of wind underestimation to inaccuracy in modelling wind speeds in unstable atmospheric conditions.

Pressure is a large-scale variable which is likely to be better represented by a global model than a limited-area model. However, the BARRA-R estimates of point-scale surface pressure are more accurate in topographically complex regions and coastlines (see Fig. S1), where the estimates from the coarser reanalyses are less representative.

3.1.2 Comparison with gridded analysis of observed 2 m temperature

The reanalyses are compared against a gridded daily 0.05× 0.05 analysis of station maximum and minimum 2 m temperature data from the Australian Water Availability Project (AWAP; Jones et al., 2009). The AWAP grids are generated using an optimized Barnes successive-correction method that applies weighted averaging to the station data. Topographical information is included by using anomalies from long-term (monthly) averages in the analysis process. The AWAP analysis errors for maximum temperature are larger near the coast around northwest Australia and around the Nullarbor Plain due to strong temperature gradients between the coast and inland deserts and a relatively sparse network (Jones et al., 2009). The coast of western Australia and parts of Northern Territory are likely to share this analysis issue. The analysis errors are larger for minimum temperature, especially over western Australia and the Nullarbor Plain.

Figure 5 shows the differences for 2007–2016 averages in daily maximum and minimum temperature from AWAP, ERA-Interim, MERRA-2, and BARRA-R. The daily statistics are derived from 3-hourly forecast fields of ERA-Interim and hourly fields of MERRA-2 and BARRA-R. While inherent biases due to sampling are expected, this comparison also highlights the advantage of higher-frequency data when examining lower and upper tail statistics. BARRA-R shows cold and warm biases (relative to AWAP) of around 1 K in daily maximum and minimum temperature, respectively, particularly over the eastern region. MERRA-2 also shows similar levels of biases but with different signs and variability. BARRA-R and MERRA-2 agree better with AWAP than ERA-Interim, which reports differences (in mean) up to 5 K in magnitude. The reduced amplitude of the diurnal cycle of temperature is a long-standing problem in the UM; experiments have shown that changes to the representation of the land surface (e.g. reductions in the amount of bare soil and changes to the scalar roughness and albedo of vegetated tiles) reduce clear-sky biases (Bush et al., 2019).

Figure 5Mean differences in (row i) daily maximum (TMax) and (ii) minimum (TMin) 2 m temperature (K) for 2007–2016 between (a) ERA-Interim and AWAP, (b) MERRA-2 and AWAP, and (c) BARRA-R and AWAP. The spatial means of the differences are reported in the text.


Figure 6 shows the monthly means of the differences in daily maximum and minimum temperature between the reanalyses and AWAP averaged across Australia. Here the OSTIA SST anomaly time series is also included, and it does not show a visible discontinuity at 2006–2007 (Sect. 2.1.3). The maximum temperature in BARRA-R appears cooler than AWAP after a strong La Niña event in 2010–2011, while the global reanalyses also show cooler trends in biases after 2010. BARRA-R and ERA-Interim show smaller levels of temporal variability than MERRA-2. The minimum temperature in BARRA-R does not show an obvious trend but is warmer during 2010–2011 when ERA-Interim and MERRA-2 are cooler. These changes do not coincide with the change in soil moisture initialization in 2014–2015 (Sect. 2.1.2) or OSTIA SST.

Figure 6Monthly mean differences in daily (a) maximum (TMax) and (b) minimum (TMin) 2 m temperature (K) averaged over Australia between (row i) BARRA-R and AWAP, (ii) ERA-Interim and AWAP, and (iii) MERRA-2 and AWAP. Black curves are shaded around the 14-year means. Green curves plot the monthly anomalies, from 2003–2016 monthly averages, of the OSTIA sea surface temperature averaged over 46–4 S and 94–174 E.


3.2 Pressure levels

To assess BARRA-R in the atmosphere, we compare the t0+6 h forecasts on pressure levels with radiosonde and pilot wind observations at 00:00 and 12:00 UTC on standard pressure levels ranging from 1000 to 10 hPa, using the harmonized dataset produced by Ramella Pralungo et al. (2014a, b). The pressure-level fields of BARRA-R and ERA-Interim analyses at time t0 are also compared, even though they are not independent from the observations; such comparisons only provide baselines to interpret the relative quality of the BARRA-R forecasts. Similar comparisons with ERA-Interim's twice-daily forecasts at these observation times are also not possible because they start from 00:00 and 12:00 UTC. The model data are interpolated horizontally to the sonde and pilot launch locations via bilinear interpolation, and RMSD is calculated at each location and pressure level. The resulting box plots of RMSD are shown in Fig. 7. Depending on the pressure level and parameter evaluated, between 54 and 203 sites were available. There is a marked variability in RMSD with the pressure levels, particularly for wind speed, due to factors such as variations in the number of observing sites, increasing sonde drift error on ascent, and differences in the dynamic range of the fields with height. A markedly higher RMSD for wind speed occurs at 200 hPa, a height a which the jet stream can be located.

Figure 7Box plots showing the RMSD distribution of BARRA-R t0+6 forecast and t0 analysis, as well as ERA-Interim analysis for (a) temperature and (b) wind speed at multiple sites in the BARRA-R domain. RMSD is calculated for temperature and wind speed at pressure levels of 10, 50, 100, 200, 400, 500, 700, 850, 925, and 1000 hPa against pilot balloon and radiosonde observations at 00:00 and 12:00 UTC. The numbers of sites are indicated in brackets.


It is difficult to discern the differences between the two analyses, suggesting that they perform similarly from assimilating the same observations. Assimilation at a coarser resolution of 0.33 (cf. 0.11 of the forecast model) in BARRA-R does not drastically improve 0.75 representations of temperature and wind at these pressure levels and at point scales. There are also small differences between the analyses and BARRA-R background, indicating that the 0.11 forecast model does not degrade from the lower-resolution analysis of BARRA-R but also does not improve upon the ERA-Interim 0.75 representation of these fields at the observation locations.

Figure 8 compares BARRA-R's 00:00 UTC analysis of air temperature at 850, 700, and 500 hPa against the analyses from ERA-Interim and MERRA-2 (M2I3NPASM). BARRA-R is cooler at 500 hPa across the domain and warmer at 850 hPa in the tropics than the global reanalyses, and the monthly differences in the zonal mean are of the order of 1 K. BARRA-R also shows a cooling shift at 700 and 500 hPa in the tropics and a warming shift south of 40 S after 2010. But when compared against MERRA-2, in the tropics, BARRA-R is warmer at 700 hPa, and the apparent shift in BARRA-R is also seen in MERRA-2 (relative to ERA-Interim) at these levels.

Figure 8Hovmöller plots of the monthly difference in zonal mean air temperature (K) at 00:00 UTC and three pressure levels, (row i) 850, (ii) 700, and (iii) 500 hPa, between (a) BARRA-R and ERA-Interim and between (b) BARRA-R and MERRA-2.


3.3 Precipitation

We consider three reference gridded datasets to compare with the reanalyses. First is the 0.05× 0.05 rain gauge analysis of daily accumulation over Australia from AWAP, produced using the Barnes method in which the ratio of observed rainfall to monthly average is used in the analysis process (Jones et al., 2009). There is a north–south gradient in the AWAP analysis errors with larger analysis errors in the northern tropical regions, where the length scales of convective rainfall events are shorter and more variable (Jones et al., 2009). Second is the 1× 1 (full data daily) rain gauge analysis over the domain from the Global Precipitation Climatology Centre (GPCC version 2018; Ziese et al., 2018), created using an empirical weighting-based interpolation method described in Becker et al. (2013). As with AWAP, GPCC is less accurate in regions where station scarcity and high precipitation variability coexist. For instance, different GPCC interpolation methods can yield very different analyses over the south Asia region (Becker et al., 2013). The third reference is the 0.25× 0.25 satellite-based analysis of 3-hourly rain rates from the Tropical Rainfall Measuring Mission (TRMM) multi-satellite precipitation analysis (TMPA 3B42 version 7; Huffman et al., 2006). TMPA 3B42 combines precipitation estimates from various satellite systems and rain gauge monthly analysis. Satellite-derived estimates of convective precipitation are largely accurate in the low latitudes (Ebert et al., 2007; Chen et al., 2013), but the TMPA product is less accurate over the ocean due to the absence of local observations used for gauge adjustments (Sapiano and Arkin, 2009) and south of 40 S due to limited local cross-sensor calibration (Huffman et al., 2006). TRMM often underestimates precipitation in high-latitude regions with significant topography due to difficulties of satellite retrievals over snow-covered surfaces and/or due to the high elevations (Barros et al., 2006; Matthews et al., 2013). TRMM is also known to underestimate light rainfall and drizzle over subtropical and high-latitude oceans (Berg et al., 2010). In addition to these considerations, there are inherent limitations in comparing the reanalyses with AWAP, GPCC, and TMPA. Specifically, products with coarser grids tend to over-represent low-threshold events occurring at spatial scales smaller than their grid sizes and under-represent high-threshold events. Further evaluation of BARRA-R precipitation estimates against point gauge observations and AWAP are reported in Acharya et al. (2019).

Neither BARRA-R nor ERA-Interim assimilated rainfall observations. Precipitation estimation from their forecast models is constrained by other observation types. Following Sect. 2.1, in BARRA-R, the microphysics scheme based on Wilson and Ballard (1999) parameterizes the atmospheric processes that transfer water between the four modelled states of water (vapour, liquid droplets, ice, and raindrops) to remove moisture resolved on the grid scale. As the 12 km model is not “storm resolving”, BARRA-R uses the mass flux convective parameterization scheme of Gregory and Rowntree (1990) with the CAPE closure to model sub-grid-scale precipitating and non-precipitating convection using an ensemble of cumulus clouds as a single entraining–detraining plume. Such a scheme prevents the unstable growth of cloudy structures on the grid, which is otherwise required for explicit vertical circulations to develop (Clark et al., 2016). The modelled convection also works independently at each grid point, and the model can only predict the area-average rainfall instead of the spectrum of rainfall rates. Consequently, BARRA-R's precipitation estimates from sub-grid convection will be more erroneous than those for large-scale precipitation. In other words, the accuracy of BARRA-R is expected to worsen during the warm season and at low latitudes and to improve during the cooler season and at high latitudes where non-convective precipitation is dominant. To allow the UM to spin up from the analysis increments, we examine the quality of the precipitation accumulation between t0+3 h and t0+9 h by comparing against gridded datasets. This also addresses the issue that the UM yields excess precipitation at analysis time (t0−3 h) due to a temporary imbalance in the moisture fields by allowing time for the model to adjust and remove the excess. For ERA-Interim, we used its first 12 h accumulation, which is considered the most accurate (Kallberg, 2011).

3.3.1 Mean annual precipitation and frequency of rain days

Figure 9i compares the 10-year (2007–2016) annual mean precipitation estimated from the five datasets. A close-up over Australia can be found in Fig. S2. BARRA-R provides a realistic depiction when compared with TMPA across the domain but shows higher precipitation over the tropics and over the Tasman Sea and Southern Ocean. BARRA-R agrees very well with AWAP and GPCC over Australian land areas, reflecting the markedly higher precipitation in the northern tropics, and western Tasmania. It also agrees with GPCC over New Zealand. BARRA-R also shows better agreement with AWAP, GPCC, and TMPA in some dry areas such as western Australia.

Figure 9(i) Mean annual precipitation (mm), (ii) fractions of light rain days with 1–10 mm of precipitation, (iii) heavy precipitation days with 10–50 mm, and (iv) very heavy precipitation days with >50 mm over 2007–2016 from (a) AWAP, (b) GPCC, (c) TMPA, (d) ERA-Interim, and (e) BARRA-R. Regions with more than 10 % missing values in AWAP are masked. Close-ups of the plots over Australia are provided in the Supplement (Fig. S2).


The frequency of days with three intensity regimes is examined next. In Fig. 9ii, we examine the frequency of light rain days with amounts between 1 and 10 mm. The 1 mm threshold is chosen to account for the tendency of the model to create light “drizzle” events with very low rain rates. Even so, the two reanalyses show significantly more rain days in the tropics than TMPA and GPCC and more rain days than TMPA over the Southern Ocean. TRMM is known to miss light rainfall events over subtropical and high-latitude oceans (Berg et al., 2010), while simulated precipitation over the Southern Ocean overestimates drizzle compared with satellite observations (Franklin et al., 2013; Wang et al., 2015). Some of these differences from TMPA are not mirrored by AWAP over Australia, suggesting a possible underestimation of rain days in TMPA over land (e.g. eastern seaboard, southwest Australia) where the gauge network is relatively dense. Despite these considerations, BARRA-R overestimates the frequency of light rain days compared with AWAP, notably in the northern and central regions of Australia and Tasmania. The UM's parameterized convection scheme assumes that there are many clouds per grid box, which is marginal at the BARRA-R resolution, and thus produces a bias towards widespread precipitation and provides little indication of the areas which could expect larger rain rates (Clark et al., 2016).

For heavy precipitation days, with amounts of 10 to 50 mm, Fig. 9iii shows greater similarities between BARRA-R, AWAP, and GPCC over land regions, such as the southeast coast of Australia and Tasmania, than for ERA-Interim. BARRA-R shows differences from AWAP and GPCC over Australia north of 30 S where the gauge analyses are poorer. Over the tropical ocean, the two reanalyses show more heavy precipitation days than TMPA.

Lastly, for the very heavy precipitation days (≥50 mm) in Fig. 9iv, it is obvious that ERA-Interim does not fully capture the frequency over land in northern Australia and Southeast Asia, whereas BARRA-R is more comparable with the three reference datasets. This agrees with the findings of Jermey and Renshaw (2016) that higher-resolution regional reanalyses show improvement in representing high-threshold events at these spatial scales. Over the ocean, BARRA-R also shows greater rainfall intensity in the tropics than ERA-Interim, but both reanalyses show lower intensity compared to TMPA. These results reflect the deficiency of the parameterized convection scheme in BARRA-R for estimating convective precipitation amounts in this region.

3.3.2 Comparison of monthly totals

Figures 10 and 11 compare differences in domain-averaged monthly totals between the reanalyses (BARRA-R and ERA-Interim) and reference data (TMPA and GPCC) over five separate sub-domains between 80 and 180 E. Precipitation over land and ocean are distinguished. Over the tropical ocean at ±10 N (Fig. 10i), the two reanalyses show different shifts in overall differences from TMPA at around 2010, and these shifts are not apparent in the other sub-domains. Across the sub-domains, the variances of the differences are similar between the two reanalyses.

Figure 10Differences in monthly precipitation total (mm) averaged over the ocean in five sub-domains (rows i–v) between (a) BARRA-R and TMPA and between (b) ERA-Interim and TMPA. Black curves are shaded around the 14-year means. Green curves plot the monthly anomalies, from 2003–2016 monthly averages, of the OSTIA sea surface temperature averaged over respective sub-domains.


Figure 11As with Fig. 10a and b, but over land. Additional comparisons are made between (c) BARRA-R and GPCC and between (d) ERA-Interim and GPCC.


Over tropical land regions, BARRA-R shows much higher totals than others (Fig. 11i) due to higher precipitation occurring in mountainous terrains in Papua New Guinea (PNG), Indonesia and Sumatra, and relatively small Indonesian islands (see Fig. S3). Other reanalyses and other gridded precipitation products disagree greatly at these locations with few observations and mountainous terrains (e.g. over PNG in Smith et al., 2013). BARRA-R (and GPCC) also shows markedly higher monthly totals below 39.2 S (Fig. 11v) than TMPA and ERA-Interim. This is due to higher BARRA-R precipitation estimates on the west coast and in the Southern Alps of New Zealand, where precipitation is likely underestimated in TMPA.

The UM can produce grid-localized high precipitation in BARRA-R, especially in unstable atmospheric conditions over steep orographic slopes. This issue is not unique to the UM but, for instance, also occurs in the Weather Research and Forecasting model (Gustafson Jr. et al., 2014). When the convective parameterization in non-convective-resolving models does not stabilize the air column, meteorological events can develop at the smallest resolvable scales in the model, producing unrealistically strong vertical velocities and precipitation; this is known as a “grid-point storm” (Scinocca and McFarlane, 2004; Williamson, 2013; Chan et al., 2014). In our cases, the model only produces isolated excessively intense rainfall over steep topography. Such storms occur more readily in models with higher horizontal resolutions (Williamson, 2013). As the resolution increases, resolved motions can produce moisture convergence and increase CAPE very rapidly, and the rate at which column instability is produced depends on the scale of moisture and heat convergence. This also tends to occur over tropical land areas, over steep topography, and during the warm seasons, when the atmosphere is unstable and there is a sufficient warm moisture supply at the surface. These considerations do not lend themselves to completely explain the observed bias in BARRA-R.

By contrast, BARRA-R shows better agreement with GPCC and TMPA in other sub-domains between 39.2 and 10.0 S (Fig. 11ii–iv). Over the land between 23 and 10 S, BARRA-R simulates wetter summer events than observed in TMPA and GPCC from 2011, when Australia was recovering from drought conditions with the onset of La Niña. Between 39 and 23 S, BARRA-R also simulated wetter events over Mt. Kosciuszko, Tasmania, and the North Island of New Zealand than TMPA after 2014. This overestimation is, however, less apparent when BARRA-R is compared with GPCC.

4 Discussion and outlook

The recent development of global and regional reanalyses addresses the need for high-quality, increasingly higher-resolution, and longer-term reanalyses, accompanied by estimates of uncertainty, within the research and broader user communities. BARRA is the first regional reanalysis that focuses on the Australasian section of the Southern Hemisphere. It is developed with significant co-investment from state-level emergency service agencies across Australia due to the advantages of a deeper understanding of past weather, including extreme events, especially in areas that have been poorly served by observation networks. The 29-year BARRA reanalysis, which is expected to be completed in 2019, will ultimately represent a collection of high-resolution gridded meteorological datasets with 12 and 1.5 km horizontal resolution and 10 min to hourly time resolution.

In this paper, we describe the BARRA 12 km regional reanalysis, BARRA-R, which is closely related to the Bureau's regional NWP system but with an updated UM, 4D-Var, variational bias correction, and automated station blacklisting systems. BARRA-R covers a significant region of the globe, including parts of Southeast Asia, the eastern Indian Ocean, the southwest Pacific, Australia, and New Zealand, and assimilates a wide range of conventional and satellite observations that have proven to improve the skill of NWP.

BARRA-R produces a credible reproduction of the meteorology at and near the surface over land as diagnosed by the selected variables. BARRA-R improves upon its global driving model, ERA-Interim, showing better agreement with point-scale observations of 2 m temperature, 10 m wind speed, and surface pressure. Results are similar when BARRA-R is compared with MERRA-2. Daily maximum and minimum statistics for 2 m temperature at 5 km resolution are captured in BARRA-R with smaller biases than ERA-Interim. There appear to be shifts in biases relative to land observation analyses over Australia amongst all the reanalyses, mirroring changes in SST. This behaviour, however, does not coincide with known changes to the forcing data (soil moisture and SST) used in BARRA-R and requires further analysis to be better understood. BARRA-R's 10 m wind fields show lower biases than ERA-Interim and MERRA-2, but the negative bias during strong winds, which is common amongst other reanalyses, remains significant. Altogether, BARRA-R provides good representation of near-surface extremes, which has implications for its uses for energy management, fire risk, and storm damage. The bias could be addressed via post-processing using methods such as those of Glahn and Lowry (1972) and Rose and Apt (2016). More generally, a variety of post-processing methods can further improve the accuracy of BARRA-R data (e.g. Berg et al., 2012; Frank et al., 2018). Our study did not discern clear merits in BARRA-R analysis and forecast relative to ERA-Interim analysis for pressure-level temperature and wind. Further, there is no conclusive explanation for the shifts in 500, 700, and 850 hPa air temperature occurring at 2010, as comparisons with ERA-Interim and MERRA-2 yield mixed results. Other evaluations of the UM GA6 configuration, including tropical cyclones, precipitation, clouds, and large-scale flow, are reported in Walters et al. (2017a, b) but in global models at coarser spatial resolutions.

Precipitation fields from BARRA-R show similarities with AWAP and GPCC rain gauge analyses over Australia, where it reflects more similar frequency statistics for heavy rain events and annual mean than ERA-Interim. While this is expected from comparing grids with different resolutions, BARRA-R contains more information pertaining to rain events at local scales. The frequency statistics (of both light and heavy rain days) of the two reanalyses are markedly different from TMPA over regions exterior to Australia. BARRA-R is likely to be positively biased over land in the regions north of 10 S and New Zealand due to higher precipitation estimates concentrated in regions with high or steep topography. This is partly due to the presence of grid-point storms that occur in non-convective-resolving models. Alas, the likely underestimation in observations associated with high elevations poses difficulties to quantify the wet bias. The characteristics of grid-point storms in terms of superficial spatial localization, precipitation amount, and vertical wind speed could be detected and screened out via post-processing. This is important as this model artefact affects the analyses of rainfall averages and extremes.

The disagreement with TMPA is also apparent over the oceans, but consensus between satellite-based products generally degrades over higher latitudes, especially over the Southern Ocean (Behrangi et al., 2014). Over the 2003–2016 period, the variability of the monthly precipitation totals is similar amongst the reanalyses, TMPA, and GPCC across the domain. Notable exceptions are a dry shift occurring in BARRA-R during 2010 over the tropical ocean and wetter summer events over land in northern and southeast Australia, as well as the North Island of New Zealand, after 2014. These coincident shifts in daily maximum 2 m temperature (over Australia), upper-air temperature (across the BARRA-R domain), and tropical precipitation in all the reanalyses suggest larger differences in large-scale synoptic patterns between them after 2010. Given all the above considerations, local evaluation of BARRA-R reanalysis before application is recommended.

Higher-resolution models used to downscale BARRA-R could alleviate the observed shortcomings by resolving sharp topographical features, resolving sub-grid processes (e.g. convection), and using science configurations more suited for a given climatic region. Assessment of the UM's first Regional Atmosphere (RA1) science configurations for convective-permitting models, recently concluded in December 2017, distinguishes two different science configurations for mid-latitude and tropical regions (RA1-M and RA1-T, respectively). Developments in RA1 have produced improvements to 2 m temperature, 10 m wind speed, and precipitation (Bush et al., 2019). Further, it is known that BARRA-R's convection scheme, involving instantaneous adjustment of cloud fields to changes in forcing (e.g. solar heating, land–sea temperature differences), can lead to unrealistic behaviour at places such as coasts and in time (e.g. incorrect diurnal cycle) (Lean et al., 2008; Clark et al., 2016). A companion article will examine the merits of downscaling BARRA-R with convective-scale models.

Finally, BARRA represents an important step in supporting the Bureau's ability to prepare for future reanalysis-related activities such as data rescue and reprocessing of observational data. Future reanalyses could use higher-resolution models and ensemble-based forecast and assimilation systems to quantify uncertainties. They will also benefit from international efforts in reprocessing historical conventional and satellite observations with enhanced quality and/or more accurate uncertainty estimates.

Code availability

All code, including the UM (version 10.2), VAR (version 2016.03.0), JULES (version 3.0), OPS (version 2016.03.0), and SURF (version 2016.07.0) systems, used to produce BARRA is version-controlled under the Met Office Science Repository Service. Readers are referred to (last access: 23 May 2019; UK Met Office, 2019) for access information.

Data availability

The first releases of the BARRA-R dataset for the period 2003–2016 are available for academic use, with subsequent releases planned for mid-2019. Readers are referred to (last access: 23 May 2019; Bureau of Meteorology, 2019) for information on available parameters and access.


The supplement related to this article is available online at:

Author contributions

PS, DJ, PFH, and CJW conceived and/or designed BARRA. CHS, NE, and PS developed the BARRA system with inputs from SR, CF, ID, and HZ. CHS and NE performed the production and evaluation. CHS prepared the paper with contributions from all co-authors.

Competing interests

The authors declare that they have no conflict of interest.


Funding for this work was provided by emergency service agencies (New South Wales Rural Fire Service, Western Australia Department of Fire and Emergency Services, South Australia Country Fire Service, South Australia Department of Environment, Water and National Resources) and research institutions (Antarctic Climate and Ecosystems Cooperative Research Centre (ACE CRC) and University of Tasmania). Funding from Tasmania is supported by the Tasmanian Government and the Australian Government, provided under the Tasmanian Bushfire Mitigation Grants Program.

BARRA-R is set up with assistance from the UKMO reanalysis team (Richard Renshaw, Peter Jermey, Jemma Davie) and colleagues (Adam Maycock, David Walters, Ian Boutle), as well as many colleagues at the Bureau of Meteorology (Tan Le, Ilia Bermous, Lawrie Rikus, Fiona Smith, Charles Sanders, Jin Lee, Gary Dietachmayer, John Le Marshall, Xudong Sun, Jim Fraser, Greg Kociuba, Chris Tingwell, Huqiang Zhang), the Commonwealth Scientific and Industrial Research Organisation (CSIRO; Martin Dix), and National Computational Infrastructure (NCI; Dale Roberts). We thank the following people: Richard Renshaw for providing the observational data from the UKMO and ECMWF archives; Stuart Moore and Trevor Carey-Smith at the National Institute of Water and Atmospheric Research (NIWA) for providing additional local observations over New Zealand; Robert Smalley and David Jones for their advice on AWAP; and Peter May, Elizabeth Ebert, Andrew Dowdy, and Tony Hirst for their feedback on early drafts of paper. BARRA-R uses the ERA-Interim data provided through the ARC Centre of Excellence for Climate System Science (Paola Petrelli) at NCI. Many of the observational datasets were provided by the ECMWF, UKMO, and NIWA. The radio occultation data (CDR v1.0 beta release) were provided by the Radio Occultation Meteorology Satellite Application Facility (ROM SAF; through K. B. Lauritsen and H. Gleisner, Danish Meteorological Institute), which is a decentralized operational RO processing centre under EUMETSAT. RO data are available at (last access: 23 May 2019). The BARRA project was undertaken with the assistance of resources and services from NCI, which is supported by the Australian Government.

ERA-Interim can be retrieved from the ECMWF at AWAP data can be requested from (last access: 23 May 2019), TMPA v7 data are retrieved via the NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC) at (last access: 23 May 2019), and GPCC v2018 data are retrieved from the Deutscher Wetterdienst.

Review statement

This paper was edited by Jatin Kala and reviewed by four anonymous referees.


Acharya, S. C., Nathan, R., Wang, Q. J., Su, C.-H., and Eizenberg, N.: An evaluation of daily precipitation from atmospheric reanalyses over Australia, Hydrol. Earth Syst. Sci. Discuss.,, in review, 2019. 

Arakawa, A. and Lamb, V. R.: Computational design of the basic dynamical processes of the UCLA general circulation model, Methods of Comp. Phys.: Adv. Res. Appl., 17, 173–265,, 1977. 

Barros, A. P., Chiao, S., Lang, T. J., Burbank, D., and Putkonen, J.: From weather to climate – Seasonal and interannual variability of storms and implications for erosion process in the Himalaya, Geol. Soc. Am. Spat. Paper 398, Penrose Conference Series, 17–38, 2006. 

Becker, A., Finger, P., Meyer-Christoffer, A., Rudolf, B., Schamm, K., Schneider, U., and Ziese, M.: A description of the global land-surface precipitation data products of the Global Precipitation Climatology Centre with sample applications including centennial (trend) analysis from 1901–present, Earth Syst. Sci. Data, 5, 71–99,, 2013. 

Behrangi, A., Stephens, G., Adler, R. F., Huffman, G. J., Lambrigtsen, B., and Lebsock, M.: An update on the oceanic precipitation rate and its zonal distribution in light of advanced observations from space, J. Climate, 27, 3957–3965,, 2014. 

Berg, W., L'Ecuyer, T., and Haynes, J. M.: The distribution of rainfall over oceans from spaceborne radars, J. Appl. Meteor. Climatol., 49, 535–543,, 2010. 

Berg, P., Feldmann, H., and Panitz, H.-J.: Bias correction of high resolution regional climate model data, J. Hydrol., 448–449, 80–92,, 2012. 

Best, M. J., Pryor, M., Clark, D. B., Rooney, G. G., Essery, R. L. H., Ménard, C. B., Edwards, J. M., Hendry, M. A., Porson, A., Gedney, N., Mercado, L. M., Sitch, S., Blyth, E., Boucher, O., Cox, P. M., Grimmond, C. S. B., and Harding, R. J.: The Joint UK Land Environment Simulator (JULES), model description – Part 1: Energy and water fluxes, Geosci. Model Dev., 4, 677–699,, 2011. 

Bollmeyer, C., Keller, J. D., Ohlwein, C., Wahl, S., Crewell, S., Friederichs, P., Hense, A., Keune, J., Kneifel, S., Pscheidt, I., Redl, S., and Steinke, S.: Towards a high-resolution regional reanalysis for the European CORDEX domain, Q. J. Roy. Meteorol. Soc., 141, 1–15,, 2015. 

Borsche, M., Kaiser-Weiss, A. K., Unden, P., and Kaspar, F.: Methodologies to characterize uncertainties in regional reanalyses, Adv. Sci. Res., 12, 207–218,, 2015. 

Borsche, M., Kaiser-Weiss, A. K., and Kaspar, F.: Wind speed variability between 10 and 116 m height from the regional reanalysis COSMO-REA6 compared to wind mast measurements over Northern Germany and the Netherlands, Adv. Sci. Res., 13, 151–161,, 2016. 

Bromwich, D., Wilson, A., Bai, L., Liu, Z., Barlage, M., Shih, C., Maldonado, S., Hines, K., Wang, S.-H., Woollen, J., Kuo, B., Lin, H., Wee, T., Serreze, M., and Walsh, J.: The Arctic System Reanalysis Version 2, B. Am. Meteorol. Soc., 99, 805–828,, 2018. 

Brown, A., Milton, S., Cullen, M., Golding, B., Mitchell, J., and Shelly, A.: Unified modeling and prediction of weather and climate: A 25-Year journey, B. Am. Meteorol. Soc., 93, 1865–1877,, 2012. 

Bureau of Meteorology: Operational implementation of the ACCESS numerical weather prediction systems, NMOC Op. Bull. No. 83, available at: (last access: 17 May 2019), 2010. 

Bureau of Meteorology: APS1 upgrade of the ACCESS-R numerical weather prediction system, NMOC Op. Bull. No. 98, available at: (last access: 17 May 2019), 2013. 

Bureau of Meteorology: APS2 upgrade to the ACCESS-G numerical weather prediction system, BNOC Op. Bull. No. 105, available at: (last access: 17 May 2019), 2016. 

Bureau of Meteorology: Atmospheric high-resolution regional reanalysis for Australia, available at:, last access: 23 May 2019. 

Bush, M., Allen, T., Bain, C., Boutle, I., Edwards, J., Finnenkoetter, A., Franklin, F., Hanley, K., Lean, H., Lock, A., Manners, J., Mittermaier, M., Morcrette, C., North, R., Petch, J., Short, C., Vosper, S., Walters, D., Webster, S., Weeks, M., Wilkinson, J., Wood, N., and Zerroukat, M.: The Met Office Unified Model/JULES Regional Atmosphere and Land configurations (RAL) – 1st release, submitted, 2019 

Carvalho, D., Rocha, A., Gomez-Gesteira, M., and Santos, C. S.: WRF wind simulation and wind energy production estimates forced by different reanalyses: comparison with observed data for Portugal, Appl. Energy, 117, 116–126,, 2014. 

Chan, S. C., Kendon, E. J., Fowler, H. J., Blenkinsop, S., Roberts, N. M., and Ferro, C. A.: The value of high-resolution Met Office regional climate models in the simulation of multihourly precipitation extremes, J. Climate, 27, 6155–6174,, 2014. 

Charney, J. G. and Phillips, N. A.: Numerical integration of the quasi-geostrophic equations for barotropic and simple baroclinic flows, J. Meteorol., 10, 71–99,<0071:NIOTQG>2.0.CO;2, 1953. 

Chen, Y., Ebert, E. E., Walsh, K. J. E., and Davidson, N. E.: Evaluation of TRMM 3B42 precipitation estimates of tropical cyclone rainfall using PACRAIN data, J. Geophys. Res.-Atmos., 118, 2184–2196,, 2013. 

Clark, P., Roberts, N., Lean, H., Ballard, S. P., and Charlton-Perez, C.: Review: Convection-permitting models: a step-change in rainfall forecasting, Meteor. App., 23, 165–181,, 2016. 

CliFlo: NIWA's National Climate Database on the Web, available at: last access: 17 February 2017. 

Davies, T., Cullen, M. J. P., Malcolm, A. J., Mawson, M. H., Staniforth, A., White, A. A., and Wood, N.: A new dynamical core for the Met Office's global and regional modelling of the atmosphere, Q. J. Roy. Meteorol. Soc., 131, 1759–1782,, 2005. 

Dee, D. P. and Uppala, S.: Variational bias correction of satellite radiance data in the ERA-Interim reanalysis, Q. J. Roy. Meteorol. Soc., 135, 1830–1841,, 2009. 

Dee, D. P., Källén, E., Simmons, A. J., and Haimberger, L.: Comments on “Reanalyses suitable for characterizing long-term trends.”, B. Am. Meteorol. Soc., 92, 65–70,, 2011. 

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol. C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Holm, E. V., Isaksen, L., Kallberg, P., Kohler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J. J., Park, B. K., Peubey, C., de Rosnay, P., Tavolato, C., Thepaut, J. N., and Vitart, F.: The Era-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597,, 2011. 

Dee, D. P., Balmaseda, M., Balsamo, G., Engelen, R., Simmons, A. J., and Thepaut, J.-N.: Towards a consistent reanalysis of the climate system, B. Am. Meteorol. Soc., 95, 1235–1248,, 2014. 

Dharssi, I. and Vinodkumar, J.: A prototype high resolution soil moisture analysis system for Australia, Bureau of Meteorology Report No. 026, available at: (last access: 17 May 2019), 2017. 

Dharssi, I., Steinle, P., and Candy, B.: Towards a Kalman filter based land surface data assimilation scheme for ACCESS, Bureau of Meteorology CAWCR Technical Report No. 54, available at: (last access: 17 May 2019), 2012. 

Dharssi, I., Steinle, P., and Fernon, J.: Improved numerical weather predictions by using optimised urban model parameter values and satellite derived tree heights, in: MODSIM2015, 21st International Congress on Modelling and Simulation, edited by: Weber, T., McPhee, M. J., and Anderssen, R. S., Modelling and Simulation Society of Australia and New Zealand, December 2015, 1161–1167, ISBN: 978-0-9872143-5-5, available at: (last access: 17 May 2019), 2015. 

Dickinson, R. E., Errico, R. M., Giorgi, F., and Bates, G. T.: A regional climate model for the western United States, Clim. Change, 15, 383–422,, 1989. 

Donlon, C. J., Martin, M., Stark, J. D., Roberts-Jones, J., Fiedler, E., and Wimmer, W.: The Operational Sea Surface Temperature and Sea Ice analysis (OSTIA) system, Remote Sens. Environ., 116, 140–158,, 2012. 

Ebert, E. E., Janowiak, J. E., and Kidd, C.: Comparison of near-real-time precipitation estimates from satellite observations and numerical models, B. Am. Meteorol. Soc., 88, 47–64,, 2007. 

Ebita, A., Kobayashi, S., Ota, Y., Moriya, M., Kumabe, R., Onogi, K., Harada, Y., Yasui, S., Miyaoka, K., Takahashi, K., Kamahori, H., Kobayashi, C., Endo, H., Soma, M., Oikawa, Y., and Ishimizu, T.: The Japanese 55-year reanalysis JRA-55: An interim report, SOLA, 7, 149–152,, 2011. 

Edwards, J. M. and Slingo, A.: Studies with a flexible new radiation code. I: Choosing a configuration for a large-scale model, Q. J. Roy. Meteorol. Soc., 122, 689–719,, 1996. 

Evans, J. P. and McCabe, M. F.: Effect of model resolution on a regional climate model simulation over southeast Australia, Clim. Res., 56, 131–145,, 2013. 

Fall, S., Niyogi, D., Gluhovsky, A., Pielke Sr., R. A., Kalnay, E., and Rochon, G.: Impacts of land use land cover on temperature trends over the continental United States: assessment using the North American Regional Reanalysis, Int. J. Climatol., 30, 1980–1993,, 2010. 

Fowler, H. J., Blenkinshop, S., and Tebaldi, C.: Linking climate change modelling to impacts studies: recent advances in downscaling techniques for hydrological modeling, Int. J. Climatol., 27, 1547–1578,, 2007. 

Frank, C. W., Wahl, S., Keller, J. D., Pospichal, B., Hense, A., and Crewell, S.: Bias correction of a novel European reanalysis data set for solar energy applications, Solar Ener., 164, 12–24,, 2018. 

Franklin, C. N., Sun, Z., Bi, D., Dix, M., Yan, H., and Bodas-Salcedo, A.: Evaluation of clouds in ACCESS using the satellite simulator package COSP: Global, seasonal, and regional cloud properties, J. Geophys. Res.-Atmos., 118, 732–748,, 2013. 

Gauthier, P. and Thépaut, J.-N.: Impact of the digital filter as a weak constraint in the preoperational 4DVar assimilation system of Météo-France, Mon. Weather Rev., 129, 2089–2102,<2089:IOTDFA>2.0.CO;2, 2001. 

Gelaro, R., McCarty, W., Suárez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., da Silva, A. M., Gu, W., Kim, G., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454,, 2017. 

Glahn, H. R. and Lowry, D. A.: The use of model output statistics (MOS) in objective weather forecasting, J. Appl. Meteor., 11, 1203–1211,<1203:TUOMOS>2.0.CO;2, 1972. 

Gregory, D. and Rowntree, P. R.: A mass flux convection scheme with representation of cloud ensemble characteristics and stability-dependent closure, Mon. Weather Rev., 118, 1483–1506,<1483:AMFCSW>2.0.CO;2, 1990. 

Gustafson Jr., W. I., Ma, P.-L., and Singh, B.: Precipitation characteristics of CAM5 physics at mesoscale resolution during MC3E and the impact of convective timescale choice, J. Adv. Model. Earth Syst., 6, 1271–1287,, 2014. 

Harris, B. A. and Kelly, G.: A satellite radiance-bias correction scheme for data assimilation, Q. J. Roy. Meteorol. Soc., 127, 1453–1468,, 2001. 

Hartmann, D. L., Klein Tank, A. M. G., Rusticucci, M., Alexander, L. V., Brönnimann, S., Charabi, Y., Dentener, F. J., Dlugokencky, E. J., Easterling, D. R., Kaplan, A., Soden, B. J., Thorne, P. W., Wild, M., and Zhai, P. M.: Observations: Atmosphere and Surface, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 2013. 

Hersbach, H. and Dee, D.: ERA5 reanalysis is in production, ECMWF Newsletter No. 147, 7, available at: (last access: 17 May 2019), 2016. 

Holt, E. and Wang, J.: Trends in wind speed at wind turbine height of 80 m over the contiguous United States using the North American Regional Reanalysis (NARR), J. Appl. Meteor. Climatol., 51, 2188–2202,, 2012. 

Howard, T. and Clark, P.: Correction and downscaling of NWP wind speed forecasts, Meteorol. Apps., 14, 105–116,, 2007. 

Huffman, G. J., Adler, R. F., Bolvin, D. T., Gu, G., Nelkin, E. J., Bowman, K. P., Hong, Y., Stocker, E. F., and Wolff, D. B.: The TRMM multisatellite precipitation analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales, J. Hydrometeor., 8, 38–54,, 2006. 

Ingleby, N. B.: The statistical structure of forecast errors and its representation in The Met. Office Global 3-D variational data assimilation scheme, Q. J. Roy. Meteorol. Soc., 127, 209-231,, 2001. 

Jakob, D., Su, C.-H., Eizenberg, N., Kociuba, G., Steinle, P., Fox-Hughes, P., and Bettio, L.: An atmospheric high-resolution regional reanalysis for Australia, B. Aus. Meteorol. Oceanog. Soc., 30, 16–23, 2017. 

Jermey, P. M. and Renshaw, R. J.: Precipitation representation over a two-year period in regional reanalysis, Q. J. Roy. Meteorol. Soc., 142, 1300–1310,, 2016. 

Jones, D. A., Wang, W., and Fawcett, R.: High-quality spatial climate data-sets for Australia, Aust. Meteorol. Oceanogr. J., 58, 233–248, 2009. 

Kallberg, P.: Forecast drift in ERA-Interim. ERA report series 10, available at: (last access: 17 May 2019), 2011. 

Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J., Zhu, Y., Chelliah, M., Ebisuzaki, W., Higgins, W., Janowiak, J., Mo, K. C., Ropelewski, C., Wang, J., Leetmaa, A., Reynolds, R., Jenne, R., and Joseph, D.: The NCEP/NCAR 40-Year Reanalysis Project, B. Am. Meteorol. Soc., 77, 437–472,<0437:TNYRP>2.0.CO;2, 1996. 

Le Marshall, J., Xiao, Y., Norman, R., Zhang, K., Rea, A., Cucurull, L., Seecamp, R., Steinle, P., Puri, K., and Le, T.: The beneficial impact of radio occultation observations on Australian region forecasts, Aust. Meteorol. Oceanogr. J., 60, 121–125, 2010. 

Le Marshall, J., Seecamp, R., Xiao, Y., Gregory, P., Jung, J., Stienle, P., Skinner, T., Tingwell, C., and Le, T.: The Operational Generation of Continuous Winds in the Australian Region and Their Assimilation with 4DVar, Weather Forecast., 28, 504–514,, 2013. 

Lean, H. W., Clark, P. A., Dixon, M., Roberts, N. M., Fitch, A., Forbes, R., and Halliwell, C.: Characteristics of high-resolution versions of the Met Office Unified Model for forecasting convection over the United Kingdom, Mon. Weather Rev., 136, 3408–3424,, 2008. 

Lock, A. P., Brown, A. R., Bush, M. R., Martin, G. M., and Smith, R. N. B.: A new boundary layer mixing scheme. Part I: Scheme description and single-column model tests, Mon. Weather Rev., 128, 3187–3199,<3187:ANBLMS>2.0.CO;2, 2000. 

Lorenc, A. C.: Modelling of error covariances by 4D-Var data assimilation, Q. J. Roy. Meteorol. Soc., 129, 3167–3182,, 2003. 

Lorenc, A. C. and Hammon, O.: Objective quality control of observations using Bayesian methods. Theory, and a practical implementation, Q. J. Roy. Meteorol. Soc., 114, 515–543,, 1988. 

Lorenc, A. C. and Payne, T. J.: 4D-Var and the butterfly Effect: Statistical four-dimensional data assimilation for a wide range of scales, Q. J. Roy. Meteorol. Soc., 133, 607–614,, 2007. 

Loveland, T. R., Reed, B. C., Brown, J. F., Ohlen, D. O., Zhu, Z., Yang, L., and Merchant, J. W.: Development of a global land cover characteristics database and IGBP DISCover from 1 km AVHRR data, Int. J. Remote Sens., 21, 1303–1330,, 2000. 

Ma, Y., Zhou, X., Bi, D., Sun, Z., and Hirst, A. C.: Improved air-sea flux algorithms in an ocean-atmosphere coupled model for simulation of global ocean SST and its tropical pacific variability, Clim. Dynam., 44, 1473–1485,, 2015. 

Mahmood, S., Davie, J., Jermey, P., Renshaw, R., George, J. P., Rajagopal, E. N., and Rani, S. I.: Indian monsoon data assimilation and analysis regional reanalysis: Configuration and performance, Atmos. Sci. Lett., 19, e808,, 2018. 

Malloy, J. W., Krahenbuhl, D. S., Bush, C. E., Balling, R. C., Santoro, M. M., White, J. R., Elder, R. C., Pace, M. B., and Cerveny, R. S.: A surface wind extremes (“wind lulls” and “wind blows”) climatology for central North America and adjoining oceans (1979–2012), J. Appl. Meteor. Climatol., 54, 643–657,, 2015. 

Martynov, A., Laprise, R., Sushama, L., Winger, K., Separovic, L., and Dugas, B.: Reanalysis-driven climate simulation over CORDEX North America domain using the Canadian Regional Climate Model, version 5: model performance evaluation, Clim. Dynam., 41, 2973–3005,, 2013. 

Masunaga, R., Nakamura, H., Miyasaka, T., Nishii, K., and Tanimoto, Y.: Separation of climatological imprints of the Kuroshio Extension and Oyashio fronts on the wintertime atmospheric boundary layer: Their sensitivity to SST resolution prescribed for atmospheric reanalysis, J. Climate, 28, 1764–1787,, 2015. 

Masunaga, R., Nakamura, H., Kamahori, H., Onogi, K., and Okajima, S.: JRA-55CHS: An atmospheric reanalysis produced with high-resolution SST, SOLA, 14, 6–13,, 2018. 

Matthews, A. J., Pickup, G., Peatman, S. C., Clews, P., and Martin, J.: The effect of the Madden-Julian Oscillation on station rainfall and riverlevel in the Fly River System, Papua New Guinea, J. Geophys. Res.-Atmos., 118, 10926–10935,, 2013. 

Mesinger, F., DiMego, G., Kalnay, E., Mitchell, K., Shafran, P. C., Ebisuzaki, W., Jović, D., Woollen, J., Rogers, E., Berbery, E. H., Ek, M. B., Fan, Y., Grumbine, R., Higgins, W., Li, H., Lin, Y., Manikin, G., Parrish, D., and Shi, W.: North American Regional Reanalysis, B. Am. Meteorol. Soc., 87, 343–360,, 2006. 

Moore, R. J.: The PDM rainfall-runoff model, Hydrol. Earth Syst. Sci., 11, 483–499,, 2007. 

Puri, K., Dietachmayer, G., Steinle, P., Dix, M., Rikus, L., Logan, L., Naughton, M., Tingwell, C., Xiao, Y., Barras, V., Bermous, I., Bowen, R., Deschamps, L., Franklin, C., Fraser, J., Glowacki, T., Harris, B., Lee, J., Le, T., Roff, G., Sulaiman, A., Sims, H., Sun, X., Sun, Z., Zhu, H., Chattopadhyay, M. and Engel, C.: Implementation of the initial ACCESS numerical weather prediction system, Aust. Meteorol. Oceanogr. J., 63, 265–284, 2013. 

Radic, V. and Clarke, G. K. C.: Evaluation of IPCC models' performance in simulating late-twentieth-century climatologies and weather Patterns over North America, J. Climate, 24, 5257–5274,, 2011. 

Ramella Pralungo, L. and Haimberger, L.: A “Global Radiosonde and tracked-balloon Archive on Sixteen Pressure levels” (GRASP) going back to 1905 – Part 2: homogeneity adjustments for pilot balloon and radiosonde wind data, Earth Syst. Sci. Data, 6, 297–316,, 2014. 

Ramella Pralungo, L., Haimberger, L., Stickler, A., and Brönnimann, S.: A global radiosonde and tracked balloon archive on 16 pressure levels (GRASP) back to 1905 – Part 1: Merging and interpolation to 00:00 and 12:00 GMT, Earth Syst. Sci. Data, 6, 185–200,, 2014. 

Randall, D. A., Wood, R. A., Bony, S., Colman, R., Fichefet, T., Fyfe, J., Kattsov, V., Pitman, A., Shukla, J., Srinivasan, J., Stouffer, R. J., Sumi, A., and Taylor, K. E.: Climate models and their evaluation. In: climate change 2007: The physical science basis, in: Contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge University Press, Cambridge and New York, NY, 2007. 

Rawlins, F., Ballard, S. P., Bovis, K. J., Clayton, A. M., Li, D., Inverarity, G. W., Lorenc, A. C., and Payne, T. J.: The Met Office global 4-dimensional data assimilation system, Q. J. Roy. Meteorol. Soc., 133, 347–362,, 2007. 

Renshaw, R., Jermey, P., Barker, D., Maycock, A., and Oxley, S.: EURO4M regional reanalysis system. Forecasting Research Technical Report No. 583, available at: (last access: 13 February 2018), 2013. 

Ridal, M., Olsson, E., Unden, P., Zimmermann, K., and Ohlsson, A.: HARMONIE reanalysis report of results and dataset, UERRA Project Deliverable D2.7, available at: (last access: 13 February 2018), 2017. 

Roberts-Jones, J., Fiedler, E. K., and Martin, M. J.: Daily, global, high-resolution SST and sea ice reanalysis for 1985–2007 Using the OSTIA system, J. Climate, 25, 6215–6232,, 2012. 

Rose, S. and Apt, J.: Quantifying sources of uncertainty in reanalysis derived wind speed, Renew. Energy, 94, 157–165,, 2016. 

Ruiz-Barradas, A. and Nigam, S.: IPCC's twentieth-century climate simulations: Varied representations of North American hydroclimate variability, J. Climate, 19, 4041–4058,, 2006. 

Sapiano, M. R. P. and Arkin, P. A.: An intercomparison and validation of high-resolution satellite precipitation estimates with 3-hourly gauge data, J. Hydrometeorol., 10, 149–166,, 2009. 

Scinocca, J. F. and McFarlane, N. A.: The variability of modeled tropical precipitation, J. Atmos. Sci., 61, 1993–2015, 2004. 

Sheridan, P., Smith, S., Brown, A., and Vosper, S.: A simple height-based correction for temperature downscaling in complex terrain, Meteor. App., 17, 329–339,, 2010. 

Simard, M., Pinto, N., Fisher, J. B., and Baccini, A.: Mapping forest canopy height globally with spaceborne lidar, J. Geophys. Res.-Biogeosci., 116, G04021,, 2011. 

Smith, I., Moise, A., Inape, K., Murphy, B., Colman, R., Power, S., and Chung, C.: ENSO-related rainfall changes over the New Guinea region, J. Geophys. Res.-Atmos., 118, 10665–10675,, 2013. 

Thorne, P. W. and Vose, R. S.: Reanalyses suitable for characterizing long-term trends, B. Am. Meteorol. Soc., 91, 353–361,, 2010. 

UK Met Office: Met Office Science Repository Service, available at:, last access: 23 May 2019. 

Walters, D., Boutle, I., Brooks, M., Melvin, T., Stratton, R., Vosper, S., Wells, H., Williams, K., Wood, N., Allen, T., Bushell, A., Copsey, D., Earnshaw, P., Edwards, J., Gross, M., Hardiman, S., Harris, C., Heming, J., Klingaman, N., Levine, R., Manners, J., Martin, G., Milton, S., Mittermaier, M., Morcrette, C., Riddick, T., Roberts, M., Sanchez, C., Selwood, P., Stirling, A., Smith, C., Suri, D., Tennant, W., Vidale, P. L., Wilkinson, J., Willett, M., Woolnough, S., and Xavier, P.: The Met Office Unified Model Global Atmosphere 6.0/6.1 and JULES Global Land 6.0/6.1 configurations, Geosci. Model Dev., 10, 1487–1520,, 2017a. 

Walters, D., Baran, A., Boutle, I., Brooks, M., Earnshaw, P., Edwards, J., Furtado, K., Hill, P., Lock, A., Manners, J., Morcrette, C., Mulcahy, J., Sanchez, C., Smith, C., Stratton, R., Tennant, W., Tomassini, L., Van Weverberg, K., Vosper, S., Willett, M., Browse, J., Bushell, A., Dalvi, M., Essery, R., Gedney, N., Hardiman, S., Johnson, B., Johnson, C., Jones, A., Mann, G., Milton, S., Rumbold, H., Sellar, A., Ujiie, M., Whitall, M., Williams, K., and Zerroukat, M.: The Met Office Unified Model Global Atmosphere 7.0/7.1 and JULES Global Land 7.0 configurations, Geosci. Model Dev. Discuss.,, in review, 2017b. 

Wang, Z., Siems, S. T., Belusic, D., Manton, M. J., and Huang, Y.: A climatology of the precipitation over the Southern Ocean as observed at Macquarie Island, J. Appl. Meteorol. Climatol., 54, 2321–2337,, 2015. 

Williamson, D. L.: The effect of time steps and time-scales on parametrization suites, Q. J. Roy. Meteorol. Soc., 139, 548–560,, 2013. 

Wilson, D. R. and Ballard, S. P.: A microphysically based precipitation scheme for the UK Meteorological Office Unified Model, Q. J. Roy. Meteorol. Soc., 125, 1607–1636,, 1999. 

Wood, N., Staniforth, A., White, A., Allen, T., Diamantakis, M., Gross, M., Melvin, T., Smith, C., Vosper, S., Zerroukat, M., and Thuburn, J.: An inherently mass-conserving semi-implicit semi-Lagrangian discretization of the deep-atmosphere global non-hydrostatic equations, Q. J. Roy. Meteorol. Soc., 140, 1505–1520,, 2014. 

Zhao, M., Zhang, H.-Q., and Dharssi, I.: Impact of land-surface initialization on ACCESS-S1 and comparison with POAMA, Bureau of Meteorology Research Report No. 023, available at: (last access: 17 May 2019), 2017. 

Zhu, H. and Dietachmayer, G.: Improving ACCESS-C convection settings, Bureau Research Report No. 008, available at: (last access: 17 May 2019), 2015. 

Ziese, M., Rauthe-Schöch, A., Becker, A., Finger, P., Meyer-Christoffer, A., and Schneider, U.: GPCC full data daily version.2018 at 1.0: Daily land-surface precipitation from rain-gauges built on GTS-based and historic data,, 2018. 

Zick, S. E. and Matyas, C. J.: Tropical cyclones in the North American Regional Reanalysis: An assessment of spatial biases in location, intensity, and structure, J. Geophys. Res.-Atmos., 120, 1651–1669,, 2015. 

Short summary
The Bureau of Meteorology Atmospheric Regional Reanalysis for Australia (BARRA) is the first regional reanalysis for Australia, NZ, and SE Asia. It offers realistic depictions of near-surface meteorology at a scale required for emergency services, defence, and other major sectors such as energy and agriculture. It uses a consistent method of analysing the atmosphere, with a higher-resolution model over 1990 to 2018, and can provide greater understanding of past weather, including extreme events.