Temporal variability in summertime bottom hypoxia in shallow areas of Mobile Bay,Alabama

This paper addresses temporal variability in bottom hypoxia in broad shallow areas of Mobile Bay, Alabama. Time-series data collected in the summer of 2004 from one station (mean depth of 4 m) exhibit bottom dissolved oxygen (DO) variations associated with various time scales of hours to days. Despite a large velocity shear, stratification was strong enough to suppress vertical mixing most of the time. Bottom DO was closely related to the vertical salinity gradient (ΔS). Hypoxia seldom occurred when ΔS (over 2.5 m) was <2 psu and occurred almost all the time when ΔS was >8 psu in the absence of extreme events like hurricanes. Oxygen balance between vertical mixing and total oxygen demand was considered for bottom water from which oxygen demand and diffusive oxygen flux were estimated. The estimated decay rates at 20°C ranging between 0.175–0.322 d⁻¹ and the corresponding oxygen consumption as large as 7.4 g O₂ m⁻² d⁻¹ fall at the upper limit of previously reported ranges. The diffusive oxygen flux and the corresponding vertical diffusivity estimated for well mixed conditions range between 8.6–9.5 g O₂ m⁻² d⁻¹ and 2.6–2.9 m² d⁻¹, respectively. Mobile Bay hypoxia is likely to be associated with a large oxygen demand, supported by both water column and sediment oxygen demands, so that oxygen supply from surface water during destratification events would be quickly exhausted to return to hypoxic conditions within a few hours to days after destratification events are terminated.     

Diffusion in single crystal of melilite: interdiffusion of Al + Al vs. Mg + Si

Interdiffusion coefficients of Al + Al vs. Mg + Si in the gehlenite–åkermanite system of melilite were determined by coupled annealing of synthesized end-member single crystals. The observed diffusion coefficients for a couple-annealed sample vary for about 2 orders of magnitude, showing strong dependence on the gehlenite–åkermanite composition: diffusion coefficient observed at 1350 °C, for example, is 3 × 10^?13 cm² s^?1 at 5 mol% åkermanite composition (Ak₅), increases to 2 × 10^?11 cm² s^?1 at Ak₈₀, and then decreases to 1 × 10^?12 cm² s^?1 at Ak₉₅. The diffusion coefficient–temperature relation indicates high activation energy of diffusion of about 420 kJ mol^?1 for gehlenite-rich melilite. The observed diffusion coefficient–composition relation may be explained by a combination of (1) the diffusion coefficient–melting temperature relation (Flynn's rule) and (2) the feasibility of local charge compensation, which can possibly be maintained more easily in the intermediate chemical composition. The high activation energy value for gehlenitic melilite appears to correspond to the complex diffusion mechanism. The observed highly variable diffusion coefficients suggest that gehlenite–åkermanite zoning in the melilite crystals in Ca, Al-rich inclusions in the carbonaceous meteorites may provide a sensitive indicator for the thermal history of the inclusions.     

Chloride as tracer of solute transport in the aquifer–aquitard system in the Pearl River Delta,China

A 1D numerical model is constructed to investigate the impact of sedimentation and sea level changes on transport of Cl^? in the aquifer–aquitard system in the Pearl River Delta (PRD), China. The model simulates the evolution of the vertical Cl^? concentration profiles during the Holocene. Sedimentation is modeled as a moving boundary problem. Chloride concentration profiles are reconstructed for nine boreholes, covering a wide area of the PRD, from northwest to southeast. Satisfactory agreement is obtained between simulated and measured Cl^? concentration profiles. Diffusion solely is adequate to reproduce the vertical Cl^? concentration profiles, which indicates that diffusion is the regionally dominant vertical transport mechanism across the aquitards in the PRD. The estimated effective diffusion coefficients of the aquitards range from 2.0?×?10^–11 to 2.0?×?10^–10 m²/s. The effective diffusion coefficients of the aquifers range from 3.0?×?10^–11 to 4.0?×?10^–10 m²/s. Advective transport tends to underestimate Cl^? concentrations in the aquitard and overestimate Cl^? concentrations in the basal aquifer. The results of this study will help understand the mechanisms of solute transport in the PRD and other deltas with similar geological and hydrogeological characteristics.     

Helium measurements of pore fluids obtained from the San Andreas Fault Observatory at Depth (SAFOD, USA) drill cores

⁴He accumulated in fluids is a well established geochemical tracer used to study crustal fluid dynamics. Direct fluid samples are not always collectable; therefore, a method to extract rare gases from matrix fluids of whole rocks by diffusion has been adapted. Helium was measured on matrix fluids extracted from sandstones and mudstones recovered during the San Andreas Fault Observatory at Depth (SAFOD) drilling in California, USA. Samples were typically collected as subcores or from drillcore fragments. Helium concentration and isotope ratios were measured 4?C6 times on each sample, and indicate a bulk ⁴He diffusion coefficient of 3.5?±?1.3?×?10^?C8 cm²?s^?C1 at 21°C, compared to previously published diffusion coefficients of 1.2?×?10^?C18 cm²?s^?C1 (21°C) to 3.0?×?10^?C15 cm²?s^?C1 (150°C) in the sands and clays. Correcting the diffusion coefficient of ⁴He_water for matrix porosity (??3%) and tortuosity (??6?C13) produces effective diffusion coefficients of 1?×?10^?C8 cm^2?s^?C1 (21°C) and 1?×?10^?C7 (120°C), effectively isolating pore fluid ⁴He from the ⁴He contained in the rock matrix. Model calculations indicate that <6% of helium initially dissolved in pore fluids was lost during the sampling process. Complete and quantitative extraction of the pore fluids provide minimum in situ porosity values for sandstones 2.8?±?0.4% (SD, n?=?4) and mudstones 3.1?±?0.8% (SD, n?=?4).     

Computed and observed currents,elevations, and salinity in a branching estuary

A one-dimensional, hydrodynamical model of the Tamar Estuary shows good agreement with measured tidal elevations and currents. Computed currents are used to drive a one-dimensional moving-element model of the salt balance. The moving-element model overcomes the numerical difficulties associated with strong tidal advection. Axial distributions of salinity at high water, computed using the moving-element model, compare well with measurements. The modelled and observed high water salinity distributions in this macrotidal estuary show little dependence on tidal range. The major variability in salinity is due to runoff. This strong and rapid dependence on runoff is a consequence of short residence (or flushing) times. Typically, residence times are less than one day throughout the year in the upper 10 km of estuary. The residence times maximize in summer, reaching 14 d for the whole estuary. During high runoff winter periods residence times are less than 5 d. Mixing coefficients for the moving-element salinity model are deduced from salinity measurements. Dispersion coefficients at fixed locations along the estuary are deduced from solutions of the salinity model. The spatially-averaged coefficients at mean spring and neap tides are 180 and 240 m² s^?1, respectively, for average runoff. Therefore, spring-neap variations in dispersion are fairly small and show a negative correlation with tidal range. The spatially-averaged dispersion coefficients at mean tides vary from 150 to 300 m² s^?1 for typical summer and winter runoff, respectively. The increase in dispersion with runoff and the decrease with tidal range implies that buoyancy-driven currents generate an important component of the shear dispersion in this estuary.     

Flow parameters of turbidity currents in a low-sinuosity giant deep-sea channel

Flow parameters (velocity and density) for turbidity currents in the Northwest Atlantic Mid-Ocean Channel (NAMOC) have been determined based on two different approaches, channel geometry and grain-size distributions of turbidites. Channel geometry has been obtained by a quantitative morphological analysis of the NAMOC which shows three genetically different segments in the upper 2000 km: (1) an upper 350 km-long ‘equilibrium channel’, (2) a middle 700 km-long ‘modified equilibrium channel’and (3) a lower ‘basement-controlled channel’which is more than 1000 km-long. In contrast to other meandering submarine channels the NAMOC has very low sinuosities and gradients. A consistently higher right-hand levee limits mean flow velocities to 3ms^?1 and channel geometry indicates mean flow velocities of 0·86 m s^?1 that decrease within the equilibrium channel to 0·05 m s^?1. Grain-size distributions on the levees and in the channel suggest strong vertical velocity and density gradients for bank-full flows with velocities of up to 8 m s^?1 and excess densities up to 87 kg m^?3 at the base, and 0·45 m s^?1 and 4 kg m^?3 at the top. The internal shear produced by these strong vertical gradients results in a decoupling of the current head and body. Channel geometry appears to be mainly the result of the slowly moving dilute body of the current.     

Physical,chemical, and biological conditions associated with the early stages of the Lake Michigan vernal thermal front

An investigation of the thermal front in southeastern Lake Michigan during April 1988 revealed a dynamic physical, chemical, and biological environment. The front was observed approximately 4 km from the coast as a distinct gradient separating cold open-lake waters from warmer nearshore waters. Surface isotherms near the front were generally parallel to one another but skewed with respect to shore. Comparison between April 22 and 29 showed that the surface isothern pattern was modulated by wind stress. The pattern from April 29 showed signs of flow instabilities with horizontal scales of 1 km to 5 km. Surface drifter trajectories provided estimates of horizontal convergence at the fron which varied from 7×10^?6 s^?1 to 20×10^?6 s^?1. Inferred rates of downwelling, which ranged from 9.5 m d^?1 to 20.7 m d^?1, were sufficient to move a water parcel from the surface to the bottom in 2 to 6 d at the front. Convergent circulation was observed on both sampling dates despite contrasts in wind stress. Concentrations of chloride, soluble silica, and chlorophyll, which were always higher inshore, were 5% to 82% larger than offshore mean values. The aquatic environment just inshore of the thermal front was characterized by chlorphyll concentrations which exceeded 5.0 μg 1^?1 while concentrations offshore were between 1.0 μg 1^?1 and 2.0 μg 1^?1. A relatively uniform vertical structure in chlorophyll concentrations in the frontal zone was consistent with the observed convergence and inferred downwelling near the front.     

Estuary-specific variation in the air-water gas exchange coefficient for oxygen

Oxygen air-water gas exchange was measured using floating chambers in two shallow tidal estuaries of differing bathymetry and local terrain, near Waquoit Bay, Massachusetts (United States). The specific chamber design permitted measurements of gas flux in 15 min, allowing analysis of the relationship with wind speed and tidal stage. Exchange coefficients ranged from 0.5 to 2.5 g O₂·m^?2 h^?1 atm^?1 (equivalent to piston velocities of 1.5 to 7 cm h^?1) for wind speeds of 0.3 to 9 m s^?1 at 10 m elevation. While the relationships for each estuary appear linear (significant linear regressions with wind speed were shown for each estuary, and the slopes were different at the 99.5% confidence level), the range of speeds differed at the two sites and an exponential function of wind speed was consistent with the combined data from both estuaries. A power function of wind speed was not an acceptable model. The exchange coefficients for our estuaries are from 57% to as low as 9% of that predicted by previously published generic equations. Because the atmospheric correction can be significant in shallow, metabolically active coastal waters, we suggest that empirically determined relationships for gas exchange versus wind for a specific estuary are preferable to the predictions of the general equations. While the floating chamber method should be used cautiously, at low winds speeds (below 8 m s^?1) and in slowly flowing waters, it provides a convenient approach for quantifying these site-specific differences. The differences, especially those between shallow sheltered systems and the open waters best fit by some published relationships, are ecologically important and do not appear yet to be measurable by other methods.     

Surface chemistry and dissolution kinetics of glassy rocks at 25°C

The weathering rates and mechanisms of three types of glassy rocks were investigated experimentally at 25 °C, pH 1.0 to 6.2, and reaction times as much as to 3 months. Changes in major element chemistry were monitored concurrently as a function of time in the aqueous solution and within the near surface region of the glass. Leach profiles, obtained by a HF leaching technique, displayed near-surface zones depleted in major cations. These zones increased in depth with increasing time and decreasing pH of reactions. Release rates into the aqueous solution were parabolic for Na and K and linear for Si and Al. A coupled weathering model, involving surface dissolution with concurrent diffusion of Na, K, and Al, produced a mass balance between the aqueous and glass phases. Steady state conditions are reached at pH 1.0 after approximately 3 weeks of reaction. Steady-state is not reached even after 3 months at pH 6.2.An interdiffusion model describes observed changes in Na diffusion profiles for perlite at pH 1.0. The calculated Na self-diffusion coefficient of 5 × 10^?19 cm²·s^?1 at 25°C approximates coefficients extrapolated from previously reported high temperature data for obsidian. The self-diffusion coefficient for H₃O⁺, 1.2 × 10^?20 cm²·s^?1, is similar to measured rates of water diffusion during hydration of obsidian to form perlite.     

Frequency spectra and intermittency of the turbulent suspension process in a sand-bed river

Previous research suggests that the turbulence-driven suspension process in sand-bed channels is dominated by intermittent, energetic eddies with length scales of the order of channel depth. Because of the scarcity of data on the turbulent suspension process in alluvial channels, the possible variability in suspension intermittency and turbulent frequency content due to contrasts in flow depth, velocity or bedforms remains unclear however. The present study analyses eddy correlation suspension signals from seven deployments in varied flow conditions around a sandy meander bend. Deployment depths at near-bankfull flow stages varied from 2 to 5.5 m, velocities at 0.75–1 m height from 0.6 to 0.9 m s^?1 and local mean suspended sand concentrations ranged from 30 to 150 mg L^?1 in the intermittence and spectral content of sand suspension between the various deployments are analysed and results are compared with previously published findings. Study data suggest that the dominant eddy sizes involved in sediment mixing across the sensor level are consistently of the order of 1–5 times flow depth and lie within the ‘energy-bearing’ turbulent range. When sand suspension is analysed in the time domain in the various deployments, energetic, burst-like suspension events occupying only 1-5% of the record duration account for 20-90% of the suspension work. The degree of intermittence in the suspension process was observed to increase in deeper flows, where mixing events contributing extreme vertical sediment fluxes appear to be relatively more frequent.     

Tropical cyclone activity and western North Atlantic stratification over the last millennium: a comparative review with viable connections

Tropical cyclones (TC) are recognized to modify the thermal structure of the upper ocean through the process of vertical mixing. Assessing the role this mixing plays in the overall stratification of the upper ocean is difficult, due to the relatively short and incomplete instrumental record. Proxy records for both TC landfalls and oceanographic stratification are preserved within the geological record and provide insight for how past changes in TC‐induced mixing have potentially affected water column structure prior to the instrumental record. Here we provide the first comparison between previously published paleo‐reconstructions of vertical ocean density and tropical cyclone activity from the western North Atlantic. A prominent lull in TC activity has been observed prior to approximately 1700 CE that extends back several centuries. This interval of low TC activity is shown to be concurrent with the timing of increased ocean stratification near Great Bahama Bank, potentially due in part to reduced TC‐induced mixing. To test whether this relationship is feasible, we present numerical results from a coarse‐resolution ocean general circulation model experiment isolating the effect of TC surface wind forcing on the upper ocean. An anomaly of roughly 0.12 kg m^?3 in vertical stratification occurs above and below the mixed layer for model runs with and without TC mixing. This anomaly is roughly 25% of the entire paleo‐density signal observed just prior to 1700 CE. These results suggest that TC mixing alone cannot completely explain the density anomaly observed prior to 1700 CE, but support TC variability as an important contributor to enhancing oceanic stratification during this interval. Copyright © 2011 John Wiley & Sons, Ltd.     

Mixing of metals during stripping of galactic gaseous halos

The mixing of metals in the intergalactic gas when a galaxy with a metal-rich envelope moves through the intergalactic medium is analyzed. Two simple models for the initial distribution of metals are considered. In the first case, the metals are concentrated in a fairly thin envelope with thickness ΔR _s =1 kpc, outer radius R _s =31 kpc, and metallicity Z=10^?3. In the second case, material with the same metallicity uniformly fills an entire spherical region of radius R _s . After 2.85 Gyr, the metals are distributed over a fairly extended volume with a typical size of ?200 kpc in the direction of the motion of the intergalactic gas, with a mean metallicity of ?4.6×10^?4 in metal-enriched regions. However, the distribution of metals remains extremely nonuniform, so that the main contribution to the overall metallicity is provided by metal-rich islands Z?6×10^?4 that occupy only ～10% of the total mixing volume. Moreover, metal-free regions remain in this volume.     

DIffusion of Eu and Gd in basalt and obsidian

Tracer diffusion coefficients of ¹⁵³Gd and ¹⁵²Eu in olivine tholeiite have been determined at temperatures between 1150 and 1440°C. The results are identical for both tracers within experimental error. Between 1440 and 1320°C the diffusion coefficients are given by D(Eu, Gd) = 0.058 exp(?40,600/ RT). Between 1320 and 1210°C, the diffusion coefficients are constant at D = (1.4 ± 0.4) × 10^?7 cm²s^?1 and between 1210 and 1150°C, the D values drop irregularly to 4 × 10^?9 cm²s^?1. The liquidus temperature (1270°C) lies within the region of constant D. Such anomalous behavior has not been encountered in previous studies of Ca, Sr, Ba and Co diffusion in basalt. To explain the constant D value near the liquidus, we speculate that the structure of the melt changes as a function of temperature in such a way that the normal temperature dependence of the diffusivity is compensated. For example, the rare earth ions may be displaced from their (high temperature) octahedral coordination sites to other sites where they are more readily dissociated and therefore become progressively more mobile. The behavior below 1210°C may be the result of relatively stable complexes or molecules in the melt or of the formation of a REE bearing crystalline phase that has so far escaped detection. Preliminary results for Eu diffusion in obsidian are D (Eu, 800°C) = 5 × 10^?13 cm² s^?1 and D (Eu, 950°C) = 1.5 × 10^?11 cm² s^?1. These data are consistent with an activation energy of 59 Kcal mole^?1. These low diffusivities indicate that the partitioning of REE in crystallizing intermediate and acidic melts may be controlled by diffusion in the melt rather than equilibrium between the crystal surface and the bulk melt.The diffusion data are applied to partial melting in the mantle, in an attempt to explain how LREE enriched tholeiites may be derived from a LREE depleted mantle source. In this model LREE diffuse from garnet bearing regions that have small melt fractions into garnet free regions that have relatively large melt fractions. REE diffusion is so slow that this process is quantitatively significant only in small partially molten bodies (diameter ～1 km or less) or in larger, but strongly flattened bodies. Internal convective motion during diapiric rise would also increase the efficiency of the process.     

Characterization of a Malaysian mangrove estuary

The Sungai Merbok estuary, in wet tropical Peninsular Malaysia, borders the Straits of Malacca. Tide, current, and salinity data are used to describe the salient hydrographic features of the mangrove-fringed system. The Sungai Merbok estuary is characterized by a 1.7 m semidiurnal tide with a 0.16 form number, peak currents of 1.3 m s^?1, and mean freshwater discharge of 20 m³ s^?1. The system is classified as 2a/2b estuary (Hansen and Rattray 1966) or 1a/1b during periods of low runoff. Gravitational circulation is highly variable (but coincides with the neap stratification) and vertical stratification varies from 10^?2 to 1. The estuary displays a pronounced fortnightly neap-spring stratification-destratification cycle. The effective longitudinal dispersion coefficient is approximately 100 m² s^?1.     

Delamination and ultra‐deep subduction of continental crust: constraints from elastic wave velocity and density measurement in ultrahigh‐pressure metamorphic rocks

Thirty‐three samples, including 22 eclogites, collected from the Dabie ultrahigh‐pressure (UHP) metamorphic belt in eastern China, have been studied for seismic properties. Compressional (V_p) and shear wave (V_s) velocities in three mutually perpendicular directions under hydrostatic pressures up to 1.0 GPa were measured for each sample. At 1.0 GPa, V_p (7.5–8.4 km s^?1), V_s (4.2–4.8 km s^?1), and densities (3.2–3.6 g cm^?3) in the UHP eclogites are higher than those of UHP orthopyroxenite (7.3–7.5 km s^?1, 4.1–4.3 km s^?1, 3.2–3.3 g cm^?3, respectively) and HP eclogites (7.1–7.9 km s^?1, 4.0–4.5 km s^?1, 3.1–3.5 g cm^?3, respectively). Kyanitites (with 99.5% kyanite) show extremely high velocities and density (9.37 km s^?1, 5.437 km s^?1, 3.581 g cm^?3, respectively). The eclogites show variation of V_p‐ and V_s‐anisotropy up to 9.70% and 9.17%, respectively. Poisson’s ratio (σ) ranges from 0.218 to 0.278 (with a mean of 0.255) for eclogites, 0.281–0.298 for granulites and 0.248 to 0.255 for amphibolites. The σ values for serpentinite (0.341) and marble (0.321) are higher than for other lithologies. The elastic moduli K, G, E of kyanitite were obtained as 163, 102 and 253 GPa, respectively. The V_p and density of representative UHP metamorphic rocks (eclogite & kyanitite) were extrapolated to mantle depth (15 GPa) following a reasonable geotherm, and compared to the one dimension mantle velocity and density model. The comparison shows that V_p and density in eclogite and kyanitite are greater than those of the ambient mantle, with differences of up to ΔV_p > 0.3 km s^?1 and Δρ > 0.3–0.4 g cm^?3, respectively. This result favours the density‐induced delamination model and also provides evidence in support of distinguishing subducted high velocity materials in the upper mantle by means of seismic tomography. Such ultra‐deep subduction and delamination processes have been recognized by seismic tomography and geochemical tracing in the postcollisional magmatism in the Dabie region.     

Dust emission and transport in Mali, West Africa

Vertical dust fluxes were measured in the Inland Delta region of the Niger River, Mali, West Africa, during April-June, 1989 and 1990. Measurements of dust flux represent, for the most part, non-dust storm conditions or ‘dust haze’ periods. The observed concentration versus height relationships are similar to data presented by other investigators. The relationship between wind shear velocity (u_*, m s^?1) and vertical dust flux (F, μg m^?2 s^?1) can be described by a relationship in which F is proportional to u_*⁴. However, there is considerable scatter within the data set which is attributable to textural controls and surface conditions. The vertical dust fluxes measured in Mali are compared to dust fluxes measured in Texas, USA, and Yukon Territories, Canada. The significantly different values for the constant of proportionality (a) in the F α a u_*⁴ relationship for these geographically diverse areas is a function of surficial controls on the release of sediments to the air stream. Dust concentrations measured in Mali were found to be uniformly high and in general exceed WHO health standards for acceptable total suspended particulate loadings. As a result background dust may be considered a long term stress on health for the people of this region of Mali.     

Amino acid composition of suspended particles,sediment-trap material,and benthic sediment in the Potomac Estuary

Sediment trap deployments in estuaries provide a method for estimating the amount of organic material transported to the sediments from the euphotic zone. The amino acid composition of suspended particles, benthic sediment, and sediment-trap material collected at 2.4 m, 5.8 m, and 7.9 m depths in the Potomac Estuary was determined in stratified summer waters, and in well-mixed oxygenated waters (DO) in late fall. The total vertical flow, or flux, of material into the top traps ranged from 3 g m^?2 d^?1 in August to 4.9 g m^?2 d^?1 in October. The carbon and nitrogen fluxes increased in the deepest traps relative to the surface traps during both sampling periods, along with that of the total material flux (up to 47.3 g m^?2 d^?1 in the deepest trap), although the actual weight percent of organic carbon and organic nitrogen decreased with depth. Amino acid concentrations ranged from 129 mg g^?1 in surface water particulate material to 22 mg g^?1 in particulate material in 9-m-deep waters and in the benthic sediment. Amino acid concentrations from 2.4-mg-depth sediment traps averaged 104±29 mg g^?1 in stratified waters and 164±81 mg g^?1 in well-mixed waters. The deep trap samples averaed, 77.3±4.8 mg g^?1 amino acids in summer waters and 37±16 mg g^?1 in oxygenated fall waters. Amino acids comprised 13% to 39% of the organic carbon and 12% to 89% of the orgnaic nitrogen in these samples. Analysis of the flux results suggest that resuspension combined with lateral advection from adjacent slopes can account for up to 27% of the material in the deep traps when the estuary was well-mixed and unstratified. When the estuary was stratified in late summer, the amino acid carbon produced by primary productivity in the euphotic zone decreased by 85% (86% for total organic carbon) at the pycnocline at 6 m depth, leaving up to 15% of the vertical organic flux available for benthic sediment deposition.     

The influence of wind speed on surface layer stability and turbulent fluxes over southern Indian peninsula station

Surface to atmosphere exchange has received much attention in numerical weather prediction models. This exchange is defined by turbulent parameters such as frictional velocity, drag coefficient and heat fluxes, which have to be derived experimentally from high-frequency observations. High-frequency measurements of wind speed, air temperature and water vapour mixing ratio (eddy covariance measurements), were made during the Integrated Ground Observation Campaign (IGOC) of Cloud Aerosol Interaction and Precipitation Enhancement Experiment (CAIPEEX) at Mahabubnagar, India (16^°44^′N, 77^°59^′E) in the south-west monsoon season. Using these observations, an attempt was made to investigate the behaviour of the turbulent parameters, mentioned above, with respect to wind speed. We found that the surface layer stability derived from the Monin–Obukhov length scale, is well depicted by the magnitude of wind speed, i.e., the atmospheric boundary layer was under unstable regime for wind speeds >4 m s^?1; under stable regime for wind speeds <2 m s^?1 and under neutral regime for wind speeds in the range of 2–3 m s^?1. All the three stability regimes were mixed for wind speeds 3–4 m s^?1. The drag coefficient shows scatter variation with wind speed in stable as well as unstable conditions.     

A chemical survey of the Mississippi estuary

A “snap shot” survey of the Mississippi estuary was made during a period of low river discharge, when the estuarine mixing zone was within the deltaic channels. Concentrations of H⁺, Ca²⁺, inorganic phosphorus and inorganic carbon suggest that the waters of the river and the low salinity (<5‰) portion of the estuary are near saturation with respect to calcite and sedimentary calcium phosphate. An input of oxidized nitrogen species and N₂O was observed in the estuary between 0 and 4‰ salinity. The concentrations of dissolved NH₄ ⁺ and O₂, over most of the estuary, appeared to be influenced by decomposition of terrestrial organic matter in bottom sediments. The estuarine bottom also appears to be a source of CH₄ which has been suggested to originate from petroleum shipping and refining operations. Estuarine mixing with offshore Gulf waters was the dominant influence on distributions of dissolved species over most of the estuary (i.e., from salinities >5‰). The phytoplankton abundance (measured as chlorophylla) increased as the depth of the mixed layer decreased in a manner consistent with that expected for a light-limited ecosystem. Fluxes of NO₃ ^?+NO₂ ^? and soluble inorganic phosphorus to the Gulf of Mexico were estimated to be 3.4±0.2×10³ g N s^?1 and 1.9±0.2 g P s^?1 respectively, at the time of this study.     

Groundwater flow dynamics of weathered hard-rock aquifers under climate-change conditions: an illustrative example of numerical modeling through the equivalent porous media approach in the north-western Pyrenees (France)

A numerical groundwater model of the weathered crystalline aquifer of Ursuya (a major water source for the north-western Pyrenees region, south-western France) has been computed based on monitoring of hydrological, hydrodynamic and meteorological parameters over 3 years. The equivalent porous media model was used to simulate groundwater flow in the different layers of the weathered profile: from surface to depth, the weathered layer (5?·?10^?8?≤?K?≤?5?·? 10^?7 m s^?1), the transition layer (7?·?10^?8?≤?K?≤?1?·? 10^?5 m s^?1, the highest values being along major discontinuities), two fissured layers (3.5?·?10^?8?≤?K?≤?5?·?? 10^?4 m s^?1, depending on weathering profile conditions and on the existence of active fractures), and the hard-rock basement simulated with a negligible hydraulic conductivity (K = 1 10 ^?9 ). Hydrodynamic properties of these five calculation layers demonstrate both the impact of the weathering degree and of the discontinuities on the groundwater flow. The great agreement between simulated and observed hydraulic conditions allowed for validation of the methodology and its proposed use for application on analogous aquifers. With the aim of long-term management of this strategic aquifer, the model was then used to evaluate the impact of climate change on the groundwater resource. The simulations performed according to the most pessimistic climatic scenario until 2050 show a low sensitivity of the aquifer. The decreasing trend of the natural discharge is estimated at about ?360 m³ y^?1 for recharge decreasing at about ?5.6 mm y^?1 (0.8 % of annual recharge).