We discuss the panchromatic properties of 99,088 galaxies selected from the SDSS Data Release 1 spectroscopic sample (a flux-limited sample for 1360 deg^2). These galaxies are positionally matched to sources detected by ROSAT, GALEX, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS. We find strong correlations between the detection fraction at other wavelengths and optical properties such as flux, colors, and emission-line strengths. Using GALEX, SDSS, and 2MASS data, we construct the UV-IR broad-band spectral energy distributions for various types of galaxies, and find that they form a nearly one-parameter family. For example, based on SDSS u- and r-band data, supplemented with redshift, the K-band 2MASS magnitudes can be "predicted" with an rms scatter of only 0.2 mag. When a dust content estimate determined from SDSS data by Kauffmann et al. (2003) is also utilized, this scatter decreases to 0.1 mag. We demonstrate that this dust content is indeed higher for galaxies detected by IRAS and that it can be used to "predict" measured IRAS 60 micron flux density within a factor of two using only SDSS data. We also show that the position of a galaxy in the emission-line-based Baldwin-Phillips-Terlevich diagram is correlated with the optical light concentration index and u-r color determined from the SDSS broad-band imaging data, and discuss changes in the morphology of this diagram induced by requiring detections at other wavelengths. We study the IR-radio correlation and find evidence that its slope may be different for AGN and star-forming galaxies and related to the H_alpha/H_beta line strength ratio.

We present a study of variable faint optical sources discovered by comparing the Sloan Digital Sky Survey (SDSS) and the Palomar Observatory Sky Survey (POSS) catalogs. We use SDSS measurements to photometrically recalibrate several publicly available POSS catalogs; a piecewise recalibration in 100 arcmin2 patches generally results in an improvement of photometric accuracy (rms) by nearly a factor of two, compared to the original data. The POSS I magnitudes can be improved to ~0.15 mag accuracy, and POSS II magnitudes to \~0.10 mag accuracy. We use the recalibrated catalogs for the ~2,000 deg2 of sky in the SDSS Data Release 1 to construct a catalog of ~60,000 sources variable on time scales 10-50 years. A series of statistical tests based on the morphology of SDSS color-color diagrams, as well as visual comparison of images and comparison with repeated SDSS observations, demonstrate the robustness of the selection methods. We quantify the distribution of variable sources in the SDSS color-color diagrams, and the variability characteristics of quasars. We detect a turn-over in quasar structure function which suggests that the characteristic time scale for quasar variability is of the order one year. The long-term (>1 year) quasar variability decreases with luminosity and rest-frame wavelength similarly to the short-term (<1 year) behavior. We also demonstrate that candidate RR Lyrae stars trace the same halo structures, such as the Sgr dwarf tidal stream, that were discovered using repeated SDSS observations. We utilize the POSS-SDSS selected candidates to constrain the halo structure in the parts of sky for which repeated SDSS observations do not exist. (abridged)

We positionally match the 6 cm GB6, 20 cm FIRST and NVSS, and 92 cm WENSS radio catalogs and find 16,500 matches in ~3,000 deg2 of sky. Using this unified radio database, we construct radio "color-magnitude-morphology" diagrams and find that they display a clear structure, rather than a random scatter. We propose a simple, yet powerful, method for morphological classification of radio sources based on FIRST and NVSS measurements. For a subset of matched sources, we find optical identifications using the SDSS Data Release 1 catalogs, and separate them into quasars and galaxies. Compact radio sources with flat radio spectra are dominated by quasars, while compact sources with steep spectra, and resolved radio sources, contain substantial numbers of both quasars and galaxies.

Recent claims by Ivezic et al. (2002) that the distribution of the radio-to-optical flux ratio, R, for quasars is bimodal (the so-called quasar radio dichotomy) were questioned on statistical grounds by Cirasuolo et al. (2003). We apply the approach suggested by Cirasuolo et al. to a sample of \~10,000 objects detected by SDSS and FIRST, and find support for the quasar radio dichotomy. The discrepancy between the claims by Cirasuolo et al. and the results presented here is most likely because 1) the ~100 times larger sample based on two homogeneous surveys that is used here allows a direct determination of the R distribution, rather than relying on indirect inferences based on Monte Carlo simulations of several heterogeneous surveys 2) the accurate SDSS colors and redshift information allow robust determination of the K-correction for R, which, if unaccounted for, introduces significant scatter that masks the intrinsic properties of the quasar R distribution.