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.

Triple collocation (TC) is routinely used to resolve approximated linear relationships between different measurements (or representations) of a geophysical variable that are subject to errors. It has been utilized in the context of calibration, validation, bias correction, and error characterization to allow comparisons of diverse data records from various direct and indirect measurement techniques including in situ remote sensing and model‐based approaches. However, successful applications of TC require sufficiently large numbers of coincident data points from three independent time series and, within the analysis period, homogeneity of their linear relationships and error structures. These conditions are difficult to realize in practice due to infrequent spatiotemporal sampling of satellite and ground‐based sensors. TC can, however, be generalized within the framework ofinstrumental variable(IV) regression theory to address some of the conceptual constraints of TC. We review the theoretics of IV and consider one possible strategy to circumvent the three‐data constraint by use of lagged variables (LV) as instruments. This particular implementation of IV is suitable for circumstances where multiple data records are limited and the geophysical variable of interest is sampled at time intervals shorter than its temporal correlation length. As a demonstration of utility, the LV method is applied to microwave satellite soil moisture data sets to recover their errors over Australia and to estimate temporal properties of their relationships with in situ and model data. These results are compared against standard two‐data linear estimators and the TC estimator as benchmark.

Global soil moisture records are essential for studying the role of hydrologic processes within the larger earth system. Various studies have shown the benefit of assimilating satellite‐based soil moisture data into water balance models or merging multisource soil moisture retrievals into a unified data set. However, this requires an appropriate parameterization of the error structures of the underlying data sets. While triple collocation (TC) analysis has been widely recognized as a powerful tool for estimating random error variances of coarse‐resolution soil moisture data sets, the estimation of error cross covariances remains an unresolved challenge. Here we propose a method—referred to as extended collocation (EC) analysis—for estimating error cross‐correlations by generalizing the TC method to an arbitrary number of data sets and relaxing the therein made assumption of zero error cross‐correlation for certain data set combinations. A synthetic experiment shows that EC analysis is able to reliably recover true error cross‐correlation levels. Applied to real soil moisture retrievals from Advanced Microwave Scanning Radiometer‐EOS (AMSR‐E) C‐band and X‐band observations together with advanced scatterometer (ASCAT) retrievals, modeled data from Global Land Data Assimilation System (GLDAS)‐Noah and in situ measurements drawn from the International Soil Moisture Network, EC yields reasonable and strong nonzero error cross‐correlations between the two AMSR‐E products. Against expectation, nonzero error cross‐correlations are also found between ASCAT and AMSR‐E. We conclude that the proposed EC method represents an important step toward a fully parameterized error covariance matrix for coarse‐resolution soil moisture data sets, which is vital for any rigorous data assimilation framework or data merging scheme.

Assimilation of remotely sensed soil moisture data (SM-DA) to correct soil water stores of rainfall-runoff models has shown skill in improving streamflow prediction. In the case of large and sparsely monitored catchments, SM-DA is a particularly attractive tool. Within this context, we assimilate satellite soil moisture (SM) retrievals from the Advanced Microwave Scanning Radiometer (AMSR-E), the Advanced Scatterometer (ASCAT) and the Soil Moisture and Ocean Salinity (SMOS) instrument, using an Ensemble Kalman filter to improve operational flood prediction within a large (> 40 000 km2) semi-arid catchment in Australia. We assess the importance of accounting for channel routing and the spatial distribution of forcing data by applying SM-DA to a lumped and a semi-distributed scheme of the probability distributed model (PDM). Our scheme also accounts for model error representation by explicitly correcting bias in soil moisture and streamflow in the ensemble generation process, and for seasonal biases and errors in the satellite data. Before assimilation, the semi-distributed model provided a more accurate streamflow prediction (Nash–Sutcliffe efficiency, NSE = 0.77) than the lumped model (NSE = 0.67) at the catchment outlet. However, this did not ensure good performance at the "ungauged" inner catchments (two of them with NSE below 0.3). After SM-DA, the streamflow ensemble prediction at the outlet was improved in both the lumped and the semi-distributed schemes: the root mean square error of the ensemble was reduced by 22 and 24%, respectively; the false alarm ratio was reduced by 9% in both cases; the peak volume error was reduced by 58 and 1%, respectively; the ensemble skill was improved (evidenced by 12 and 13% reductions in the continuous ranked probability scores, respectively); and the ensemble reliability was increased in both cases (expressed by flatter rank histograms). SM-DA did not improve NSE. Our findings imply that even when rainfall is the main driver of flooding in semi-arid catchments, adequately processed satellite SM can be used to reduce errors in the model soil moisture, which in turn provides better streamflow ensemble prediction. We demonstrate that SM-DA efficacy is enhanced when the spatial distribution in forcing data and routing processes are accounted for. At ungauged locations, SM-DA is effective at improving some characteristics of the streamflow ensemble prediction; however, the updated prediction is still poor since SM-DA does not address the systematic errors found in the model prior to assimilation.