Nutrient exchange across the sediment-water interface in the Potomac River estuary

The flux of ammonia, phosphate, silica and radon-222 from Potomac tidal river and estuary sediments is controlled by processes occurring at the sediment-water interface and within surficial sediment. Calculated diffusive fluxes range between 0·6 and 6·5 mmol m^?2 day^?1 for ammonia, 0·020 and 0·30 mmol m^?2 day^?1 for phosphate, and 1·3 and 3·8 mmol m^?2 day^?1 for silica. Measured in situ fluxes range between 1 and 21 mmol m^?2 day^?1 for ammonia, 0·1 and 2·0 mmol m^?2 day^?1 for phosphate, and 2 and 19 mmol m^?2 day^?1 for silica. The ratio of in situ fluxes to diffusive fluxes (flux enhancement) varied between 1·6 and 5·2 in the tidal river, between 2·0 and 20 in the transition zone, and from 1·3 to 5·1 in the lower estuary. The large flux enhancements from transition zone sediments are attributed to macrofaunal irrigation. Nutrient flux enhancements are correlated with radon flux enhancements, suggesting that fluxes may originate from a common region and that nutrients are regenerated within the upper 10–20 cm of the sediment column.The low fluxes of phosphate from tidal viver sediments reflect the control benthic sediment exerts on phosphorus through sorption by sedimentary iron oxyhydroxides. In the tidal river, benthic fluxes of ammonia and phosphate equal one-half and one-third of the nutrient input of the Blue Plains sewage treatment plant. In the tidal Potomac River, benthic sediment regeneration supplies a significant fraction of the nutrients utilized by primary producers in the water column during the summer months.     

Error analyses of resonant orbits for geodesy

An error analysis of resonant orbits for geodesy indicates that attempts to use resonance to recover high order geopotential coefficients may be seriously hampered by errors in the geopotential. This effect, plus the very high correlations (up to .99) of the resonant coefficients with each other and the orbital period in single satellite solutions, makesindividual resonant orbits of limited value for geodesy. Multiple-satellite, single-plane solutions are only a slight improvement over the single satellite case. Accurate determination of high order coefficients from low altitude resonant satellites requires multiple orbit planes and small drift-periods to reduce correlations and effects of errors of non-resonant geopotential terms. Also, the effects of gravity model errors on low-altitude resonant satellites make the use of tracking arcs exceeding two to three weeks of doubtful validity. Because high-altitude resonant orbits are less affected by non-resonant terms in the geopotential, much longer tracking arcs can be used for them.     

A comparison of the photospheric and interplanetary magnetic fields during the flight of Mariner IV

The polarity of the interplanetary magnetic field observed by the Mariner IV spacecraft is compared to the polarity of the photospheric magnetic field observed with the solar magnetograph at Mt. Wilson Observatory. Unlike the results obtained from observations during the flight of IMP-I, these polarities are not well correlated when the photospheric polarity is determined from data along a narrow latitudinal strip. It is suggested that the structure of the interplanetary field is often related to major features in the photospheric field that are observed over a broad range of solar latitudes.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

XUV Photometer System (XPS): Improved Solar Irradiance Algorithm Using CHIANTI Spectral Models

Solar soft X-ray (XUV) radiation is highly variable on all time scales and strongly affects Earth’s ionosphere and upper atmosphere; consequently, the solar XUV irradiance is important for atmospheric studies and for space weather applications. Although there have been several recent measurements of the solar XUV irradiance, detailed understanding of the solar XUV irradiance, especially its variability during flares, has been hampered by the broad bands measured in the XUV range. In particular, the simple conversion of the XUV photometer signal into irradiance, in which a static solar spectrum is assumed, overestimates the flare variations by more than a factor of two as compared to the atmospheric response to the flares. To address this deficiency in the simple conversion, an improved algorithm using CHIANTI spectral models has been developed to process the XUV Photometer System (XPS) measurements with its broadband photometers. Model spectra representative of quiet Sun, active region, and flares are combined to match the signals from the XPS and produce spectra from 0.1 to 40 nm in 0.1-nm intervals for the XPS Level 4 data product. The two XPS instruments are aboard NASA’s Solar Radiation and Climate Experiment (SORCE) and Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics (TIMED) satellites. In addition, the XPS responsivities have been updated for the latest XPS data processing version. The new XPS results are consistent with daily variations from the previous simple conversion technique used for XPS and are also consistent with spectral measurements made at wavelengths longer than 27 nm. Most importantly, the XPS flare variations are reduced by factors of 2 – 4 at wavelengths shorter than 14 nm and are more consistent, for the first time, with atmospheric response to solar flares. Along with the details of the new XPS algorithm, several comparisons to dayglow and photoelectron measurements and model results are also presented to help verify the accuracy of the new XUV irradiance spectra.     

Reanalysis of the Historic AFTAC Bolide Infrasound Database

We have recently digitized and partially reanalyzed the historic bolide infrasonic database. These 10 events were originally detected by the U.S. Air Force Technical Applications Center (AFTAC) from ∼1960 to 1974. In this paper we present the first preliminary reanalysis results for two of the 10 bolide events, namely the Revelstoke bolide of 3/31/1965 as well as the Prince Edward Islands (P.E.I). S. African bolide of 8/03/1963, which were among the largest bolides detected during this time period. These bolides have been investigated initially since they are most likely to have had a significant effect on the computed global influx rate of ReVelle (Global Infrasonic Monitoring of Large Bolides, pp 483–490, 2001) as indicated in Brown et al. (Nature, 420:314–316, 2002). We are in the process of recomputing all relevant infrasonic propagation quantities such as plane wave back azimuth, signal velocities, power spectra, spectrograms, as well as energy estimates using multiple techniques. In a future paper we will present a complete digital reanalysis of the AFTAC bolide infrasonic data and its final resulting global bolide influx implications.     

Acoustic-Gravity Waves from Bolide Sources

We have developed a new approach to modeling the acoustic-gravity wave (AGW) radiation from bolide sources. This first effort involves entry modeling of bolide sources that have available satellite data through procedures developed in ReVelle (Earth Moon Planets 95, 441–476, 2004a; in: A. Milani, G. Valsecchi, D. Vokrouhlicky (eds) NEO Fireball Diversity: Energetics-based Entry Modeling and Analysis Techniques, Near-earth Objects: Our Celestial Neighbors (IAU S236), 2007b). Results from the entry modeling are directly coupled to AGW production through line source blast wave theory for the initial wave amplitude and period at (at 10 blast wave radii and perpendicular to the trajectory). The second effort involves the prediction of the formation and or dominance of the propagation of the atmospheric Lamb, edge-wave composite mode in a viscous fluid (Pierce, J. Acoust. Soc. Amer. 35, 1798–1807, 1963) as a function of the source energy, horizontal range and source altitude using the Lamb wave frequency that was deduced directly during the entry modeling and that is used as a surrogate for the source energy. We have also determined that Lamb wave production by bolides at close range decreases dramatically as either the source energy decreases or the source altitude increases. Finally using procedures in Gill (Atmospheric-Ocean Dynamics, 1982) and in Tolstoy (Wave Propagation, 1973), we have analyzed two simple dispersion relationships and have calculated the expected dispersion for the Lamb edge-wave mode and for the excited, propagating internal acoustic waves. Finally, we have used the above formalism to fully evaluate these techniques for four large bolides, namely: the Tunguska bolide of June 30, 1908; the Revelstoke bolide of March 31, 1965; the Crete bolide of June 6, 2002 and the Antarctic bolide of September 3, 2004. Due to page limitations, we will only present results in detail for the Revelstoke bolide.