Thresholds of hydrologic flow regime of a river and investigation of climate change impact—the case of the Lower Brahmaputra river Basin

The sustainability of social-ecological systems depends on river flows being maintained within a range to which those systems are adapted. In order to determine the extent of this natural range of variation, we assess ecological flow thresholds and the occurrence of potentially damaging flood events to society in the context of the Lower Brahmaputra river basin. The ecological flow threshold was calculated using twenty-two ‘Range of Variability (RVA)’ parameters, considering the range between?±?1 standard deviation from the mean of the natural flow. Damaging flood events were calculated using flood frequency analysis of Annual Maxima series and using the flood classification of Bangladesh. The climate change impacts on future river flow were calculated by using a weighted ensemble analysis of twelve global circulation models (GCMs) outputs driving a large-scale hydrologic model. The simulated climate change induced altered flow regime of the Lower Brahmaputra River Basin was then investigated and compared with the calculated threshold flows. The results demonstrate that various parameters including the monthly mean of low flow (January, February and March) and high flow (June, July and August) periods, the 7-day average minimum flow, and the yearly maximum flow will exceed the threshold conditions by 1956–1995 under the business-as-usual A1B and A2 future scenarios. The results have a number of policy level implications for government agencies of the Lower Brahmaputra River Basin, specifically for Bangladesh. The calculated thresholds may be used as a good basis for negotiations with other riparian countries of the basin. The methodological approach presented in this study can be applied to other river basins and provide a useful basis for transboundary water resources management.     

Triggering of El Niño by westerly wind events in a coupled general circulation model

Two ten-members ensemble experiments using a coupled ocean-atmosphere general circulation model are performed to study the dynamical response to a strong westerly wind event (WWE) when the tropical Pacific has initial conditions favourable to the development of a warm event. In the reference ensemble (CREF), no wind perturbation is introduced, whereas a strong westerly wind event anomaly is introduced in boreal winter over the western Pacific in the perturbed ensemble (CWWE). Our results demonstrate that an intense WWE is capable of establishing the conditions under which a strong El Niño event can occur. First, it generates a strong downwelling Kelvin wave that generates a positive sea surface temperature (SST) anomaly in the central-eastern Pacific amplified through a coupled ocean-atmosphere interaction. This anomaly can be as large as 2.5°C 60 days after the WWE. Secondly, this WWE also initiates an eastward displacement of the warm-pool that promotes the occurrence of subsequent WWEs in the following months. These events reinforce the initial warming through the generation of additional Kelvin waves and generate intense surface jets at the eastern edge of the warm-pool that act to further shift warm waters eastward. The use of a ten-members ensemble however reveals substantial differences in the coupled response to a WWE. Whereas four members of CWWE ensemble develop into intense El Niño warming as described above, four others display a moderate warming and two remains in neutral conditions. This diversity between the members appears to be due to the internal atmospheric variability during and following the inserted WWE. In the four moderate warm cases, the warm-pool is initially shifted eastward following the inserted WWE, but the subsequent weak WWE activity (when compared to the strong warming cases) prevents to further shift the warm-pool eastwards. The seasonal strengthening of trade winds in June–July can therefore act to shift warm waters back into the western Pacific, reducing the central-eastern Pacific warming. This strong sensitivity of the coupled response to WWEs may therefore limit the predictability of El Niño events, as the high frequency wind variability over the warm pool region remains largely unpredictable even at short time lead.     

Will climate change exacerbate water stress in Central Asia?

Millions of people in the geopolitically important region of Central Asia depend on water from snow- and glacier-melt driven international rivers, most of all the Syr Darya and Amu Darya. The riparian countries of these rivers have experienced recurring water allocation conflicts ever since the Soviet Union collapsed. Will climate change exacerbate water stress and thus conflicts? We have developed a coupled climate, land-ice and rainfall-runoff model for the Syr Darya to quantify impacts and show that climatic changes are likely to have consequences on runoff seasonality due to earlier snow-melt. This will increase water stress in unregulated catchments because less water will be available for irrigation in the summer months. Threats from geohazards, above all glacier lake outbursts, are likely to increase as well. The area at highest risk is the densely populated, agriculturally productive, and politically unstable Fergana Valley. Targeted infrastructural developments will be required in the region. If the current mismanagement of water and energy resources can be replaced with more effective resource allocation mechanisms through the strengthening of transboundary institutions, Central Asia will be able to successfully address these future climate-related challenges.     

A procedure for automated quality control and homogenization of historical daily temperature and precipitation data (APACH): part 1: quality control and application to the Argentine weather service stations

The present paper describes the quality-control component of an automatic procedure (APACH: A Procedure for Automated Quality Control and Homogenization of Weather Station Data) developed to control quality and homogenize the historical daily temperature and precipitation data from meteorological stations. The quality-control method is based on a set of decision-tree algorithms analyzing separately precipitation and minimum and maximum temperature. All our tests are non-parametric and therefore are potentially useful in regions or countries presenting different climates as those observed in Argentina. The method is applied to the 1959–2005 historical daily database of the Argentine National Weather Service. Our results are coherent with the history of the Weather Service and more specifically with the history of implementation of systematized quality control processes. In temperature, our method detects a larger number of suspect values before 1967 (when there was no quality control) and after 1997 (when only real-time quality control had been applied). In precipitation, the detection of error in extreme precipitations is complex, but our method clearly detected a strong decrease in the number of potential outliers after 1976 when the National Weather Service was militarized, and the network was strongly reduced, focusing more on airport weather stations. Also in precipitation, we analyze in detail the long dry sequences and are able to identify potential long erroneous sequences. This is important for the use of the data for hydrological or agricultural impact studies. Finally, all the data are flagged with codes representing the path followed by the record in our decision-tree algorithms. While each code is associated to one of the categories (“Useful”, “Need-Check”, “Doubtful” or “Suspect”), the final user is free to redefine such category-assignment.     

Frequency of Boundary-Layer-Top Fluctuations in Convective and Stable Conditions Using Laser Remote Sensing

The planetary boundary-layer (PBL) height is determined with high temporal and altitude resolution from lidar backscatter profiles. Then, the frequencies of daytime thermal updrafts and downdrafts and of nighttime gravity waves are obtained applying a fast Fourier transform on the temporal fluctuation of the PBL height. The principal frequency components of each spectrum are related to the dominant processes occurring at the daytime and nighttime PBL top. Two groups of cases are selected for the study: one group combines daytime cases, measured in weak horizontal wind conditions and dominated by convection. The cases show higher updraft and downdraft frequencies for the shallow, convective boundary layer and lower frequencies for a deep PBL. For cases characterized by strong horizontal winds, the frequencies directly depend on the wind speed. The temporal variation of the PBL height is determined also in the likely presence of lee waves. For nighttime cases, the main frequency components in the spectra do not show a real correlation with the nocturnal PBL height. Altitude fluctuations of the top of the nocturnal boundary layer are observed even though the boundary layer is statically stable. These oscillations are associated with the wind shear effect and with buoyancy waves at the PBL top.     

Twentieth century ENSO characteristics in the IPCC database

In this paper, we assess and compare to observations the spatial characteristics of the twentieth Century ENSO SST variability simulated by 23 models of the IPCC-AR4/CMIP3 database. The analysis is confined to the SST anomalies along the equatorial Pacific and is based on the use of a non-linear neural classification algorithm, the Self-Organizing Maps. Systematic biases include a larger than observed proportion for modelled ENSO maximum variability occurring in the Western Pacific. No clear relationship is found between this bias and the characteristics of the modelled mean state bias in the equatorial Pacific. This bias is mainly related to a misrepresentation of both El Niño and La Niña termination phases for most of the models. In contrast, the onset phase is quite well simulated. Modelled El Niño and La Niña peak phases display an asymmetric bias. Whereas the main bias of the modelled El Niño peak is to exhibit a maximum in the western Pacific, the simulated La Niña bias mainly occurs in the central Pacific. In addition, some models are able to capture the observed El Niño peak characteristics while none of them realistically simulate La Niña peaks. It also arises that the models closest to the observations score unevenly in reproducing the different phases, preventing an accurate classification of the models quality to reproduce the overall ENSO-like variability.     

One-step solution to local tie vector determination at co-located GNSS/VLBI sites

The local tie vector, which connects the different space geodetic techniques at a co-located site, plays an important role in the realization of the International Terrestrial Reference Frame (ITRF). This paper presents a new method to determine the tie vector between the GNSS and very long baseline interferometry tracking points. The parameters of the local tie vector and the axes offsets are introduced into constraint equations. The parameters are then resolved using the 3D constrained least squares adjustment. With the surveying data collected at two different sites (Kunming and Urumqi) in China, the proposed method can precisely determine the local tie vectors in a geocentric frame. The root mean square error (RMSE) is (1.2, 2.3 and 1.5 mm) and (1.0, 1.5 and 1.4 mm) for the three coordinate components at the sites in Kunming and Urumqi, respectively. The offset between the primary and secondary axes of the VLBI telescopes is estimated to be 7.5 mm in Kunming’s site and 4.0 mm in Urumqi’s site, and the corresponding RMSE is 1.8 mm and 2.0 mm for the two sites, respectively.     

Seismic motion in urban sites consisting of blocks in welded contact with a soft layer overlying a hard half-space

We address the problem of the response to a seismic wave of an urban site consisting of   N _b   blocks overlying a soft layer underlain by a hard substratum. The results of a theoretical analysis, appealing to a space–frequency mode-matching (MM) technique, are compared to those obtained by a space–time finite-element (FE) technique. The two methods are shown to give rise to the same prediction of the seismic response for   N _b = 1  , 2 and 40 blocks. The mechanism of the interaction between blocks and the ground, as well as that of the mutual interaction between blocks, are studied. It is shown, in the first part of this paper, that the presence of a small number of blocks modifies the seismic disturbance in a manner which evokes qualitatively, but not quantitatively, what was observed during the 1985 Michoacan earthquake in Mexico City. Anomalous earthquake response at a much greater level, in terms of duration, peak and cumulative amplitude of motion, is shown, by a theoretical and numerical analysis in the second part of this paper, to be induced by the presence of a large (≥10) number of identical equi-spaced blocks that are present in certain districts of many cities.     

Water-removed spectra increase the retrieval accuracy when estimating savanna grass nitrogen and phosphorus concentrations

Information about the distribution of grass foliar nitrogen (N) and phosphorus (P) is important for understanding rangeland vitality and for facilitating the effective management of wildlife and livestock. Water absorption effects in the near-infrared (NIR) and shortwave-infrared (SWIR) regions pose a challenge for nutrient estimation using remote sensing. The aim of this study was to test the utility of water-removed (WR) spectra in combination with partial least-squares regression (PLSR) and stepwise multiple linear regression (SMLR) to estimate foliar N and P, compared to spectral transformation techniques such as first derivative, continuum removal and log-transformed (Log(1/R)) spectra. The study was based on a greenhouse experiment with a savanna grass species (Digitaria eriantha). Spectral measurements were made using a spectrometer. The D. eriantha was cut, dried and chemically analyzed for foliar N and P concentrations. WR spectra were determined by calculating the residual from the modelled leaf water spectra using a nonlinear spectral matching technique and observed leaf spectra. Results indicated that the WR spectra yielded a higher N retrieval accuracy than a traditional first derivative transformation (R²=0.84, RMSE = 0.28) compared to R²=0.59, RMSE = 0.45 for PLSR. Similar trends were observed for SMLR. The highest P retrieval accuracy was derived from WR spectra using SMLR (R²=0.64, RMSE = 0.067), while the traditional first derivative and continuum removal resulted in lower accuracy. Only when using PLSR did the first derivative result in a higher P retrieval accuracy (R²=0.47, RMSE = 0.07) than the WR spectra (R²=0.43, RMSE = 0.070). It was concluded that the water removal technique is a promising technique to minimize the perturbing effect of foliar water content when estimating grass nutrient concentrations.