Using snakes for the registration of topographic road database objects to ALS features

For historical reasons many national mapping agencies store their topographic data in a dual system consisting of a Digital Landscape Model (DLM) and a Digital Terrain Model (DTM). The DLM contains 2D vector data representing objects on the Earth’s surface, such as roads and rivers, whereas the DTM is a 2.5D representation of the related height information, often acquired by Airborne Laser Scanning (ALS). Today, many applications require reliable 3D topographic data. Therefore, it is advantageous to convert the dual system into a 3D DLM. However, as a result of different methods of acquisition, processing, and modelling, the registration of the two data sets often presents difficulties. Thus, a straightforward integration of the DTM and DLM might lead to inaccurate and semantically incorrect 3D objects.In this paper we propose a new method for the fusion of the two data sets that exploits parametric active contours (also called snakes), focusing on road networks. For that purpose, the roads from a DLM initialise the snakes, defining their topology and their internal energy, whereas ALS features exert external forces to the snake via the image energy. After the optimisation process the shape and position of the snakes should coincide with the ALS features. With respect to the robustness of the method several known modifications of snakes are combined in a consistent framework for DLM road network adaptation. One important modification redefines the standard internal energy and thus the geometrical model of the snake in order to prevent changes in shape or position not caused by significant features in the image energy. For this purpose, the initial shape is utilized creating template-like snakes with the ability of local adaptation. This is one crucial point towards the applicability of the entire method considering the strongly varying significance of the ALS features. Other concepts related to snakes are integrated which enable our method to model network and ribbon-like characteristics simultaneously. Additionally, besides ALS road features information about context objects, such as bridges and buildings, is introduced as part of the image energy to support the optimisation process. Meaningful examples are presented that emphasize and evaluate the applicability of the proposed method.     

Palaeovariations in the East-Asian Monsoon Regime Geochemically Recorded in Varved Sediments of Lake Sihailongwan (Northeast China, Jilin Province). Part 2: a 200-Year Record of Atmospheric Lead-210 Flux Variations and its Palaeoclimatic Implications

Palaeovariations of the atmospheric ²¹⁰Pb flux in Northeast China (Long Gang area, Jilin province), quantified by high resolution ²¹⁰Pb measurements on seasonally laminated sediments of Lake Sihailongwan are presented on decadal scale. The mean flux of unsupported ²¹⁰Pb between 1790 and 1970 for the centre of the lake basin is 517 Bq m⁻² yr⁻¹ with maximum deviations between −23% and +27% of this value. Flux rates above this average were found between 1783 and 1813, around 1836, and between 1860 and 1901. The mean ²¹⁰Pb^exc flux rates derived were clearly lower around 1821, 1908, 1930, and 1953. The ²¹⁰Pb^exc flux reached its minimum of 344 Bq m⁻² yr⁻¹ during the period 1977–1982. The atmospheric flux of unsupported ²¹⁰Pb is correlated with the precipitation frequency during the summer monsoon and shows coinciding variability with geochemical proxies that document the groundwater inflow into the lake. Al₂O₃-rich dust of remote provenance scavenged by wet-deposition in the rainy season is the major carrier of the atmospheric ²¹⁰Pb^exc flux.     

Turbidity-discharge hysteresis in a meso-scale catchment: The importance of intermediate scale events

In-situ sensors for riverine water quality monitoring are a powerful tool to describe temporal variations when efficient and informative analyses are applied to the large quantities of data collected. Concentration-discharge hysteresis patterns observed during storm events give insights into headwater catchment processes. However, the applicability of this approach to larger catchments is less well known. Here, we evaluate the potential for high-frequency turbidity-discharge (Q) hysteresis patterns to give insights into processes operating in a meso-scale (722 km²) northern mixed land use catchment. As existing event identification methods did not work, we developed a new, objective method based on hydrograph characteristics and identified 76 events for further analysis. Qualitative event analysis identified three recurring patterns. Events with low mean Q (≤ 2 m³/s) often showed short-term, quasi-periodic turbidity variation, to a large extent disconnected from Q variation. High max Q events (≥15 m³/s) were often associated with spring flood or snowmelt, and showed a disconnection between turbidity and Q. Intermediate Q events (mean Q: 2–11 m³/s) were the most informative when applying hysteresis indexes, since changes in turbidity and Q were actually connected. Hysteresis indexes could be calculated on a subset of 60 events, which showed heterogeneous responses: 38% had a clockwise response, 12% anticlockwise, 12% figure eight (clockwise–anticlockwise), 10% reverse figure eight (anticlockwise–clockwise) and 28% showed a complex response. Clockwise hysteresis responses were associated with the wetter winter and spring seasons. Generally, changes in Q and turbidity were small during anticlockwise hysteresis events. Precipitation often influenced figure-eight patterns, while complex patterns often occurred during summer low flows. Analysis of intermediate Q events can improve process understanding of meso-scale catchments and possibly aid in choosing appropriate management actions for targeting a specific observed pattern.     

Methodology for construction, calibration and validation of a national hydrological model for Denmark

An integrated groundwater/surface water hydrological model with a 1 km² grid has been constructed for Denmark covering 43,000 km². The model is composed of a relatively simple root zone component for estimating the net precipitation, a comprehensive three-dimensional groundwater component for estimating recharge to and hydraulic heads in different geological layers, and a river component for streamflow routing and calculating stream–aquifer interaction. The model was constructed on the basis of the MIKE SHE code and by utilising comprehensive national databases on geology, soil, topography, river systems, climate and hydrology. The present paper describes the modelling process for the 7330 km² island of Sjælland with emphasis on the problems experienced in combining the classical paradigms of groundwater modelling, such as inverse modelling of steady-state conditions, and catchment modelling, focussing on dynamic conditions and discharge simulation. Three model versions with different assumptions on input data and parameter values were required until the performance of the final, according to pre-defined accuracy criteria, model was evaluated as being satisfactory. The paper highlights the methodological issues related to establishment of performance criteria, parameterisation and assessment of parameter values from field data, calibration and validation test schemes. Most of the parameter values were assessed directly from field data, while about 10 ‘free’ parameters were subject to calibration using a combination of inverse steady-state groundwater modelling and manual trial-and-error dynamic groundwater/surface water modelling. Emphasising the importance of tests against independent data, the validation schemes included combinations of split-sample tests (another period) and proxy-basin tests (another area).     

Bioassessment Techniques for Monitoring of Eutrophication and Nutrient Limitation in Coastal Ecosystems

Bioassesment by the use of the macroalga, Ulva lactuca L., was carried out in the Limfjord, Denmark, to assess the significance of nitrogen and phosphorus as limiting factors for primary production during 1985, 1993, 1994 and 1995 and for the detection of changes in eutrophication levels.

Minimum and critical tissue concentrations for nitrogen and phosphorus in macroalgae were identified. The concentrations of nitrogen were generally below the critical concentration level in June–October in 1985, 1993, and 1995 but in 1994 nitrogen was only limiting for primary production in short periods. Only in early spring in 1985 and 1993 were the tissue concentrations of phosphorus below the critical concentration level, whereas in 1994 up to 3–4 months showed phosphorus limited growth, indicating that significant changes in limitation patterns can occur between different years.

It was concluded that the use of biomonitoring techniques is well suited as a bioassessment method for direct detection and for providing a time-integrated measure of nutrient availability in coastal waters, and thus for assessing ecosystem health with regard to eutrophication. It is recommended that biomonitors and the concept of critical tissue concentrations should be used in environmental management and incorporated in future monitoring programmes.     

Mechanical compaction and ultrasonic velocity of sands with different texture and mineralogical composition

This study presents the results of experimental compaction while measuring ultrasonic velocities of sands with different grain size, shape, sorting and mineralogy. Uniaxial mechanical compaction tests up to a maximum of 50 MPa effective stress were performed on 29 dry sand aggregates derived from eight different sands to measure the rock properties. A good agreement was found between the Gassmann saturated bulk moduli of dry and brine saturated tests of selected sands. Sand samples with poor sorting showed low initial porosity while sands with high grain angularity had high initial porosity. The sand compaction tests showed that at a given stress well‐sorted, coarse‐grained sands were more compressible and had higher velocities (Vp and Vs) than fine‐grained sands when the mineralogy was similar. This can be attributed to grain crushing, where coarser grains lead to high compressibility and large grain‐to‐grain contact areas result in high velocities. At medium to high stresses the angular coarse to medium grained sands (both sorted sands and un‐sorted whole sands) showed high compaction and velocities (Vp and Vs). The small grain‐to‐grain contact areas promote higher deformation at grain contacts, more crushing and increased porosity loss resulting in high velocities. Compaction and velocities (Vp and Vs) increased with decreasing sorting in sands. However, at the same porosity, the velocities in whole sands were slightly lower than in the well‐sorted sands indicating the presence of loose smaller grains in‐between the framework grains. Quartz‐poor sands (containing less than 55% quartz) showed higher velocities (Vp and Vs) compared to that of quartz‐rich sands. This could be the result of sintering and enlargement of grain contacts of ductile mineral grains in the quartz‐poor sands increasing the effective bulk and shear stiffness. Tests both from wet measurements and Gassmann brine substitution showed a decreasing Vp/Vs ratio with increasing effective stress. The quartz‐rich sands separated out towards the higher side of the Vp/Vs range. The Gassmann brine substituted Vp and Vs plotted against effective stress provide a measure of the expected velocity range to be found in these and similar sands during mechanical compaction. Deviations of actual well log data from experimental data may indicate uplift, the presence of hydrocarbon, overpressure and/or cementation. Data from this study may help to model velocity‐depth trends and to improve the characterization of reservoir sands from well log data in a low temperature (<80–100^o C) zone where compaction of sands is mostly mechanical.     

Predicting Soil-Water and Soil-Air Transport Properties and Their Effects on Soil-Vapor Extraction Efficiency

Accurate prediction of water and air Iran sport parameters in variably saturated soil is necessary for modeling of soil-vapor extraction (SVE) at soil sites contaminated with volatile organic chemicals (VOCs). An expression for predicting saturated water permeability (k_l,s) in undisturbed soils from the soil total porosity and the field capacity soil-water content was developed by fitting a tortuous-tube fluid flow model to measured water permeability and gas diffusivity data. The new k_l,s expression gave accurate predictions when tested against independent k_l,s data. The k_l,s expression was implemented in the Campbell relative water permeability model to yield a predictive model for water permeability in variably saturated, undisturbed soil. The water permeability model, together with recently developed predictive equations for gas permeability and gas diffusivity, was used in a two-dimensional numerical SVE model that also included non-equilibrium mass transfer of VOC from a separate phase (nonaqueous phase liquid [NAPL]) to the air phase. SVE: calculations showed that gas permeability is likely the most important factor controlling VOC migration and vapor extraction efficiency. Water permeability and gas diffusivity effects became significant at water contents near and above field capacity. The NAPL-air mass transfer coefficient also had large impacts on simulated vapor extraction efficiency. The calculations suggest that realistic SVE models need to include predictive expressions for both conveciive, diffusive. and phase-partitioning processes in natural, undisturbed soils.     

Influence of instruments on the H/V spectral ratios of ambient vibrations

For an optimal analysis of the H/V curve, it appears necessary to check the instrument signal to noise ratio in the studied frequency band, to ensure that the signal from the ground noise is well above the internal noise. We assess the reliability and accuracy of various digitizers, sensors and/or digitizer-sensor couples. Although this study is of general interest for any kind of seismological study, we emphasize the influence of equipment on H/V analysis results. To display the impact of the instrumental part on the H/V behavior, some series of tests have been carried out following a step-by-step procedure: first, the digitizers have been tested in the lab (sensitivity, internal noise...), then the three components sensors, still in the lab, and finally the usual user digitizers-sensors couple in lab and outdoors. In general, the digitizer characteristics, verified during this test, correspond well to the manufacturer specifications, however, depending on the digitizer, the quality of the digitized waveform can be very good to very poor, with variation from a channel to another channel (gain, time difference etc.). It appears very clearly that digitizers need a warming up time before the recording to avoid problems in the low-frequency range. Regarding the sensors, we recommend strongly to avoid the use of “classical” accelerometers (i.e., usual force balance technology). The majority of tested seismometers (broadband and short period, even 4.5 Hz) can be used without problems from 0.4 to 25 Hz. In all cases, the instrumentation should be checked first to verify that it works well for the defined study aim, but also to define its limit of use (frequency, sensitivity...).