MARGINAL AGRICULTURAL LAND CLASSIFICATION: A NEW APPROACH

Land Economics/Use,

Simulating Star Formation and Feedback in Galactic Disk Models

We use a high-resolution grid-based hydrodynamics method to simulate the multi-phase interstellar medium in a Milky Way-size quiescent disk galaxy. The models are global and three-dimensional, and include a treatment of star formation and feedback. We examine the formation of gravitational instabilities and show that a form of the Toomre instability criterion can successfully predict where star formation will occur. Two common prescriptions for star formation are investigated. The first is based on cosmological simulations and has a relatively low threshold for star formation, but also enforces a comparatively low efficiency. The second only permits star formation above a number density of 1000 cm^-3 but adopts a high efficiency. We show that both methods can reproduce the observed slope of the relationship between star formation and gas surface density (although at too high a rate for our adopted parameters). A run which includes feedback from type II supernovae is successful at driving gas out of the plane, most of which falls back onto the disk. This feedback also substantially reduces the star formation rate. Finally, we examine the density and pressure distribution of the ISM, and show that there is a rough pressure equilibrium in the disk, but with a wide range of pressures at a given location (and even wider for the case including feedbackComment: 14 pages, 12 figures, accepted to Astrophysical Journa

Design of \u27Iris\u27, a Small Autonomous Surveillance UAV

This paper documents the design process used for a small autonomous surveillance UAV. The most significant requirements for the plane were size (man-packable), endurance (about 1 hour) and cost (essentially disposable). The plane that resulted, named Iris , is a tailless plane with a 45 cm wing span and a total mass of less than 200g. During flight tests, it achieved an endurance of 52 minutes. The estimated cost to manufacture the planes was $343, excluding the autopilot

Gas Stripping in Simulated Galaxies with a Multiphase ISM

Cluster galaxies moving through the intracluster medium (ICM) are expected to lose some of their interstellar medium (ISM) through ISM-ICM interactions. We perform high resolution (40 pc) three-dimensional hydrodynamical simulations of a galaxy undergoing ram pressure stripping including radiative cooling in order to investigate stripping of a multiphase medium. The clumpy, multiphase ISM is self-consistently produced by the inclusion of radiative cooling, and spans six orders of magnitude in gas density. We find no large variations in the amount of gas lost whether or not cooling is involved, although the gas in the multiphase galaxy is stripped more quickly and to a smaller radius. We also see significant differences in the morphology of the stripped disks. This occurs because the multiphase medium naturally includes high density clouds set inside regions of lower density. We find that the lower density gas is stripped quickly from any radius of the galaxy, and the higher density gas can then be ablated. If high density clouds survive, through interaction with the ICM they lose enough angular momentum to drift towards the center of the galaxy where they are no longer stripped. Finally, we find that low ram pressure values compress gas into high density clouds that could lead to enhanced star formation, while high ram pressure leads to a smaller amount of high-density gas.Comment: 17 pages, 12 figures, accepted in Ap

A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data

We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l, followed by a particular deterministic rescaling operation of the sky signal. The acceptance probability for this joint move depends on the sky map only through the difference of chi-squared between the original and proposed sky sample, which is close to unity in the low signal-to-noise regime. The algorithm is completed by alternating this move with a standard Gibbs move. Together, these two proposals constitute a computationally efficient algorithm for mapping out the full joint CMB posterior, both in the high and low signal-to-noise regimes.Comment: Submitted to Ap

A re-analysis of the three-year WMAP temperature power spectrum and likelihood

We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do find two small but potentially important discrepancies: On large angular scales there is a small power excess in the WMAP spectrum (5-10% at l<~30) primarily due to likelihood approximation issues between 13 <= l <~30. On small angular scales there is a systematic difference between the V- and W-band spectra (few percent at l>~300). Recently, the latter discrepancy was explained by Huffenberger et al. (2006) in terms of over-subtraction of unresolved point sources. As far as the low-l bias is concerned, most parameters are affected by a few tenths of a sigma. The most important effect is seen in n_s. For the combination of WMAP, Acbar and BOOMERanG, the significance of n_s =/ 1 drops from ~2.7 sigma to ~2.3 sigma when correcting for this bias. We propose a few simple improvements to the low-l WMAP likelihood code, and introduce two important extensions to the Gibbs sampling method that allows for proper sampling of the low signal-to-noise regime. Finally, we make the products from the Gibbs sampling analysis publically available, thereby providing a fast and simple route to the exact likelihood without the need of expensive matrix inversions.Comment: 14 pages, 7 figures. Accepted for publication in ApJ. Numerical results unchanged, but interpretation sharpened: Likelihood approximation issues at l=13-30 far more important than potential foreground issues at l <= 12. Gibbs products (spectrum and sky samples, and "easy-to-use" likelihood module) available from http://www.astro.uio.no/~hke/ under "Research

CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator

We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky coverage and uniform noise, and an excellent approximation for sky cuts and scanning patterns relevant for modern satellite experiments such as WMAP and Planck. A single evaluation of this estimator between l=2 and 200 takes ~0.2 CPU milliseconds, while for comparison, a single pixel space likelihood evaluation between l=2 and 30 for a map with ~2500 pixels requires ~20 seconds. We apply this tool to the 5-year WMAP temperature data, and re-estimate the angular temperature power spectrum, $C_{\ell}$, and likelihood, L(C_l), for l<=200, and derive new cosmological parameters for the standard six-parameter LambdaCDM model. Our spectrum is in excellent agreement with the official WMAP spectrum, but we find slight differences in the derived cosmological parameters. Most importantly, the spectral index of scalar perturbations is n_s=0.973 +/- 0.014, 1.9 sigma away from unity and 0.6 sigma higher than the official WMAP result, n_s = 0.965 +/- 0.014. This suggests that an exact likelihood treatment is required to higher l's than previously believed, reinforcing and extending our conclusions from the 3-year WMAP analysis. In that case, we found that the sub-optimal likelihood approximation adopted between l=12 and 30 by the WMAP team biased n_s low by 0.4 sigma, while here we find that the same approximation between l=30 and 200 introduces a bias of 0.6 sigma in n_s.Comment: 10 pages, 7 figures, submitted to Ap

Small-angle X-ray Scattering Studies of the Oligomeric State and Quaternary Structure of the Trifunctional Proline Utilization A (PutA) Flavoprotein from \u3ci\u3eEscherichia coli\u3c/i\u3e

Background: Trifunctional proline utilization A (PutA) proteins are multifunctional flavoproteins that catalyze two reactions and repress transcription of the put regulon. Results: PutA from Escherichia coli is a V-shaped dimer, with the DNA-binding domain mediating dimerization. Conclusion: Oligomeric state and quaternary structures are not conserved by PutAs. Significance: The first three-dimensional structural information for any trifunctional PutA is reported