Modelling biomass of mountainous grasslands by including a species composition map

High mountain grasslands offer multiple goods and services to society but are severely threatened by improper land use practices such as abandonment or rapid intensification. In order to reduce abandonment and strengthen the common extensive agricultural practice a sustainable land use management of high mountain grasslands is needed. A spatially detailed yield assessment helps to identify possible meadows or, on the contrary, areas with a low carrying capacity in a region, making it easier to manage these sites. Such assessments are rarely available for remote and inaccessible areas. Remotely sensed vegetation indices are able to provide valuable information on grassland properties. These indices tend, however, to saturate for high biomass. This affects their applicability to assessments of high-yield grasslands.The main aim of this study was to model a spatially explicit grassland yield map and to test whether saturation issues can be tackled by consideration of plant species composition in the modelling process. The high mountain grassland of the subalpine belt (1800 – 2500 m a.s.l.) in the Kazbegi region, Greater Caucasus, Georgia, was chosen as test site for its strong species composition and yield gradients.We first modelled the species composition of the grassland described as metrically scaled gradients in the form of ordination axes by random forest regression. We then derived vegetation indices from Rapid Eye imagery, and topographic variables from a digital elevation model, which we used together with the multispectral bands as predictive variables. For comparison, we performed two yield models, one excluding the species composition maps and one including the species composition map as predictors. Moreover, we performed a third individual model, with species composition as predictors and a split dataset, to produce the final yield map.Three main grassland types were found in the vegetation analysis: Hordeum violaceum-meadows, Gentianella caucasea-grassland and Astragalus captiosus-grassland. The three random forest regression models for the ordination axes explained 64%, 33% and 46% of the variance in species composition. Independent validation of modelled ordination scores against a validation data set resulted in an R² of 0.64, 0.32 and 0.46 for the first, second and third axes, respectively. The model based on species composition resulted in a R² = 0.55, whereas the benchmark model showed weaker relationships between yield and the multispectral reflectance, vegetation indices, and topographical parameters (R² = 0.42). The final random forest yield model used to derive the yield map resulted in 62% variance explained and an R² = 0.64 between predicted and observed biomass. The results further indicate that high yields are generally difficult to predict with both models.The benefit of including a species composition map as a predictor variable for grassland yield lies in the preservation of ecologically meaningful features, especially the occurrence of high yielding vegetation type of Hordeum violaceum meadows is depicted accurately in the map. Even though we used a gradient based design, sharp boundaries or immediate changes in productivity were visible, especially in small structures such as arable fields or roads (Fig. 6b), making it a valuable tool for sustainable land use management. The saturation effect however, was mitigated by using species composition as predictor variables but is still present at high yields.     

An Algorithm for Calculation of n-Dimensional Nonlinear Dependencies – Behaviour Tests on Rats

The evaluation of the data obtained during the behaviour tests always leads to the problem of multiple correlation, very often with non-linear dependencies on the target. All mathematical and statistical procedures that have been used so far are based on the assumption of an equation for the desired correlation for which parameters and related statistical equivalents are determined eventually. The MODAK system applied here (MODAK = algorithms of modelling for the calculation of multi-dimensional non-linear mathematical models) breaks down a complex correlation into individual dependencies in a mathematical and statistical way and selects suitable equations for each of them independently and determines the corresponding parameters. The numerical example evaluates data of behaviour tests on rats. First results obtained on the correlations of various behaviour tests indicate both the possibility of selecting suitable tests independent of each other and a better interpretation of the observed patterns of behaviour taking into account the interrelations between the tests. In addition, MODAK is a method which can be applied as a matter of course in a general way to all cases which call for the reduction and analysis of data occurring in process and system analysis and in the evaluation of test results requiring statistical modelling. So far, MODAK applications range from engineering sciences to medicine.     

Estimating woody and herbaceous vegetation cover from time series satellite observations

In this paper we test a method to estimate the tree and grass vegetation cover over Australia from satellite-derived normalized difference vegetation index (NDVI) time series (monthly 1981–91, ≈5 km pixels) observations. The evergreen cover is assumed to track along the base of the NDVI time series, which is assumed to be equivalent to the woody vegetation cover. The base of the NDVI time series is estimated using modifications to a classical econometric model (i.e. time series is the sum of trend, seasonal and random components). Estimates of the average evergreen component during 1982–85 and 1986–89 were generally consistent with known vegetation distributions. Changes in evergreen cover were largely restricted to the south-west and south-east of Australia. Those changes were largely the result of differences in rainfall between the two periods. The proposed method for estimating woody vegetation cover is found to be generally robust. However, there are some regions where the grass (or pasture) is mostly evergreen. Some possible refinements are proposed to handle such cases.     

Estimation of uncertainties in X-ray refinement results by use of perturbed structures

The uncertainties in the refined parameters for a 1.5-A X-ray structure of carbon-monoxy (FeII) myoglobin are estimated by combining energy minimization with least-squares refinement against the X-ray data. The energy minimizations, done without reference to the X-ray data, provide perturbed structures which are used to restart conventional X-ray refinement. The resulting refined structures have the same, or better, R-factor and stereochemical parameters as the original X-ray structure, but deviate from it by 0.13 A rms for the backbone atoms and 0.31 A rms for the sidechain atoms. Atoms interacting with a disordered sidechain, Arg 45 CD3, are observed to have larger positional uncertainties. The uncertainty in the B-factors, within the isotropic harmonic motion approximation, is estimated to be 15%. The resulting X-ray structures are more consistent with the energy parameters used in simulations.     

Diurnal water storage in the stems of Picea sitchensis (Bong.) Carr.

Abstract. Two models of the relationship between diurnal variation in shoot water potential and transpiration in 14-year-old Picea sitchensis (Bong.) Carr. were compared. The first model was a physiologically based resistance-capacitance (R-C) analogue with its associated differential equations. The second was a non-physiological discrete-difference (D-D) or stochastic transfer function model. The RC model included only the effect of water storage in the phloem and bark while the D-D model implicity included all storage mechanisms. The R-C and D-D models explained similar fractions (62% and 68% respectively) of the variation in shoot water potential due to diurnal changes in transpiration rate. However, the D-D model had fewer parameters than the R-C model. The results from the D-D model showed that the resistance to flow from soil to shoots along the trunk, (RT), was 5 × 10³ MPa kg^-1s and the capacitance of the phloem and bark treated as a single store, (C_s), was 1.6 kg MPa^-1. It is suggested that the resistance to flow into storage (R_s) is much greater than R_T and can be disregarded. A non-linear version of the D-D model suggested [hat resistance to flow in the trunk increases with increasing transpiration rate.     

Hydro-sedimentology of the Johnstone River estuary

This paper examines the physical consequences of increased catchment sediment yields on the sediment budget and the hydrodynamic setting of the South Johnstone River estuary in North Queensland. A combined study involving hydrological monitoring, assessment of sediment sources, estimation of riverine sediment budget and hydro-sedimentological numerical modelling for estuarine sediment transport is currently underway. Initial field and laboratory observations indicate that the sediment delivery from highly erosion-prone sugar cane cultivations in the tropical catchment has increased dramatically during the last 10 years. This has subsequently given rise to elevated flood levels in both the lower and upper catchment areas, as well as significant modifications to the river bed morphology.     

Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison

A method for the prediction of hydrogen positions in proteins is presented. The method is based on the knowledge of the heavy atom positions obtained, for instance, from X-ray crystallography. It employs an energy minimization limited to the environment of the hydrogen atoms bound to a common heavy atom or to a single water molecule. The method is not restricted to proteins and can be applied without modification to nonpolar hydrogens and to nucleic acids. The method has been applied to the neutron diffraction structures of trypsin ribonuclease A, and bovine pancreatic trypsin inhibitor. A comparison of the constructed and the observed hydrogen positions shows few deviations except in situations in which several energetically similar conformations are possible. Analysis of the potential energy of rotation of Lys amino and Ser, Thr, Tyr hydroxyl groups reveals that the conformations of lowest intrinsic torsion energies are statistically favored in both the crystal and the constructed structures.