Development of an operational rainfall data quality-control scheme based on radar-raingauge co-kriging analysis

Abstract

Automatic raingauge data often serve as an important input to hydrological and weather warning operations. They are not only fundamental in quantitative rainfall analysis, but also act as the ground truth in warning operation and forecast validation. Quality control is required before the data can be used quantitatively due to systematic and random errors. Extremely large random errors and unreasonably small or false zero values can hamper effective monitoring of heavy rain. Yet both are difficult to detect in real-time by objective means. In an attempt to address these problems, a rainfall data quality-control scheme based on radar-raingauge co-kriging analysis was developed. The important threshold values required in the data quality control of 60-min raingauge rainfall were determined from a detailed analysis of the distributions of rainfall residuals defined as the arithmetic difference and the logarithm of the ratio between a raingauge measurement and its co-kriging estimate. The scheme has been developed and is in real-time use in Hong Kong, a coastal city of about 1100 km² area with more than 150 raingauges installed. Geographically, it is located in the subtropics and dominated by heavy convective rainfall in the summer. As a basis of the quality-control scheme, the co-kriging rainfall analysis was shown through a verification exercise to be superior to those obtained by the Barnes analysis and ordinary kriging of raingauge data. The performance of the quality-control algorithm was assessed using selected cases and controlled tests, and was found to be satisfactory, with a high error detection rate for the two targeted types of error. Limitations and operational issues identified during a real-time trial of the quality-control scheme are also discussed.

Citation Yeung, H.Y., Man, C., Chan, S.T., and Seed, A., 2014. Development of an operational rainfall data quality-control scheme based on radar-raingauge co-kriging analysis. Hydrological Sciences Journal, 59 (7), 1285–1299. http://dx.doi.org/10.1080/02626667.2013.839873     

Implications of subsurface thermal–hydraulic–mechanical processes associated with glaciation on Shield flow system evolution and performance assessment

A deep geologic repository (DGR) situated on the Canadian Shield will be subject to long-term climate change that will markedly alter surface conditions as a result of glaciation and permafrost penetration. Systematic, two-dimensional and three-dimensional coupled thermal–hydraulic–mechanical finite-element simulations with varying degrees of coupling, including depth-dependent salinity (represented as a change in groundwater density) and temperature-dependent density and viscosity, were undertaken to address the implications of glaciation on groundwater flow system dynamics as it could affect DGR performance. The modelling domain consisted of a 1.6-km deep sub-regional scale (≈100 km²) fractured Shield flow system. Initial and transient thermal, hydraulic and mechanical boundary conditions were developed from two realizations of the University of Toronto Glacial Systems Model of the last Laurentide glaciation. Results indicate that during the glacial loading/unloading cycle, for this particular conceptual model, there is limited penetration of glacial meltwaters to depth and small residual anomalous hydraulic head. During glacial coverage, the mechanical factor of safety increases in the moderately fractured and sparsely fractured rock mass, but principal effective stress reorientation also occurs. Given the assumed nonglacial in situ state of stress and mechanical properties, the fracture zones were predicted to be less stable under glacial conditions.

Quantifying potential propagule pressure of aquatic invasive species from the commercial shipping industry in Canada

We quantify and compare different measures of potential propagule pressure (PPP) of aquatic invasive species (AIS) from commercial vessels in Canada. We used ship arrivals and ballast water discharge volumes as proxies for PPP from ballast water organisms, and wetted surface area (WSA) as a proxy for hull fouling PPP, to determine their relative contributions to total PPP. For three regions studied, PPP proxies correlated significantly across ports and some vessel categories. Relative contributions of ship arrivals, ballast discharge, and WSAs to PPP, evidenced by non-significant correlations across these measures, varied across regions, ports, vessel types, and seasons. Flow-through (dominant on east and west coasts) and empty-refill (in Great Lakes-St. Lawrence region) were the major ballast water exchange methods employed by the vessels surveyed. These methods have different biological efficacy for AIS removal, influencing PPP. Our study illustrates benefits and limitations of using different PPP proxies to estimate invasion risk.     

Industrial structural changes in Hong Kong, China under one country, two systems framework

Hong Kong has undergone substantial economic transformations and developed into a sophisticated business and financial center in the Asia-Pacific region,since the return of sovereignty to China as a Special Administrative Region under the One Country Two Systems (OCTS) in 1997.This paper discusses and analyzes the industrial structural changes of Hong Kong in recent decades as well as its future challenges and opportunities.The data and finding reveal that even though Hong Kong will face fierce competition from the Mainland’s cities as the rise of China,the important role as a bridge between China and the rest of the world will brace Hong Kong itself under the OCTS for developing into a service hub for business and trade in the Asia-Pacific region.     

Projections of severe droughts in future climate in Southeast Brazil: a case study in Southern Minas Gerais State,Brazil

Theoretical and Applied Climatology - South of Minas Gerais state, in Southeast Brazil, is known for the coffee crop production (more than 30% of country’s production) and hydroelectricity...     

High-resolution subtropical summer precipitation derived from dynamical downscaling of the NCEP/DOE reanalysis: how much small-scale information is added by a regional model?

This study assesses the regional-scale summer precipitation produced by the dynamical downscaling of analyzed large-scale fields. The main goal of this study is to investigate how much the regional model adds smaller scale precipitation information that the large-scale fields do not resolve. The modeling region for this study covers the southeastern United States (Florida, Georgia, Alabama, South Carolina, and North Carolina) where the summer climate is subtropical in nature, with a heavy influence of regional-scale convection. The coarse resolution (2.5° latitude/longitude) large-scale atmospheric variables from the National Center for Environmental Prediction (NCEP)/DOE reanalysis (R2) are downscaled using the NCEP/Environmental Climate Prediction Center regional spectral model (RSM) to produce precipitation at 20?km resolution for 16 summer seasons (1990?C2005). The RSM produces realistic details in the regional summer precipitation at 20?km resolution. Compared to R2, the RSM-produced monthly precipitation shows better agreement with observations. There is a reduced wet bias and a more realistic spatial pattern of the precipitation climatology compared with the interpolated R2 values. The root mean square errors of the monthly R2 precipitation are reduced over 93% (1,697) of all the grid points in the five states (1,821). The temporal correlation also improves over 92% (1,675) of all grid points such that the domain-averaged correlation increases from 0.38 (R2) to 0.55 (RSM). The RSM accurately reproduces the first two observed eigenmodes, compared with the R2 product for which the second mode is not properly reproduced. The spatial patterns for wet versus dry summer years are also successfully simulated in RSM. For shorter time scales, the RSM resolves heavy rainfall events and their frequency better than R2. Correlation and categorical classification (above/near/below average) for the monthly frequency of heavy precipitation days is also significantly improved by the RSM.     

Numerical simulation of nonlinear dispersive waves propagating over a submerged bar by IB–VOF model

It is a good test for a numerical model to simulate progressive waves propagating over a submerged bar with a relatively high ratio of slopes. In this paper, the combined IB–VOF model is used to predict nonlinear dispersive waves propagating over a submerged bar with both slopes of 1:2. The predicted free surface elevations are compared with the experimental data and numerical results presented by other researchers. The comparison shows that the IB–VOF model is able to provide satisfactory wave profiles in the shallow water with strong nonlinear effects and in the wave transmitted region with strong wave dispersion in particular. Moreover, the wave evolution behind the submerged bar is described in detail, including the spatial wave profile modulation, spectral analysis of the time-series waves, flow velocity and pressure fields, and kinetic energy distribution. The effect of fluid viscosity on the numerical simulations is also studied, and it is found that the effect on the wave evolution considered in this paper is not significant. Finally, the hydrodynamic force acting on the bar is calculated using the IB–VOF model.     

U–Pb and Hf isotopic study of zircons from migmatised amphibolites in the Cathaysia Block: Implications for the early Paleozoic peak tectonothermal event in Southeastern China

Amphibolites occur in a number of localities in the Cathaysia Block, some of them have been migmatised and their protoliths represent basaltic magmas erupted in various tectonic settings. Four migmatised amphibolites were collected from Jiangxi and Fujian Provinces. Cathodo-luminescence images of zircons extracted from the representative amphibolites show unzoned or sector-zoned structure. LA-ICP-MS analysis indicates that most zircons have high Th/U ratios and yield U–Pb zircon ages of 446 ± 5, 435 ± 2, 434 ± 4 and 423 ± 2 Ma, respectively. Lu–Hf isotopic analysis on these zircons gives Hf model ages ranging from 900 to 1200 Ma. Based on lithological observations and previously published geochronological data, we interpret that these U–Pb ages record an important tectonothermal event that led to the migmatization. This early Paleozoic (Caledonian) tectonothermal event in the Southeastern China has a great tectonic implication for the evolutionary history of the Cathaysia Block.     

An empirical model of subcritical bedform migration

The ability to predict bedform migration in rivers is critical for estimating bed material load, yet there is no relation for predicting bedform migration (downstream translation) that covers the full range of conditions under which subcritical bedforms develop. Here, the relation between bedform migration rates and transport stage is explored using a field and several flume data sets. Transport stage is defined as the non‐dimensional Shields stress divided by its value at the threshold for sediment entrainment. Statistically significant positive correlations between both ripple and dune migration rates and transport stage are found. Stratification of the data by the flow depth to grain‐size ratio improved the amount of variability in migration rates that was explained by transport stage to ca 70%. As transport stage increases for a given depth to grain‐size ratio, migration rates increase. For a given transport stage, the migration rate increases as the flow depth to grain‐size ratio gets smaller. In coarser sediment, bedforms move faster than in finer sediment at the same transport stage. Normalization of dune migration rates by the settling velocity of bed sediment partially collapses the data. Given the large amount of variability that arises from combining data sets from different sources, using different equipment, the partial collapse is remarkable and warrants further testing in the laboratory and field.