Numerical experiments on the stability of preplanetary disks

Gravitational stability of gaseous protostellar disks is relevant to theories of planetary formation. Stable gas disks favor formation of planetesimals by the accumulation of solid material; unstable disks allow the possibility of direct condensation of gaseous protoplanets. We present the results of numerical experiments designed to test the stability of thin disks against large-scale, self-gravitational disruption. The disks are represented by a distribution of about 6 × 10⁴ point masses on a two-dimensional (r, φ) grid. The motions of the particles in the self-consistent gravity field are calculated, and the evolving density distributions are examined for instabilities. Two parameters that have major influences on stability are varied: the initial temperature of the disk (represented by an imposed velocity dispersion), and the mass of the protostar relative to that of the disk. It is found that a disk as massive as 1M_⊙, surrounding a 1M_⊙ protostar, can be stable against long-wavelength gravitational disruption if its temperature is about 300°K or greater. Stability of a cooler disk requires that it be less massive, but even at 100°K a stable disk can have an appreciable fraction ($\text{～}\text{1}\text{3}$) of a solar mass.     

Interference in solar oscillations

We have analyzed magnetograph observations of the 5-min oscillations. We find that most of the oscillatory power is concentrated in space and frequency. Interference effects where these concentrations overlap can explain some of the variations in amplitude of the oscillation.Of the National Bureau of Standards and University of Colorado.     

A BME solution of the inverse problem for saturated groundwater flow

In most real-world hydrogeologic situations, natural heterogeneity and measurement errors introduce major sources of uncertainty in the solution of the inverse problem. The Bayesian Maximum Entropy (BME) method of modern geostatistics offers an efficient solution to the inverse problem by first assimilating various physical knowledge bases (hydrologic laws, water table elevation data, uncertain hydraulic resistivity measurements, etc.) and then producing robust estimates of the subsurface variables across space. We present specific methods for implementing the BME conceptual framework to solve an inverse problem involving Darcys law for subsurface flow. We illustrate one of these methods in the case of a synthetic one-dimensional case study concerned with the estimation of hydraulic resistivity conditioned on soft data and hydraulic head measurements. The BME framework processes the physical knowledge contained in Darcys law and generates accurate estimates of hydraulic resistivity across space. The optimal distribution of hard and soft data needed to minimize the associated estimation error at a specified sampling cost is determined. This work was supported by grants from the National Institute of Environmental Health Sciences (Grant no. 5 P42 ES05948 and P30ES10126), the National Aeronautics and Space Administration (Grant no. 60-00RFQ041), the Army Research Office (Grant no. DAAG55-98-1-0289), and the National Science Foundation under Agreement No. DMS-0112069.     

On the Dynamics of the Jovian Ionosphere and Thermosphere: III. The Modelling of Auroral Conductivity

Recent work has been concerned with calculating the three-dimensional ion concentrations and Pedersen and Hall conductivities within the auroral region of Jupiter for varying conditions of incident electron precipitation. Using the jovian ionospheric model, we present results that show the auroral ionospheric response to changing the incoming flux of precipitating electrons (for constant initial energy) and also the response to changing the initial energy (for both constant flux and constant energy flux). The results show that, for expected energy fluxes of precipitating particles, the average auroral integrated Pedersen conductivity attains values in excess of 1 mho. In addition, it is shown that electrons with an initial energy of around 60 keV are particularly effective at generating auroral conductivity: Particles of this energy penetrate most effectively to the layer of the jovian ionosphere at which the auroral conductivity is at a maximum.     

First Results from Viper: Detection of Small-scale Anisotropy at 40 GHz

Using all available major samples of Seyfert galaxies and their corresponding closely matched control samples of nonactive galaxies, we find that the bar ellipticities (or axial ratios) in Seyfert galaxies are systematically different from those in nonactive galaxies. Overall, there is a deficiency of bars with large ellipticities (i.e., "thin" or "strong" bars) in Seyfert galaxies compared to nonactive galaxies. Accompanied with a large dispersion due to small number statistics, this effect is strictly speaking at the 2 sigma level. To obtain this result, the active galaxy samples of near-infrared surface photometry were matched to those of normal galaxies in type, host galaxy ellipticity, absolute magnitude, and, to some extent, redshift. We discuss possible theoretical explanations of this phenomenon within the framework of galactic evolution, and, in particular, of radial gas redistribution in barred galaxies. Our conclusions provide further evidence that Seyfert hosts differ systematically from their nonactive counterparts on scales of a few kiloparsecs.     

A River Water Quality Model for Chlorophyll and Dissolved Oxygen that Accommodates Zebra Mussel Metabolism

The development, testing and application of a dynamic two-dimensional (longitudinal-vertical) mass balance model for dissolved oxygen (DO) and chlorophyll (Chl) for rivers is documented that for the first time accommodates both the oxygen demand and filtering effects of zebra mussels. The test system is a phytoplankton-rich section ( 15 km long) of the Seneca River, NY, that is believed to represent an upper bound of the impact of this exotic invader. Changes in common measures of water quality of the river brought about by the zebra mussel invasion are reviewed and related longitudinal patterns in DO, Chl, and Secchi disc transparency are documented. Model testing is supported by comprehensive measurements of DO, Chl, and various forcing conditions over a three-month period, and independent determinations of several model coefficients. Wide variations in the areal consumption rate of DO (ZOD; g·m^–2·d^–1) and filtering rate (m³·m^–2·d^–1) of zebra mussels, as determined through model calibration, occurred over the study period. Values of ZOD in areas with dense zebra mussel populations at times (e.g., > 50 g·m^–2·d^–1) were an order of magnitude greater than the sediment oxygen demand associated with organically enriched deposits. The value of determinations of these fluxes from model calibration procedures is evaluated within the context of the limitations of protocols presently available to support independent specification of these rates. Model analyses are conducted to evaluate the relative magnitude of source and sink processes for DO and Chl, the potential operation and implications of feedback from low DO levels on oxygen consumption by zebra mussels, and the sensitivity of model simulations to selected sources of uncertainty and variability. Model projections of oxygen resources of the river are presented in a probabilistic format in evaluating reductions in zebra mussel biomass that would be necessary to eliminate violations of standards and regain assimilative capacity.     

A spatially and temporally adaptive solution of Richards’ equation

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.     

Modeling studies of the upper ocean response to a tropical cyclone

A coupled ocean and boundary layer flux numerical modeling system is used to study the upper ocean response to surface heat and momentum fluxes associated with a major hurricane, namely, Hurricane Dennis (July 2005) in the Gulf of Mexico. A suite of experiments is run using this modeling system, constructed by coupling a Navy Coastal Ocean Model simulation of the Gulf of Mexico to an atmospheric flux model. The modeling system is forced by wind fields produced from satellite scatterometer and atmospheric model wind data, and by numerical weather prediction air temperature data. The experiments are initialized from a data assimilative hindcast model run and then forced by surface fluxes with no assimilation for the time during which Hurricane Dennis impacted the region. Four experiments are run to aid in the analysis: one is forced by heat and momentum fluxes, one by only momentum fluxes, one by only heat fluxes, and one with no surface forcing. An equation describing the change in the upper ocean hurricane heat potential due to the storm is developed. Analysis of the model results show that surface heat fluxes are primarily responsible for widespread reduction (0.5°–1.5°C) of sea surface temperature over the inner West Florida Shelf 100–300 km away from the storm center. Momentum fluxes are responsible for stronger surface cooling (2°C) near the center of the storm. The upper ocean heat loss near the storm center of more than 200 MJ/m² is primarily due to the vertical flux of thermal energy between the surface layer and deep ocean. Heat loss to the atmosphere during the storm’s passage is approximately 100–150 MJ/m². The upper ocean cooling is enhanced where the preexisting mixed layer is shallow, e.g., within a cyclonic circulation feature, although the heat flux to the atmosphere in these locations is markedly reduced.     

Enter the working forest: Discourse analysis in the Northern Forest

In 1988, four states in the northeastern USA commissioned a study to address land use changes in the Northern Forest, 26 million acres of temperate and boreal forest extending from Maine to eastern New York State. Against a backdrop of economic destabilization and concerns regarding social and ecological implications of a real estate boom, the sustained deliberative dialogue catalyzed by this study has come to rely heavily on the ambiguous concept of “working forest.” To clarify political and environmental dynamics in the region, we analyzed how people respond to and seek to capitalize on the interpretive flexibility of the term working forest. We combine an analysis of socio-political discourses of working forest based on a structured literature review with an assessment of local peoples’ definitions of working forest based on a survey conducted in a pair of contrasting New York State communities. The first study site represents an amenity-oriented community (i.e. a place where the forest supports a service economy including recreation and tourism) and the other study site represents a timber-dependent community. By linking data from community-level analysis to data derived from a general analysis of forest politics, we seek to develop a more robust perspective. By comparing discourses across differently structured communities, we investigate how local forest politics are mediated by local economic development processes. Our study empirically illustrates contested and geographically uneven processes of social construction of environment and rural development in a region confronting pressures of globalization. Results indicate that timber harvesting is a heavily privileged management objective, as a logic of ‘the forest that pays is the forest that stays’ dominates. Environmental politics in the region, and perhaps more generally, increasingly conforms to a form of pragmatism in which economic opportunities structure conservation planning and investment.