The variation and controls of mean transit times in Australian headwater catchments

Determining mean transit times in headwater catchments is critical for understanding catchment functioning and understanding their responses to changes in landuse or climate. Determining whether mean transit times (MTTs) correlate with drainage density, slope angle, area, or land cover permits a better understanding of the controls on water flow through catchments and allows first-order predictions of MTTs in other catchments to be made. This study assesses whether there are identifiable controls on MTTs determined using ³H in headwater catchments of southeast Australia. Despite MTTs at baseflow varying from a few years to >100 years, it was difficult to predict MTTs using single or groups of readily-measured catchment attributes. The lack of readily-identifiable correlations hampers the prediction of MTTs in adjacent catchments even where these have similar geology, land use, and topography. The long MTTs of the Australian headwater catchments are probably in part due to the catchments having high storage volumes in deeply-weathered regolith, combined with low recharge rates due to high evapotranspiration. However, the difficulty in estimating storage volumes at the catchment scale hampers the use of this parameter to estimate MTTs. The runoff coefficient (the fraction of rainfall exported via the stream) is probably also controlled by evapotranspiration and recharge rates. Correlations between the runoff coefficient and MTTs in individual catchments allow predictions of MTTs in nearby catchments to be made. MTTs are shorter in high rainfall periods as the catchments wet up and shallow water stores are mobilized. Despite the contribution of younger water, the major ion geochemistry in individual catchments commonly does not correlate with MTTs, probably reflecting heterogeneous reactions and varying degrees of evapotranspiration. Documenting MTTs in catchments with high storage volumes and/or low recharge rates elsewhere is important for understanding MTTs in diverse environments.     

Analysis of nonstationarity in low flow in the Loess Plateau of China

An approach for nonstationary low‐flow frequency analysis is developed and demonstrated on a dataset from the rivers on the Loess Plateau of China. Nonstationary low‐flow frequency analysis has drawn significant attention in recent years by establishing relationships between low‐flow series and explanatory variables series, but few studies have tested whether the time‐varying moments of low flow can be fully described by the time‐varying moments of the explanatory variables. In this research, the low‐flow distributions are analytically derived from the 2 basic explanatory variables—the recession duration and the recession coefficient—with the assumption that the recession duration and recession coefficient variables follow exponential and gamma distributions, respectively; the derived low‐flow distributions are applied to test whether the time‐varying moments of explanatory variables can explain the nonstationarities found in the low‐flow variable. The effects of ecosystem construction measures, that is, check dam, terrace, forest, and grassland, on the recession duration and recession coefficient are further discussed. Daily flow series from 11 hydrological stations from the Loess Plateau are used and processed with a moving average technique. Low‐flow data are extracted following the pit under threshold approach. Six of the 11 low‐flow series show significant nonstationarities at the 5% significance level, and the trend curves of the moments of low flow are in close agreement with the curves estimated from the derived distribution with time‐dependent moments of the recession duration and time‐constant moments of the recession coefficient. It is indicated that the nonstationarity in the low‐flow distribution results from the nonstationarity in the recession duration in all 6 cases, and the increase in the recession duration is resulted from large‐scale ecosystem constructions rather than climate change. The large‐scale ecosystem constructions are found to have more influence on the decrease in streamflow than on the increase in watershed storage, thus resulting in the reduction of low flow. A high return period for the initial fixed design value decreases dramatically with an increasing recession duration.     

Glacial and environmental changes over the last 60 000 years in the Polar Ural Mountains,Arctic Russia,inferred from a high‐resolution lake record and other observations from adjacent areas

Our knowledge about the glaciation history in the Russian Arctic has to a large extent been based on geomorphological mapping supplemented by studies of short stratigraphical sequences found in exposed sections. Here we present new geochronological data from the Polar Ural Mountains along with a high‐resolution sediment record from Bolshoye Shchuchye, the largest and deepest lake in the mountain range. Seismic profiles show that the lake contains a 160‐m‐thick sequence of unconsolidated lacustrine sediments. A well‐dated 24‐m‐long core from the southern end of the lake spans the last 24 cal. ka. From downward extrapolation of sedimentation rates we estimate that sedimentation started about 50–60 ka ago, most likely just after a large glacier had eroded older sediments from the basin. Terrestrial cosmogenic nuclide (TCN) exposure dating (¹⁰Be) of boulders and Optically Stimulated Luminescence (OSL) dating of sediments indicate that this part of the Ural Mountains was last covered by a coherent ice‐field complex during Marine Isotope Stage (MIS) 4. A regrowth of the glaciers took place during a late stage of MIS 3, but the central valleys remained ice free until the present. The presence of small‐ and medium‐sized glaciers during MIS 2 is reflected by a sequence of glacial varves and a high sedimentation rate in the lake basin and likewise from ¹⁰Be dating of glacial boulders. The maximum extent of the mountain glaciers during MIS 2 was attained prior to 24 cal. ka BP. Some small present‐day glaciers, which are now disappearing completely due to climate warming, were only slightly larger during the Last Glacial Maximum (LGM) as compared to AD 1953. A marked decrease in sedimentation rate around 18–17 cal. ka BP indicates that the glaciers then became smaller and probably disappeared altogether around 15–14 cal. ka BP.     

A Global Terrestrial Reference Frame from simulated VLBI and SLR data in view of GGOS

In this study, we assess the impact of two combination strategies, namely local ties (LT) and global ties (GT), on the datum realization of Global Terrestrial Reference Frames in view of the Global Geodetic Observing System requiring 1 mm-accuracy. Simulated Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data over a 7 year time span was used. The LT results show that the geodetic datum can be best transferred if the precision of the LT is at least 1 mm. Investigating different numbers of LT, the lack of co-located sites on the southern hemisphere is evidenced by differences of 9 mm in translation and rotation compared to the solution using all available LT. For the GT, the combination applying all Earth rotation parameters (ERP), such as pole coordinates and UT1-UTC, indicates that the rotation around the Z axis cannot be adequately transferred from VLBI to SLR within the combination. Applying exclusively the pole coordinates as GT, we show that the datum can be transferred with mm-accuracy within the combination. Furthermore, adding artificial stations in Tahiti and Nigeria to the current VLBI network results in an improvement in station positions by 13 and 12%, respectively, and in ERP by 17 and 11%, respectively. Extending to every day VLBI observations leads to 65% better ERP estimates compared to usual twice-weekly VLBI observations.     

On the consistency of the current conventional EOP series and the celestial and terrestrial reference frames

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.     

The extension of the parametrization of the radio source coordinates in geodetic VLBI and its impact on the time series analysis

The radio sources within the most recent celestial reference frame (CRF) catalog ICRF2 are represented by a single, time-invariant coordinate pair. The datum sources were chosen mainly according to certain statistical properties of their position time series. Yet, such statistics are not applicable unconditionally, and also ambiguous. However, ignoring systematics in the source positions of the datum sources inevitably leads to a degradation of the quality of the frame and, therefore, also of the derived quantities such as the Earth orientation parameters. One possible approach to overcome these deficiencies is to extend the parametrization of the source positions, similarly to what is done for the station positions. We decided to use the multivariate adaptive regression splines algorithm to parametrize the source coordinates. It allows a great deal of automation, by combining recursive partitioning and spline fitting in an optimal way. The algorithm finds the ideal knot positions for the splines and, thus, the best number of polynomial pieces to fit the data autonomously. With that we can correct the ICRF2 a priori coordinates for our analysis and eliminate the systematics in the position estimates. This allows us to introduce also special handling sources into the datum definition, leading to on average 30 % more sources in the datum. We find that not only the CPO can be improved by more than 10 % due to the improved geometry, but also the station positions, especially in the early years of VLBI, can benefit greatly.     

GPS/BDS short-term ISB modelling and prediction

GPS Solutions - The Chinese BeiDou Navigation Satellite System (BDS) has completed its first milestone by providing coverage of the Asia–Pacific area navigation service since December 27,...     

Performance of the WRF model to simulate the seasonal and interannual variability of hydrometeorological variables in East Africa: a case study for the Tana River basin in Kenya

This study investigates the ability of the regional climate model Weather Research and Forecasting (WRF) in simulating the seasonal and interannual variability of hydrometeorological variables in the Tana River basin (TRB) in Kenya, East Africa. The impact of two different land use classifications, i.e., the Moderate Resolution Imaging Spectroradiometer (MODIS) and the US Geological Survey (USGS) at two horizontal resolutions (50 and 25 km) is investigated. Simulated precipitation and temperature for the period 2011–2014 are compared with Tropical Rainfall Measuring Mission (TRMM), Climate Research Unit (CRU), and station data. The ability of Tropical Rainfall Measuring Mission (TRMM) and Climate Research Unit (CRU) data in reproducing in situ observation in the TRB is analyzed. All considered WRF simulations capture well the annual as well as the interannual and spatial distribution of precipitation in the TRB according to station data and the TRMM estimates. Our results demonstrate that the increase of horizontal resolution from 50 to 25 km, together with the use of the MODIS land use classification, significantly improves the precipitation results. In the case of temperature, spatial patterns and seasonal cycle are well reproduced, although there is a systematic cold bias with respect to both station and CRU data. Our results contribute to the identification of suitable and regionally adapted regional climate models (RCMs) for East Africa.