Search

A fundamental challenge for wide-field imaging surveys is obtaining follow-up spectroscopic observations: there are > $10^9$ photometrically cataloged sources, yet modern spectroscopic surveys are limited to ~few x $10^6$ targets. As we approach the Large Synoptic Survey Telescope (LSST) era, new algorithmic solutions are required to cope with the data deluge. Here we report the development of a machine-learning framework capable of inferring fundamental stellar parameters (Teff, log g, and [Fe/H]) using photometric-brightness variations and color alone. A training set is constructed from a systematic spectroscopic survey of variables with Hectospec/MMT. In sum, the training set includes ~9000 spectra, for which stellar parameters are measured using the SEGUE Stellar Parameters Pipeline (SSPP). We employed the random forest algorithm to perform a non-parametric regression that predicts Teff, log g, and [Fe/H] from photometric time-domain observations. Our final, optimized model produces a cross-validated root-mean-square error (RMSE) of 165 K, 0.39 dex, and 0.33 dex for Teff, log g, and [Fe/H], respectively. Examining the subset of sources for which the SSPP measurements are most reliable, the RMSE reduces to 125 K, 0.37 dex, and 0.27 dex, respectively, comparable to what is achievable via low-resolution spectroscopy. For variable stars this represents a ~12-20% improvement in RMSE relative to models trained with single-epoch photometric colors. As an application of our method, we estimate stellar parameters for ~54,000 known variables. We argue that this method may convert photometric time-domain surveys into pseudo-spectrographic engines, enabling the construction of extremely detailed maps of the Milky Way, its structure, and history.

We present near- and mid-infrared photometry and spectroscopy from PAIRITEL, IRTF, and Spitzer of a metallicity-unbiased sample of 117 cool, hydrogen-atmosphere white dwarfs from the Palomar-Green survey and find five with excess radiation in the infrared, translating to a 4.3+2.7-1.2% frequency of debris disks. This is slightly higher than, but consistent with the results of previous surveys. Using an initial-final mass relation, we apply this result to the progenitor stars of our sample and conclude that 1-7Msol stars have at least a 4.3% chance of hosting planets; an indirect probe of the intermediate-mass regime eluding conventional exoplanetary detection methods. Alternatively, we interpret this result as a limit on accretion timescales as a fraction of white dwarf cooling ages; white dwarfs accrete debris from several generations of disks for ~10Myr. The average total mass accreted by these stars ranges from that of 200km asteroids to Ceres-sized objects, indicating that white dwarfs accrete moons and dwarf planets as well as Solar System asteroid analogues., Comment: ApJ, in press

We present the results of a machine-learning (ML) based search for new R Coronae Borealis (RCB) stars and DY Persei-like stars (DYPers) in the Galaxy using cataloged light curves from the All-Sky Automated Survey (ASAS) Catalog of Variable Stars (ACVS). RCB stars - a rare class of hydrogen-deficient carbon-rich supergiants - are of great interest owing to the insights they can provide on the late stages of stellar evolution. DYPers are possibly the low-temperature, low-luminosity analogs to the RCB phenomenon, though additional examples are needed to fully establish this connection. While RCB stars and DYPers are traditionally identified by epochs of extreme dimming that occur without regularity, the ML search framework more fully captures the richness and diversity of their photometric behavior. We demonstrate that our ML method can use newly discovered RCB stars to identify additional candidates within the same data set. Our search yields 15 candidates that we consider likely RCB stars/DYPers: new spectroscopic observations confirm that four of these candidates are RCB stars and four are DYPers. Our discovery of four new DYPers increases the number of known Galactic DYPers from two to six; noteworthy is that one of the new DYPers has a measured parallax and is m ~ 7 mag, making it the brightest known DYPer to date. Future observations of these new DYPers should prove instrumental in establishing the RCB connection. We consider these results, derived from a machine-learned probabilistic classification catalog, as an important proof-of-concept for the efficient discovery of rare sources with time-domain surveys., Comment: 18 pages, 2 new figures, accepted for publication in ApJ

The rate of image acquisition in modern synoptic imaging surveys has already begun to outpace the feasibility of keeping astronomers in the real-time discovery and classification loop. Here we present the inner workings of a framework, based on machine-learning algorithms, that captures expert training and ground-truth knowledge about the variable and transient sky to automate 1) the process of discovery on image differences and, 2) the generation of preliminary science-type classifications of discovered sources. Since follow-up resources for extracting novel science from fast-changing transients are precious, self-calibrating classification probabilities must be couched in terms of efficiencies for discovery and purity of the samples generated. We estimate the purity and efficiency in identifying real sources with a two-epoch image-difference discovery algorithm for the Palomar Transient Factory (PTF) survey. Once given a source discovery, using machine-learned classification trained on PTF data, we distinguish between transients and variable stars with a 3.8% overall error rate (with 1.7% errors for imaging within the Sloan Digital Sky Survey footprint). At >96% classification efficiency, the samples achieve 90% purity. Initial classifications are shown to rely primarily on context-based features, determined from the data itself and external archival databases. In the ~one year since autonomous operations, this discovery and classification framework has led to several significant science results, from outbursting young stars to subluminous Type IIP supernovae to candidate tidal disruption events. We discuss future directions of this approach, including the possible roles of crowdsourcing and the scalability of machine learning to future surveys such a the Large Synoptical Survey Telescope (LSST)., Comment: 43 pages, 13 figures, submitted to Publications of the Astronomical Society of the Pacific (PASP)

In August 2008, the accreting milli-second X-ray pulsar (AMXP), IGR J00291+5934, underwent an outburst lasting ~ 100 days, the first since its discovery in 2004. We present data from the double-peaked outburst from Faulkes Telescope North, the INT, the Keck Telescope, PAIRITEL, the Westerbork Synthesis Radio Telescope and the Swift, XMM-Newton and RXTE X-ray missions. We study the outburst's evolution at various wavelengths. We study the light curve morphology, presenting the first radio-X-ray Spectral Energy Distributions (SEDs) for this source and the most detailed UV-IR SEDs for any outbursting AMXP. We show simple models that attempt to identify the emission mechanisms responsible. We analyse short-timescale optical variability, and compare a medium resolution optical spectrum with those from 2004. The outburst morphology is unusual for an AMXP, comprising two peaks, the second containing a 'plateau' of ~ 10 days at maximum brightness within 30 days of the initial activity. This has implications on duty cycles of short-period X-ray transients. The X-ray spectrum can be fitted by a single, hard power-law. We detect optical variability of ~ 0.05 magnitudes, on timescales of minutes, but find no periodic modulation. In the optical, the SEDs contain a blue component, indicative of an irradiated disc, and a transient near-infrared (NIR) excess. This excess is consistent with a simple model of an optically thick synchrotron jet (as seen in other outbursting AMXPs). The optical spectrum shows a double-peaked H alpha profile, a diagnostic of an accretion disc, but we do not clearly see other lines (e.g. He I, II) reported in 2004. Optical/IR observations of AMXPs are excellent for studying the evolution of both the outer accretion disc and the inner jet, and may eventually provide us with tight constraints to model disc-jet coupling in accreting neutron stars., Comment: 15 pages, 10 figures, accepted for publication in A&A

We present broadband (gamma-ray, X-ray, near-infrared, optical, and radio) observations of the gamma-ray burst (GRB) 090709A and its afterglow in an effort to ascertain the origin of this high-energy transient. Previous analyses suggested that GRB 090709A exhibited quasi-periodic oscillations with a period of 8.06 s, a trait unknown in long-duration GRBs but typical of flares from soft gamma-ray repeaters. When properly accounting for the underlying shape of the power-density spectrum of GRB 090709A, we find no conclusive (> 3 sigma) evidence for the reported periodicity. In conjunction with the location of the transient (far from the Galactic plane and absent any nearby host galaxy in the local universe) and the evidence for extinction in excess of the Galactic value, we consider a magnetar origin relatively unlikely. A long-duration GRB, however, can account for the majority of the observed properties of this source. GRB 090709A is distinguished from other long-duration GRBs primarily by the large amount of obscuration from its host galaxy (A_K,obs >~ 2 mag)., Comment: 12 pages; 6 figures; Submitted to AJ; Comments welcome

Swift captured for the first time a smoothly rising X-ray re-brightening of clear non-flaring origin after the steep decay in a long gamma-ray burst (GRB): GRB081028. A rising phase is likely present in all GRBs but is usually hidden by the prompt tail emission and constitutes the first manifestation of what is later to give rise to the shallow decay phase. Contemporaneous optical observations reveal a rapid evolution of the injection frequency of a fast cooling synchrotron spectrum through the optical band, which disfavours the afterglow onset (start of the forward shock emission along our line of sight when the outflow is decelerated) as the origin of the observed re-brightening. We investigate alternative scenarios and find that the observations are consistent with the predictions for a narrow jet viewed off-axis. The high on-axis energy budget implied by this interpretation suggests different physical origins of the prompt and (late) afterglow emission. Strong spectral softening takes place from the prompt to the steep decay phase: we track the evolution of the spectral peak energy from the gamma-rays to the X-rays and highlight the problems of the high latitude and adiabatic cooling interpretations. Notably, a softening of both the high and low spectral slopes with time is also observed. We discuss the low on-axis radiative efficiency of GRB081028 comparing its properties against a sample of Swift long GRBs with secure E_gamma,iso measurements., Comment: 22 pages, 11 figures; MNRAS accepted

The Palomar Transient Factory (PTF) is a fully-automated, wide-field survey aimed at a systematic exploration of the optical transient sky. The transient survey is performed using a new 8.1 square degree camera installed on the 48-inch Samuel Oschin telescope at Palomar Observatory; colors and light curves for detected transients are obtained with the automated Palomar 60-inch telescope. PTF uses eighty percent of the 1.2-m and fifty percent of the 1.5-m telescope time. With an exposure of 60-s the survey reaches a depth of approximately 21.3 in g' and 20.6 in R (5 sigma, median seeing). Four major experiments are planned for the five-year project: 1) a 5-day cadence supernova search; 2) a rapid transient search with cadences between 90 seconds and 1 day; 3) a search for eclipsing binaries and transiting planets in Orion; and 4) a 3-pi sr deep H-alpha survey. PTF provides automatic, realtime transient classification and follow up, as well as a database including every source detected in each frame. This paper summarizes the PTF project, including several months of on-sky performance tests of the new survey camera, the observing plans and the data reduction strategy. We conclude by detailing the first 51 PTF optical transient detections, found in commissioning data., Comment: 12 pages, 11 figures, 3 tables, submitted to PASP

Long-duration gamma-ray bursts (GRBs) are widely believed to be highly-collimated explosions (opening angle theta ~ 1-10 deg). As a result of this beaming factor, the true energy release from a GRB is usually several orders of magnitude smaller than the observed isotropic value. Measuring this opening angle, typically inferred from an achromatic steepening in the afterglow light curve (a "jet" break), has proven exceedingly difficult in the Swift era. Here we undertake a study of five of the brightest (in terms of the isotropic prompt gamma-ray energy release, E(gamma, iso)) GRBs in the Swift era to search for jet breaks and hence constrain the collimation-corrected energy release. We present multi-wavelength (radio through X-ray) observations of GRBs 050820A, 060418, and 080319B, and construct afterglow models to extract the opening angle and beaming-corrected energy release for all three events. Together with results from previous analyses of GRBs 050904 and 070125, we find evidence for an achromatic jet break in all five events, strongly supporting the canonical picture of GRBs as collimated explosions. The most natural explanation for the lack of observed jet breaks from most Swift GRBs is therefore selection effects. However, the opening angles for the events in our sample are larger than would be expected if all GRBs had a canonical energy release of ~ 10e51 erg. The total energy release we measure for those "hyper-energetic" (E(total) >~ 10e52 erg) events in our sample is large enough to start challenging models with a magnetar as the compact central remnant., Comment: Submitted to ApJ, comments welcome

In this work we present the first results of our imaging campaign at Keck Observatory to identify the host galaxies of "dark" gamma-ray bursts (GRBs), events with no detected optical afterglow or with detected optical flux significantly fainter than expected from the observed X-ray afterglow. We find that out of a uniform sample of 29 Swift bursts rapidly observed by the Palomar 60-inch telescope through March 2008 (14 of which we classify as dark), all events have either a detected optical afterglow, a probable optical host-galaxy detection, or both. Our results constrain the fraction of Swift GRBs coming from very high redshift (z > 7), such as the recent GRB 090423, to between 0.2-7 percent at 80% confidence. In contrast, a significant fraction of the sample requires large extinction columns (host-frame A_V > 1 mag, with several events showing A_V > 2-6 mag), identifying dust extinction as the dominant cause of the dark GRB phenomenon. We infer that a significant fraction of GRBs (and, by association, of high-mass star formation) occurs in highly obscured regions. However, the host galaxies of dark GRBs seem to have normal optical colors, suggesting that the source of obscuring dust is local to the vicinity of the GRB progenitor or highly unevenly distributed within the host galaxy., Comment: 21 pages, 10 figures. Published in AJ

We report on strong H2 and CO absorption from gas within the host galaxy of gamma-ray burst (GRB) 080607. Analysis of our Keck/LRIS afterglow spectrum reveals a very large HI column density (NHI = 10^22.70 cm^-2) and strong metal-line absorption at z_GRB = 3.0363 with a roughly solar metallicity. We detect a series of A-X bandheads from CO and estimate N(CO) = 10^16.5 cm^-2 and T_ex^CO > 100K. We argue that the high excitation temperature results from UV pumping of the CO gas by the GRB afterglow. Similarly, we observe H2 absorption via the Lyman-Werner bands and estimate N(H2) = 10^21.2 cm^-2 with T_ex^H2 = 10--300K. The afterglow photometry suggests an extinction law with R_V=4 and A_V=3.2 mag and requires the presence of a modest 2175A bump. Additionally, modeling of the Swift/XRT X-ray spectrum confirms a large column density with N(H) = 10^22.58 cm^-2. Remarkably, this molecular gas has extinction properties, metallicity, and a CO/H2 ratio comparable to those of translucent molecular clouds of the Milky Way, suggesting that star formation at high z proceeds in similar environments as today. However, the integrated dust-to-metals ratio is sub-Galactic, suggesting the dust is primarily associated with the molecular phase while the atomic gas has a much lower dust-to-gas ratio. Sightlines like GRB 080607 serve as powerful probes of nucleosynthesis and star-forming regions in the young universe and contribute to the population of "dark" GRB afterglows., Comment: 6 pages, 4 figures, 1 table. Accepted to ApJL on December 3, 2008

Aims: With this paper we want to investigate the highly variable afterglow light curve and environment of gamma-ray burst (GRB) 060526 at $z=3.221$. Methods: We present one of the largest photometric datasets ever obtained for a GRB afterglow, consisting of multi-color photometric data from the ultraviolet to the near infrared. The data set contains 412 data points in total to which we add additional data from the literature. Furthermore, we present low-resolution high signal-to-noise spectra of the afterglow. The afterglow light curve is modeled with both an analytical model using broken power law fits and with a broad-band numerical model which includes energy injections. The absorption lines detected in the spectra are used to derive column densities using a multi-ion single-component curve-of-growth analysis from which we derive the metallicity of the host of GRB 060526. Results: The temporal behaviour of the afterglow follows a double broken power law with breaks at $t=0.090\pm0.005$ and $t=2.401\pm0.061$ days. It shows deviations from the smooth set of power laws that can be modeled by additional energy injections from the central engine, although some significant microvariability remains. The broadband spectral-energy distribution of the afterglow shows no significant extinction along the line of sight. The metallicity derived from \ion{S}{II} and \ion{Fe}{II} of [S/H] = --0.57 $\pm$0.25 and [Fe/H] = --1.09$\pm$0.24 is relatively high for a galaxy at that redshift but comparable to the metallicity of other GRB hosts at similar redshifts. At the position of the afterglow, no host is detected to F775W(AB) = 28.5 mag with the HST, implying an absolute magnitude of the host M(1500 \AA{})$>$--18.3 mag which is fainter than most long-duration hosts, although the GRB may be associated with a faint galaxy at a distance of 11 kpc., Comment: 20 pages, 8 figures, A&A, in press, three additional authors, additional minor changes after second referee report

We present extensive early photometric (ultraviolet through near-infrared) and spectroscopic (optical and near-infrared) data on supernova (SN) 2008D as well as X-ray data analysis on the associated Swift/X-ray transient (XRT) 080109. Our data span a time range of 5 hours before the detection of the X-ray transient to 150 days after its detection, and detailed analysis allowed us to derive constraints on the nature of the SN and its progenitor; throughout we draw comparisons with results presented in the literature and find several key aspects that differ. We show that the X-ray spectrum of XRT 080109 can be fit equally well by an absorbed power law or a superposition of about equal parts of both power law and blackbody. Our data first established that SN 2008D is a spectroscopically normal SN Ib (i.e., showing conspicuous He lines), and show that SN 2008D had a relatively long rise time of 18 days and a modest optical peak luminosity. The early-time light curves of the SN are dominated by a cooling stellar envelope (for \Delta t~0.1- 4 day, most pronounced in the blue bands) followed by 56^Ni decay. We construct a reliable measurement of the bolometric output for this stripped-envelope SN, and, combined with estimates of E_K and M_ej from the literature, estimate the stellar radius R_star of its probable Wolf-Rayet progenitor. According to the model of Waxman et al. and of Chevalier & Fransson, we derive R_star^{W07}= 1.2+/-0.7 R_sun and R_star^{CF08}= 12+/-7 R_sun, respectively; the latter being more in line with typical WN stars. Spectra obtained at 3 and 4 months after maximum light show double-peaked oxygen lines that we associate with departures from spherical symmetry, as has been suggested for the inner ejecta of a number of SN Ib cores., Comment: Accepted to ApJ, v3 contains more data than v2 and more references, conclusions not significantly changed, 28 pages in emulateapj, 17 figures

The first gamma-ray burst (GRB) confirmed to be bright enough to be seen with the naked eye, GRB 080319B at redshift z = 0.937, allowed for exquisite follow-up observations across the electromagnetic spectrum. We present our detailed optical and infrared observations of the afterglow, consisting of over 5000 images starting 51 s after the GRB trigger, in concert with our own analysis of the Swift data. The event is extreme not only in observed properties but intrinsically: it was the most luminous event ever recorded at optical and infrared wavelengths and had an exceedingly high isotropic-equivalent energy release in gamma-rays. At early times, the afterglow evolution is broadly consistent with being reverse-shock dominated, but then is subsumed by a forward shock at around 1000 s. The overall spectral energy distribution, spanning from ultraviolet through near-infrared wavelengths, shows no evidence for a significant amount of dust extinction in the host frame. The afterglow evolution, however, is highly chromatic: starting at about 1000 s the index shifts blueward before shifting back to the red at late times. In our deepest late-time observations, we find tentative evidence for an optical jet break and a luminous supernova. Finally, we examine the detectability of such events with current and future facilities and find that such an event could be detected in gamma-rays by BAT out to z = 10.7 (8 sigma), while the nominal EXIST sensitivity would allow detection to z ~ 32. At K band, this source would have been easily detected with meter-class telescopes to z ~ 17., Comment: Accepted to ApJ (18 Sept 2008), 45 pages, 9 figures. Updated following referee comments. Added photometry from early-time (saturated) IR imaging starting at 51 s post trigger. A full machine-readable table of all photometry used in our plots, figures, and analysis is given at http://astro.berkeley.edu/~dperley/grbpub/080319b/photometry.txt

Temporal sampling does more than add another axis to the vector of observables. Instead, under the recognition that how objects change (and move) in time speaks directly to the physics underlying astronomical phenomena, next-generation wide-field synoptic surveys are poised to revolutionize our understanding of just about anything that goes bump in the night (which is just about everything at some level). Still, even the most ambitious surveys will require targeted spectroscopic follow-up to fill in the physical details of newly discovered transients. We are now building a new system intended to ingest and classify transient phenomena in near real-time from high-throughput imaging data streams. Described herein, the Transient Classification Project at Berkeley will be making use of classification techniques operating on ``features'' extracted from time series and contextual (static) information. We also highlight the need for a community adoption of a standard representation of astronomical time series data (i.e., ``VOTimeseries'')., Comment: 4 pages, refereed proceedings of "Hot-wiring the Transient Universe," June 2007 (Tuscon). Editors: A. Allan, R. Seaman, J. S. Bloom

We report on observations of a gamma-ray burst (GRB 061126) with an extremely bright (R ~ 12 mag at peak) early-time optical afterglow. The optical afterglow is already fading as a power law 22 seconds after the trigger, with no detectable prompt contribution in our first exposure, which was coincident with a large prompt-emission gamma-ray pulse. The optical--infrared photometric spectral energy distribution is an excellent fit to a power law, but it exhibits a moderate red-to-blue evolution in the spectral index at about 500 s after the burst. This color change is contemporaneous with a switch from a relatively fast decay to slower decay. The rapidly decaying early afterglow is broadly consistent with synchrotron emission from a reverse shock, but a bright forward-shock component predicted by the intermediate- to late-time X-ray observations under the assumptions of standard afterglow models is not observed. Indeed, despite its remarkable early-time brightness, this burst would qualify as a dark burst at later times on the basis of its nearly flat optical-to-X-ray spectral index. Our photometric spectral energy distribution provides no evidence of host-galaxy extinction, requiring either large quantities of grey dust in the host system (at redshift 1.1588 +/- 0.0006, based upon our late-time Keck spectroscopy) or separate physical origins for the X-ray and optical afterglows., Comment: Revised version submitted to ApJ. Contains significantly expanded discussion, an additional figure, and numerous other changes

We present lightcurves of the afterglow of GRB050502A, including very early data at t-t_{GRB} < 60s. The lightcurve is composed of unfiltered ROTSE-IIIb optical observations from 44s to 6h post-burst, R-band MDM observations from 1.6 to 8.4h post-burst, and PAIRITEL J H K_s observations from 0.6 to 2.6h post-burst. The optical lightcurve is fit by a broken power law, where t^{alpha} steepens from alpha = -1.13 +- 0.02 to alpha = -1.44 +- 0.02 at \~5700s. This steepening is consistent with the evolution expected for the passage of the cooling frequency nu_c through the optical band. Even in our earliest observation at 44s post-burst, there is no evidence that the optical flux is brighter than a backward extrapolation of the later power law would suggest. The observed decay indices and spectral index are consistent with either an ISM or a Wind fireball model, but slightly favor the ISM interpretation. The expected spectral index in the ISM interpretation is consistent within 1 sigma with the observed spectral index beta = -0.8 +- 0.1; the Wind interpretation would imply a slightly (~2 sigma) shallower spectral index than observed. A small amount of dust extinction at the source redshift could steepen an intrinsic spectrum sufficiently to account for the observed value of beta. In this picture, the early optical decay, with the peak at or below 4.7e14 Hz at 44s, requires very small electron and magnetic energy partitions from the fireball., Comment: 22 pages, including 3 tables and 1 figure, Accepted by ApJ

The explosion that results in a cosmic gamma-ray burst (GRB) is thought to produce emission from two physical processes -- the activity of the central engine gives rise to the high-energy emission of the burst through internal shocking and the subsequent interaction of the flow with the external environment produces long-wavelength afterglow. While afterglow observations continue to refine our understanding of GRB progenitors and relativistic shocks, gamma-ray observations alone have not yielded a clear picture of the origin of the prompt emission nor details of the central engine. Only one concurrent visible-light transient has been found and was associated with emission from an external shock. Here we report the discovery of infrared (IR) emission contemporaneous with a GRB, beginning 7.2 minutes after the onset of GRB 041219a. Our robotic telescope acquired 21 images during the active phase of the burst, yielding the earliest multi-colour observations of any long-wavelength emission associated with a GRB. Analysis of an initial IR pulse suggests an origin consistent with internal shocks. This opens a new possibility to study the central engine of GRBs with ground-based observations at long wavelengths., Comment: Accepted to Nature on March 1, 2005. 9 pages, 4 figures, nature12.cls and nature1.cls files included. This paper is under press embargo until print publication

The internet has brought about great change in the astronomical community, but this interconnectivity is just starting to be exploited for use in instrumentation. Utilizing the internet for communicating between distributed astronomical systems is still in its infancy, but it already shows great potential. Here we present an example of a distributed network of telescopes that performs more efficiently in synchronous operation than as individual instruments. RAPid Telescopes for Optical Response (RAPTOR) is a system of telescopes at LANL that has intelligent intercommunication, combined with wide-field optics, temporal monitoring software, and deep-field follow-up capability all working in closed-loop real-time operation. The Telescope ALert Operations Network (TALON) is a network server that allows intercommunication of alert triggers from external and internal resources and controls the distribution of these to each of the telescopes on the network. TALON is designed to grow, allowing any number of telescopes to be linked together and communicate. Coupled with an intelligent alert client at each telescope, it can analyze and respond to each distributed TALON alert based on the telescopes needs and schedule., Comment: Presentation at SPIE 2004, Glasgow, Scotland (UK)

The mining of Virtual Observatories (VOs) is becoming a powerful new method for discovery in astronomy. Here we report on the development of SkyDOT (Sky Database for Objects in the Time domain), a new Virtual Observatory, which is dedicated to the study of sky variability. The site will confederate a number of massive variability surveys and enable exploration of the time domain in astronomy. We discuss the architecture of the database and the functionality of the user interface. An important aspect of SkyDOT is that it is continuously updated in near real time so that users can access new observations in a timely manner. The site will also utilize high level machine learning tools that will allow sophisticated mining of the archive. Another key feature is the real time data stream provided by RAPTOR (RAPid Telescopes for Optical Response), a new sky monitoring experiment under construction at Los Alamos National Laboratory (LANL)., Comment: to appear in SPIE proceedings vol. 4846, 11 pages, 5 figures

The Rapid Telescope for Optical Response (RAPTOR) program consists of a network of robotic telescopes dedicated to the search for fast optical transients. The pilot project is composed of three observatories separated by approximately 38 kilometers located near Los Alamos, New Mexico. Each of these observatories is composed of a telescope, mount, enclosure, and weather station, all operating robotically to perform individual or coordinated transient searches. The telescopes employ rapidly slewing mounts capable of slewing a 250 pound load 180 degrees in under 2 seconds with arcsecond precision. Each telescope consists of wide-field cameras for transient detection and a narrow-field camera with greater resolution and sensitivity. The telescopes work together by employing a closed-loop system for transient detection and follow-up. Using the combined data from simultaneous observations, transient alerts are generated and distributed via the Internet. Each RAPTOR telescope also has the capability of rapidly responding to external transient alerts received over the Internet from a variety of ground-based and satellite sources. Each observatory may be controlled directly, remotely, or robotically while providing state-of-health and observational results to the client and the other RAPTOR observatories. We discuss the design and implementation of the spatially distributed RAPTOR system., Comment: To appear in SPIE proceedings vol. 4845

The Rapid Telescopes for Optical Response (RAPTOR) experiment is a spatially distributed system of autonomous robotic telescopes that is designed to monitor the sky for optical transients. The core of the system is composed of two telescope arrays, separated by 38 kilometers, that stereoscopically view the same 1500 square-degree field with a wide-field imaging array and a central 4 square-degree field with a more sensitive narrow-field "fovea" imager. Coupled to each telescope array is a real-time data analysis pipeline that is designed to identify interesting transients on timescales of seconds and, when a celestial transient is identified, to command the rapidly slewing robotic mounts to point the narrow-field ``fovea'' imagers at the transient. The two narrow-field telescopes then image the transient with higher spatial resolution and at a faster cadence to gather light curve information. Each "fovea" camera also images the transient through a different filter to provide color information. This stereoscopic monitoring array is supplemented by a rapidly slewing telescope with a low resolution spectrograph for follow-up observations of transients and a sky patrol telescope that nightly monitors about 10,000 square-degrees for variations, with timescales of a day or longer, to a depth about 100 times fainter. In addition to searching for fast transients, we will use the data stream from RAPTOR as a real-time sentinel for recognizing important variations in known sources. Altogether, the RAPTOR project aims to construct a new type of system for discovery in optical astronomy--one that explores the time domain by "mining the sky in real time"., Comment: 11 pages, To appear in the Proceedings of the SPIE, Volume 4845

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.

The TC4 mission in Central America during summer 2007 examined convective transport into the tropical Upper Troposphere/Lower Stratosphere (UTLS) and the evolution of cirrus clouds. The tropical tropopause layer (TTL) circulation is dominated by the Asian monsoon anticyclone and westward winds that stretch from the western Pacific into the Atlantic. During TC4, TTL westward flow over Central America was stronger than normal. Incidence of cold clouds over the Central American region was the third lowest out of 34 years sampled. The major factor was an incipient La Nina, specifically anomalously cold temperatures off the Pacific Coast of South America. Weakness in the low level Caribbean jet caused a shift in the coldest clouds from the Caribbean to the Pacific side of Central America. The character of tropopause temperature variability was that of upward propagating waves generated by local and nonlocal convection. These waves produced tropopause temperature variations of 3 K, with peak-to-peak variations of 8 K. At low levels in Central America, flow from the Sahara desert predominated; further south, the air came from the Amazon region. Convectively influenced air in the upper troposphere came from Central America, the northern Amazon region, the Atlantic ITCZ, and the North American monsoon. In the TTL, Asian and African convection affected the observed air masses. North of 10N in the Central American TTL, African and Asian convection may have contributed as much to the air masses as Central and South American convection. South of 8N, Asian and African convection had far less impact.

A model is presented that permits the simulation of the optical properties of cirrus clouds as measured with depolarization Raman lidars. It comprises a one-dimensional cirrus model with explicit microphysics and an optical module that transforms the microphysical model output to cloud and particle optical properties. The optical model takes into account scattering by randomly oriented or horizontally aligned planar and columnar monocrystals and polycrystals. Key cloud properties such as the fraction of plate-like particles and the number of basic crystals per polycrystal are parameterized in terms of the ambient temperature, the nucleation temperature, or the mass of the particles. The optical-microphysical model is used to simulate the lidar measurement of a synoptically forced cirrostratus in a first case study. It turns out that a cirrus cloud consisting of only monocrystals in random orientation is too simple a model scenario to explain the observations. However, good agreement between simulation and observation is reached when the formation of polycrystals or the horizontal alignment of monocrystals is permitted. Moreover, the model results show that plate fraction and morphological complexity are best parameterized in terms of particle mass, or ambient temperature which indicates that the ambient conditions affect cirrus optical properties more than those during particle formation. Furthermore, the modeled profiles of particle shape and size are in excellent agreement with in situ and laboratory studies, i.e., (partly oriented) polycrystalline particles with mainly planar basic crystals in the cloud bottom layer, and monocrystals above, with the fraction of columns increasing and the shape and size of the particles changing from large thin plates and long columns to small, more isometric crystals from cloud center to top. The findings of this case study corroborate the microphysical interpretation of cirrus measurements with lidar as suggested previously.

Formation and evolution of cirrostratus in response to weak, uniform and constant synoptic forcing is simulated using a one-dimensional numerical model with explicit microphysics, in which the particle size distribution in each grid box is fully resolved. A series of tests of the model response to nucleation modes (homogeneous-freezing-only/heterogeneous nucleation) and heterogeneous nucleation parameters are performed. In the case studied here, nucleation is first activated in the prescribed moist layer. A continuous cloud-top nucleation zone with a depth depending on the vertical humidity gradient and one of the nucleation parameters is developed afterward. For the heterogeneous nucleation cases, intermittent nucleation zones in the mid-upper portion of the cloud form where the relative humidity is on the rise, because existent ice crystals do not uptake excess water vapor efficiently, and ice nuclei (IN) are available. Vertical resolution as fine as 1 m is required for realistic simulation of the homogeneous-freezing-only scenario, while the model resolution requirement is more relaxed in the cases where heterogeneous nucleation dominates. Bulk microphysical and optical properties are evaluated and compared. Ice particle number flux divergence, which is due to the vertical gradient of the gravity-induced particle sedimentation, is constantly and rapidly changing the local ice number concentration, even in the nucleation zone. When the depth of the nucleation zone is shallow, particle number concentration decreases rapidly as ice particles grow and sediment away from the nucleation zone. When the depth of the nucleation zone is large, a region of high ice number concentration can be sustained. The depth of nucleation zone is an important parameter to be considered in parametric treatments of ice cloud generation.

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D, see a brief review by Tao 2003). Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research, at NOAA GFDL, at the U. K. Met. Office, at Colorado State University and at NASA Goddard Space Flight Center (Tao 2003). At Goddard, a 3D Goddard Cumulus Ensemble (GCE) model was used to simulate periods during TOGA COARE (December 19-27, 1992), GATE (September 1-7, 1974), SCSMEX (June 2-11, 1998), ARM (June 26-30, 1997) and KWAJEX (August 7-13, August 18-21, and August 29-September 12, 1999) using a 512 by 512 km domain (with 2-km resolution). The results indicate that surface precipitation and latent heating profiles are similar between the 2D and 3D GCE model simulations. However, there are difference in radiation, surface fluxes and precipitation characteristics. The 2D GCE model was used to perform a long-term integration on ARM/GCSS case 4 (22 days at the ARM Southern Great Plains site in March 2000). Preliminary results showed a large temperature bias in the upper troposphere that had not been seen in previous tropical cases. The major objectives of this paper are: (1) to determine the sensitivities to model configuration (i.e., 2D in west-east, south-north or 3D), (2) to identify the differences and similarities in the organization and entrainment rates of convection between 2D- and 3D-simulated ARM cloud systems, and (3) assess the impact of upper tropospheric forcing on tropical and ARM case 4 cases.

As part of the International H2O Project (MOP-2002), three NASNGSFC lidars acquired high-resolution clear air data of wind, CBL evolution, and water vapor mixing ratio profiles during a dryline event that occurred on 22 May 2002. Together with the anciliary suite of instruments deployed in MOP-2002, these lidar data sets offer a unique look into the mesoscale evolution and convective scale dynamics in and around a dryline.