Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Depths of Partial Crystallization of H2O-bearing MORB: Phase Equilibria Simulations of Basalts at the MAR near Ascension Island (7 11{degrees}S)

Phase equilibria simulations were performed on naturally quenchedbasaltic glasses to determine crystallization conditions priorto eruption of magmas at the Mid-Atlantic Ridge (MAR) east ofAscension Island (7–11°S). The results indicate thatmid-ocean ridge basalt (MORB) magmas beneath different segmentsof the MAR have crystallized over a wide range of pressures(100–900 MPa). However, each segment seems to have a specificcrystallization history. Nearly isobaric crystallization conditions(100–300 MPa) were obtained for the geochemically enrichedMORB magmas of the central segments, whereas normal (N)-MORBmagmas of the bounding segments are characterized by polybariccrystallization conditions (200–900 MPa). In addition,our results demonstrate close to anhydrous crystallization conditionsof N-MORBs, whereas geochemically enriched MORBs were successfullymodeled in the presence of 0·4–1 wt% H₂O in theparental melts. These estimates are in agreement with direct(Fourier transform IR) measurements of H₂O abundances in basalticglasses and melt inclusions for selected samples. Water contentsdetermined in the parental melts are in the range 0·04–0·09and 0·30–0·55 wt% H₂O for depleted and enrichedMORBs, respectively. Our results are in general agreement (within±200 MPa) with previous approaches used to evaluate pressureestimates in MORB. However, the determination of pre-eruptiveconditions of MORBs, including temperature and water contentin addition to pressure, requires the improvement of magma crystallizationmodels to simulate liquid lines of descent in the presence ofsmall amounts of water. KEY WORDS: MORB; Mid-Atlantic Ridge; depth of crystallization; water abundances; phase equilibria calculations; cotectic crystallization; pressure estimates; polybaric fractionation     

Pole-to-pole moisture conditions for the IGY

Summary A study of the mean atmospheric humidity conditions on a planetary scale during the IGY covering the calendar year 1958 is presented. The fields of mean precipitable water content and of the zonal and meridional transports of water vapor are analyzed for the entire globe. Zonally averaged values of the various quantities at several levels are presented in tabular form and compared whenever possible with previous results or indirectly are analyzed on the basis of information obtained from different sources. The structure of these fields is studied and the corresponding implications for the general circulations of the atmosphere are discussed. Finally the zonal water balance for all the globe is discussed and its implications analyzed.     

Particulate barium fluxes and export production in the northwestern Mediterranean

Biogenic barium, mostly in the barite (BaSO₄) form, has been proposed as a tracer for export production in the ocean. Here we report on biogenic barium (Ba_xs) and particulate organic carbon (POC) fluxes from sediment traps deployed at the DYFAMED site in the Northwestern Mediterranean Sea. Ba_xs fluxes display average values of 37 ± 45 and 50 ± 58 μg/m²/d at 200 and 1000 m respectively, and are linearly correlated to POC fluxes (mean values of 7.9 ± 9.3 and 6.8 ± 6.8 mg C/m²/d at 200 and 1000 m). Export production estimates, calculated using published Ba_xs- or POC-based algorithms, all fall below or close to the lower limit of potential export values proposed in the literature. This work clearly demonstrates the usefulness of Ba_xs as a tracer of oceanic export production in the Northwestern Mediterranean Sea. However, development of a quantitative export production proxy requires a clear understanding of the underlying cause(s) for the observed spatial variations in the relationship between Ba_xs and POC fluxes. The present study confirms that the processes leading to barite formation differ between margin and open-ocean sites and probably account for much of the regional variability in the POC/Ba_xs ratio.     

Swath mapping of the seafloor and its application to deep-bottom fisheries in New Caledonia

A bottom longline fishery operated in the EEZ of New Caledonia from 1988 to 1991. Fishing focused on five seamounts with summits at depths ranging from 500 to 750 m. The target species was alfonsino, Beryx splendens. As the soundings available from marine charts were not detailed enough, the fishing masters had to make their own charts in order to set the gear in the right location. A series of 11 scientific cruises devoted to a research program on alfonsino was launched in late 1991. During the first scientific cruises, several days were spent mapping the seabed to improve the existing knowledge of the topography of three seamounts by coupling the EDO echosounder depth measurements to the GPS positions. As this procedure is slow, it is applicable to limited areas or, if a wider grid is used, to large zones to locate major structures such as ridges and seamounts. The emergence of the multibeam echosounder has greatly improved seabed mapping performance. This tool seems to fit the requirement of exploratory deep-bottom fishing particularly well, as it covers large areas while providing details of the bottom. The EM 12 was used to obtain bathymetry and imagery of the SE portion of the EEZ of New Caledonia, covering an approximate area of 70,000 km² in two weeks. For depths ranging from 500 to 4,500 m, the results were impressive. They confirmed the known major features and provided greater topographical detail, revealing fine unknown structures. They also provided information on the type of substratum, information which might influence the fishing strategy. Finally, they made it possible to obtain an accurate estimate of the exploitable area which in turn led to new stock estimates.     

The Loyalty—New Hebrides Arc collision: Effects on the Loyalty Ridge and basin system,Southwest Pacific (first results of the ZoNéCo programme)

The ZoNéCo 1 and 2 cruises of Ifremer's Research Vessel L'Atalante, collected new swath bathymetry and geophysical data over the southern and northern segments of the basins and ridges forming the Loyalty system. Between the two surveyed areas, previous studies found evidence for the resistance of the Loyalty Ridge to subduction beneath the New Hebrides trench near 22°S–169°E. On the subducted plate, except for seismicity related to the downbending of the Australian plate, recorded shallow seismicity is sparse within the Loyalty system (Ridge and Basin) where reliable focal mechanism solutions are almost absent.Swath bathymetry, seismic reflection and magnetic data acquired during the ZoNéCo 1 and 2 cruises revealed a transverse asymmetric morphology in the Loyalty system, and an along-strike horst and graben structure on the discontinuous Loyalty Ridge. South of 23°50S and at 20°S, the two WSW-ENE-trending fault systems, respectively, sinistral and dextral, that crosscut the southern and northern segments of the Loyalty system, are interpreted as due to the early effects of collision with the New Hebrides Arc. A NNW-SSE trend, evident along the whole Loyalty system and on the island of New Caledonia, is interpreted as an inherited structural trend that may have been reactivated through flexure of the Australian lithospheric plate at the subduction zone.Overall then, the morphology, structure and evolution of the southern and northern segments of the Loyalty system probably result from the combined effects of the Australian plate lithospheric bulge, the active Loyalty-New Hebrides collision and the overthrust of the New Caledonian ophiolite.     

Oxidative stress in Laeonereis acuta (Polychaeta, Nereididae): environmental and seasonal effects

Laeonereis acuta was seasonally collected in an industrially polluted site (P) and in an unpolluted site (UP) at the Patos Lagoon estuary (southern Brazil). Glutathione S-transferase (GST), superoxide dismutase (SOD) and catalase (CAT) activity (U/mg protein) was determined in five groups of worms from each site. Metallothionein (MT - mol GSH/g ww) and lipid peroxides content (LPO - nmoles of cumene hydroperoxide/g ww) were also measured. Annual mean values for CAT (UP=3.7+/-0.3; P=5.7+/-0.6), GST (UP=0.034+/-0.003; P=0.045+/-0.004) and MT (UP=0.15+/-0.02; P=0.23+/-0.03) were higher (p<0.05) in worms from the P site. In autumn, CAT activity was higher (p<0.05) in worms from the P site (7.6 +/- 1.3) than in those from the UP site (3.6 +/- 0.4). In summer, MT concentration was higher in worms from the P site (0.37 +/- 0.03) than in those from the UP site (0.19 +/- 0.01). No significant difference (p>0.05) in the LPO content was observed in worms from the different sites or collected in different seasons. These results indicate that worms from the polluted site showed higher antioxidant responses than those from the unpolluted site, sufficient to prevent oxidative damage in terms of LPO.     

Use of unit response functions for management of regional multilayered aquifers: application to the North Aquitaine Tertiary system (France)

Although there are many numerical applications used in aquifer management, few methodologies seem to be sufficiently adapted to provide solutions for regional-scale problems. Their applications for multilayered aquifers are also too limited. A new groundwater management tool adapted to such large physical systems has been developed, using the unit response function (URF) approach. This tool is based on the link between a numerical groundwater flow model and a mathematical optimisation tool. Pre- and post processing are based on a geographical information system (GIS). The developed tool has been applied to a real management case—the regional Aquitaine multilayered aquifer (France). A representative hydrodynamic model was built using MODFLOW 2000 code. First tests show that the linearity condition needed for the application of URF methodology is respected despite the high complexity of the system. Two management scenarios were tested. Optimal solutions thus obtained allow for viable exploitation alternatives to be proposed, which integrate complex environmental constraints.

---

Résumé Malgré un nombre important d’applications existant dans le domaine de la gestion des aquifères, peu de méthodologies semblent adaptées pour proposer des solutions à des problèmes réels de dimension régionale. Leurs application pour des systèmes multicouches sont souvent trop limitées. Un outil de gestion adapté à ces systèmes de large taille a été développé, basé sur la théorie des Fonctions de Réponses Unitaires (FRU). Cet outil est basé sur le couplage d’un modèle numérique d’écoulement et d’un code d’optimisation mathématique. Les entrées et sorties du modèle sont réalisées à l’aide d’un Système d’Information Géographique (SIG). L’outil développé est appliqué à un cas de gestion réel—le système aquifère régional multicouche aquitain (France). Un modèle hydrodynamique représentatif a été construit, basé sur le code MODFLOW 2000. Les premiers tests montrent que les conditions de linéarité nécessaires à l’application de la méthode FRU sont respectées, malgré l’importante complexité du système. Deux scénarii de gestion sont testés. Les solutions optimales ainsi obtenues permettent de proposer des solutions d’exploitation alternatives viables, intégrant des contraintes environnementales complexes.

---

Resumen Aunque hay muchas aplicaciones numéricas usadas en la gestión de acuíferos, unas pocas metodologías parecen estar adaptadas suficientemente para dar las soluciones a los problemas de escala regional. También son limitadas sus aplicaciones para los acuíferos multicapa. Se ha desarrollado una nueva herramienta para el manejo del agua subterránea, adaptada a tales sistemas físicos grandes, usando el avance de la Función de Respuesta Unitaria (URF). Esta herramienta se basa en la combinación entre un modelo numérico de flujo de agua subterránea y una herramienta de optimización matemática. El pre-proceso y el post-proceso se basan en un Sistema de Información Geográfico (GIS). La herramienta desarrollada se ha aplicado a un caso real de manejo en el acuífero regional multicapa de Aquitania (Francia). Se construyó un modelo hidrodinámico representativo, que usó el código MODFLOW 2000. Las primeras pruebas muestran que la condición de linearidad, necesaria para la aplicación de la metodología URF, se respeta a pesar de la complejidad alta del sistema. Se probaron dos escenarios de manejo. Las soluciones óptimas así obtenidas permiten la proyección de alternativas viables de explotación, las cuales integran restricciones medioambientales complejas.