Holocene sediments of Wistari Reef: towards a global quantification of coral reef related neritic sedimentation in the Holocene

Wistari Reef. within the southern Great Barrier Reef. is a shallow coral reef platform featuring a very clearly defined leeward accretionary wedge of carbonate sediments. The total global area of shallowly submerged coral reef has been quantified as 255 000 km(2). The question then becomes: What additional area of sediment of significant thickness is associated with the measured shallow reef areas T At Wistari Reef, the leeward sedimentary wedge has an area and a thickness that are roughly equal to the Holocene sediments that have accumulated on the platform. Several important observations can be made from these data. Firstly. the area of significant neritic carbonate sedimentation ( > 1 m/ka) associated with coral reefs is near 500000 km(2). Secondly, the production rate of neritic carbonates at Wistari Reef is almost 50%, less than the accumulation rate needed to obtain the volume of Holocene reef sediments observed. This implies that both production and accumulation neritic carbonate must have been more than a factor of two higher in the early to mid Holocene. (C) 2001 Elsevier Science B.V. All rights reserved.

Integrating high and moderate spatial resolution image data to estimate forest age structure

The collection of spatial information to quantify changes to the state and condition of the environment is a fundamental component of conservation or sustainable utilization of tropical and subtropical forests, Age is an important structural attribute of old-growth forests influencing biological diversity in Australia eucalypt forests. Aerial photograph interpretation has traditionally been used for mapping the age and structure of forest stands. However this method is subjective and is not able to accurately capture fine to landscape scale variation necessary for ecological studies. Identification and mapping of fine to landscape scale vegetative structural attributes will allow the compilation of information associated with Montreal Process indicators lb and ld, which seek to determine linkages between age structure and the diversity and abundance of forest fauna populations. This project integrated measurements of structural attributes derived from a canopy-height elevation model with results from a geometrical-optical/spectral mixture analysis model to map forest age structure at a landscape scale. The availability of multiple-scale data allows the transfer of high-resolution attributes to landscape scale monitoring. Multispectral image data were obtained from a DMSV (Digital Multi-Spectral Video) sensor over St Mary's State Forest in Southeast Queensland, Australia. Local scene variance levels for different forest tapes calculated from the DMSV data were used to optimize the tree density and canopy size output in a geometric-optical model applied to a Landsat Thematic Mapper (TU) data set. Airborne laser scanner data obtained over the project area were used to calibrate a digital filter to extract tree heights from a digital elevation model that was derived from scanned colour stereopairs. The modelled estimates of tree height, crown size, and tree density were used to produce a decision-tree classification of forest successional stage at a landscape scale. The results obtained (72% accuracy), were limited in validation, but demonstrate potential for using the multi-scale methodology to provide spatial information for forestry policy objectives (ie., monitoring forest age structure).

Sediment yield and impacts from river catchments to the Great Barrier Reef lagoon

Land use intensification is estimated to result in an overall increase in sediment delivery to the Great Barrier Reef lagoon by a factor of approximately four. Modelling suggests that, following land use intensification, croplands cause the greatest increase of sediment yield and sediment concentration, whereas erosion of grazing land is the main contemporary source of sediments, primarily owing to the large spatial extent of this land use. The spatial pattern of sediment yield to the coast after land use intensification is strongly correlated with the pattern under natural conditions, although the greatest increase is estimated to have occurred in the wet-dry catchments. Sediment transport and resuspension processes have led to the development of a strongly sediment-partitioned shelf, with modern mud-rich sediments almost exclusively restricted to the inner and inner-middle shelf, northward-facing embayments and in the lee of headlands. Elevated sediment concentrations increase the potential transport rates of nutrients and other pollutants. Whether increased sediment supply to the coastal zone has impacted on reefs remains a point of contention. More sediment load data need to be collected and analysed in order to make detailed estimates of catchment yields and establish the possible sediment impact on the Great Barrier Reef.

High-frequency climatic oscillations recorded in a Holocene coral reef at Leizhou Peninsula, South China sea

A detailed study of the Goniopora reef profile at Dengloujiao, Xuwen County, Leizhou Peninsula, the northern coast of the South China Sea suggests that a series of high-frequency, large-amplitude and abrupt cold events occurred during the Holocene Hypsithermal, an unusual phenomenon termed Leizhou Events in this paper. This period (corresponding to C-14 age of 6.2 -6.7 kaBP or calendar age of 6.7-7.2 kaBP), when the climatic conditions were ideal for coral. reefs to develop, can be divided into at least nine stages. Each stage (or called a climate optimum), lasting about 20 to 50 a, was terminated by an abrupt cold nap and (or) a sea-level lowering event in winter, leading to widespread emergence and death of the Goniopora corals, and growth discontinuities on the coral surface. Such a cyclic process resulted in the creation of a > 4m thick Goniopora reef flat. During this period, the crust subsided periodically but the sea level was rising. The reef profile provides valuable archives for the study of decadal-scale mid-Holocene climatic oscillations in the tropical area of South China. Our results provide new evidence for high-frequency climate instability in the Holocene Hypsithermal, and challenge the traditional understanding of Holocene climate.

DMSP in corals and benthic algae from the Great Barrier Reef

In this study the first measurements of DMSP in six species of corals and ten species of benthic algae collected from four coral reefs in the Great Barrier Reef are reported, together with DMSP measurements made on cultured zooxanthellae. Concentrations ranged from 21 to 3831 (mean=743) fmol DMSP zooxanthellae(-1) in corals, 0.16 to 2.96 nmol DMSP cm(-2) (mean=90) for benthic macroalgae, and 48-285 fmol DMSP zooxanthellae(-1) (mean=153) for cultured zooxanthellae. The highest concentrations of DMSP in corals occurred in Acropora formosa (mean= 371 fmol DMSP zooxanthellae(-1)) and Acropora palifera (mean=3341 fmol DMSP zooxanthellae(-1)) with concentrations in A. palifera the highest DMSP concentrations reported in corals examined to date. As well as inter-specific differences in DMSP, intra-specific variation was also observed. Adjacent colonies of A. formosa that are known to have different thermal bleaching thresholds and morphologically distinct zooxanthellae, were also observed to have different DMSP concentrations, with the zooxanthellae in the colony that bleached containing DMSP at an average concentration of 436 finol zooxanthellae(-1), whilst the non-bleaching colony contained DMSP at an average concentration of 171 finol zooxanthellae(-1). The results of the present study have been used to calculate the area normalized DMSP concentrations in benthic algae (mean=0.015 mmol m(-2)) and corals (mean=2.22 mmol m(-2)) from the GBR. This data indicates that benthic algae and corals are a significant reservoir of DMSP in GBR waters. (C) 2002 Published by Elsevier Science Ltd.

A quantitative state and transition model for the Mitchell grasslands of central western Queensland

Concerns of reduced productivity and land degradation in the Mitchell grasslands of central western Queensland were addressed through a range monitoring program to interpret condition and trend. Botanical and eclaphic parameters were recorded along piosphere and grazing gradients, and across fenceline impact areas, to maximise changes resulting from grazing. The Degradation Gradient Method was used in conjunction with State and Transition Models to develop models of rangeland dynamics and condition. States were found to be ordered along a degradation gradient, indicator species developed according to rainfall trends and transitions determined from field data and available literature. Astrebla spp. abundance declined with declining range condition and increasing grazing pressure, while annual grasses and forbs increased in dominance under poor range condition. Soil erosion increased and litter decreased with decreasing range condition. An approach to quantitatively define states within a variable rainfall environment based upon a time-series ordination analysis is described. The derived model could provide the interpretive framework necessary to integrate on-ground monitoring, remote sensing and geographic information systems to trace states and transitions at the paddock scale. However, further work is needed to determine the full catalogue of states and transitions and to refine the model for application at the paddock scale.

Photosynthetic responses of the coral Montipora digitata to cold temperature stress

Coral bleaching events have become more frequent and widespread, largely due to elevated sea surface temperatures. Global climate change could lead to increased variability of sea surface temperatures, through influences on climate systems, e.g. El Nino Southern Oscillation (ENSO). Field observations in 1999, following a strong ENSO, revealed that corals bleached in winter after unusually cold weather. To explore the basis for these observations, the photosynthetic responses of the coral species Montipora digitata Studer were investigated in a series of temperature and light experiments. Small replicate coral colonies were exposed to ecologically relevant lower temperatures for varying durations and under light regimes that ranged from darkness to full sunlight. Photosynthetic efficiency was analyzed using a pulse amplitude modulated (PAM) fluorometer (F-0, F-m, F-v/F-m), and chlorophyll a (chl a) content and symbiotic dinoflagellate density were analyzed with spectrophotometry and microscopy, respectively. Cold temperature stress had a negative impact on M digitata colonies indicated by decreased photosynthetic efficiency (F-v/F-m), loss of symbiotic dinoflagellates and changes in photosynthetic pigment concentrations. Corals in higher light regimes were more susceptible to cold temperature stress, Moderate cold stress resulted in photoacclimatory responses, but severe cold stress resulted in photodamage, bleaching and increased mortality. Responses to cold temperature stress of M digitata appeared similar to that observed in corals exposed to warmer than normal temperatures, suggesting a common mechanism. The results of this study suggest that corals and coral reefs may also be impacted by exposure to cold as well as warm temperature extremes as climate change occurs.

Using different satellite imagery and classification techniques to assess the contribution of trees outside forests in the municipality of Maringá, Brazil

Forest cover of the Maringá municipality, located in northern Parana State, was mapped in this study. Mapping was carried out by using high-resolution HRC sensor imagery and medium resolution CCD sensor imagery from the CBERS satellite. Images were georeferenced and forest vegetation patches (TOFs - trees outside forests) were classified using two methods of digital classification: reflectance-based or the digital number of each pixel, and object-oriented. The areas of each polygon were calculated, which allowed each polygon to be segregated into size classes. Thematic maps were built from the resulting polygon size classes and summary statistics generated from each size class for each area. It was found that most forest fragments in Maringá were smaller than 500 m². There was also a difference of 58.44% in the amount of vegetation between the high-resolution imagery and medium resolution imagery due to the distinct spatial resolution of the sensors. It was concluded that high-resolution geotechnology is essential to provide reliable information on urban greens and forest cover under highly human-perturbed landscapes.

Dinossáurios Eocretácicos de Lagosteiros

According to an ancient folkloric legend, Our Lady, stepping down from the sea, would have rided on a mule to the platform above the cliffs named Pedra da Mua at Lagosteiros'bay, near Espichel cape. Mule's footprints, regarded by fishermen as evidence, would be clearly recognizable on exposed surfaces of the rocks. Indeed there are footprints but from Dinosaurs of latest Jurassic, Portlandian age, this spectacular locality being specially rich in giant Sauropod tracks (that have seldom been found elsewhere in Europe). As we proceeded to its study, another locality with Dinosaur footprints, Lower Cretaceous (Hauterivian) in age, was found on the northern cliffs at Lagosteiros. It is probably the richest one in european Lower Cretaceous and the only of this age known in Portugal, so we decided to give priority to its study. Dinosaur tracks have been printed on calciclastic sands in a lagoonal environment protected by fringing coral reefs. There have been emersion episodes; beaches were frequented by Dinosaurs. Later on, the marine barremian ingression restablished a gulf and such animals could not come here any more. Under a paleogeographical viewpoint, the evidence of a marine regression near the end of Hauterivian is to be remarked. Five types of tracks and footprints have been recognized: - Neosauropus lagosteirensis, new morphogenus and species, tracks from a giant Sauropod, perhaps from Camarasaurus; with its proportions the total length of the author would be about 15,5 m. These are the only Sauropod tracks known till now in Europe's Lower Cretaceous. - tracks from a not so big quadruped, maybe a Sauropod (young individual?); however it is not impossible that they were produced by Stegosaurians or Ankylosaurians. -Megalosauropus (?Eutynichnium) gomesi new morphospecies, four Theropod tracks most probably produced by megalosaurs. - Iguanodon sp., represented by some footprints and specially by a set corresponding to the feet and tail from an individual standing in a rest position. - problematical, quite small-sized biped (maybe an Ornithopod related to Camptosaurus). Evidence points to a richer fauna than that known in barremian "Dinosaur sandstones" from a nearby locality, Boca do Chapim. Lagosteiros' association clearly indicates the predominance of herbivores, which required large amounts of vegetable food in the neighbourhood. This is an indirect evidence of the vegetal wealth, also suggested by associations of plant macrofossils, polen and spores found in early Cretaceous sediments at the same region. The relatively high proportion of Theropoda is related to the wealth of the whole fauna, which comprised a lot of the prey needed by such powerful flesh-eaters. The evidence, as a whole, points out to a warm and moist climate. All the tracks whose direction could be measured are directed to the southern quadrants, this being confirmed by the approximative direction of other footprints. Massive displacements (migration?) could take place during a brief emersion episode. This may result from the ingression of barremian seas, flooding the region and restablishing here a small gulf. Even if the arrival of the waters damaged certain footprints it has not destroyed them completely, thus allowing the preservation of such evidence from a remote past.

Land use/ land cover change and its impact on soil erosion process in Begnas Tal Rupa Tal watershed using geospatial tools, Kaski district, Nepal

Dissertation submitted in partial fulfilment of the requirements for the Degree of Master of Science in Geospatial Technologies

Remote experimentation: integrating research, education, and industrial application

This paper presents a low-cost scaled model of a silo for drying and airing cereal grains. It allows the control and monitor of several parameters associated to the silo's operation, through a remote accessible infrastructure. The scaled model consists of a 2.50 m wide × 2.10 m long plant with all control and monitor capacities provided by micro-Web servers. An application running on the micro-Web servers enables storing all parameters in a data basis for later analysis. The implemented model aims to support a remote experimentation facility for technological education, research-oriented tutorials, and industrial applications. Given the low-cost requirement, this remote facility can be easily replicated in other institutions to support a network of remote labs, which encompasses the concurrent access of several users (e.g. students).

GPU implementation of a hyperspectral coded aperture algorithm for compressive sensing

This paper presents a new parallel implementation of a previously hyperspectral coded aperture (HYCA) algorithm for compressive sensing on graphics processing units (GPUs). HYCA method combines the ideas of spectral unmixing and compressive sensing exploiting the high spatial correlation that can be observed in the data and the generally low number of endmembers needed in order to explain the data. The proposed implementation exploits the GPU architecture at low level, thus taking full advantage of the computational power of GPUs using shared memory and coalesced accesses to memory. The proposed algorithm is evaluated not only in terms of reconstruction error but also in terms of computational performance using two different GPU architectures by NVIDIA: GeForce GTX 590 and GeForce GTX TITAN. Experimental results using real data reveals signficant speedups up with regards to serial implementation.

Parallel hyperspectral compressive sensing method on GPU

Remote hyperspectral sensors collect large amounts of data per flight usually with low spatial resolution. It is known that the bandwidth connection between the satellite/airborne platform and the ground station is reduced, thus a compression onboard method is desirable to reduce the amount of data to be transmitted. This paper presents a parallel implementation of an compressive sensing method, called parallel hyperspectral coded aperture (P-HYCA), for graphics processing units (GPU) using the compute unified device architecture (CUDA). This method takes into account two main properties of hyperspectral dataset, namely the high correlation existing among the spectral bands and the generally low number of endmembers needed to explain the data, which largely reduces the number of measurements necessary to correctly reconstruct the original data. Experimental results conducted using synthetic and real hyperspectral datasets on two different GPU architectures by NVIDIA: GeForce GTX 590 and GeForce GTX TITAN, reveal that the use of GPUs can provide real-time compressive sensing performance. The achieved speedup is up to 20 times when compared with the processing time of HYCA running on one core of the Intel i7-2600 CPU (3.4GHz), with 16 Gbyte memory.

Unmixing hyperspectral data: independent and dependent component analysis

The development of high spatial resolution airborne and spaceborne sensors has improved the capability of ground-based data collection in the fields of agriculture, geography, geology, mineral identification, detection [2, 3], and classification [4–8]. The signal read by the sensor from a given spatial element of resolution and at a given spectral band is a mixing of components originated by the constituent substances, termed endmembers, located at that element of resolution. This chapter addresses hyperspectral unmixing, which is the decomposition of the pixel spectra into a collection of constituent spectra, or spectral signatures, and their corresponding fractional abundances indicating the proportion of each endmember present in the pixel [9, 10]. Depending on the mixing scales at each pixel, the observed mixture is either linear or nonlinear [11, 12]. The linear mixing model holds when the mixing scale is macroscopic [13]. The nonlinear model holds when the mixing scale is microscopic (i.e., intimate mixtures) [14, 15]. The linear model assumes negligible interaction among distinct endmembers [16, 17]. The nonlinear model assumes that incident solar radiation is scattered by the scene through multiple bounces involving several endmembers [18]. Under the linear mixing model and assuming that the number of endmembers and their spectral signatures are known, hyperspectral unmixing is a linear problem, which can be addressed, for example, under the maximum likelihood setup [19], the constrained least-squares approach [20], the spectral signature matching [21], the spectral angle mapper [22], and the subspace projection methods [20, 23, 24]. Orthogonal subspace projection [23] reduces the data dimensionality, suppresses undesired spectral signatures, and detects the presence of a spectral signature of interest. The basic concept is to project each pixel onto a subspace that is orthogonal to the undesired signatures. As shown in Settle [19], the orthogonal subspace projection technique is equivalent to the maximum likelihood estimator. This projection technique was extended by three unconstrained least-squares approaches [24] (signature space orthogonal projection, oblique subspace projection, target signature space orthogonal projection). Other works using maximum a posteriori probability (MAP) framework [25] and projection pursuit [26, 27] have also been applied to hyperspectral data. In most cases the number of endmembers and their signatures are not known. Independent component analysis (ICA) is an unsupervised source separation process that has been applied with success to blind source separation, to feature extraction, and to unsupervised recognition [28, 29]. ICA consists in finding a linear decomposition of observed data yielding statistically independent components. Given that hyperspectral data are, in given circumstances, linear mixtures, ICA comes to mind as a possible tool to unmix this class of data. In fact, the application of ICA to hyperspectral data has been proposed in reference 30, where endmember signatures are treated as sources and the mixing matrix is composed by the abundance fractions, and in references 9, 25, and 31–38, where sources are the abundance fractions of each endmember. In the first approach, we face two problems: (1) The number of samples are limited to the number of channels and (2) the process of pixel selection, playing the role of mixed sources, is not straightforward. In the second approach, ICA is based on the assumption of mutually independent sources, which is not the case of hyperspectral data, since the sum of the abundance fractions is constant, implying dependence among abundances. This dependence compromises ICA applicability to hyperspectral images. In addition, hyperspectral data are immersed in noise, which degrades the ICA performance. IFA [39] was introduced as a method for recovering independent hidden sources from their observed noisy mixtures. IFA implements two steps. First, source densities and noise covariance are estimated from the observed data by maximum likelihood. Second, sources are reconstructed by an optimal nonlinear estimator. Although IFA is a well-suited technique to unmix independent sources under noisy observations, the dependence among abundance fractions in hyperspectral imagery compromises, as in the ICA case, the IFA performance. Considering the linear mixing model, hyperspectral observations are in a simplex whose vertices correspond to the endmembers. Several approaches [40–43] have exploited this geometric feature of hyperspectral mixtures [42]. Minimum volume transform (MVT) algorithm [43] determines the simplex of minimum volume containing the data. The MVT-type approaches are complex from the computational point of view. Usually, these algorithms first find the convex hull defined by the observed data and then fit a minimum volume simplex to it. Aiming at a lower computational complexity, some algorithms such as the vertex component analysis (VCA) [44], the pixel purity index (PPI) [42], and the N-FINDR [45] still find the minimum volume simplex containing the data cloud, but they assume the presence in the data of at least one pure pixel of each endmember. This is a strong requisite that may not hold in some data sets. In any case, these algorithms find the set of most pure pixels in the data. Hyperspectral sensors collects spatial images over many narrow contiguous bands, yielding large amounts of data. For this reason, very often, the processing of hyperspectral data, included unmixing, is preceded by a dimensionality reduction step to reduce computational complexity and to improve the signal-to-noise ratio (SNR). Principal component analysis (PCA) [46], maximum noise fraction (MNF) [47], and singular value decomposition (SVD) [48] are three well-known projection techniques widely used in remote sensing in general and in unmixing in particular. The newly introduced method [49] exploits the structure of hyperspectral mixtures, namely the fact that spectral vectors are nonnegative. The computational complexity associated with these techniques is an obstacle to real-time implementations. To overcome this problem, band selection [50] and non-statistical [51] algorithms have been introduced. This chapter addresses hyperspectral data source dependence and its impact on ICA and IFA performances. The study consider simulated and real data and is based on mutual information minimization. Hyperspectral observations are described by a generative model. This model takes into account the degradation mechanisms normally found in hyperspectral applications—namely, signature variability [52–54], abundance constraints, topography modulation, and system noise. The computation of mutual information is based on fitting mixtures of Gaussians (MOG) to data. The MOG parameters (number of components, means, covariances, and weights) are inferred using the minimum description length (MDL) based algorithm [55]. We study the behavior of the mutual information as a function of the unmixing matrix. The conclusion is that the unmixing matrix minimizing the mutual information might be very far from the true one. Nevertheless, some abundance fractions might be well separated, mainly in the presence of strong signature variability, a large number of endmembers, and high SNR. We end this chapter by sketching a new methodology to blindly unmix hyperspectral data, where abundance fractions are modeled as a mixture of Dirichlet sources. This model enforces positivity and constant sum sources (full additivity) constraints. The mixing matrix is inferred by an expectation-maximization (EM)-type algorithm. This approach is in the vein of references 39 and 56, replacing independent sources represented by MOG with mixture of Dirichlet sources. Compared with the geometric-based approaches, the advantage of this model is that there is no need to have pure pixels in the observations. The chapter is organized as follows. Section 6.2 presents a spectral radiance model and formulates the spectral unmixing as a linear problem accounting for abundance constraints, signature variability, topography modulation, and system noise. Section 6.3 presents a brief resume of ICA and IFA algorithms. Section 6.4 illustrates the performance of IFA and of some well-known ICA algorithms with experimental data. Section 6.5 studies the ICA and IFA limitations in unmixing hyperspectral data. Section 6.6 presents results of ICA based on real data. Section 6.7 describes the new blind unmixing scheme and some illustrative examples. Section 6.8 concludes with some remarks.