The fit of the continents around the North Atlantic Ocean

A new reassembly of the continents around the North Atlantic Ocean is presented. The first criterion used for this reassembly is the identification of the structural framework related to the opening which consists of marginal fracture zones generated by offsets of the Rift. The Africa—North America, Eurasia—Greenland, Greenland—North America and Eurasia—North America adjustments are successively discussed. It is argued that the adjustments are best made at the 3000-meters isobath between Africa and North America and at the 2000-meters isobath for the younger rifts. The difference is attributed to subsidence and modification of continental margins with time. The importance of the Late Paleozoic tectonic phase in determining the subsequent pattern of Mesozoic rifting is emphasized.     

Regression‐based approach to low flow prediction in the Maryland Piedmont region under joint climate and land use change

We examine the low flow records for six urbanized watersheds in the Maryland Piedmont region and develop regression equations to predict annual minimum low flow events. The effects of both future climate (based on precipitation and temperature projections from two climate models: Hadley and the Canadian Climate Centre (CCC)) and land use change are incorporated to illustrate possible future trends in low flows. A regression modelling approach is pursued to predict the minimum annual 7‐day low flow estimates for the proposed future scenarios. A regional regression model was calibrated with between 10 and 50 years of daily precipitation, daily average temperature, annual imperviousness, and the daily observed flow time‐series across six watersheds. Future simulations based on a 55 km² urbanizing watershed just north of Washington, DC, were performed. When land use and climate change were employed singly, the former predicted no trends in low flows and the latter predicted significant increasing trends under Hadley and no trends under CCC. When employed jointly, however, low flows were predicted to decrease significantly under CCC, whereas Hadley predicted no significant trends in low flows. Antecedent precipitation was the most influential predictor on low flows, followed by urbanization. Copyright © 2006 John Wiley & Sons, Ltd.     

INTEGRAL/XMM views on the MeV source GRO J1411-64

The COMPTEL unidentified source GRO J 1411-64 was observed by INTEGRAL and XMM-Newton in 2005. The Circinus Galaxy is the only source detected within the 4σ location error of GRO J1411-64, but in here excluded as the possible counterpart. At soft X-rays, 22 reliable and statistically significant sources (likelihood >10) were extracted and analyzed from XMM-Newton data. Only one of these sources, XMMU J141255.6-635932, is spectrally compatible with GRO J1411-64 although the fact the soft X-ray observations do not cover the full extent of the COMPTEL source position uncertainty make an association hard to quantify and thus risky. At the best location of the source, detections at hard X-rays show only upper limits, which, together with MeV results obtained by COMPTEL suggest the existence of a peak in power output located somewhere between 300–700 keV for the so-called low state. Such a spectrum resembles those in blazars or microquasars, and might suggest at work by the models accordingly. However, an analysis using a microquasar model consisting on a magnetized conical jet filled with relativistic electrons, shows that it is hard to comply with all observational constrains. This fact and the non-detection at hard X-rays introduce an a-posteriori question mark upon the physical reality of this source, what is discussed here.     

Calibrating a glaciological model of the Greenland ice sheet from the Last Glacial Maximum to present-day using field observations of relative sea level and ice extent

We constrain a three-dimensional thermomechanical model of Greenland ice sheet (GrIS) evolution from the Last Glacial Maximum (LGM, 21 ka BP) to the present-day using, primarily, observations of relative sea level (RSL) as well as field data on past ice extent. Our new model (Huy2) fits a majority of the observations and is characterised by a number of key features: (i) the ice sheet had an excess volume (relative to present) of 4.1 m ice-equivalent sea level at the LGM, which increased to reach a maximum value of 4.6 m at 16.5 ka BP; (ii) retreat from the continental shelf was not continuous around the entire margin, as there was a Younger Dryas readvance in some areas. The final episode of marine retreat was rapid and relatively late (c. 12 ka BP), leaving the ice sheet land based by 10 ka BP; (iii) in response to the Holocene Thermal Maximum (HTM) the ice margin retreated behind its present-day position by up to 80 km in the southwest, 20 km in the south and 80 km in a small area of the northeast. As a result of this retreat the modelled ice sheet reaches a minimum extent between 5 and 4 ka BP, which corresponds to a deficit volume (relative to present) of 0.17 m ice-equivalent sea level. Our results suggest that remaining discrepancies between the model and the observations are likely associated with non-Greenland ice load, differences between modelled and observed present-day ice elevation around the margin, lateral variations in Earth structure and/or the pattern of ice margin retreat.     

Arsenic in an alkaline AMD treatment sludge: Characterization and stability under prolonged anoxic conditions

Lime treatment of acid mine drainage (AMD) generates large volumes of neutralization sludge that are often stored under water covers. The sludge consists mainly of calcite, gypsum and a widespread ferrihydrite-like Fe phase with several associated species of metal(loid) contaminants. The long-term stability of metal(loid)s in this chemically ill-defined material remains unknown. In this study, the stability and speciation of As in AMD sludge subjected to prolonged anoxic conditions is determined. The total As concentration in the sludge is 300 mg kg⁻¹. In the laboratory, three distinct water cover treatments were imposed on the sludge to induce different redox conditions (100%N₂, 100%N₂ + glucose, 95%N₂:5%H₂). These treatments were compared against a control of oxidized, water-saturated sludge. Electron micro-probe (EMP) analysis and spatially resolved synchrotron X-ray fluorescence (SXRF) results indicate that As is dominantly associated with Fe in the sludge. In all treatments and throughout the experiment, measured concentrations of dissolved As were less than 5 μg L⁻¹. Dissolved Mn concentration in the N₂ + glucose treatment increased significantly compared to other treatments. Manganese and As K-edge X-ray absorption near edge structure spectroscopy (XANES) analyses showed that Mn was the redox-active element in the solid-phase, while As was stable. Arsenic(V) was still the dominant species in all water-covered sludges after 9 months of anoxic treatments. In contrast, Mn(IV) in the original sludge was partially reduced into Mn(II) in all water-covered sludges. The effect was most pronounced in the N₂ + glucose treatment, suggesting microbial reduction. Micro-scale SXRF and XANES analysis of the treated sludge showed that Mn(II) accumulated in areas already enriched in Fe and As. Overall, the study shows that AMD sludges remain stable under prolonged anoxic conditions. External sources of chemical reductants or soluble C were needed to induce lower redox state in the systems, and even under these imposed treatments, only weakly reducing conditions (Mn threshold) developed. The results suggest that As(V) in AMD sludge will remain stable under prolonged anoxic conditions as long as Mn(IV) is present and organic matter accumulation is negligible.     

Tagging-the-triangle algorithm for partitioning features with inconsistent boundaries

In spatial data sets, gaps or overlaps among features are frequently found in spatial tessellations due to the non-abutting edges with adjacent features. These non-abutting edges in loose tessellations are also called inconsistent boundaries or slivers; polygons containing at least one inconsistent boundary are called inconsistent polygons or sliver polygons. The existing algorithms to solve topological inconsistencies in sliver polygons suffer from one or more of three major issues, namely determination of tolerances, excessive CPU processing time for large data sets and loss of vertex history. In this article, we introduce a new algorithm that mitigates these three issues. Our algorithm efficiently searches the features with inconsistent polygons in a given spatial data set and logically partitions them among adjacent features. The proposed algorithm employs the constrained Delaunay triangulation technique to generate labelled triangles from which inconsistent polygons with gaps and overlaps are identified using label counts. These inconsistent polygons are then partitioned using the straight skeleton method. Moreover, each of these partitioned gaps or overlaps is distributed among the adjacent features to improve the topological consistency of the spatial data sets. We experimentally verified our algorithm using the real land cadastre data set. The comparison results show that the proposed algorithm is four times faster than the existing algorithm for data sets with 200,000 edges.     

Geophysical exploration on the subsurface geology of La Garrotxa monogenetic volcanic field (NE Iberian Peninsula)

We applied self-potential (SP) and electrical resistivity tomography (ERT) to the exploration of the uppermost part of the substrate geology and shallow structure of La Garrotxa monogenetic volcanic field, part of the European Neogene–Quaternary volcanic province. The aim of the study was to improve knowledge of the shallowest part of the feeding system of these monogenetic volcanoes and of its relationship with the subsurface geology. This study complements previous geophysical studies carried out at a less detailed scale and aimed at identifying deeper structures, and together will constitute the basis to establish volcanic susceptibility in La Garrotxa. SP study complemented previous smaller-scale studies and targeted key areas where ERT could be conducted. The main new results include the generation of resistivity models identifying dykes and faults associated with several monogenetic cones. The combined results confirm that shallow tectonics controlling the distribution of the foci of eruptive activity in this volcanic zone mainly correspond to NNW–SSE and accessorily by NNE–SSW Neogene extensional fissures and faults and concretely show the associated magmatic intrusions. These structures coincide with the deeper ones identified in previous studies, and show that previous Alpine tectonic structures played no apparent role in controlling the loci of this volcanism. Moreover, the results obtained show that the changes in eruption dynamics occurring at different vents located at relatively short distances in this volcanic area are controlled by shallow stratigraphical, structural and hydrogeological differences underneath these monogenetic volcanoes.     

Can a Training Image Be a Substitute for a Random Field Model?

In most multiple-point simulation algorithms, all statistical features are provided by one or several training images (TI) that serve as a substitute for a random field model. However, because in practice the TI is always of finite size, the stochastic nature of multiple-point simulation is questionable. This issue is addressed by considering the case of a sequential simulation algorithm applied to a binary TI that is a genuine realization of an underlying random field. At each step, the algorithm uses templates containing the current target point as well as all previously simulated points. The simulation is validated by checking that all statistical features of the random field (supported by the simulation domain) are retrieved as an average over a large number of outcomes. The results are as follows. It is demonstrated that multiple-point simulation performs well whenever the TI is a complete (infinitely large) realization of a stationary, ergodic random field. As soon as the TI is restricted to a limited domain, the statistical features cannot be obtained exactly, but integral range techniques make it possible to predict how much the TI should be extended to approximate them up to a prespecified precision. Moreover, one can take advantage of extending the TI to reduce the number of disruptions in the execution of the algorithm, which arise when no conditioning template can be found in the TI.     

High-frequency radar observations of the June 2013 US East Coast meteotsunami

We report here on the observation and offline detection of the meteotsunami off the New Jersey coast on June 13, 2013, using coastal radar systems and tide gauges. This work extends the previous observations of tsunamis originating in Japan and Indonesia. The radars observed the meteotsunami 23 km offshore, 47 min before it arrived at the coast. Subsequent observations showed it moving onshore. The neighboring tide gauge height reading provides confirmation of the radar observations near the shore.