Assessment of Transverse Isotropy in Clinical-Level CT Images of Trabecular Bone Using the Gradient Structure Tensor

The aim of this study was to develop a GST-based methodology for accurately measuring the degree of transverse isotropy in trabecular bone. Using femoral sub-regions scanned in high-resolution peripheral QCT (HR-pQCT) and clinical-level-resolution QCT, trabecular orientation was evaluated using the mean intercept length (MIL) and the gradient structure tensor (GST) on the HR-pQCT and QCT data, respectively. The influence of local degree of transverse isotropy (DTI) and bone mineral density (BMD) was incorporated into the investigation. In addition, a power based model was derived, rendering a 1:1 relationship between GST and MIL eigenvalues. A specific DTI threshold (DTI thres) was found for each investigated size of region of interest (ROI), above which the estimate of major trabecular direction of the GST deviated no more than 30° from the gold standard MIL in 95% of the remaining ROIs (mean error: 16°). An inverse relationship between ROI size and DTI thres was found for discrete ranges of BMD. A novel methodology has been developed, where transversal isotropic measures of trabecular bone can be obtained from clinical QCT images for a given ROI size, DTI thres and power coefficient. Including DTI may improve future clinical QCT finite-element predictions of bone strength and diagnoses of bone disease.

Measurement of structural anisotropy in femoral trabecular bone using clinical-resolution CT images

Discrepancies in finite-element model predictions of bone strength may be attributed to the simplified modeling of bone as an isotropic structure due to the resolution limitations of clinical-level Computed Tomography (CT) data. The aim of this study is to calculate the preferential orientations of bone (the principal directions) and the extent to which bone is deposited more in one direction compared to another (degree of anisotropy). Using 100 femoral trabecular samples, the principal directions and degree of anisotropy were calculated with a Gradient Structure Tensor (GST) and a Sobel Structure Tensor (SST) using clinical-level CT. The results were compared against those calculated with the gold standard Mean-Intercept-Length (MIL) fabric tensor using micro-CT. There was no significant difference between the GST and SST in the calculation of the main principal direction (median error=28°), and the error was inversely correlated to the degree of transverse isotropy (r=−0.34, p<0.01). The degree of anisotropy measured using the structure tensors was weakly correlated with the MIL-based measurements (r=0.2, p<0.001). Combining the principal directions with the degree of anisotropy resulted in a significant increase in the correlation of the tensor distributions (r=0.79, p<0.001). Both structure tensors were robust against simulated noise, kernel sizes, and bone volume fraction. We recommend the use of the GST because of its computational efficiency and ease of implementation. This methodology has the promise to predict the structural anisotropy of bone in areas with a high degree of anisotropy, and may improve the in vivo characterization of bone.

The importance of biotic interactions and local adaptation for plant response to environmental changes: field evidence along an elevational gradient

Predicting the response of species to environmental changes is a great and on-going challenge for ecologists, and this requires a more in-depth understanding of the importance of biotic interactions and the population structuration in the landscape. Using a reciprocal transplantation experiment, we tested the response of five species to an elevational gradient. This was combined to a neighbour removal treatment to test the importance of local adaptation and biotic interactions. The trait studied was performance measured as survival and biomass. Species response varied along the elevational gradient, but with no consistent pattern. Performance of species was influenced by environmental conditions occurring locally at each site, as well as by positive or negative effects of the surrounding vegetation. Indeed, we observed a shift from competition for biomass to facilitation for survival as a response to the increase in environmental stress occurring in the different sites. Unlike previous studies pointing out an increase of stress along the elevation gradient, our results supported a stress gradient related to water availability, which was not strictly parallel to the elevational gradient. For three of our species, we observed a greater biomass production for the population coming from the site where the species was dominant (central population) compared to population sampled at the limit of the distribution (marginal population). Nevertheless, we did not observe any pattern of local adaptation that could indicate adaptation of populations to a particular habitat. Altogether, our results highlighted the great ability of plant species to cope with environmental changes, with no local adaptation and great variability in response to local conditions. Our study confirms the importance of taking into account biotic interactions and population structure occurring at local scale in the prediction of communities’ responses to global environmental changes.

Improved methods for measuring forest landscape structure: LiDAR complements field-based habitat assessment

Conservation and monitoring of forest biodiversity requires reliable information about forest structure and composition at multiple spatial scales. However, detailed data about forest habitat characteristics across large areas are often incomplete due to difficulties associated with field sampling methods. To overcome this limitation we employed a nationally available light detection and ranging (LiDAR) remote sensing dataset to develop variables describing forest landscape structure across a large environmental gradient in Switzerland. Using a model species indicative of structurally rich mountain forests (hazel grouse Bonasa bonasia), we tested the potential of such variables to predict species occurrence and evaluated the additional benefit of LiDAR data when used in combination with traditional, sample plot-based field variables. We calibrated boosted regression trees (BRT) models for both variable sets separately and in combination, and compared the models’ accuracies. While both field-based and LiDAR models performed well, combining the two data sources improved the accuracy of the species’ habitat model. The variables retained from the two datasets held different types of information: field variables mostly quantified food resources and cover in the field and shrub layer, LiDAR variables characterized heterogeneity of vegetation structure which correlated with field variables describing the understory and ground vegetation. When combined with data on forest vegetation composition from field surveys, LiDAR provides valuable complementary information for encompassing species niches more comprehensively. Thus, LiDAR bridges the gap between precise, locally restricted field-data and coarse digital land cover information by reliably identifying habitat structure and quality across large areas.

Vertical and Horizontal Vegetation Structure across Natural and Modified Habitat Types at Mount Kilimanjaro

In most habitats, vegetation provides the main structure of the environment. This complexity can facilitate biodiversity and ecosystem services. Therefore, measures of vegetation structure can serve as indicators in ecosystem management. However, many structural measures are laborious and require expert knowledge. Here, we used consistent and convenient measures to assess vegetation structure over an exceptionally broad elevation gradient of 866–4550m above sea level at Mount Kilimanjaro, Tanzania. Additionally, we compared (human)-modified habitats, including maize fields, traditionally managed home gardens, grasslands, commercial coffee farms and logged and burned forests with natural habitats along this elevation gradient. We distinguished vertical and horizontal vegetation structure to account for habitat complexity and heterogeneity. Vertical vegetation structure (assessed as number, width and density of vegetation layers, maximum canopy height, leaf area index and vegetation cover) displayed a unimodal elevation pattern, peaking at intermediate elevations in montane forests, whereas horizontal structure (assessed as coefficient of variation of number, width and density of vegetation layers, maximum canopy height, leaf area index and vegetation cover) was lowest at intermediate altitudes. Overall, vertical structure was consistently lower in modified than in natural habitat types, whereas horizontal structure was inconsistently different in modified than in natural habitat types, depending on the specific structural measure and habitat type. Our study shows how vertical and horizontal vegetation structure can be assessed efficiently in various habitat types in tropical mountain regions, and we suggest to apply this as a tool for informing future biodiversity and ecosystem service studies.