Exploring uncertainty in remotely sensed data with parallel coordinate plots

The existence of uncertainty in classified remotely sensed data necessitates the application of enhanced techniques for identifying and visualizing the various degrees of uncertainty. This paper, therefore, applies the multidimensional graphical data analysis technique of parallel coordinate plots (PCP) to visualize the uncertainty in Landsat Thematic Mapper (TM) data classified by the Maximum Likelihood Classifier (MLC) and Fuzzy C-Means (FCM). The Landsat TM data are from the Yellow River Delta, Shandong Province, China. Image classification with MLC and FCM provides the probability vector and fuzzy membership vector of each pixel. Based on these vectors, the Shannon's entropy (S.E.) of each pixel is calculated. PCPs are then produced for each classification output. The PCP axes denote the posterior probability vector and fuzzy membership vector and two additional axes represent S.E. and the associated degree of uncertainty. The PCPs highlight the distribution of probability values of different land cover types for each pixel, and also reflect the status of pixels with different degrees of uncertainty. Brushing functionality is then added to PCP visualization in order to highlight selected pixels of interest. This not only reduces the visualization uncertainty, but also provides invaluable information on the positional and spectral characteristics of targeted pixels.     

Large area forest inventory using Landsat ETM+: A geostatistical approach

Large area forest inventory is important for understanding and managing forest resources and ecosystems. Remote sensing, the Global Positioning System (GPS), and geographic information systems (GIS) provide new opportunities for forest inventory. This paper develops a new systematic geostatistical approach for predicting forest parameters, using integrated Landsat 7 Enhanced Thematic Mapper Plus (ETM+) images, GPS, and GIS. Forest parameters, such as basal area, height, health conditions, biomass, or carbon, can be incorporated as a response variable, and the geostatistical approach can be used to predict parameter values for uninventoried points. Using basal area as the response and Landsat ETM+ images of pine stands in Georgia as auxiliary data, this approach includes univariate kriging (ordinary kriging and universal kriging) and multivariable kriging (co-kriging and regression kriging). The combination of bands 4, 3, and 2, as well as the combination of bands 5, 4, and 3, normalized difference vegetation index (NDVI), and principal components (PCs) were used in this study with co-kriging and regression kriging. Validation based on 200 randomly sampling points withheld field inventory was computed to evaluate the kriging performance and demonstrated that band combination 543 performed better than band combination 432, NDVI, and PCs. Regression kriging resulted in the smallest errors and the highest R-squared indicating the best geostatistical method for spatial predictions of pine basal area.     

Altering intertidal sediment topography: effects on biodiversity and ecosystem functioning

Sedimentary marine systems are often highly productive and perform important nutrient regeneration functions as they efficiently decompose organic material. In recent years the role of habitat effects and of species composition in ecosystem functioning has become of interest. Estuarine environments are frequently subject to considerable anthropogenic pressures whilst supporting a variety of habitats ranging from well sorted soft muds through biogenically stable sediments to highly mobile coarse sands. There is therefore considerable spatial complexity in habitat type and faunal composition. This study set out to observe the effects of altering the topographical habitat features of an estuarine mudflat on a range of porewater nutrient concentrations (NH₄⁺, NO₃⁻, NO₂⁻, PO₄³⁻ and SiO) collected from four depths (3, 5, 9, and 12 cm) and on faunal composition. Two treatments (Shelled Nets and Net Controls) were used to alter the topography from simple mud to a mussel shell crumble and were compared to un-manipulated Control areas. Sediment granulometry and organic matter content analyses alongside biological traits analysis of the fauna were also conducted.
Differences were observed in porewater nutrient concentrations between the Control and both netted treatments at 5 cm depth only; the species diversity and abundance were also different in the netted treatments compared to the Controls, although no difference between the two manipulated treatments were observed. The changes in faunal composition were attributed solely to the altered topography and the observed nutrient changes were attributed to the faunal alteration rather than the topographical manipulation.     

Second-order wave maker theory using force-feedback control. Part I: A new theory for regular wave generation

Second-order wave maker theory has long been established; the most extensive and detailed approach given by Schäffer [1996. Second-order wave maker theory for irregular waves. Ocean Engineering 23, 47–88]. However, all existing theories assume the wave paddle is driven by a position-feedback motion controller. Early research in the wave power field led to the design of a force-controlled absorbing wave machine [Salter, S., 1982. Absorbing wave-makers and wide tanks. In: Directional Wave Spectra Applications, pp. 185–200]. In addition to operating as an excellent absorber, this machine seemed to introduce very little spurious harmonic content when driven with a first-order command signal. The present paper provides a mathematical model for the operation of wave makers using force-feedback control and seeks to explain this apparent advantage. The model is developed to second-order so that a command signal compensating for the remaining spurious wave is also provided. Due to the complexity of the problem, the model has been limited to flap-type wave machines and the generation of regular waves. A variety of numerical tests in force-control mode have been conducted, indicating that the spurious wave content is greatly reduced when compared to the position-control mode. A separate experimental study validating the theory is presented in a part II paper by the same authors.     

Wintertime rates of net primary production and nitrate and ammonium uptake in the southern Benguela upwelling system

The elevated levels of primary productivity associated with eastern boundary currents are driven by nutrient- rich waters upwelled from depth, such that these regions are typically characterised by high rates of nitrate-fuelled phytoplankton growth. Production studies from the southern Benguela upwelling system (SBUS) tend to be biased towards the summer upwelling season, yet winter data are required to compute annual budgets and understand seasonal variability. Net primary production (NPP) and nitrate and ammonium uptake were measured concurrently at six stations in the SBUS in early winter. While euphotic zone NPP was highest at the stations nearest to the coast and declined with distance from the shore, a greater proportion was potentially exportable from open-ocean surface waters, as indicated by the higher specific nitrate uptake rates and f-ratios (ratio of nitrate uptake to total nitrogen consumption) at the stations located off the continental shelf. Near the coast, phytoplankton growth was predominantly supported by ammonium despite the high ambient nitrate concentrations. Along with ammonium concentrations as high as 3.6 µmol l^–1, this strongly suggests that nitrate uptake in the inshore SBUS, and by extension carbon drawdown, is inhibited by ammonium, at least in winter, although this has also been hypothesised for the summer.     

Preliminary assessment of the ATHENA/WFI non-X-ray background

We present a preliminary assessment of the non-X-ray background for the WFI on board ATHENA conducted at IAAT in the context of the collaborative background and radiation damage working group activities. Our main result is that in the baseline configuration originally assumed for the camera the requirement on the level of non-X-ray background could not be met. In light of the results of Geant4 simulations we propose and discuss a possible optimization of the camera design and pinpoint some open issues to be addressed in the next phase of investigation. One of these concerns the possible contribution to the non-X-ray background from soft protons and ions funneled to the focal plane through the optics. This is not quantified at this stage, here we just briefly report on our ongoing activities aimed at validating the mechanisms of proton scattering at grazing incidence.     

Variation in Coronal Activity from Solar Cycle 24 Minimum to Maximum Using Three-Dimensional Reconstructions of the Coronal Electron Density from STEREO/COR1

Three-dimensional electron density distributions in the solar corona are reconstructed for 100 Carrington rotations (CR 2054?–?2153) during 2007/03?–?2014/08 using the spherically symmetric method from polarized white-light observations with the inner coronagraph (COR1) onboard the twin Solar Terrestrial Relations Observatory (STEREO). These three-dimensional electron density distributions are validated by comparison with similar density models derived using other methods such as tomography and a magnetohydrodynamics (MHD) model as well as using data from the Solar and Heliospheric Observatory (SOHO)/Large Angle and Spectrometric Coronagraph (LASCO)-C2. Uncertainties in the estimated total mass of the global corona are analyzed based on differences between the density distributions for COR1-A and -B. Long-term variations of coronal activity in terms of the global and hemispheric average electron densities (equivalent to the total coronal mass) reveal a hemispheric asymmetry during the rising phase of Solar Cycle 24, with the northern hemisphere leading the southern hemisphere by a phase shift of 7?–?9 months. Using 14 CR (\(\approx13\)-month) running averages, the amplitudes of the variation in average electron density between Cycle 24 maximum and Cycle 23/24 minimum (called the modulation factors) are found to be in the range of 1.6?–?4.3. These modulation factors are latitudinally dependent, being largest in polar regions and smallest in the equatorial region. These modulation factors also show a hemispheric asymmetry: they are somewhat larger in the southern hemisphere. The wavelet analysis shows that the short-term quasi-periodic oscillations during the rising and maximum phases of Cycle 24 have a dominant period of 7?–?8 months. In addition, it is found that the radial distribution of the mean electron density for streamers at Cycle 24 maximum is only slightly larger (by \(\approx30\%\)) than at cycle minimum.     

Joint modelling of wave spectral parameters for extreme sea states

Characterising the dependence between extremes of wave spectral parameters such as significant wave height (H_S) and spectral peak period (T_P) is important in understanding extreme ocean environments and in the design and assessment of marine structures. For example, it is known that mean values of wave periods tend to increase with increasing storm intensity. Here we seek to characterise joint dependence in a straightforward manner, accessible to the ocean engineering community, using a statistically sound approach.Many methods of multivariate extreme value analyses are based on models which assume implicitly that in some joint tail region each parameter is either independent of or asymptotically dependent on other parameters; yet in reality the dependence structure in general is neither of these. The underpinning assumption of multivariate regular variation restricts these methods to estimation of joint regions in which all parameters are extreme; but regions where only a subset of parameters are extreme can be equally important for design. The conditional approach of Heffernan and Tawn (2004), similar in spirit to that of Haver (1985) but with better theoretical foundation, overcomes these difficulties.We use the conditional approach to characterise the dependence structure of H_S and T_P. The key elements of the procedure are: (1) marginal modelling for all parameters, (2) transformation of data to a common standard Gumbel marginal form, (3) modelling dependence between data for extremes of pairs of parameters using a form of regression, (4) simulation of long return periods to estimate joint extremes. We demonstrate the approach in application to measured and hindcast data from the Northern North Sea, the Gulf of Mexico and the North West Shelf of Australia. We also illustrate the use of data re-sampling techniques such as bootstrapping to estimate the uncertainty in marginal and dependence models and accommodate this uncertainty in extreme quantile estimation.We discuss the current approach in the context of other approaches to multivariate extreme value estimation popular in the ocean engineering community.     

The Eltanin asteroid impact: possible South Pacific palaeomegatsunami footprint and potential implications for the Pliocene–Pleistocene transition

Large asteroid impacts are rare, and those into the deep ocean are rarer still. The Eltanin asteroid impact around 2.51 ± 0.07 Ma occurred at a time of great climatic and geological change associated with the Pliocene–Pleistocene boundary. Numerical models of the event indicate that a megatsunami was generated, although there is debate concerning its magnitude and the region‐wide extent of its influence. We summarise the existing evidence for possible Eltanin megatsunami deposits in Antarctica, Chile and New Zealand, while also examining other potential sites from several locations, mainly around the South Pacific region. In reviewing these data we note that these events were unfolding at the same time as those associated with the Pliocene–Pleistocene boundary and, as such, most of the geological evidence from that time has a climatic interpretation. The potential climatic and geological ramifications of the Eltanin asteroid impact, however, have failed to be considered by most researchers studying this time period. Although we are not advocating that all geological activity at that time is connected with the Eltanin asteroid impact, it raises interesting questions about the role potentially played by such catastrophic events in contributing to or even triggering epochal transitions. Copyright © 2012 John Wiley & Sons, Ltd.