Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.     

Determining the maximum degree of harmonic coefficients in geopotential models by Monte Carlo methods

Random errors for the harmonic coefficients of a geopotential model are generated from the matrix of normal equations by a parallel computer applying the Gibbs sampler. This leads to random values for the harmonic coefficients. They are transformed by nonlinear, quadratic transformations to random values for the square roots of degree variances, of mean squares of geoid undulations and gravity anomalies. The expected values of these quantities are not equal to the values of these quantities computed by the estimated harmonic coefficients, due to correlations and errors in the estimation. By hypothesis tests estimated harmonic coefficients distorted by correlations and errors are detected. Applying the tests to the geopotential model ITG-CHAMP01 of the Institute of Theoretical Geodesy in Bonn it is concluded that above the degree 62 the harmonic coefficients cannot add any information to the geopotential model.     

Effects of 12 years' operation of a sewage treatment plant on trace metal occurrence within a Mediterranean commercial sponge (Spongia officinalis, Demospongiae)

The present field study uses Spongia officinalis for assessing trace metals occurrence in time and space within Mediterranean rocky communities. Nine sites were selected in the Marseille area for studying spatial trends in 12 metal concentrations. Long term changes in 8 metal concentrations were assessed at sites that had been sampled before and 12 years after the opening of a treatment plant. Spongia officinalis highly concentrated all the trace metal surveyed excepted Hg and Cd. The overall contamination level registered provided a classification of the study sites which is congruent with that given by other studies on pollutant accumulation in neighbouring sandy-bottoms or benthic assemblages. Among the metals studied, Fe, Pb, Cr are those that best highlighted a pollution gradient. In the present study, only Cd concentration did not vary in space. Except for Ni, all pollutant concentrations clearly decreased between 1984 and 1999. This very impressive decrease in heavy metal concentrations within the Marseille area represents an indisputable evidence of the improvement of the seawater quality resulting from 12 years' operation of the Marseille sewage plant. Moreover, the significant decrease also recorded in the reference population at Port-Cros might reflect an overall improvement in the seawater quality of the NW Mediterranean.     

Magnetostratigraphy of the lower Triassic volcanics from deep drill SG6 in western Siberia: evidence for long-lasting Permo–Triassic volcanic activity

The deep drill hole SG6 in western Siberia (66°N, 78.5°E) penetrated 1.1  km of lower Triassic basalts, which are possibly an extension of the central Siberian Permo– Triassic flood basalt province. About 300 samples of these basalts were progressively demagnetized and measured. Principal component analysis often shows multiple magnetizations carried by haematite and magnetite. The corrected mean inclinations are +77° and −77° for the haematite component. A magnetostratigraphic scale was derived and showed a N–R–N–R–N succession. This is quite different from the Noril'sk and Taimyr typical polarity scale, R–N.
  The basalts found in the SG6 deep drill hole are slightly younger than those of central Siberia and Taimyr. They correspond to middle–upper Induan age, whereas the Noril'sk and Taimyr sections correspond to an uppermost Permian and lower Induan age. Altogether they indicate that, after a high output rate of volcanic material near the Permo–Triassic boundary, this activity slowed down drastically on the Siberian platform and Taimyr, but persisted for several million years in western Siberia.     

Wave train model for knickpoint migration

Rivers respond to a drop in their base level by incising the topography. The upstream propagation of an incision, as usually depicted by a knickpoint migration, is thought to depend on several parameters such as the drainage area, lithology, and the amplitude of the base level drop. We first investigate the case of the Messinian Salinity Crisis that was characterized by the extreme base level fall (1500 m) of the Mediterranean Sea at the end of the Miocene. The response of drainage areas of three orders of magnitude (10³ to 10⁶ km²) highlights the dominant role of the drainage area (with a square root relationship) in controlling the knickpoint migration after a base level fall. A compilation of mean rates of knickpoint propagation for time durations ranging from 10² to 10⁷ years displays a similar relationship indicating that successive wave trains of knickpoint can migrate in a river: first, wave trains linked to the release of the alluvial cover and then, wave trains related to the bedrock incision, which correspond to the real time response of rivers. Wave trains with very low retreat rates (long lived knickpoints > 1 My) rather correspond to the response time of regional landscape.     

Estimating the error in vertically interpolated forecast‐observation differences: U‐ and v‐wind components

Abstract

Data assimilation in numerical weather forecasting corrects current forecast values by subtracting a portion of interpolated forecast‐minus‐observation differences at the points of a three‐dimensional grid. Deviations used in updating a forecast data field are forecast errors obtained or derived from observations available at update time. When observations are missing at mandatory levels, construction of full vertical soundings by interpolation introduces extraneous errors. The present paper is concerned with determination of the error in vertical extrapolations of surface winds, and of aircraft and satellite cloud‐tracked winds. In addition it examines the effect on accuracy of using location‐specific statistics compared to averaged statistics as the basis for the interpolation weighting scheme and compares errors of one‐ and two‐variable interpolations.

Interpolation accuracy tests demonstrate the influence of the interpolation scheme on the quality of interpolated information used in forecast updating. The results show that the level of accuracy exceeds the benchmark provided by monthly mean forecast error values only with bivariate interpolation of wind components from off‐level data sources.     

Climate and climate change in western canada as simulated by the Canadian regional climate model

Abstract

A þrst climate simulation performed with the novel Canadian Regional Climate Model (CRCM) is presented. The CRCM is based on fully elastic non‐hydrostatic þeld equations, which are solved with an efþcient semi‐implicit semi‐Lagrangian (SISL) marching algorithm, and on the parametrization package of subgrid‐scale physical effects of the second‐generation Canadian Global Climate Model (GCMII). Two 5‐year integrations of the CRCM nested with GCMII simulated data as lateral boundary conditions are made for conditions corresponding to current and doubled CO₂ scenarios. For these simulations the CRCM used a grid size of 45 km on a polar‐stereographic projection, 20 scaled‐height levels and a time step of 15 min; the nesting GCMII has a spectral truncation of T32, 10 hybrid‐pressure levels and a time step of 20 min. These simulations serve to document: (1) the suitability of the SISL numerical scheme for regional climate modelling, (2) the use of GCMII physics at much higher resolution than in the nesting model, (3) the ability of the CRCM to add realistic regional‐scale climate information to global model simulations, and (4) the climate of the CRCM compared to that of GCMII under two greenhouse gases (GHG) scenarios.     

Ozone production due to emissions from vegetation burning

Ozone has been observed in elevated concentrations by satellites over areas previously believed to be background. There is meteorological evidence, that these ozone plumes found over the Atlantic Ocean originate from vegetation fires on the African continent.In a previous study (DECAFE-88), we have investigated ozone and assumed precursor compounds over African tropical forest regions. Our measurements revealed large photosmog layers at altitudes from 1.5 to 4 km. Both chemical and meteorological evidence point to savanna fires up to several thousand km upwind as sources.Here we describe ozone mixing ratios observed over western Africa and compare ozone production ratios from different field measurement campaigns related to vegetation burning. We find that air masses containing photosmog ingredients require several days to develop their oxidation potential, similar to what is known from air polluted by emissions from fossil fuel burning. Finally, we estimate the global ozone production due to vegetation fires and conclude that this source is comparable in strength to the stratospheric input.     

Carboxylic monoacids in the air of mayombe forest (Congo): Role of the forest as a source or sink

In the tropical rain forests of the Congo during the dry season, from June to September 1987, carboxylic acid partial pressures (P _gas) in the air above the canopy, at ground level, and at the boundary layer, were estimated from water samples such as fog and rainwater. The concentrations of these acids were also measured in the sap of tree leaves. Tree leaves act as a sink or as a source if the acid P _gas is greater of lower than the acid concentrations in molecular form in sap. For each of these soluble gases, there is a value of P _gas where the exchange is nul. This is called the compensation point. Values of the compensation point for some tree leaves were evaluated according to Henry's law. Henry's law coefficients at ppm levels were redetermined for formic (HCOOH), acetic (CH₃COOH), propionic (CH₃CH₂COOH), and isobutyric (CH₃CH(CH₃)COOH) acids.By comparison of P _gas and compensation points, it is concluded that the forest was a potential source for these acids. The soil-or the litter-acts as a significant source of a carboxylic acid of C₃ or C₄ atoms in the aliphatic chain. This carboxylic acid, not yet fully characterized, could play an important role in the rain acidity in forested zones of the equatorial areas.The direct emission of these carboxylic acids by vegetation was the main source in the boundary layer above the forest. The average emissions were 3.1×10⁹, 7.8×10⁹, and 8.4×10⁹ molecules cm^-2 s^-1 for HCOOH, CH₃COOH, and CH₃CH₂COOH, respectively. The savanna is an exogenous source of HCOOH and CH₃CH₂COOH during moderately rainy days for 30% of the time. The ozonolysis of isoprene seems to be a small source of HCOOH.