Improving the Elevation Accuracy of Digital Elevation Models: A Comparison of Some Error Detection Procedures

The widespread availability of powerful desktop computers, easy‐to‐use software tools and geographic datasets have raised the quality problem of input data to be a crucial one. Even though accuracy has been a concern in every serious application, there are no general tools for its improvement. Some particular ones exist, however, and some results are presented here for a particular case of quantitative raster data: Digital Elevation Models (DEM). Two procedures designed to detect anomalous values (also named gross errors, outliers or blunders) in DEMs, but valid also for other quantitative raster datasets, were tested. A DEM with elevations varying from 181 to 1044 m derived from SPOT data has been used as a contaminated sample, while a manually derived DEM obtained from aerial photogrammetry was regarded as the ground truth to allow a direct performance comparison for the methods with real errors. It is assumed that a “better” value can be measured or obtained through some methodology once an outlier location is suggested. The options are different depending upon the user (DEM producers might go to the original data and make another reading, while end users might use interpolation). Both choices were considered in this experiment. Preliminary results show that for the available dataset, the accuracy might be improved to some extent with very little effort. Effort is defined here as the percentage of points suggested by the methodology in relation with its total number: thus 100 per cent effort implies that all points have been checked. The method proposed by López (1997) gave poor results, because it has been designed for errors with low spatial autocorrelation (which is not the case here). A modified version was then designed and compared with the method proposed by Felicísimo (1994). The three procedures can be applied both for error detection during DEM generation and by end users, and they might be of use for other quantitative raster data. The choice of the best methodology is different depending on the effort involved. The conclusions have been derived for a photogrammetrically obtained DEM; other production procedures might lead to different results.     

Tectonics of an extinct ridge-transform intersection,Drake Passage (Antarctica)

New swath bathymetric, multichannel seismic and magnetic data reveal the complexity of the intersection between the extinct West Scotia Ridge (WSR) and the Shackleton Fracture Zone (SFZ), a first-order NW-SE trending high-relief ridge cutting across the Drake Passage. The SFZ is composed of shallow, ridge segments and depressions, largely parallel to the fracture zone with an `en echelon' pattern in plan view. These features are bounded by tectonic lineaments, interpreted as faults. The axial valley of the spreading center intersects the fracture zone in a complex area of deformation, where N120° E lineaments and E–W faults anastomose on both sides of the intersection. The fracture zone developed within an extensional regime, which facilitated the formation of oceanic transverse ridges parallel to the fracture zone and depressions attributed to pull-apart basins, bounded by normal and strike-slip faults.On the multichannel seismic (MCS) profiles, the igneous crust is well stratified, with numerous discontinuous high-amplitude reflectors and many irregular diffractions at the top, and a thicker layer below. The latter has sparse and weak reflectors, although it locally contains strong, dipping reflections. A bright, slightly undulating reflector observed below the spreading center axial valley at about 0.75 s (twt) depth in the igneous crust is interpreted as an indication of the relict axial magma chamber. Deep, high-amplitude subhorizontal and slightly dipping reflections are observed between 1.8 and 3.2 s (twt) below sea floor, but are preferentially located at about 2.8–3.0 s (twt) depth. Where these reflections are more continuous they may represent the Mohorovicic seismic discontinuity. More locally, short (2–3 km long), very high-amplitude reflections observed at 3.6 and 4.3 s (twt) depth below sea floor are attributed to an interlayered upper mantle transition zone. The MCS profiles also show a pattern of regularly spaced, steep-inclined reflectors, which cut across layers 2 and 3 of the oceanic crust. These reflectors are attributed to deformation under a transpressional regime that developed along the SFZ, shortly after spreading ceased at the WSR. Magnetic anomalies 5 to 5 E may be confidently identified on the flanks of the WSR. Our spreading model assumes slow rates (ca. 10–20 mm/yr), with slight asymmetries favoring the southeastern flank between 5C and 5, and the northwestern flank between 5 and extinction. The spreading rate asymmetry means that accretion was slower during formation of the steeper, shallower, southeastern flank than of the northwestern flank.     

Galaxy destruction and diffuse light in clusters

Deep images of the Centaurus and Coma clusters reveal two spectacular arcs of diffuse light that stretch for over 100 kpc, yet are just a few kiloparsecs wide. At a surface brightness of m _b ∼27–28 mag arcsec⁻² , the Centaurus arc is the most striking example known of structure in the diffuse light component of a rich galaxy cluster. We use numerical simulations to show that the Centaurus feature can be reproduced by the tidal debris of a spiral galaxy that has been tidally disrupted by the gravitational potential of NGC 4709. The surface brightness and narrow dimensions of the diffuse light suggest that the disc was corotating with its orbital path past pericentre. Features this prominent in clusters will be relatively rare, although at fainter surface brightness levels the diffuse light will reveal a wealth of structure. Deeper imaging surveys may be able to trace this feature for several times its presently observed extent, and somewhere along the tidal debris, a fraction of the original stellar component of the disc will remain bound, but transformed into a faint spheroidal galaxy. It should be possible to confirm the galactic origin of the Centaurus arc by observing planetary nebulae along its length with redshifts close to that of NGC 4709.     

Lithium Abundances in the Atmospheres of SLR C-Giants WZ Cas and WX Cyg from Resonance and Subordinate Li I Lines

Lithium abundances in the atmospheres of the super Li-rich C-giants WZ Cas and WX Cyg are derived by the spectral synthesis technique using the Li I resonance line at λ670.8 nm and three subordinate lines at λλ 812.6, 610.4 and 497.2 nm. The differences between the Li abundances derived from the λ670.8 nm line and the λλ 497.2, 812.6 nm lines do not exceed ±0.5 dex. The lithium line at λ610.4 nm provides typically lower abundances than the resonance line (by ≈ 1 dex). The mean LTE and NLTE Li abundances from three Li I lines (excluding λ610.4 nm) are 4.7, 4.9 for WZ Cas, and 4.6, 4.8 for WX Cyg, respectively. This revised version was published online in July 2006 with corrections to the Cover Date.     

Soil liquefaction during the Arequipa Mw 8.4, June 23, 2001 earthquake, southern coastal Peru

The Arequipa June 23, 2001, earthquake with a moment magnitude of Mw 8.4 struck southern Peru, northern Chile and western Bolivia. This shallow (29 km deep) interplate event, occurring in the coupled zone of the Nazca subduction next to the southeast of the subducting Nazca ridge, triggered very localized but widely outspread soil liquefaction. Although sand blows and lateral spreading of river banks and road bridge abutments were observed 390 km away from the epicenter in the southeast direction (nearing the town of Tacna, close to the Chile border), liquefaction features were only observed in major river valleys and delta and coastal plains in the meizoseismal area. This was strongly controlled by the aridity along the coastal strip of Southern Peru. From the sand blow distribution along the coastal area, a first relationship of isolated sand blow diameter versus epicentral distance for a single event is ever proposed. The most significant outcome from this liquefaction field reconnaissance is that energy propagation during the main June 23, 2001, event is further supported by the distribution and size of the isolated sand blows in the meizoseismal area. The sand blows are larger to the southeast of the epicenter than its northwestern equivalents. This can be stated in other words as well. The area affected by liquefaction to the northwest is less spread out than to the southeast. Implications of these results in future paleoliquefaction investigations for earthquake magnitude and epicentral determinations are extremely important. In cases of highly asymmetrical distribution of liquefaction features such as this one, where rupture propagation tends to be mono-directional, it can be reliably determined an epicentral distance (between earthquake and liquefaction evidence) and an earthquake magnitude only if the largest sand blow is found. Therefore, magnitude estimation using this uneven liquefaction occurrence will surely lead to underrating if only the shortest side of the meizoseismal area is unluckily studied, which can eventually be the only part exhibiting liquefaction evidence, depending on the earthquake location and the distribution of liquefaction-prone environments.     

Phenological Changes of Blooming Diatoms Promoted by Compound Bottom-Up and Top-Down Controls

Understanding phytoplankton species-specific responses to multiple biotic and abiotic stressors is fundamental to assess phenological and structural shifts at the community level. Here, we present the case of Thalassiosira curviseriata, a winter-blooming diatom in the Bahía Blanca Estuary, Argentina, which displayed a noticeable decrease in the past decade along with conspicuous changes in phenology. We compiled interannual field data to assess compound effects of environmental variations and grazing by the invasive copepod Eurytemora americana. The two species displayed opposite trends over the period examined. The diatom decreased toward the last years, mainly during the winters, and remained relatively constant over the other seasons, while the copepod increased toward the last years, with an occurrence restricted to winter and early spring. A quantitative assessment by structural equation modeling unveiled that the observed long-term trend of T. curviseriata resulted from the synergistic effects of environmental changes driven by water temperature, salinity, and grazing. These results suggest that the shift in the abundance distribution of T. curviseriata toward higher annual ranges of temperature and salinity—as displayed by habitat association curves—constitutes a functional response to avoid seasonal overlapping with its predator in late winters. The observed changes in the timing and abundance of the blooming species resulted in conspicuous shifts in primary production pulses. Our results provide insights on mechanistic processes shaping the phenology and structure of phytoplankton blooms.     

The new national climate change documents of Mexico: what do the regional climate change scenarios represent?

This paper presents a review of the methodology applied for generating the regional climate change scenarios utilized in important National Documents of Mexico, such as the Fourth National Communication to the United Nations Framework Convention on Climate Change, the Fourth National Report to the Convention on Biological Diversity and The Economics of Climate Change in Mexico. It is shown that these regional climate change scenarios, which are one of the main inputs to support the assessments presented in these documents, are an example of the erroneous use of statistical downscaling techniques. The arguments presented here imply that the work based on such scenarios should be revised and therefore, these documents are inadequate for supporting national decision- making.     

Low calcification rates and calcium carbonate production in Porites panamensis at its northernmost geographic distribution

Porites panamensis is a hermatypic coral present in the eastern Pacific Ocean. Skeletal growth parameters have been reported, but studies of the relationship between annual calcification rates and environmental controls are scarce. In this study, we investigated three aspects of the annual calcification rates of P. panamensis: growth parameters among three P. panamensis populations; the sea surface temperature as a calcification rate control spanning a latitudinal gradient; and calcium carbonate production among three sites. Growth parameters varied among the sites due to the colony growth form. Massive colonies in the north showed a higher calcification rate than encrusting colonies in the south (mean: 1.22–0.49 g CaCO₃ · cm^?2 · yr^?1), where variations in calcification rates were related to growth rate (0.91–0.38 cm · yr^?1) rather than to skeletal density differences (overall mean ± SD, 1.31 ± 0.04 g CaCO₃ · cm^?3). Our results showed a positive linear relationship between annual calcification rates and sea surface temperatures within these P. panamensis populations. Differences were related to distinct oceanographic environments (within and at the entrance of the Gulf of California) with different sea surface temperature regimes and other chemical properties. Different populations calcified under different environmental conditions. Calcium carbonate production was dependent upon the calcification rate and coral cover and so carbonate production was higher in the north (coral cover 12%) than in the south (coral cover 3.5). Thus, the studied sites showed low calcium carbonate production (0.25–0.43 kg CaCO₃ · m^?2 · yr^?1). Our results showed reduced calcification rates, regional temperature regime control over calcification rates, different growth forms, low coral cover and low calcium carbonate production rates in P. panamensis.