A numerical simulation of residual circulation in Tampa Bay. Part I: Low-frequency temporal variations

The residual (time-average) salinity and circulation in a numerical ocean model of the Tampa Bay estuary are shown to experience significant temporal variation under realistic forcing conditions. A version of the Estuarine Coastal Ocean Model developed for Tampa Bay with 70 by 100 horizontal grid points and 11 sigma levels is examined for the years 2001–2003. Model output variables are averaged over the entire time of the simulation to generate long-term residual fields. The residual axial current is found to be dominated by the buoyancy-driven baroclinic circulation with an outflow (southwestward) at the surface and to the sides of the shipping channel, and an inflow (northeastward) usually occurring subsurface within or above the shipping channel. Averages over 30 d are used to examine variations in the residual fields. During the simulation the average surface salinity near the head of Tampa Bay varies with the freshwater inflow, from 12‰ to 33%. At the bay mouth salinity varies from 30%. to 36%.. A localized measure of the baroclinic circulation in the shipping channel indicates the residual circulation can vary strongly, attaining a magnitude triple the long-term mean value. The baroclinic circulation can be disrupted, going to near zero or even reversing, when the buoyancy-driven flow is weak and the surface winds are to the northeast. Three time periods, representing different environmental conditions, are chosen to examine these results in detail. A scaling argument indicates the relative strength of buoyancy versus wind as ΔρgH²(LC _Dω²)⁻¹, where δρ is head-to-mouth density difference across the bay,g is gravitational acceleration,H is depth,L is bay length,C _D is the surface wind drag coefficient, andw is wind speed. Tampa Bay is usually in the buoyancy dominated regime. The importance of winds in the weak-buoyancy case is demonstrated in an additional simulation without wind stress.     

Recolonization patterns of benthic assemblages after relict sand dredging in the central Tyrrhenian sea

Recolonization patterns of benthic assemblages after relict sand dredging were investigated in an area offshore of Capo d’Anzio (central Tyrrhenian Sea). Surveys were performed before, during, and after dredging in 2007. Sampling stations were located within and outside the dredging site in water depths between 40 and 56 m. Defaunation due to sand extraction was the main impact observed in the inner stations and in stations located in proximity of the dredging site. At the inner stations, 9 months after dredging, a significant increase in species richness and diversity was observed: the sediment removal led to an increase of the sandy sediment fraction, favouring the settlement of sabulicolous species. A decrease in number of individuals and species was also observed in most of the outer stations, probably due to fine sediment redeposition. Recolonization of macrobenthic assemblages was essentially achieved at inner stations both in terms of abundance and species richness, while at outer stations it was still in progress. Results of this study, providing a picture of recovery times and processes in the Tyrrhenian Sea, could be used to develop predictions of the effects of future offshore sand dredging projects in comparable areas. Monitoring programs and research have a key role to provide a more detailed overview of biological recovery processes and times in different regions and with different dredging intensity.     

A surge of Skobreen, Svalbard

Surging glaciers are common in Svalbard yet relatively few glaciers have been observed during a surge. This paper presents observations of the currently surging glacier Skobreen, in southern Spitsbergen. The study is based on examinations of new and archival photographs and maps. Skobreen, an 18 km² valley glacier terminating into the lower part of the glacier Paulabreen, has not been registered previously as a surging glacier. Skobreen experienced a build-up in its upper part, while there has been a lowering of the surface in the terminal region. Photographs from 1990 show incipient crevassing in the upper part. Photographs from 2003 show a slight advance of the terminus and marginal crevassing, indicating an initiation period of about 15 years for a surge of this glacier. In June 2005 transverse crevassing appeared in the upper part of the glacier, while the middle section moved as a block with strong shear margins and a pronounced drawdown of the ice surface. No traces of a surge front could be seen in the crevasse pattern. However, the crevasse pattern indicates an initiation area in the transition zone between the transverse crevassing in the upper part and the block of ice in the middle region.     

Procida volcanic history: new insights into the evolution of the Phlegraean Volcanic District (Campania region,Italy)

New stratigraphic, major- and trace-element, and Sr-, Nd- and Pb- isotopic data on volcanic deposits older than 14 ka from the island of Procida, Italy, are presented and compared with published analyses from the rest of the Phlegraean Volcanic District (PVD). Procida rocks range in composition from basalt to shoshonite and trachyte and show ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb ratios varying within the ranges 0.70523–0.70678, 0.512539–0.512630, 18.99–19.29, 15.67–15.69 and 39.10–39.39, respectively. The presence of a compositional gap in the range SiO₂=54–59 wt % is evidence of magma bimodality, and suggests that the feeding magmatic system was formed by at least two different reservoirs located at different depths. Geochemical and isotopic variations with increasing differentiation can be explained by fractional crystallization mechanisms, that in some cases were associated with crustal contamination that occurred in both deeper and shallower reservoirs; the most evolved magmas formed in the shallower one. Mixing/mingling processes also occurred. The variation of isotopic composition through time observed both for Procida and for Campi Flegrei and Ischia rocks is evidence of strong affinity between magmas that erupted on the entire PVD until about 40 ka. This indicates that they share a common origin and a common plumbing system. Most of the PVD eruptive centers active until about 40 ka lie within a NE-SW-oriented volcano-tectonic belt extending from the southeastern part of Ischia, through Procida and Torregaveta volcano to the northeastern sector of the present Campi Flegrei caldera. This not only indicates the existence of a link between regional structures and volcanism in the area, but also suggests that PVD magma genesis and evolution were strongly regulated by extensional tectonics. In the last 40 ka the mafic rocks erupted along this extensional structure – from Torregaveta and the islands of Ischia and Procida – indicate that it still represents an important crustal discontinuity that focuses mantle-derived magmas. Procida trachybasalts are enriched in large ion lithophile elements (LILE) and light and middle rare earth elements (LREE and MREE), and show slight negative anomalies in the high field strength elements (HFSE) relative to average MORB. A slight depletion in HREE is present. Trace element and isotope systematics can be referred to a lithospheric mantle source. The lithospheric mantle source carries intra-plate and slab-derived components, the latter probably inherited from a previous subduction event.     

Joint Inversion of Body-Wave Arrival Times and Surface-Wave Dispersion for Three-Dimensional Seismic Structure Around SAFOD

We incorporate body-wave arrival time and surface-wave dispersion data into a joint inversion for three-dimensional P-wave and S-wave velocity structure of the crust surrounding the site of the San Andreas Fault Observatory at Depth. The contributions of the two data types to the inversion are controlled by the relative weighting of the respective equations. We find that the trade-off between fitting the two data types, controlled by the weighting, defines a clear optimal solution. Varying the weighting away from the optimal point leads to sharp increases in misfit for one data type with only modest reduction in misfit for the other data type. All the acceptable solutions yield structures with similar primary features, but the smaller-scale features change substantially. When there is a lower relative weight on the surface-wave data, it appears that the solution over-fits the body-wave data, leading to a relatively rough V _s model, whereas for the optimal weighting, we obtain a relatively smooth model that is able to fit both the body-wave and surface-wave observations adequately.     

Nursery Habitat Shifts in an Estuarine Ecosystem<Emphasis Type="Bold">:</Emphasis> Patterns of Use by Sympatric Catfish Species

The seasonal and spatial distribution (density and biomass) of five size classes of two catfish species (Cathorops spixii and Cathorops agassizii) were studied along an estuarine ecocline to test the relative importance of the nursery function of each habitat. Seasonal vs. area interactions were significant for all size classes of both species. During the early rainy season, the middle estuary is an important nursery habitat for juveniles of both species. When environmental conditions change during the late rainy season, the C. spixii primary nursery habitat shifts to the lower estuary. During this period, juveniles of C. agassizii remain in the middle estuary. Another important ecological area is the upper estuary, which becomes a breeding, spawning and hatchery area during the late dry season for both species. The nursery function of habitats shifts according to the seasonal fluctuation of salinity and dissolved oxygen, and each species responds differently to this change.     

Scenarios of Earthquake-Generated Tsunamis for the Italian Coast of the Adriatic Sea

We calculated the expected impact on the Italian coast of the Adriatic Sea of a large set of tsunamis resulting from potential earthquakes generated by major fault zones. Our approach merges updated knowledge on the regional tectonics and scenario-like calculations of expected tsunami impact. We selected six elongated potential source zones. For each of them we determined a Maximum Credible Earthquake and the associated Typical Fault, described by its size, geometry and kinematics. We then let the Typical Fault float along strike of its parent source zone and simulated all tsunamis it could generate. Simulations are based on the solution of the nonlinear shallow water equations through a finite-difference technique. For each run we calculated the wavefields at specified simulation times and the maximum water height field (above mean sea level), then generated travel-time maps and maximum water height profiles along the target coastline. Maxima were also classified in a three-level code of expected tsunami threat. We found that the southern portion of Apulia facing Albania and the Gargano promontory are especially prone to the tsunami threat. We also found that some bathymetric features are crucial in determining the focalization-defocalization of tsunami energy. We suggest that our results be taken into account in the design of early-warning strategies.     

Determination of Carbon,Hydrogen, Nitrogen and Sulfur in Geological Materials Using Elemental Analysers

This article describes a series of methods developed for the determination of total carbon (C_Total), organic carbon (C_org), hydrogen, nitrogen and sulfur. The following elemental analysers were used: LECO model RC‐412 for the determination of organic carbon, total carbon and hydrogen; LECO model CS‐200 for the determination of total carbon and sulfur; LECO model TN‐400 for the determination of nitrogen; and LECO model TruSpec CHNS for the determination of organic carbon, total carbon, hydrogen, nitrogen and sulfur. Uncertainty and limits of detection and quantification were calculated for each method, as well as the running costs to define the most effective instrument for each material and each analyte. Accuracy was checked by the application of the Sutarno–Steger test. Finally, a compilation of the results obtained in the determination of C_Total, C_org, H, N and S in forty‐nine reference materials is presented.     

Effects of uncertainty of lithofacies, conductivity and porosity distributions on stochastic interpretations of a field scale tracer test

We investigate the importance of selecting two different methodologies for the determination of hydraulic conductivity from available grain-size distributions on the stochastic modeling of the depth-averaged breakthrough curve observed during a forced-gradient tracer test experiment. The latter was performed in the Lauswiesen alluvial aquifer, located near the city of Tübingen, Germany, by injecting NaBr into a well at a distance of about 50 m from a pumping well. We also examine the joint effect of the choice of the transport model adopted to describe solute transport at the site and the way the spatial distribution of porosity is assessed. In the absence of direct measurements of porosity, we consider: (a) the model used by Riva et al. (J Contam Hydrol 88:92–118, 2006; J Contam Hydrol 101:1–13, 2008), which relates the natural logarithms of effective porosity and conductivity through an empirical, experimentally-based, linear relationship derived for a nearby experimental site; and (b) a model based on a commonly used relationship linking the total porosity to the coefficient of uniformity of grain size distributions. Transport is described in terms of a purely advective process and/or by including mass exchange processes between mobile and immobile regions. Modeling of flow and transport is performed within a Monte Carlo framework, upon conceptualizing the aquifer as a random composite medium. Our results indicate that the model adopted to describe the correlation between conductivity and porosity and the way grain-sieve information are incorporated to depict the heterogeneous distribution of hydraulic conductivity can have relevant effects in the interpretation of the data at the site. All the conceptual models employed to describe the structural heterogeneity of the system and transport features can reasonably reproduce the global characteristics of the experimental depth-averaged breakthrough curve. Specific details, such as the peak concentration and the time of first arrival, can be better reproduced by a double porosity transport model when a correlation between conductivity and porosity based on grain size information at the site is considered. The best prediction of the late-time behavior of the measured breakthrough curves, in terms of the observed heavy tailing, is offered by directly linking porosity distribution to the spatial variability of particle size information.