We have statistically analyzed a set of 115 low frequency (Deca- Hectometer wavelengths range) type II and type III bursts associated with major Solar Energetic Particle (SEP: Ep > 10 MeV) events and their solar causes such as solar flares and coronal mass ejections (CMEs) observed from 1997 to 2014. We classified them into two sets of events based on the duration of the associated solar flares:75 impulsive flares (duration < 60 min) and 40 gradual flares (duration > 60 min).The impulsive flare-associated SEP events (Rt = 989.23 min: 2.86 days) are short lived and they quickly reach their peak intensity (shorter rise time) when compared with gradual flares associated events (Rt =1275.45 min: 3.34 days). We found a good correlation between the logarithmic peak intensity of all SEPs and properties of CMEs (space speed: cc = 0.52, SEcc = 0.083), and solar flares (log integrated flux: cc = 0.44, SEcc = 0.083). This particular result gives no clear cut distinction between flare-related and CME-related SEP events for this set of major SEP events. We derived the peak intensity, integrated intensity, duration and slope of these bursts from the radio dynamic spectra observed by Wind/WAVES. Most of the properties (peak intensity, integrated intensity and starting frequency) of DH type II bursts associated with impulsive and gradual flare events are found to be similar in magnitudes. In addition, we also found a significant correlation between the properties of SEPs and key parameters of DH type III bursts. This result shows a closer association of peak intensity of the SEPs with the properties of DH type III radio bursts than with the properties DH type II radio bursts, at least for this set of 115 major SEP events.

It is generally accepted that extreme space weather events tend to be related to strong flares and fast halo coronal mass ejections CMEs. In the present paper, we carefully identify the chain of events from the Sun to the Earth induced by all 12 X-class flares that occurred in 2002. In this small sample, we find an unusual high rate (58\%) of solar sources with a longitude larger than 74 degrees. Yet, all 12 X-class flares are associated with at least one CME. The fast halo CMEs (50\% ) are related to interplanetary CMEs (ICMEs) at L1 and weak Dst minimum values ($> -51\;$nT); while 5 (41\%) of the 12 X-class flares are related to solar proton events (SPE). We conclude that: (i) All twelve analyzed solar events, even those associated with fast halo CMEs originating from the central disk region, and those ICMEs and SPEs were not very geo-effective. This unexpected result demonstrates that the suggested events in the chain (fast halo CME, X-class flares, central disk region, ICME, SPE) are not infallible proxies for geo-effectiveness. (ii) The low value of integrated and normalized southward component of the IMF ($B^*_z$) may explain the low geo-effectiveness for this small sample. In fact, $B^*_z$ is well correlated to the weak Dst and low auroral electrojet (AE) activity. Hence, the only space weather impact at Earth in 2002 we can explain is based on $B^*_z$ at L1.

The Alfv\'{e}n speed and plasma beta in photospheric bright points observed by the Broadband Filter Imager (BFI) of the Solar Optical Telescope (SOT) on board the \textit{Hinode} satellite, are estimated seismologically. The diagnostics is based on the theory of slow magnetoacoustic waves in a non-isothermally stratified photosphere with a uniform vertical magnetic field. We identify and track bright points in a G-band movie by using the 3D region growing method, and align them with blue continuum images to derive their brightness temperatures. From the Fourier power spectra of 118 continuum light curves made in the bright points, we find that light curves of 91 bright points have oscillations with properties which are significantly different from oscillation in quiet regions, with the periods ranging 2.2--16.2~min. We find that the model gives a moderate value of the plasma beta when $\gamma$ lies at around 5/3. The calculated Alfv\'{e}n speed is 9.68$\pm$2.02~km~s$^{-1}$, ranging in 6.3--17.4~km~s$^{-1}$. The plasma beta is estimated to be of 0.93$\pm$0.36, ranging in 0.2--1.9.

Impulsive solar energetic particle (SEP) events originate from the energy dissipation process in small solar flares. Anomalous abundances in impulsive SEP events provide an evidence on unique, yet unclear, acceleration mechanism. The pattern of heavy-ion enhancements indicates that the temperature of the source plasma that is accelerated is low and not flare-like. We examine the solar source of the $^3$He-rich SEP event of 2012 November 20 using Solar Dynamics Observatory (SDO)/ Atmospheric Imaging Assembly (AIA) images and investigate its thermal variation. The examined event is associated with recurrent coronal jets. The Differential Emission Measure (DEM) analysis is applied to study the temperature evolution/distribution of the source regions. Preliminary results show that the temperature of the associated solar source is ranged between 1.2-3.1 MK.

It was recently discovered that the amplitudes of pulsating red giants and supergiants vary significantly on time scales of 20-30 pulsation periods. Here, we analyze the amplitude variability in 29 pulsating yellow supergiants (5 RVa, 4 RVb, 9 SRd, 7 long-period Cepheid, and 4 yellow hypergiant stars), using visual observations from the AAVSO International Database, and Fourier and wavelet analysis using the AAVSO's VSTAR package. We find that these stars vary in amplitude by factors of up to 10 or more (but more typically 3-5), on a mean time scale (L) of 33 +/- 4 pulsation periods (P). Each of the five sub-types shows this same behavior, which is very similar to that of the pulsating red giants, for which the median L/P was 31. For the RVb stars, the lengths of the cycles of amplitude variability are the same as the long secondary periods, to within the uncertainty of each.

To measure the magnetic field strength in the solar corona, we examined 10 fast (> 1000 km/s) limb CMEs which show clear shock structures in SOHO/LASCO images. By applying piston-shock relationship to the observed CME's standoff distance and electron density compression ratio, we estimated the Mach number, Alfven speed, and magnetic field strength in the height range 3 to 15 solar radii (Rs). Main results from this study are: (1) the standoff distance observed in solar corona is consistent with those from a magnetohydrodynamic (MHD) model and near-Earth observations; (2) the Mach number as a shock strength is in the range 1.49 to 3.43 from the standoff distance ratio, but when we use the density compression ratio, the Mach number is in the range 1.47 to 1.90, implying that the measured density compression ratio is likely to be underestimated due to observational limits; (3) the Alfven speed ranges from 259 to 982 km/s and the magnetic field strength is in the range 6 to 105 mG when the standoff distance is used; (4) if we multiply the density compression ratio by a factor of 2, the Alfven speeds and the magnetic field strengths are consistent in both methods; (5) the magnetic field strengths derived from the shock parameters are similar to those of empirical models and previous estimates.

We present here results on the composite luminosity functions of galaxies in the clusters of galaxies selected from the SDSS Cut and Enhance cluster catalog (CE; Goto et al. 2001). We construct the composite luminosity function in the five SDSS bands, u,g,r,i and z, using 204 CE clusters ranging from z=0.02 to z=0.25. We use photometric redshifts to construct composite luminosity functions. The robustness of the weighting scheme was tested using Monte Carlo simulation. We find the slope of composite LFs become flatter toward redder color band. Comparing with field LFs of SDSS (Blanton et al. 2001), cluster LFs have brighter characteristic magnitude and flatter slopes in g, r, i and z bands. These results are consistent with the hypothesis that the cluster LF has two distinct underlying populations i.e. the bright end of the LF is dominated by bright early types that follow a gaussian--like luminosity distribution, while the faint--end of the cluster LF is a steep power--law like function dominated by star--forming (bluer) galaxies. We also study the composite luminosity functions for early--type and late--type galaxies using profile fits, a concentration parameter and u-r color to classify galaxy morphology. The strong dependence of LF on galaxy morphology is found.The faint end slope of the LF is always flatter for early--type galaxies than late--type regardless of passband and methodology. These results are consistent with the hypothesis that the cluster region are dominated by bright elliptical galaxies.

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

We present a catalog of 799 clusters of galaxies in the redshift range z_est = 0.05 - 0.3 selected from ~400 deg^2 of early SDSS commissioning data along the celestial equator. The catalog is based on merging two independent selection methods -- a color-magnitude red-sequence maxBCG technique (B), and a Hybrid Matched-Filter method (H). The BH catalog includes clusters with richness \Lambda >= 40 (Matched-Filter) and N_gal >= 13 (maxBCG), corresponding to typical velocity dispersion of \sigma_v >~ 400 km s^{-1} and mass (within 0.6 h^{-1) Mpc radius) >~ 5*10^{13} h^{-1} M_sun. This threshold is below Abell richness class 0 clusters. The average space density of these clusters is 2*10^{-5} h^3 Mpc^{-3}. All NORAS X-ray clusters and 53 of the 58 Abell clusters in the survey region are detected in the catalog; the 5 additional Abell clusters are detected below the BH catalog cuts. The cluster richness function is determined and found to exhibit a steeply decreasing cluster abundance with increasing richness. We derive observational scaling relations between cluster richness and observed cluster luminosity and cluster velocity dispersion; these scaling relations provide important physical calibrations for the clusters. The catalog can be used for studies of individual clusters, for comparisons with other sources such as X-ray clusters and AGNs, and, with proper correction for the relevant selection functions, also for statistical analyses of clusters.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

We describe an automated method, the Cut & Enhance method (CE) for detecting clusters of galaxies in multi-color optical imaging surveys. This method uses simple color cuts, combined with a density enhancement algorithm, to up-weight pairs of galaxies that are close in both angular separation and color. The method is semi-parametric since it uses minimal assumptions about cluster properties in order to minimize possible biases. No assumptions are made about the shape of clusters, their radial profile or their luminosity function. The method is successful in finding systems ranging from poor to rich clusters of galaxies, of both regular and irregular shape. We determine the selection function of the CE method via extensive Monte Carlo simulations which use both the real, observed background of galaxies and a randomized background of galaxies. We use position shuffled and color shuffled data to perform the false positive test. We have also visually checked all the clusters detected by the CE method. We apply the CE method to the 350 deg^2 of the SDSS (Sloan Digital Sky Survey) commissioning data and construct a SDSS CE galaxy cluster catalog with an estimated redshift and richness for each cluster. The CE method is compared with other cluster selection methods used on SDSS data such as the Matched Filter (Postman et al. 1996, Kim et al. 2001), maxBCG technique (Annis et al. 2001) and Voronoi Tessellation (Kim et al. 2001). The CE method can be adopted for cluster selection in any multi-color imaging surveys.

We identify close pairs of galaxies from 278 deg^2 of Sloan Digital Sky Survey commissioning imaging data. The pairs are drawn from a sample of 330,041 galaxies with 18 < r^* < 20. We determine the angular correlation function of galaxy pairs, and find it to be stronger than the correlation function of single galaxies by a factor of 2.9 +/- 0.4. The two correlation functions have the same logarithmic slope of 0.77. We invert Limber's equation to estimate the three-dimensional correlation functions; we find clustering lengths of r_0= 4.2 +/- 0.4 h^{-1} Mpc for galaxies and 7.8 +/- 0.7 h^{-1} Mpc for galaxy pairs. These results agree well with the global richness dependence of the correlation functions of galaxy systems.

We positionally match sources observed by the Sloan Digital Sky Survey (SDSS), the Two Micron All Sky Survey (2MASS), and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. Practically all 2MASS sources are matched to an SDSS source within 2 arcsec; ~11% of them are optically resolved galaxies and the rest are dominated by stars. About 1/3 of FIRST sources are matched to an SDSS source within 2 arcsec; ~80% of these are galaxies and the rest are dominated by quasars. Based on these results, we project that by the completion of these surveys the matched samples will include about 10^7 stars and 10^6 galaxies observed by both SDSS and 2MASS, and about 250,000 galaxies and 50,000 quasars observed by both SDSS and FIRST. Here we present a preliminary analysis of the optical, infrared and radio properties for the extragalactic sources from the matched samples. In particular, we find that the fraction of quasars with stellar colors missed by the SDSS spectroscopic survey is probably not larger than ~10%, and that the optical colors of radio-loud quasars are ~0.05 mag. redder (with 4-sigma significance) than the colors of radio-quiet quasars.

One of the most vital observational clues for unraveling the origin of Brightest Cluster Galaxies (BCG) is the observed alignment of the BCGs with their host cluster and its surroundings. We have examined the BCG-cluster alignment effect, using clusters of galaxies detected from the Sloan Digital Sky Survey (SDSS). We find that the BCGs are preferentially aligned with the principal axis of their hosts, to a much higher redshift (z >~ 0.3) than probed by previous studies (z <~ 0.1). The alignment effect strongly depends on the magnitude difference of the BCG and the second and third brightest cluster members: we find a strong alignment effect for the dominant BCGs, while less dominant BCGs do not show any departure from random alignment with respect to the cluster. We therefore claim that the alignment process originates from the same process that makes the BCG grow dominant, be it direct mergers in the early stage of cluster formation, or a later process that resembles the galactic cannibalism scenario. We do not find strong evidence for (or against) redshift evolution between 0<z<0.45, largely due to the insufficient sample size (< 200 clusters). However, we have developed a framework by which we can examine many more clusters in an automated fashion for the upcoming SDSS cluster catalogs, which will provide us with better statistics for systematic investigations of the alignment with redshift, richness and morphology of both the cluster and the BCG.

We present a comparison of three cluster finding algorithms from imaging data using Monte Carlo simulations of clusters embedded in a 25 deg^2 region of Sloan Digital Sky Survey (SDSS) imaging data: the Matched Filter (MF; Postman et al. 1996), the Adaptive Matched Filter (AMF; Kepner et al. 1999) and a color-magnitude filtered Voronoi Tessellation Technique (VTT). Among the two matched filters, we find that the MF is more efficient in detecting faint clusters, whereas the AMF evaluates the redshifts and richnesses more accurately, therefore suggesting a hybrid method (HMF) that combines the two. The HMF outperforms the VTT when using a background that is uniform, but it is more sensitive to the presence of a non-uniform galaxy background than is the VTT; this is due to the assumption of a uniform background in the HMF model. We thus find that for the detection thresholds we determine to be appropriate for the SDSS data, the performance of both algorithms are similar; we present the selection function for each method evaluated with these thresholds as a function of redshift and richness. For simulated clusters generated with a Schechter luminosity function (M_r^* = -21.5 and alpha = -1.1) both algorithms are complete for Abell richness >= 1 clusters up to z ~ 0.4 for a sample magnitude limited to r = 21. While the cluster parameter evaluation shows a mild correlation with the local background density, the detection efficiency is not significantly affected by the background fluctuations, unlike previous shallower surveys.

We present 13 Broad Absorption Line (BAL) quasars, including 12 new objects, which were identified in the Sloan Digital Sky Survey (SDSS) and matched within 2'' to sources in the FIRST radio survey catalog. The surface density of this sample of radio-detected BAL quasars is 4.5 +- 1.2 per 100 deg^2, i.e. approximately 4 times larger than previously found by the shallower FIRST Bright Quasar Survey (FBQS). A majority of these radio-detected BAL quasars are moderately radio-loud objects. The fraction of BAL quasars in the entire radio quasar sample, 4.8 +- 1.3 %, is comparable to the fraction of BAL quasars among the SDSS optical quasar sample (ignoring selection effects). We estimate that the true fraction of BAL quasars (mostly HiBALs) in the radio sample is 9.2 +- 2.6 % once selection effects are accounted for. We caution that the absorption troughs of 4 of the 13 radio-detected quasars considered do not strictly satisfy the standard BALnicity criterion. One or possibly two of the new radio-detected BAL quasars are of the rare ``iron LoBAL'' type. BAL quasars are generally redder than the median SDSS quasar at the same redshift.

We present bright galaxy number counts in five broad bands ($u', g', r', i', z'$) from imaging data taken during the commissioning phase of the Sloan Digital Sky Survey (SDSS). The counts are derived from two independent stripes of imaging scans along the Celestial Equator, one each toward the North and the South Galactic cap, covering about 230 and 210 square degrees, respectively. A careful study is made to verify the reliability of the photometric catalog. For galaxies brighter than $r^* = 16$, the catalog produced by automated software is examined against eye inspection of all objects. Statistically meaningful results on the galaxy counts are obtained in the magnitude range $12 \le r^* \le 21$, using a sample of 900,000 galaxies. The counts from the two stripes differ by about 30% at magnitudes brighter than $r^*= 15.5$, consistent with a local $2\sigma$ fluctuation due to large scale structure in the galaxy distribution. The shape of the number counts-magnitude relation brighter than $r^* = 16$ is well characterized by $N \propto 10^{0.6m}$, the relation expected for a homogeneous galaxy distribution in a ``Euclidean'' universe. In the magnitude range $ 16 < r^* < 21$, the galaxy counts from both stripes agree very well, and follow the prediction of the no-evolution model, although the data do not exclude a small amount of evolution. We use empirically determined color transformations to derive the galaxy number counts in the $B$ and $I_{814}$ bands. We compute the luminosity density of the universe at zero redshift in the five SDSS bands and in the $B$ band. We find ${\cal L}_{B} = 2.4 \pm 0.4 \times 10^8L_\odot h $Mpc$^{-3}$, for a reasonably wide range of parameters of the Schechter luminosity function in the $B$ band.

We report on optical spectroscopic observations of the X-ray Transient XTE J1118+480 covering the period from April 7, 2000 to July 4, 2000. The spectrum is characterized by weak, broad, double-peaked Balmer and He lines on top of a blue continuum of slope p ~ 1/3, as expected for an optically-thick accretion disk. The weak Bowen blend seen in our spectra may indicate a low intrinsic metallicity for the source. The presence of a partial S-wave pattern in the HeII lambda-4686 line appears consistent with the reported photometric orbital period Porb = 4.1 hr for XTE J1118+480. By using a combination of Doppler mapping and various theoretical arguments, we constrain plausible orbital parameters for the system: a mass ratio 0.02 <~ q <~ 0.1, an inclination i >~ 70 deg for a neutron star primary, or 30 deg <~ i <~ 50 deg for a black hole primary with a mass between 4 and 10 Msun. CaII absorption features observed at very high resolution constrain the interstellar hydrogen absorption column log[N_HI (cm^{-2})] ~ 20.45 +/- 0.2 and the identification of three absorbing clouds indicate a distance to the source <~ 1 kpc, assuming the line-of-sight to XTE J1118+480 has average high-latitude properties. These results are discussed in the context of previous multiwavelength observations of this unusual system.

We outline here the next generation of cluster-finding algorithms. We show how advances in Computer Science and Statistics have helped develop robust, fast algorithms for finding clusters of galaxies in large multi-dimensional astronomical databases like the Sloan Digital Sky Survey (SDSS). Specifically, this paper presents four new advances: (1) A new semi-parametric algorithm - nicknamed ``C4'' - for jointly finding clusters of galaxies in the SDSS and ROSAT All-Sky Survey databases; (2) The introduction of the False Discovery Rate into Astronomy; (3) The role of kernel shape in optimizing cluster detection; (4) A new determination of the X-ray Cluster Luminosity Function which has bearing on the existence of a ``deficit'' of high redshift, high luminosity clusters. This research is part of our ``Computational AstroStatistics'' collaboration (see Nichol et al. 2000) and the algorithms and techniques discussed herein will form part of the ``Virtual Observatory'' analysis toolkit.

Clusters of galaxies are the most massive objects in the Universe and mapping their location is an important astronomical problem. This paper describes an algorithm (based on statistical signal processing methods), a software architecture (based on a hybrid layered approach) and a parallelization scheme (based on a client/server model) for finding clusters of galaxies in large astronomical databases. The Adaptive Matched Filter (AMF) algorithm presented here identifies clusters by finding the peaks in a cluster likelihood map generated by convolving a galaxy survey with a filter based on a cluster model and a background model. The method has proved successful in identifying clusters in real and simulated data. The implementation is flexible and readily executed in parallel on a network of workstations.

The Sloan Digital Sky Survey has obtained 450 square degrees of photometric scan data, in five bands (u',g',r',i',z'), which we use to identify clusters of galaxies. We illustrate how we do star-galaxy separation, and present a simple and elegant method of detecting overdensities in the galaxy distribution, using the Voronoi Tessellation.

To probe the weakly non-linear regime, past the point where simple linear theory is sufficient to describe the statistics of the density distribution, we measure the skewness (S_3) and kurtosis (S_4) of the Count Probability Distribution Function (CPDF) of the IRAS 1.2 Jy sample obtained from counts in cells. These quantities are free parameters in a maximum likelihood fit of an Edgeworth expansion convolved with a Poissonian to the observed CPDF. This method, applicable on scales greater than 5 Mpc, is appreciably less sensitive to the tail of the distribution than are measurements of S_3 and S_4 from moments of the CPDF. We measure S_3 and S_4 to l~50 h^{-1} Mpc; the data are consistent with scale invariance, yielding averages of <S_3> = 2.83 +/- 0.09, and <S_4> = 6.89 +/- 0.68. These values are higher than those found using the moments method on the same data set, <S_3> = 1.5 +/- 0.5 and <S_4> = 4.4 +/- 3.7, due to lack of correction in the latter work for finite-volume effects. Unlike the moments method, our results are quite robust to the fact that IRAS galaxies are under-represented in cluster cores. We use N-body simulations to show that our method yields unbiased results.