The All-Sky Automated Survey for Supernovae (ASAS-SN) has recently used over 2000 days of data to identify more than 50,000 variable stars, automatically classify these, determine periods and amplitudes for those that are periodic -- part of a remarkable project to classify 412,000 known variable stars, and determine their basic properties. This information about the newly-discovered variables, along with the photometric data is freely available on-line. In this paper, we analyze ASAS-SN data on two small random samples of pulsating red giants (PRGs) in detail, and compare our results with those found by ASAS-SN. For the majority of a sample of 29 mostly semi-regular (SR) PRGs, the ASAS-SN results are incorrect or incomplete: either the ASAS-SN periods are 2, 3, or 4 times the actual period, or the ASAS-SN period is a "long secondary period" with a shorter pulsation present, or the star is multi-periodic or otherwise complex, or the star's data are contaminated by instrumental effects. For almost all of a sample of 20 of the longest-period Mira stars (period 640 days or more), the ASAS-SN period is actually 2 or more times the actual period. Our results are not surprising, given the very complex behavior of PRGs.

Photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN), along with the AAVSO VSTAR time-series analysis package has been used to study the pulsational properties of 31 R Coronae Borealis (RCB) stars. Periods have been derived for many of the stars, but the variability tends to be small and complex, occasionally multiperiodic, with noticeable variation in amplitude. As with other RCB stars, the periods are a few weeks.

We have measured the times of onset of recent fadings in four R CrB stars -- V854 Cen, RY Sgr, R CrB, and S Aps. These times continue to be locked to the stars' pulsation periods, though with some scatter. In RY Sgr, the onsets of fading tend to occur at or a few days after pulsation maximum. We have studied the pulsation properties of RY Sgr through its recent long maximum using (O-C) analysis and wavelet analysis. The period "wanders" by a few percent. This wandering can be modelled by random cycle-to-cycle period fluctuations, as in some other types of pulsating stars. The pulsation amplitude varies between 0.05 and 0.25 in visual light, non-periodically but on a time scale of about 20 pulsation periods.

We have used visual and Johnson V observations from the American Association of Variable Star Observers (AAVSO) International Database, and the AAVSO VSTAR time-series package, and (O-C) analysis to investigate possible changes in the pulsation period and amplitude of the pulsating red giants U Del, EU Del, X Her, and Y Lyn around the cycle of their long secondary periods (LSPs). We find no such changes in period or amplitude. This suggests (weakly) that the process -- whatever it is -- that causes the LSPs does not significantly change the physical properties of the visible hemisphere of the stars as the LSP cycle proceeds.

We have used many decades of visual observations from the AAVSO International Database, and the AAVSO time-series analysis package VSTAR, to study the long-term changes in period, amplitude, and mean magnitude in about 30 normal pulsating red giants (PRGs) i.e. those without large secular changes in period, as well as a few of the rare PRGs which do have such secular period changes. The periods of the typical PRGs "wander" on time scales of about 40 pulsation periods -- significantly longer than the time scales of amplitude variation which are 20-35 with a mean of 27 pulsation periods. We have also studied the range and time scale of the long-term changes in pulsation amplitude and mean magnitude, as well as period, and looked for correlations between these. Changes in mean magnitude are larger in stars in stars with larger mean amplitude, but correlate negatively with changes in amplitude. There is a weak positive correlation between the long-term period changes and amplitude changes. The causes of these three kinds of long-term changes are still not clear. We note, from the presence of harmonics in the Fourier spectra, that the longest-period PRGs have distinctly non-sinusoidal phase curves. For studying PRGs, we demonstrate the advantage of studying stars with minimal seasonal gaps in the observations, such as those near the celestial poles. We studied Z UMi, misclassified as a possible Mira star, but actually an R CrB star. We determined times of onset of 16 fadings. We were not able to determine a coherent pulsation period for this star at maximum, with an amplitude greater than 0.05. We did, however, find that the times of onset of fadings were "locked" to a 41.98-day period -- a typical pulsation period for an R CrB star.

This paper reports on time-series analysis of 156 pulsating red giants (21 SRa, 52 SRb, 33 SR, 50 Lb) in the AAVSO (American Association of Variable Star Observers) observing program for which there are no more than 150-250 observations in total. Some results were obtained for 68 of these stars: 17 SRa, 14 SRb, 20 SR, and 17 Lb. These results generally include only an average period and amplitude. Many, if not most of the stars are undoubtedly more complex; pulsating red giants are known to have wandering periods, variable amplitudes, and often multiple periods including "long secondary periods" of unknown origin. These results (or lack thereof) raise the question of how the AAVSO should best manage the observation of these and other sparsely-observed pulsating red giants.

In order to extend our previous studies of the unexplained phenomenon of cyclic amplitude variations in pulsating red giants, we have used the AAVSO time-series analysis package VSTAR to analyze long-term AAVSO visual observations of 50 such stars, mostly Mira stars. The relative amount of the variation, typically a factor of 1.5, and the time scale of the variation, typically 20-35 pulsation periods, are not significantly different in longer-period, shorter-period, and carbon stars in our sample, and they also occur in stars whose period is changing secularly, perhaps due to a thermal pulse. The time scale of the variations is similar to that in smaller-amplitude SR variables, but the relative amount of the variation appears to be larger in smaller-amplitude stars, and is therefore more conspicuous. The cause of the amplitude variations remains unknown.

In this review, I reflect on four decades of my experience in linking astronomy research and education by supervising variable-star research projects by undergraduates, and by outstanding senior high school students. I describe the evolution of my experience, the students I have supervised, the nature of their projects, the educational contexts of the projects, the need for "best practices", the journals in which we publish, and the special role of the American Association of Variable Star Observers (AAVSO). I then describe our recent research on pulsating red giants and related objects, including three astrophysical mysteries that we have uncovered. Finally, I suggest how my projects might be scaled up or extended by others who supervise student research.

We have used AAVSO visual and photoelectric V data, and the AAVSO time-series package VSTAR and the Lomb-Scargle time-series algorithm to determine improved pulsation periods, "long secondary periods" (LSPs). and their amplitudes in 51 shorter-period pulsating red giants in the AAVSO photoelectric photometry program, and the AAVSO LPV (long period variable) binocular program. As is well known, radial pulsation becomes detectable in red giants at about spectral type M0, with periods of about 20 days. We find that the LSP phenomenon is also first detectable at about M0. Pulsation and LSP amplitudes increase from near zero to about 0.1 at pulsation periods of about 100 days. At longer periods, the pulsation amplitudes continue to increase, but the LSP amplitudes are generally between 0.1 and 0.2 on average. The ratios of LSP to pulsation period cluster around 5 and 10, presumably depending on whether the pulsation period is the fundamental or first overtone. The pulsation and LSP phase curves are generally close to sinusoidal, except when the amplitude is small, in which case they may be distorted by observational scatter or, in the case of the LSP amplitude, by the pulsational variability. As with longer-period stars, the LSP amplitude increases and decreases by a factor of two or more, for unknown reasons, on a median time scale of about 20 LSPs. The LSP phenomenon is thus present and simiklar in radially pulsating red giants of all periods. Its cause remains unknown.

We review ways in which observations of the changing periods of pulsating variable stars can be used to detect and directly measure their evolution. We briefly describe the two main techniques of analysis -- (O-C) analysis and wavelet analysis -- and results for pulsating variable star types which are reasonably periodic: type I and II Cepheids, RR Lyrae stars, beta Cephei stars, and Mira stars. We comment briefly on delta Scuti stars and pulsating white dwarfs. For some of these variable star types, observations agree approximately with the predictions of evolutionary models, but there still exist significant areas of disagreement that challenge future models of stellar evolution. There may be a need, for instance, to include processes such as rotation, mass loss, and magnetic fields. There may also be non-evolutionary processes which are contributing to the period changes.

We have used systematic, sustained visual observations from the AAVSO International Database, and the AAVSO time-series analysis package VSTAR to study the unexplained "long secondary periods" (LSPs) in 27 pulsating red giants. In our sample, the LSPs range from 479 to 2967 days, and are on average 8.1 +/- 1.3 times the pulsation period. There is no evidence for more than one LSP in each star. In stars with both the fundamental and first overtone radial period present, the LSP is more often about 10 times the latter. The visual amplitudes of the LSPs are typically 0.1 magnitude and do not correlate with the LSP. The phase curves tend to be sinusoidal, but at least two are sawtooth. The LSPs are stable, within their errors, over the timespan of our data, which is typically 25,000 days. The amplitudes, however, vary by up to a factor of two or more on a timescale of roughly 20-30 LSPs. There is no obvious difference between the behavior of the carbon (C) stars and the normal oxygen (M) stars. Previous multicolor photoelectric observations showed that the LSP color variations are similar to those of the pulsation period, and of the LSPs in the Magellanic Clouds, and not like those of eclipsing stars. We note that the LSPs are similar to the estimated rotation periods of the stars, though the latter have large uncertainties. This suggests that the LSP phenomenon may be a form of modulated rotational variability.

This project began as a search for triple-mode pulsating red giants (PRGs) among known double-mode PRGs,, using AAVSO visual data. It then expanded to include a check on the reality and significance of some double-mode PRGs. We identified several "double-mode" PRGs in which one of the periods might simply be an alias of the other. We also comment on the possibility that the additional periods might be harmonics, rather than overtones. As for the triple-mode pulsators: in our sample, we find no convincing cases in which there is a significant, non-alias third period which might be an additional radial period. We note, however, that a few triple-mode pulsators have been found in a very large sample of PRGs in the LMC/ We discuss the results in the context of recent theoretical models of pulsating red giants.

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.

Several classes of pulsating stars are now known to undergo slow changes in amplitude; these include pulsating red giants and supergiants, and yellow supergiants. We have used visual observations from the AAVSO International Database, and wavelet analysis of 39 red giants, 7 red supergiants, and 3 yellow supergiants, to test the hypothesis that an increase in amplitude would result in an increase in period, because of non-linear effects in the pulsation. For most of the stars, the results are complex and/or indeterminate, due to the limitations of the data, the small amplitude or amplitude variation, or other processes such as random cycle-to-cycle period fluctuations. For the dozen stars which have substantial amplitude variation, and reasonably simple behavior, there is a 75-80% tendency to show a positive correlation between amplitude and period.

We have used long-term AAVSO visual observations and Fourier and wavelet analysis to identify periods and study long-term amplitude variations in 44 red supergiants. Of these, 12 stars had data which were too sparse and/or had low amplitude and/or were without conspicuous peaks in the Fourier spectrum; 6 stars had only long (2500-4000 days) periods without significant amplitude variation. The other 26 stars had one or two periods, either "short" (hundreds of days) or "long" (thousands of days), whose amplitudes varied by up to a factor of 8, but more typically 2-4. The median timescale of the amplitude variation was 18 periods. We interpret the shorter periods as due to pulsation, and the longer periods as analogous to the "long secondary periods" found in pulsating red giants.

MV Sgr is a hot, hydrogen-deficient star which has undergone R CrB fadings. We have used self-correlation analysis and Fourier analysis of CCD V photometry in the AAVSO International Database to identify a period of 8.0 days in this star; the amplitude is about 0.03 mag. The variability is most likely due to pulsation.

We describe a study of period changes in 59 RR Lyrae stars, using times of maximum brightness from the GEOS database. The work was carried out by outstanding senior high school students in the University of Toronto Mentorship Program. This paper is written in such a way that high school or undergraduate physics and astronomy students could use it as a guide and template for carrying out original research, by studying period changes in these and other types of variable stars.

A collection of active photometric observations over the last half decade, archival data from the past 120 years, radial velocity observations from 1984, and recent monitoring through a pro-am collaboration reveal that the 9th magnitude F9 Ib supergiant HDE 344787 is a double-mode Cepheid variable of extremely small amplitude. It displays remarkably similar, but much more extreme, properties to the exotic Cepheid Polaris, including a rapidly-increasing period and sinusoidal light variations of decreasing amplitude suggesting that pulsational stability may occur as early as 2045. Unlike Polaris, HDE 344787 displays sinusoidal light variations at periods of both 5.4 and 3.8 days, corresponding to canonical fundamental mode and overtone pulsation. But it may be similar to Polaris in helping to define a small subgroup of Cepheids that display characteristics consistent with a first crossing of the instability strip. An update of 2010 observations of this remarkable star is presented.

Random changes in pulsation period occur in cool pulsating Mira variables, Type A, B, and C semiregular variables, RV Tauri variables, and in most classical Cepheids. The physical processes responsible for such fluctuations are uncertain, but presumably originate in temporal modifications of the envelope convection in such stars. Such fluctuations are seemingly random over a few pulsation cycles of the stars, but are dominated by the regularity of the primary pulsation over the long term. The magnitude of stochasticity in pulsating stars appears to be linked directly to their dimensions, although not in simple fashion. It is relatively larger in M supergiants, for example, than in short-period Cepheids, but is common enough that it can be detected in visual observations for many types of pulsating stars. Although chaos was discovered in such stars 80 years ago, detection of its general presence in the group has only been possible in recent studies.

P Cygni is a prototype of the Luminous Blue Variables (or S Doradus variables), and the star displays photometric and emission line variability on a timescale of years (known as the "short S Doradus phase" variations). Here we present new high resolution H-alpha spectroscopy of P Cyg that we combine with earlier spectra and concurrent V-band photometry to document the emission and continuum flux variations over a 24 y time span. We show that the emission and continuum fluxes vary in concert on timescales of 1.6 y and longer, but differ on shorter timescales. The H-alpha profile shape also varies on the photometric timescales, and we describe the observed co-variations of the emission peak and absorption trough properties. We argue that the episodes of photometric and emission brightening are caused by increases in the size of the emission region that are related to variations in wind mass loss rate and outflow speed. We find evidence of blueward accelerating, Discrete Absorption Components (DACs) in the absorption trough of the H-alpha profile, and these features have slower accelerations and longer durations than those observed in other lines. The DAC strengths also appear to vary on the photometric timescales, and we suggest that the propagation of the DAC-related wind structures is closely related to changes in the overall wind mass loss rate and velocity.

Random changes in pulsation period are well established in cool pulsating stars, in particular the red giant variables: Miras, semi-regulars of types A and B, and RV Tau variables. Such effects are also observed in a handful of Cepheids, the SX Phe variable XX Cyg, and, most recently, the red supergiant variable, BC Cyg, a type C semi-regular. The nature of such fluctuations is seemingly random over a few pulsation cycles of the stars, yet the regularity of the primary pulsation mechanism dominates over the long term. The degree of stochasticity is linked to the dimensions of the stars, the randomness parameter 'e' appearing to correlate closely with mean stellar radius through the period 'P', with an average value of e/P = 0.0136+-0.0005. The physical processes responsible for such fluctuations are uncertain, but presumably originate in temporal modifications of envelope convection in such stars.

We provide a quantitative analysis of time-variable phenomena in the photospheric, near-star, and outflow regions of the late-B supergiant (SG) HD 199478. The analysis is based primarily on optical spectroscopic datasets secured between 1999 and 2000 from the Bulgarian NAO, Tartu, and Ritter Observatories. The temporal behaviour of HD 199478 is characterised by three key empirical properties: (i) systematic central velocity shifts in the photospheric absorption lines, including C II and He I, over a characteristic time-scale of abou 20 days; (ii) extremely strong, variable H alpha emission with no clear modulation signal, and (iii) the occurrence in 2000 of a (rare) high-velocity absorption (HVA) event in H alpha, which evolved over about 60 days, showing the clear signature of mass infall and outflows. In these properties HD 199478 resembles few other late-B SGs with peculiar emission and HVAs in H alpha (HD 91619, HD 34085, HD 96919). Non-LTE line synthesis modelling is conducted using FASTWIND for these late-B SGs to constrain and compare their fundamental parameters within the context of extreme behaviour in the H alpha lines. Our analysis indicate that at the cooler temperature edge of B SGs, there are objects whose wind properties, as traced by H alpha, are inconsistent with the predictions of the smooth, spherically symmetric wind approximation. This discordance is still not fully understood and may highlight the role of a non-spherical, disk-like, geometry, which may result from magnetically-driven equatorial compression of the gas. Ordered dipole magnetic fields may also lead to confined plasma held above the stellar surface, which ultimately gives rise to transient HVA events.

The preliminary results of a 1996 multisite campaign on the pulsating lambda Bootis star 29 Cyg are presented. This campaign, initiated by the Central Asian Network (CAN), allowed to collect data during 48 photometric and 2 spectroscopic nights at observatories of Ukraine, Kazakhstan, Turkey, Hungary, Canada and USA. The multiplicity of excited modes is confirmed in both spectroscopy and photometry. The analysis of all 1995 and 1996 data, combined with new CAN observations made in July-October 1997, is under way.