Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

Recovering the initial conditions of our local Universe from NOG and PSCz catalogues

We apply the ztrace algorithm to the optical NOG and infrared PSC z galaxy catalogues to reconstruct the pattern of primordial fluctuations that have generated our local Universe. We check that the density fields traced by the two catalogues are well correlated, and consistent with a linear relation [either in δ or in  log (1 +δ)  ] with relative bias (of NOG with respect to PSC z )   b _rel= 1.1 ± 0.1  . The relative bias relation is used to fill the optical zone of avoidance at  | b | < 20°  using the PSC z galaxy density field.
We perform extensive testing on simulated galaxy catalogues to optimize the reconstruction. The quality of the reconstruction is predicted to be good at large scales, up to a limiting wavenumber   k _lim≃ 0.4 h Mpc⁻¹  beyond which all information is lost. We find that the improvement arising from the denser sampling of the optical catalogue is compensated by the uncertainties connected to the larger zone of avoidance.
The initial conditions reconstructed from the NOG catalogue are found (analogously to those from the PSC z ) to be consistent with a Gaussian paradigm. We use the reconstructions to produce sets of initial conditions ready to be used for constrained simulations of our local Universe.     

Polycyclic aromatic hydrocarbons in water, seston and copepods in a harbour area in the Western Mediterranean (Ligurian Sea)

Polycyclic aromatic hydrocarbons (PAHs) are widespread carcinogenic, mutagenic and teratogenic pollutants that can significantly affect ecosystems and exert an influence on marine planktonic and benthic communities as well as on higher organisms such as fish and mammals, including humans. In this paper the concentrations of six PAHs in water, seston and copepods were examined in a harbour zone of the Genova Gulf (Ligurian Sea) from November 1997 to June 1998. The importance of freshwater inputs coming from a polluted stream and the influence of seasonal and atmospheric factors were also considered. Highest PAH values were recorded during spring. Fluoranthene, Benzo[b]fluoranthene and Benzo[k]fluoranthene were the most abundant PAHs either in water, seston or copepods. The concentration of total PAHs in seawater was often lower than the detection limit: recorded values ranged between 0.001 and 0.06 ppb. In seston, total PAHs ranged between 0.033 and 6.26 ng·g^?1 and in copepods between 0.01 and 2.34 ng·g^?1. With this study, PAH contamination in the harbour zone of Genova Sestri‐Ponente and the importance of inputs of organic pollutants coming from land and freshwaters have been emphasized.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Evidence of recent climatic anomalies at Trieste (Italy)

A recent noticeable climatic anomaly in some meteomarine quantities is analysed based on observations during the period March 1, 1971-February 28, 1993 taken in Trieste (Italy). Marked deviations from the normal values of the coupled quantities atmospheric pressure-solar irradiance and atmospheric pressure-sea level elevation, during the last 22 years, are analysed. The procedures used bring out a clear and persistent configuration of high pressure, high irradiance and low sea level from 1989 on, previously unmatched.     

Gravity field determination using boundary element methods

The Boundary Element Method (BEM), a numerical technique for solving boundary integral equations, is introduced to determine the earth's gravity field. After a short survey on its main principles, we apply this method to the fixed gravimetric boundary value problem (BVP), i.e. the determination of the earth's gravitational potential from measurements of the intensity of the gravity field in points on the earth's surface. We show how to linearize this nonlinear BVP using an implicit function theorem and how to transform the linearized BVP into a boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation into a linear system of equations. We discuss the major problems of this approach for setting up and solving the linear system. The BVP is numerically solved for a bounded part of the earth's surface using a high resolution reference gravity model, measured gravity values of high density, and a 50 50 m² digital terrain model to describe the earth's surface. We obtain a gravity field resolution of 1 1 km² with an accuracy of the order 10^–3 to 10^–4 in about 1 CPU-hour on a Siemens/Fujitsu SIMD vector pipeline machine using highly sophisticated numerical integration techniques and fast equation solvers. We conclude that BEM is a powerful numerical tool for solving boundary value problems and may be an alternative to classical geodetic techniques.     

Intercomparison of measurements of stratospheric hydrogen fluoride

Observations of the vertical profile of hydrogen fluoride (HF) vapor in the stratosphere and of the vertical column amounts of HF above certain altitudes were made using a variety of spectroscopic instruments in the 1982 and 1983 Balloon Intercomparison Campaigns. Both emission instruments working in the far infrared spectral region and absorption instruments using solar occultation in the 2.5m region were employed. No systematic differences were seen in results from the two spectral regions. A mean profile from 20–45 km is presented, with uncertainties ranging from 20% to 50%. Total columns measured from ground and from 12 km are consistent with the profile if the mixing ratio for HF is small in the tropophere and low stratosphere.     

Late Holocene sea-level changes and isostasy in western Denmark

Cores and exposed cliff sections in salt marshes around Ho Bugt, a tidal embayment in the northernmost part of the Danish Wadden Sea, were subjected to ¹⁴C dating and litho- and biostratigraphical analyses to reconstruct paleoenvironmental changes and to establish a late Holocene relative sea-level history. Four stages in the late Holocene development of Ho Bugt can be identified: (1) groundwater-table rise and growth of basal peat (from at least 2300 BC to AD 0); (2) salt-marsh formation (0 to AD 250); (3) a freshening phase (AD 250 to AD 1600?), culminating in the drying out of the marshes and producing a distinct black horizon followed by an aeolian phase with sand deposition; and (4) renewed salt-marsh deposition (AD 1600? to present). From 16 calibrated AMS radiocarbon ages on fossil plant fragments and 4 calibrated conventional radiocarbon ages on peat, we reconstructed a local relative sea-level history that shows a steady sea-level rise of 4 m since 4000 cal yr BP. Contrary to suggestions made in the literature, the relative sea-level record of Ho Bugt does not contain a late Holocene highstand. Relative sea-level changes at Ho Bugt are controlled by glacio-isostatic subsidence and can be duplicated by a glacial isostatic adjustment model in which no water is added to the world's oceans after ca. 5000 cal yr BP.