Elemental mass-balance of material carried by major world rivers

An estimate of average river particulate matter (RPM) composition was bàsed on analyses of more than 40 elements in the Amazon, Congo, Ganges, Magdalena, Mekong, Parana and Orinoco rivers, to which were added literature data for 13 other major world rivers, covering the whole spectrum of morphoclimatic features. Geographic variations of major elements in RPM are mostly linked to weathering types and to the balance between weathering rate and river transport. As a result of chemical erosion, Al, Fe and Ti are enriched in RPM with respect to the average parent rock, while Na, Ca, Mg and Sr are strongly depleted. These figures are directly related to the relative importance of dissolved and particulate transport in rivers; this has been computed for each of 40 elements. In order to study weathering on a global scale, the total observed elemental fluxes (dissolved + particulate) have been computed and compared to theoretical ones. The latter were derived from the elemental content in the average parent rock and the total quantity of weathered material, computed from the Al ratio in RPM and in parent rock. Observed and theoretical fluxes are balanced for the less mobilized elements (rare earths, Co, Cr, Cs, Fe, Mn, Rb, Si, Th, Ti, U and V) for which no enrichment relative to Al is noted in RPM, and for B, Ba, Ca, K, Mg, Na, Sr which are relatively depleted in RPM due to their high dissolved transport. Additional fluxes have been found for Br, Sb, Pb, Cu, Mo, Zn and are possible also for Ni and P. This is reflected by marked enrichments in RPM relative to Al for the poorly or moderately dissolved transports (Pb, Cu, Zn). Several hypotheses involving either the natural origin (volcanic dust, marine aerosols, geochemical fractionation) or the artificial origin (worldwide pollution) are discussed to explain these discrepancies, assuming river transport and weathering either to be in a steady state on a global scale or not. However, none of them can fully account for these additional fluxes. It is most likely that these excesses have multiple origins, anthropogenic or natural or both. The comparison between RPM and deep-sea clay compositions emphasizes the prime influence of river input on oceanic sedimentation of Si, Al, Fe, Ti, lanthanides, Sc, Rb, V, etc. A few elements such as Zn, Sb, occur in excess in RPM as compared to deep-sea clays; in order to balance this excess, a remobilization of these elements out of the sediment can be considered. Finally, the enrichment of Co, Cu, Mn and Ni in deep-sea clays compared to RPM is discussed and attributed to several sources and processes.     

Variability of the axial morphology and of the gravity structure along the Central Spreading Ridge (North Fiji Basin): evidence for contrasting thermal regimes

The Central Spreading Ridge (CSR) is located in the central part of the North Fiji Basin, a complex back-arc basin created 12 Ma ago between the Pacific and Indo-Australian plates. The 3.5 Ma old CSR is the best developed, for both structure and magmatism, of all the spreading centers identified in the basin, and may be one of the largest spreading systems of the west Pacific back-arc basins. It is more than 800 km long and 50–60 km wide, and has been intensively explored during the French-Japanese STARMER project (1987–1991).The CSR is segmented into three first order segments named, from north to south, N160°, N15° and N-S according to their orientation. This segmentation pattern is similar to that found at mid-ocean ridges. The calculated spreading rate is intermediate and ranges from 83 mm/yr at 20°30 S to 50 mm/yr at 17°S. In addition, there is a change in the axial ridge morphology and gravity structure between the northern and southern sections of the CSR. The axial morphology changes from a deep rift valley (N160° segment), to a dome split by an axial graben (N15° segment) and to a rectangular flat top high (N-S segment). The Mantle Bouguer Anomalies obtained on the northern part of the CSR (N160°/N15° segments) show bull's eye structures associated with mantle upwelling at the 16°50S triple junction and also in the middle of the segments. The Mantle Bouguer Anomalies of the southern part of the ridge (N-S segment) are more homogeneous and consistent with the observed smooth topography associated with axial isostatic compensation.At these intermediate spreading rates the contrast in bathymetry and gravity structure between the segments may reflect differences in heat supply. We suggest that the N160° and N15° segments are cold with respect to the hot N-S segment. We use a non-steady-state thermal model to test this hypothesis. In this model, the accretion is simulated as a nearly steady-state seafloor spreading upon which are superimposed periodic thermal inputs. With the measured spreading rate of 50 mm/yr, a cooling cycle of 200,000 yr develops a thermal state that permits to explain the axial morphology and gravity structure observed on the N160° segment. A spreading rate of 83 mm/yr and a cooling cycle of 120,000 yr would generate the optimal thermal structure to explain the characteristics of the N-S segment. The boundaries between the hot N-S segment and its cold bounding segments are the 18°10 S and 20°30 S propagating rifts. A heat propagation event along the N-S segment at the expense of the adjacent colder failing segments, can explain the sharp changes in the observed morphology and structure between the segments.     

Outflow of trace metals into the Laptev Sea by the Lena River

Freshwater concentrations confirm the pristine character of the Lena River environment as already pointed-out in a previous study with a limited set of data (Martin et al., 1993). Total dissolved concentrations of the freshwater are 13.8 ± 1.6 nM, Cu, 4.4 ± 0.1 nM, Ni, 0.054 ± 0.047 nM, Cd, 642 ± 208 nM, Fe, 0.2–0.3 nM Pb and 1.2 ± 1.0 nM, Zn. For Zn and Pb, a simple mixing of the Lena River waters with the Arctic waters is observed. Relationships with salinity suggest that for Cu, Ni and Cd, there is a mobilization of the dissolved fraction from the suspended matter, with an increase of the dissolved concentration of 1.5, 3 and 6 times, respectively. For Fe, the total dissolved concentrations follow an exponential decrease in the mixing zone and 80% of the total “dissolved” Fe is removed from the solution. For Cu, Ni, Cd and Fe, the riverine end-members are 20 nM, 12 nM, 0.3 nM and 47 nM, respectively. When considering the input of total dissolved metals to the Arctic Ocean, the fraction attributed to the freshwaters from the Arctic rivers appears to be small (4% of the input of dissolved metal to the Arctic Ocean for Cd, 27% for Cu, 11 % for Ni and 2% for Zn). Metal concentrations in the Laptev Sea and Arctic Ocean are very similar, indicating a generally homogeneous distribution in the areas sampled.     

Modelling the effect of temperature on unsaturated soil behaviour

A simple thermohydromechanical (THM) constitutive model for unsaturated soils is described. The effective stress concept is extended to unsaturated soils with the introduction of a capillary stress. This capillary stress is based on a microstructural model and calculated from attraction forces due to water menisci. The effect of desaturation and the thermal softening phenomenon are modelled with a minimal number of material parameters and based on existing models. THM process is qualitatively and quantitatively modelled by using experimental data and previous work to show the application of the model, including a drying path under mechanical stress with transition between saturated and unsaturated states, a heating path under constant suction and a deviatoric path with imposed suction and temperature. The results show that the present model can simulate the THM behaviour in unsaturated soils in a satisfactory way.     

Mid-infrared detection of large longitudinal asymmetries in Io's SO2 atmosphere

We have observed about 16 absorption lines of the ν₂ SO₂ vibrational band on Io, in disk-integrated 19-μm spectra taken with the TEXES high spectral resolution mid-infrared spectrograph at the NASA Infrared Telescope Facility in November 2001, December 2002, and January 2004. These are the first ground-based infrared observations of Io's sunlit atmosphere, and provide a new window on the atmosphere that allows better longitudinal and temporal monitoring than previous techniques. Dramatic variations in band strength with longitude are seen that are stable over at least a 2 year period. The depth of the strongest feature, a blend of lines centered at 530.42 cm⁻¹, varies from about 7% near longitude 180° to about 1% near longitude 315° W, as measured at a spectral resolution of 57,000. Interpretation of the spectra requires modeling of surface temperatures and atmospheric density across Io's disk, and the variation in non-LTE ν₂ vibrational temperature with altitude, and depends on the assumed atmospheric and surface temperature structure. About half of Io's 19-μm radiation comes from the Sun-heated surface, and half from volcanic hot spots with temperatures primarily between 150 and 200 K, which occupy about 8% of the surface. The observations are thus weighted towards the atmosphere over these low-temperature hot spots. If we assume that the atmosphere over the hot spots is representative of the atmosphere elsewhere, and that the atmospheric density is a function of latitude, the most plausible interpretation of the data is that the equatorial atmospheric column density varies from about 1.5×¹⁰¹⁷ cm⁻² near longitude 180° W to about 1.5×¹⁰¹⁶ cm⁻² near longitude 300° W, roughly consistent with HST UV spectroscopy and Lyman-α imaging. The inferred atmospheric kinetic temperature is less than about 150 K, at least on the anti-Jupiter hemisphere where the bands are strongest, somewhat colder than inferred from HST UV spectroscopy and millimeter-wavelength spectroscopy. This longitudinal variability in atmospheric density correlates with the longitudinal variability in the abundance of optically thick, near-UV bright SO₂ frost. However it is not clear whether the correlation results from volcanic control (regions of large frost abundance result from greater condensation of atmospheric gases supported by more vigorous volcanic activity in these regions) or sublimation control (regions of large frost abundance produce a more extensive atmosphere due to more extensive sublimation). Comparison of data taken in 2001, 2002, and 2004 shows that with the possible exception of longitudes near 180° W between 2001 and 2002, Io's atmospheric density does not appear to decrease as Io recedes from the Sun, as would be expected if the atmosphere were supported by the sublimation of surface frost, suggesting that the atmosphere is dominantly supported by direct volcanic supply rather than by frost sublimation. However, other evidence such as the smooth variation in atmospheric abundance with latitude, and atmospheric changes during eclipse, suggest that sublimation support is more important than volcanic support, leaving the question of the dominant atmospheric support mechanism still unresolved.     

Cloud filling of ocean colour and sea surface temperature remote sensing products over the Southern North Sea by the Data Interpolating Empirical Orthogonal Functions methodology

Optical remote sensing data is now being used systematically for marine ecosystem applications, such as the forcing of biological models and the operational detection of harmful algae blooms. However, applications are hampered by the incompleteness of imagery and by some quality problems. The Data Interpolating Empirical Orthogonal Functions methodology (DINEOF) allows calculation of missing data in geophysical datasets without requiring a priori knowledge about statistics of the full dataset and has previously been applied to SST reconstructions. This study demonstrates the reconstruction of complete space–time information for 4 years of surface chlorophyll a (CHL), total suspended matter (TSM) and sea surface temperature (SST) over the Southern North Sea (SNS) and English Channel (EC). Optimal reconstructions were obtained when synthesising the original signal into 8 modes for MERIS CHL and into 18 modes for MERIS TSM. Despite the very high proportion of missing data (70%), the variability of original signals explained by the EOF synthesis reached 93.5% for CHL and 97.2% for TSM. For the MODIS TSM dataset, 97.5% of the original variability of the signal was synthesised into 14 modes. The MODIS SST dataset could be synthesised into 13 modes explaining 98% of the input signal variability. Validation of the method is achieved for 3 dates below 2 artificial clouds, by comparing reconstructed data with excluded input information. Complete weekly and monthly averaged climatologies, suitable for use with ecosystem models, were derived from regular daily reconstructions. Error maps associated with every reconstruction were produced according to Beckers et al. (2006). Embedded in this error calculation scheme, a methodology was implemented to produce maps of outliers, allowing identification of unusual or suspicious data points compared to the global dynamics of the dataset. Various algorithm artefacts were associated with high values in the outlier maps (undetected cloud edges, haze areas, contrails, and cloud shadows). With the production of outlier maps, the data reconstruction technique becomes also a very efficient tool for quality control of optical remote sensing data and for change detection within large databases.     

Magnetic changes observed in the formation of two filaments in a complex active region: TRACE and MSDP observations

This paper is focused on the formation of two filaments in a complex center of decaying active regions (AR 8329 and AR 8326), located in the northern hemisphere. The observations were obtained in Hα by the Multi-channel Subtractive Double Pass spectrograph (MSDP mounted on the German telescope VTT in Tenerife) and EUV lines with TRACE (Transition Region And Corona Explorer). High Doppler shifts are found to be related to the ends of filament segments where canceling magnetic fields are also located (as seen on magnetograms from Big Bear Solar Observatory). At these locations, velocities along the line of sight, derived by using a cloud model method reach −20 km s⁻¹, the segments of filaments merge and frequently a time-related sub-flare is observed by TRACE. The chirality of the filament segments has been determined by different methods: the segments of dextral chirality join together and form a long dextral filament, and a single filament of sinistral chirality forms end to end with the dextral filament but does not merge with it. Assuming a model of twisted flux tube for filament material, we suggest that the dextral filament has negative helicity and a relationship between its formation and the close by sunspot with the same sign of helicity.     

Variational interpolation of high-frequency radar surface currents using DIVAnd

Ocean Dynamics - DIVAnd (Data-Interpolating Variational Analysis, in n-dimensions) is a tool to interpolate observations on a regular grid using the variational inverse method. We have extended...