Long-Term Changes in Water Quality and Productivity in the Patuxent River Estuary: 1985 to 2003

We conducted a quantitative assessment of estuarine ecosystem responses to reduced phosphorus and nitrogen loading from sewage treatment facilities and to variability in freshwater flow and nonpoint nutrient inputs to the Patuxent River estuary. We analyzed a 19-year dataset of water quality conditions, nutrient loading, and climatic forcing for three estuarine regions and also computed monthly rates of net production of dissolved O₂ and physical transport of dissolved inorganic nitrogen (DIN) and phosphorus (DIP) using a salt- and water-balance model. Point-source loading of DIN and DIP to the estuary declined by 40–60% following upgrades to sewage treatment plants and correlated with parallel decreases in DIN and DIP concentrations throughout the Patuxent. Reduced point-source nutrient loading and concentration resulted in declines in phytoplankton chlorophyll-a (chl-a) and light-saturated carbon fixation, as well as in bottom-layer O₂ consumption for upper regions of the estuary. Despite significant reductions in seaward N transport from the middle to lower estuary, chl-a, turbidity, and surface-layer net O₂ production increased in the lower estuary, especially during summer. This degradation of water quality in the lower estuary appears to be linked to a trend of increasing net inputs of DIN into the estuary from Chesapeake Bay and to above-average river flow during the mid-1990s. In addition, increased abundance of Mnemiopsis leidyi significantly reduced copepod abundance during summer from 1990 to 2002, which favored increases in chl-a and allowed a shift in total N partitioning from DIN to particulate organic nitrogen. These analyses illustrate (1) the value of long-term monitoring data, (2) the need for regional scale nutrient management that includes integrated estuarine systems, and (3) the potential water quality impacts of altered coastal food webs.     

Coastal Habitat Use and Residency of Juvenile Atlantic Sharpnose Sharks (Rhizoprionodon terraenovae)

Coastal habitat use and residency of a coastal bay by juvenile Atlantic sharpnose sharks, Rhizoprionodon terraenovae, were examined by acoustic monitoring, gillnet sampling, and conventional tag–recapture. Acoustic monitoring data were used to define the residency and movement patterns of sharks within Crooked Island Sound, Florida. Over 3 years, sharks were monitored for periods of 1–37 days, with individuals regularly moving in and out of the study site. Individual sharks were continuously present within the study site for periods of 1–35 days. Patterns of movement could not be correlated with time of day. Home range sizes were typically small (average?=?1.29 km²) and did not vary on a yearly basis. Gillnet sampling revealed that juvenile Atlantic sharpnose sharks were present in all habitat types found within Crooked Island Sound, and peaks in abundance varied depending on month within a year. Although telemetry data showed that most individuals remained within the study site for short periods of time before emigrating, conventional tag–recapture data indicates some individuals return to Crooked Island Sound after extended absences (maximum length?=?1,352 days). Although conventional shark nursery theory suggests small sharks remain in shallow coastal waters to avoid predation, juvenile Atlantic sharpnose sharks frequently exited from protected areas and appear to move through deeper waters to adjacent coastal bays and estuaries. Given the high productivity exhibited by this species, the benefit gained through a nursery that reduces predation may be limited for this species.     

Groundwater flow modeling for effective implementation of landslide stability enhancement measures

This paper deals with groundwater hydrology at a prominent fracture zone landslide slope (Nuta–Yone landslides) in Japan with an objective to explore an efficient method for the application of landslide stability enhancement measures. The correlation analyses between the hydrological parameters and ground surface movement data at this landslide resulted in low correlation values indicating that the geological formation of the area is extremely complex. For the purpose of understanding the groundwater flow behavior in the landslide area, a three-dimensional transient groundwater flow model was prepared for a part of the landslide slope, where the levels of effectiveness of applied landslide stability enhancement measures (in the form of multilayered deep horizontal drains) are different, and was calibrated against the measured water surface elevations at different piezometer locations. The parameter distributions in the calibrated model and the general directions of the groundwater flow in terms of flow vectors and the results of particle tracking at the model site were interpreted to understand the reasons for variations in effectiveness of existing landslide stability enhancement measures and to find potentially better locations for the implementation of future landslide stability enhancement measures. From the modeling results, it was also understood that groundwater flow model can be effectively used in better planning and locating the landslide stability enhancement measures.     

A New Algorithm for Wind-Vector Retrieval From Scatterometers

A new efficient algorithm for retrieving wind-vector solutions from scatterometers is developed based on a criterion of minimum normalized standard deviation (NSD) of wind speed derived from backscatter measurements using a geophysical model function (GMF). Its performance has been evaluated through simulations using QSCAT-1 GMF and the QuikSCAT observational geometry. The present algorithm, named the NSD algorithm, is found to be computationally more efficient (two to three times) besides being at par with the maximum-likelihood estimator (MLE) algorithm in terms of retrieval skill, retrieval errors, and distribution of solutions, on the basis of simulations as well as comparison of limited QuikSCAT-data-derived winds with National Centers for Environmental Prediction and European Centre for Medium-Range Weather Forecasts model winds. Simulation results and analysis of sample QuikSCAT data are presented.      

Developing a SPASE Query Language

The advent of the Virtual Observatory has begun an evolution in the space physics data environment. A number of nascent and discipline specific Virtual Observatories have started to emerge with an emphasis on data search and retrieval. As this new data environment takes shape an emphasis will be placed on interdisciplinary communication in attempts to address large scale and global problems. To this end we formulate the development of a query language to facilitate Virtual Observatory to Virtual Observatory communication. Furthermore, we outline the goals of such a language, how it would work and how existing community efforts can be leveraged to speed the development of this query language.

The gravitational field of topographic-isostatic masses and the hypothesis of mass condensation II-the topographic-isostatic geoid

In Part I we focussed on a convergent representation of the gravitational potential generated bytopographic masses on top of the equipotential surface atMean Sea Level, thegeoid, and by those masses which compensate topography. Topographic masses have also been condensated, namely represented by a single layer. Part II extends the computation of the gravitational field of topographic-isostatic masses by a detailed analysis of itsforce field in terms ofvector-spherical harmonic functions. In addition, the discontinuous mass-condensated topographic gravitational force vector (head force) is given. Once we identify theMoho discontinuity asone interface of isostatically compensated topographical masses, we have computed the topographic potential and the gravitational potential which is generated by isostatically compensated masses atMean Sea Level, the geoid, and illustrated by various figures of geoidal undulations. In comparison to a data oriented global geoid computation ofJ. Engels (1991) the conclusion can be made that the assumption of aconstant crustal mass density, the basic condition for isostatic modeling, does not apply. Insteaddensity variations in the crust, e.g. betweenoceanic and continental crust densities, have to be introduced in order to match the global real geoid and its topographic-isostatic model. The performed analysis documents that thestandard isostatic models based upon aconstant crustal density areunreal.     

The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

Testing for quantized redshifts. I. The project

A project intended to examine the long-standing claims that extragalactic redshifts are periodic or quantized was initiated some years ago at the Royal Observatory, Edinburgh. The approach taken is outlined, and the main conclusions to date are summarized. The existence of a galactocentric redshift quantization is confirmed at a high confidence level.     

The qualitative zoning record of minerals. A method for determining the duration of metamorphic events?

Summary The efficacy of ion exchange in petrological systems, for example the Fe-Mg exchange between garnet and biotite, is a function of grainsize, temperature, rates of temperature change and diffusion parameters. The combination of these variables determines the final zoning profile of minerals. Therefore, zoning profiles may be used to derive one of these variables if the others are known. For example, a mineral grain that experienced a short thermal event may still preserve a zoning profile characteristic of the heating path. In contrast, if that grain experienced a long thermal event it may develop a zoning profile typical of the cooling path. Conversely, for the same temperature-time cycle, large grains may not completely equilibrate at the metamorphic peak, and in smaller grains any previous record may have been erased. This is commonly observed in natural rocks where different grains within one thin section preserve often qualitatively different zoning profiles. Thus, a critical grainsize,l _crit, may exist that separates grains with qualitatively different zoning profiles so that grains of sizel > l _crit still retain zoning information about the heating path and grains of sizel <l _crit contain only information about the cooling path. If the critical grainsize can be measured and an independent estimate for the peak metamorphic temperature exists, the duration of the thermal event may, in principle, be estimated.The applicability of this method to natural garnets is hampered by a range of uncertainties. However, a parameterisation of the critical grain size may be used to illustrate the relative importance of grain size, temperature and event duration to the equilibration of minerals. In this paper, the critical grainsize is parameterized for various temperature-time cycles and its dependence on diffusion parameters is discussed. It is shown that, for Barrovian conditions,l _crit between 0.1 and 1 mm separates garnets retaining prograde information from garnets retaining cooling path information. For illustration, we compare the results with the critical grainsize of garnets from two metamorphic terranes, the Prydz Bay region (Antarctica) and the Koralm complex (Eastern Alps). Despite the large range of uncertainties attached to the method, it is shown that the critical grainsize of garnets in both terranes is consistent with a very short duration of the last thermal event that affected the two regions.

---

Die qualitative Form von Zonierungsprofilen in Mineralen. Information über die Dauer metamorpher Ereignisse?

Zusammenfassung Zonierungsprofile von Paragenesen im Ionenaustausch, zum Beispiel das Fe-Mg Austauschsgleichgewicht zwischen Granat und Biotit, sind eine Funktion von (i) Korngröße, (ii) Temperatur, (iii) Kühl- und Heizgeschwindigkeit sowie (iv) den Diffusionskonstanen. Es sollte daher möglich sein, Zonierungsprofile dazu zu verwenden, einen dieser Parameter zu bestimmen, wenn die anderen bekannt sind. Zonierungsprofile von Kristallen die ein kurzes thermisches Ereignis erfahren haben, mögen daher prograde Information beinhalten, wogegen derselbe Kristall ein retrogrades Profil aufweisen mag, wenn er einem langlebigen thermischen Ereignis unterlag (prograd und retrograd ist hier als: bei ansteigender Temperatur und bei abfallender Temperatur definiert). In Kristallen die dem gleichen Temperatur-Zeitpfad unterworfen waren, könnten kleine Körner nur die Kühlgeschichte dokumentieren, wogegen große noch Zonierungsinformation vom Heizpfad aufweisen. In natürlichen Gesteinen wird das oft dadurch beobachtet, daß verschiedene Körner auseinem Dünnschliff qualitativ verschiedene Zonierungsprofile aufweisen. Es ist daher möglich eine kritische Korngöße,l _crit, zu definieren, die Korngrößen mit qualitativ verschiedenen Zonierungsprofilen voneinander trennt. Körner mit einem Durchmesserl > l _crit haben, zumindest teilweise, noch prograde Profile, wogegen Körner mitl <l _crit nur retrograde Information dokumentieren. Wenn man diese kritische Korngröße messen kann, sollte sie dazu benutzt werden können, etwas über die Dauer des thermischen Ereignisses auszusagen.Die Anwendbarkeit dieser Methode ist durch eine Reihe von Fehlern limitiert. Nichtsdestotrotz, ist eine Parameterisierung der kritischen Korngröße nützlich um die relative Wichtigkeit von Korngröße, Temperatur und Zeit, für die Entwicklung von Zonierungsprofilen, zu illustrieren. Unser Modell kann dazu benutzt werden, um die Größenordnung von Granaten abzuschätzen, die noch prograde Information dokumentieren können. Die Abhängigkeit der kritischen Korngröße von Diffusionsparametern und der Form des Temperatur-Zeit Pfades wird ebenfalls diskutiert. Es wird gezeigt, daß für Metamorphose in der mittleren Grünschiefer und Amphibolit fazies,l _crit zwischen 0.1 and 1 mm prograd zonierte von retrograd zonierten Granaten trennt. Um das Modell zu illustrieren, werden die Ergebnisse auf die Prydz Bay Region (Antarktis) und den Koralm Komplex (Ostalpen) angewendet. Trotz der großen Fehlergrenzen der Methode glauben wir zeigen zu können, daß die kritische Korngröße beider Terrains andeutet, daß das letzte thermische Ereignis in beiden Gebieten kurzlebig war.

---

---

With 5 Figures     

LREE distribution in perovskite,apatite and titanite from South West Ugandan xenoliths and kamafugite lavas

Summary In-situ microprobe LREE analyses of perovskite and titanite (La, Ce, Nd), and apatite (La, Ce), from SW Ugandan clinopyroxenite xenoliths and kamafugite lavas indicate that LREE distribution in these minerals is determined by a number of factors related to their different parageneses: In particular LREE content is affected by whether the LREE-bearing minerals have crystallised from metasomatic carbonate or from silicate (i.e. metasomatic or magmatic) melts in the mantle. In this situation LREE partition favours carbonate over silicate melts. Distribution of LREE in perovskite and apatite crystallised from magmatic mantle melts or mantle-derived lavas is chiefly determined by preference of LREE for perovskite > apatite > titanite. LREE zoning in perovskite is influenced by changes in melt structure: increasing melt polymerisation enhancing mineral_LREE/melt_LREE partition into perovskite rims in magmatic xenoliths; decreasing melt polymerisation depleting LREE in lava perovskite rims. This zoning is reinforced by perovskite competition with apatite for LREE: perovskite (cores/rims) co-crystallising with apatite is reduced in LREE. There are 37 instances of perovskitewith Ce below detection while La and Nd levels are normal. These occur in both xenoliths and lavas; in grain zones or whole grains. Likewise Ce alone of the LREE is below detection in six out of ten titanite analyses. These observations are interpreted as evidence for increased f_{O 2}, Ce⁴ + being excluded from these mineral structures. Recognition of these various processes can elucidate the interpretation of bulk rock and bulk mineral LREE signatures in kamafugite volcanism.

---

LREE Verteilung in Perovskit, Apatit und Titanit aus Xenolithen und kamafugitischen Laven Südwest-Ugandas

Zusammenfassung In-situ LREE Analysen von Perovskit und Titanit (La, Ce, Nd) und Apatit (La, Ce) aus Klinopyroxenit-Xenolithen und kamafugitischen Laven Südwest-Ugandas zeigen, daß die LREE Verteilung in diesen Mineralen durch eine Vielzahl von Faktoren, die mit Unterschieden in den Paragenesen zusammenhängen, bestimmt wird: Der LREE-Gehalt wird im besonderen davon bestimmt, ob die LREE-führenden Minerale aus metasomatischen Karbonat- oder aus (metasomatischen oder magmatischen) Silikatschmelzen im Mantel auskristallisierten. Dabei erfolgt die LREE Fraktionierung zu Gunsten der Karbonatschmelzen. Die LREE-Verteilung von Perovskit und Apatit, die aus magmatischen Mantelschmelzen oder -laven kristallisierten, wird vorrangig durch den bevorzugten Einbau der LREE in Perovskit > Apatit > Titanit kontrolliert. Der LREE Zonarbau von Perovskit wird durch die Änderungen der Schmelzstruktur beinflußt: Verstärkte Schmelzpolymerisation führt zu verstärkter Mineral_LFEE/Schmelze_LREE Fraktionierung in den Perovskiträndern magmatischer Xenolithe, eine Abnahme der Schmelzpolymerisation hingegen resultiert in einer Abreicherung der LREE in den Perovskiträndern. Diese Art der Zonierung wird durch den Wettbewerb von Perovskit mit Apatit um die LREE verstärkt. Perovskit (Kerne/Ränder), der mit Apatit gemeinsam auskristallisierte, ist ärmer an LREE. 37 Fälle, in denenCe nicht nachweisbar war, La und Nd aber in normaler Konzentration auftreten, wurden sowohl in den Xenolithen als auch in den Laven gefunden; und zwar entweder in Kornbereichen oder in ganzen Körnern. Vergleichsweise liegt Ce nur in sechs von zehn Titanitproben unterhalb der Nachweisgrenze. Diese Beobachtungen werden als Hinweise auf erhöhte SauerstoffFugazitäten, bei denen Ce^4– aus der Mineralstruktur ausgeschlossen wird, angesehen.Ein Verständnis dieser verschiedenen Prozesse kann zur besseren Interpretation von LREE Gesamtgesteins- und Gesamtmineral-Signaturen in Kamafugiten beitragen.

---

---

With 3 Figures