Constraining groundwater modeling with magnetic resonance soundings

Magnetic resonance sounding (MRS) is a noninvasive geophysical method that allows estimating the free water content and transmissivity of aquifers. In this article, the ability of MRS to improve the reliability of a numerical groundwater model is assessed. Thirty-five sites were investigated by MRS over a ～5000 km(2) domain of the sedimentary Continental Terminal aquifer in SW Niger. Time domain electromagnetic soundings were jointly carried out to estimate the aquifer thickness. A groundwater model was previously built for this section of the aquifer and forced by the outputs from a distributed surface hydrology model, to simulate the observed long-term (1992 to 2003) rise in the water table. Uncertainty analysis had shown that independent estimates of the free water content and transmissivity values of the aquifer would facilitate cross-evaluation of the surface-water and groundwater models. MRS results indicate ranges for permeability (K = 1 × 10(-5) to 3 × 10(-4) m/s) and for free water content (w = 5% to 23% m(3) /m(3) ) narrowed by two orders of magnitude (K) and by ～50% (w), respectively, compared to the ranges of permeability and specific yield values previously considered. These shorter parameter ranges result in a reduction in the model's equifinality (whereby multiple combinations of model's parameters are able to represent the same observed piezometric levels), allowing a better constrained estimate to be derived for net aquifer recharge (～22 mm/year).     

The mesoscale variability in the Caribbean Sea. Part II: Energy sources

The processes which drive the production and the growth of the strong mesoscale eddy field in the Caribbean Sea are examined using a general circulation model. Diagnostics of the simulations suggest that:(1) The mean currents in the Caribbean Sea are intrinsically unstable. The nature of the instability and its strength vary spatially due to strong differences of current structure among basins.(2) The greatest and most energetic eddies of the Caribbean Sea originate in the Venezuela Basin by mixed barotropic-baroclinic instability of an intense jet, formed with waters mostly from the surface return flow of the Meridional Overturning Circulation and the North Equatorial Current which converge and accelerate through the Grenada Passage. The vertical shear of this inflow is enhanced by an eastward undercurrent, which flows along the south American Coast between 100 and 250 m depth. The shallow eddies (less than 200 m depth) formed in the vicinity of the Grenada Passage get rapidly deeper (down to 1000 m depth) and stronger by their interaction with the deep interior flow of the Subtropical Gyre, which enters through passages north of St. Lucia. These main eastern Caribbean inflows merge and form the southern Caribbean Current, whose baroclinic instability is responsible for the westward growth and strengthening of these eddies from the Venezuela to the Colombia Basin.(3) Eddies of lesser strength are produced in other regions of the Caribbean Sea. Their generation and growth is also linked with instability of the local currents. First, cyclones are formed in the cyclonic shear of the northern Caribbean Current, but appear to be rapidly dissipated or absorbed by the large anticyclones coming from the southern Caribbean. Second, eddies in the Cayman Sea, which impact the Yucatan region, are locally produced and enhanced by barotropic instability of the deep Cayman Current.(4) The role of the North Brazil Current (NBC) rings is mostly to act as a finite perturbation for the instability of the mean flow. Their presence near the Lesser Antilles is ubiquitous and they appear to be linked with most of the Caribbean eddies. There are some evidences that the frequency at which they form near the Grenada Passage is influenced by the frequency at which the NBC rings impinge the Lesser Antilles. But large Caribbean eddies also form without a close influence of any ring, and comparison between simulations shows that mean eddy kinetic energy and eddy population in the Caribbean Sea are not substantially different in absence or presence of NBC rings: their presence is not a necessary condition for the generation and growth of the Caribbean eddies.     

Changes in water quality following tidal inundation of coastal lowland acid sulfate soil landscapes

This study examines the remediation of surface water quality in a severely degraded coastal acid sulfate soil landscape. The remediation strategy consisted of partial restoration of marine tidal exchange within estuarine creeks and incremental tidal inundation of acidified soils, plus strategic liming of drainage waters. Time-series water quality and climatic data collected over 5 years were analysed to assess changes in water quality due to this remediation strategy. A time-weighted rainfall function (TWR) was generated from daily rainfall data to integrate the effects of antecedent rainfall on shallow groundwater levels in a way that was relevant to acid export dynamics. Significant increases in mean pH were evident over time at multiple monitoring sites. Regression analysis at multiple sites revealed a temporal progression of change in significant relationships between mean daily electrical conductivity (EC) vs. mean daily pH, and TWR vs. mean daily pH. These data demonstrate a substantial decrease over time in the magnitude of creek acidification per given quantity of antecedent rainfall. Data also show considerable increase in soil pH (2–3 units) in formerly acidified areas subject to tidal inundation. This coincides with a decrease in soil pe, indicating stronger reducing conditions. These observations suggest a fundamental shift has occurred in sediment geochemistry in favour of proton-consuming reductive processes. Combined, these data highlight the potential effectiveness of marine tidal inundation as a landscape-scale acid sulfate soil remediation strategy.     

Post-glacial persistence of turbiditic activity within the Rhône deep-sea turbidite system (Gulf of Lions,Western Mediterranean): Linking the outer shelf and the basin sedimentary records

Emplacement of post-glacial turbidites is commonly controlled by rapid changes in sea level or by seismicity. On the continental rise of the Gulf of Lions (Western Mediterranean), an aseismic area, we identified turbiditic beds deposited during the rising stage and highstand of sea level. Swath bathymetry, sediment cores, in situ Cone Penetrating Tests (CPTU), heavy mineral associations and radiocarbon dating determined the source, composition, distribution and age of the turbiditic beds. Turbidites are composed of homogeneous to positively graded silts to medium sand with quartz (up to 90%), shell debris and shelfal benthic faunas. Their distribution on the sea floor is very patchy and controlled by abundant inherited erosional bedforms. Their source is found in relict regressive sands at the outershelf. Their deposition occurred just after the onset of the post-glacial sea level rise and the concomitant sediment starvation of the Rhône deep sea turbiditic system until recently. Whilst canyons are fed with sand by strong seasonal hydro-sedimentary dynamics on the outershelf, the emplacement of post-glacial turbidites is not controlled by sea level changes but probably by the periodic flushing of the canyons. Our study revealed that this low energy aseismic margin undergoes significant transport of sand, down to the base of slope, during the sea-level rise and the Holocene highstand.     

Scapolite stoichiometry

Over the range Ca/(Ca + Na) from zero to 0.75, scapolites vary linearly between Na₄Al₃Si₉O₂₄Cl and NaCa₃Al₅Si₇O₂₄CO₃, by means of a coupled replacement of Na₃Si₂Cl by Ca₃Al₂CO₃. Ca includes minor Sr, Na includes K, and CO₃ may include substantial SO₄. From Ca/(Ca + Na) = 0.75 to 1.00, variation is between NaCa₃Al₅Si₇O₂₄CO₃ and Ca₄Al₆Si₆O₂₄CO₃, by means of a coupled replacement of NaSi by CaAl, as in plagioclase. Anion substitution is complete at the 0.75 point. Recalculation of scapolite analyses on the basis of Al + Si atoms = 12.0, as suggested by computed unit cell contents, offers a reasonable basis for normalization of atomic proportions. On this basis, Ca + Sr + Na + K + Fe atoms total very close to 4.0. New limits on the minor element content of scapolite are proposed. In view of the inflection at 0.75, subdivision of the scapolites at Ca/(Ca + Na) ratios of 0.25, 0.5, and 0.75 seems more logical than the 0.2, 0.5, and 0.8 limits used hitherto for marialite, dipyre, mizzonite, and meionite.     

Predicting road system speeds using spatial structure variables and network characteristics

Spatial regression is applied to GPS floating car measurements to build a predictive model of road system speed as a function of link type, time period, and spatial structure. The models correct for correlated spatial errors and autocorrelation of speeds. Correlation neighborhoods are based on either Euclidean or network distance. Econometric and statistical methods are used to choose the best model form and statistical neighborhood. Models of different types have different coefficient estimates and fit quality, which might affect inferences. Speed predictions are validated against a holdout sample to illustrate the usefulness of spatial regression in road system speed monitoring.      

Alkali-calcic and alkaline post-orogenic (PO) granite magmatism: petrologic constraints and geodynamic settings

The end of an orogenic Wilson cycle corresponds to amalgamation of terranes into a Pangaea and is marked by widespread magmatism dominated by granitoids. The post-collision event starts with magmatic processes still influenced by subducted crustal materials. The dominantly calc-alkaline suites show a shift from normal to high-K to very high-K associations. Source regions are composed of depleted and later enriched orogenic subcontinental lithospheric mantle, affected by dehydration melting and generating more and more K- and LILE-rich magmas. In the vicinity of intra-crustal magma chambers, anatexis by incongruent melting of hydrous minerals may generate peraluminous granitoids bearing mafic enclaves. The post-collision event ends with emplacement of bimodal post-orogenic (PO) suites along transcurrent fault zones. Two suites are defined, (i) the alkali-calcic monzonite–monzogranite–syenogranite–alkali feldspar granite association characterised by [biotite+plagioclase] fractionation and moderate [LILE+HFSE] enrichments and (ii) the alkaline monzonite–syenite–alkali feldspar granite association characterised by [amphibole+alkali feldspar] fractionation and displaying two evolutionary trends, one peralkaline with sodic mafic mineralogy and higher enrichments in HFSE than in LILE, and the other aluminous biotite-bearing marked by HFSE depletion relative to LILE due to accessory mineral precipitation. Alkali-calcic and alkaline suites differ essentially in the amounts of water present within intra-crustal magma chambers, promoting crystallisation of various mineral assemblages. The ultimate enriched and not depleted mantle source is identical for the two PO suites. The more primitive LILE and HFSE-rich source rapidly replaces the older orogenic mantle source during lithosphere delamination and becomes progressively the thermal boundary layer of the new lithosphere. Present rock compositions are a mixture of major mantle contribution and various crustal components carried by F-rich aqueous fluids circulating within convective cells created around magma chambers. In favourable areas, PO suites pre-date a new orogenic Wilson cycle.     

Horizontal and Vertical Co2 Advection In A Sloping Forest

A system measuring the horizontal and vertical advection was devised and installed in a sloping forest at the Vielsalm site, Belgium. The measurements showed that under stable conditions a flow regime established below the canopy: air flowed horizontally along the slope and entrained the air above the canopy vertically. This movement occurs during stable nights characterised by strongly negative net radiation. It creates negative air concentration gradients in both the vertical and horizontal directions. The advection fluxes associated with these movements are opposite and of a similar order of magnitude. This implies that the horizontal advection cannot be ignored in the carbon budget equation at night. Unfortunately, the large variability of, and considerable uncertainty about, advection fluxes does not enable one to produce estimates of the source term from these equations. Advection measurement systems should be improved in order to enable such estimates to be made. Particular attention should be paid to the estimation of the vertical velocity above the canopy and to the vertical profiles of the horizontal velocity and horizontal CO₂ gradient below the canopy.