Retrieval of Large-Scale Hydrological Signals in Africa from GRACE Time-Variable Gravity Fields

Since its launch in April 2002, the Gravity Recovery and Climate Experiment (GRACE) mission is recording the Earth’s time-variable gravity field with temporal and spatial resolutions of typically 7–30?days and a few hundreds of kilometers, allowing the monitoring of continental water storage variations from both continental and river-basin scales. We investigate here large scale hydrological variations in Africa using different GRACE spherical harmonic solutions, using different processing strategies (constrained and unconstrained solutions). We compare our GRACE estimates to different global hydrology models, with different land-surface schemes and also precipitation forcing. We validate GRACE observations through two different techniques: first by studying desert areas, providing an estimate of the precision. Then we compare GRACE recovered mass variations of main lakes to volume changes derived from radar altimetry measurements. We also study the differences between different publicly available precipitation datasets from both space measurements and ground rain gauges, and their impact on soil-moisture estimates.     

Hydrological responses to rainfall variability and dam construction: a case study of the upper Senegal River basin

The understanding of the spatial and temporal dynamic of river systems is essential for developing sustainable water resource management plan. For the Senegal River, this subject is very complex according to the context of (1) transboundary basin, (2) several contrasted climatic zones (Guinea, South Sudanian, North Sudanian and Sahelian) with high rainfall variability and (3) high human pressures (dam construction and water uses). From 1954 to 2000, 80% (mean value) of the Senegal River flows recorded downstream part of the basin are provided by three majors tributaries (Bafing, Bakoye and Faléme) located in the upstream part. Then, in our study, this upper Senegal River basin was chosen in order to investigate the hydrological responses to rainfall variability and dam construction. Two nonparametric statistical methods, Mann–Kendall and Hubert test, were used to detect the long-term changes in the time series of precipitation and water discharge (1954–2000) at the annual and seasonal scales. The continuous wavelet transform (Morlet Wavelet) was employed to characterize the different mode in the water discharge variability. Flow duration curve and cumulative curve methods were used to assess the impact of dams on the hydrological regime of the Senegal River. Results showed that the Senegal River flows have been changing under the influence of both rainfall variation and dam construction. The long-term evolution of water discharge depend on long-term rainfall variability: The wet periods of the 1950s and 1960s correspond to periods of higher river flows, while the droughts of the 1970s and 1980s led to unprecedented river flows deficits. The new period, since 1994, show a high inter-annual variability of rainfall and discharge without clear trend. At seasonal scale, the results showed also a strong relationship between rainfall and runoff (R ² > 0.8) resulting from alternating wet and dry seasons and rapid hydrological responses according to annual rainfall. Nevertheless, the observed flows during dry seasons highlighted the influence of water storage and restitution of infiltrated waters in soils and surficial formations during wet seasons. In the dry seasons, the water budget of the three upstream tributaries showed a water deficit at the downstream gauging station. This deficit was characterized by water loss to underlying aquifers and highlighted the influence of geological setting on water balance. However, in this context, water restitution during the dry season remained dependent on climatic zone and on the total annual rainfall volume during the previous wet season. The results have highlighted an impact of the Manantali dam previously obscured: The dam has no effect on the regulation of high river flows. That is what explains that since its construction in 1988, flooding of coastal cities, like Saint-louis, by seasonal river floods has not ceased. The flooding risk in coastal cities is not avoided, and the dams caused hyper-salinization of the Senegal lower estuary. The breach created in the coastal barrier of the Langue of Barbary in October 2003 promotes direct export of excess floodwater to the sea and reduces this risk of flooding in the delta area. But, this solution led to considerable loss of potential water resources, and the authors recommend a new water management plan with a global focus. However, this study shows the positives impacts of the two dams. They allow the availability of freshwater in order to support agricultural irrigation in the valley and delta zone, in particular during low flows periods.     

Influence of interaction between small asperities on various types of slow earthquakes in a 3-D simulation for a subduction plate boundary

Recently, the occurrence of slow earthquakes such as low-frequency earthquakes and very low-frequency earthquakes have been recognized at depths of about 30 km in southwest Japan and Cascadia. These slow earthquakes occur sometimes in isolation and sometimes break into chain-reaction, producing tremor that migrates at a speed of about 5–15 km/day and suggesting a strong interaction among nearby small asperities. In this study, we formulate a 3-D subduction plate boundary model with two types of small asperities chained along the trench at the depth of 30 km. Our simulation succeeds in representing various types of slow earthquakes including low-frequency earthquakes and rapid slip velocity in the same asperity, and indicates that interaction between asperities may cause the very low-frequency earthquakes. Our simulation also shows chain reaction along trench with propagation speed that can be made consistent with observations by adjusting model parameters, which suggests that the interactions also explain the observed migration of slow earthquakes.     

Hydrodynamics of suspended particulate matter in the tidal freshwater zone of a macrotidal estuary (the Seine Estuary, France)

The effects of fortnightly, semidiurnal, and quaterdiurnal lunar tidal cycles on suspended particle concentrations in the tidal freshwater zone of the Seine macrotidal estuary were studied during periods of medium to low freshwater flow. Long-term records of turbidity show semidiurnal and spring-neap erosion-sedimentation cycles. During spring tide, the rise in low tide levels in the upper estuary leads to storage of water in the upper estuary. This increases residence time of water and suspended particulate matter (SPM). During spring tide periods, significant tidal pumping, measured by flux calculations, prevents SPM transit to the middle estuary which is characterized by the turbidity maximum zone. On a long-term basis, this tidal pumping allows marine particles to move upstream for several tens of kilometers into the upper estuary. At the end of the spring tide period, when the concentrations of suspended particulate matter are at their peak values and the low-tide level drops, the transport of suspended particulate matter to the middle estuary reaches its highest point. This period of maximum turbidity is of short duration because a significant amount of the SPM settles during neap tide. The particles, which settle under these conditions, are trapped in the upper estuary and cannot be moved to the zone of maximum turbidity until the next spring tide. From the upper estuary to the zone of maximum turbidity, particulate transport is generated by pulses at the start of the spring-neap tide transition period.     

Turbulence structure in temperature inversions and in convection fields as observed by doppler sodar

The structure of turbulence in an inversion layer and in an homogeneous convective field of the planetary boundary layer is described. In the first part of the paper, we validate the sodar estimates of turbulent dissipation, by using measurements with an hot-wire anemometric system in situ. Limitations of an ε measurement technique using structure function calculations are given, taking account of atmospheric properties and acoustic Doppler instrumental effects. By comparison between isopleths of backscattering intensity and of turbulent dissipation rates, we observe that in the early morning, turbulence is advected by mechanical turbulence generated by wind shear. The same mechanism seems to be operating in the case of an inversion layer capping thermal instability, when the convective activity is not too greatly developed. A turbulent kinetic energy budget is examined using aircraft, sodar, and tower measurements. This indicates a constant turbulent dissipation profile through a deep convective layer.     

Bidirectional Reflectance of Gaussian Random Surfaces and Its Scaling Properties

Using rigorous probabilistic techniques we define and compute the shadowing factor, the conditional mean slope and the conditional standard deviation of the slope for a random rough surface described as an isotropic Gaussian stochastic process with Gaussian autocorrelation function. We use these quantities in order to obtain the bidirectional reflectance distribution function for incoherent light scattering by rough surfaces. The calculated quantities depend on reduced, affine invariant parameters, such that changing roughness is equivalent to rescaling the slopes of incident and emergent direction. We discuss some possible applications of these scaling properties to remote sensing.     

Weathering, atmospheric deposition and vegetation uptake: role for ecosystem sensitivity to acid deposition and critical load

Critical loads of acidity represent the maximum acceptable atmospheric deposition for an ecosystem type. Two hundred and forty-one ecosystem types have been defined in France using pedologic, geologic and vegetation data. Weathering rate plays the most important part in soil buffering capacity, but for poor weatherable soils, non-marine atmospheric deposition represents up to 80% of base-cation inputs. Base-cation vegetation uptake decreases significantly the buffering capacity in case of high-productivity forests. Ecosystems combining low weathering rate and low non-marine base-cation deposition with high biomass productivity are the most sensitive to acidification. To cite this article: D. Moncoulon et al., C. R. Geoscience 336 (2004).

Résumé

Les charges critiques d'acidité représentent le dépôt atmosphérique maximal admissible pour un écosystème. Deux cent quarante et un types d'écosystèmes ont été définis en France à partir de données pédologiques, géologiques et de végétation. L'altération joue un rôle prépondérant contre l'acidification, mais pour les sols faiblement altérables, les dépôts atmosphériques non marins peuvent représenter jusqu'à 80 % des apports de cations basiques. Le prélèvement de cations par la végétation contribue significativement à diminuer le pouvoir tampon des sols pour les forêts à forte productivité. Les écosystèmes combinant faible altération et faibles dépôts de cations non marins ainsi qu'une forte productivité sont les plus sensibles à l'acidification. Pour citer cet article : D. Moncoulon et al., C. R. Geoscience 336 (2004).     

The distribution of geothermal fields along the East Pacific Rise from 13°10′ N to 8°20′ N: Implications for deep seated origins

In 1983 a combined SeaMARC I, Sea Beam swath mapping expedition traversed the East Pacific Rise from 13°20 N to 9°50 N, including most of the Clipperton Transform Fault at 10°15 N, and a chain of seamounts at 9°50 N which runs obliquely to both the ridge axis and transform fault trends. We collected temperature, salinity and magnetic data along the same track. These data, combined with Deep-Tow data and French hydrocasts, are used to construct a thermal section of the rise axis from 13°10 N to 8°20 N.Thermal data collected out to 25 km from the rise axis and along the Clipperton Transform Fault indicate that temperatures above the rise axis are uniformly warmer by 0.065°C than bottom water temperatures at equal depths off the axis. The rise axis thermal structure is punctuated by four distinct thermal fields with an average spacing of 155 km. All four of these fields are located on morphologic highs. Three fields are characterized by lenses of warmed water 20 km in length and 300 m thick. Additional clues to hydrothermal activity are provided in two cases by high concentrations of CH₄, dissolved Mn and ³He in the water column and in another case by concentrations of benthic animals commonly associated with hydrothermal regions.We use three methods to estimate large-scale heat loss. Heat flow estimates range from 1250 MW to 5600 MW for one thermal field 25 km in length. Total convective heat loss for the four major fields is estimated to lie between 2100 MW and 9450 MW. If we add the amount of heat it takes to warm the rest of the rise axis (489 km in length) by 0.065.°C, then the calculated axial heat loss is from 12,275 to 38,525 MW (19–61% of the total heat theoretically emitted from crust between 0 and 1 m.y. in age).