Debate over the existence of branches in the stellar activity-rotation diagrams continues. Application of modern time series analysis tools to study the mean cycle periods in chromospheric activity index is lacking. We develop such models, based on Gaussian processes, for one-dimensional time series and apply it to the extended Mount Wilson Ca H&K sample. Our main aim is to study how the previously commonly used assumption of strict harmonicity of the stellar cycles as well as handling of the linear trends affects the results. We introduce three methods of different complexity, starting with the simple Bayesian harmonic model and followed by Gaussian Process models with periodic and quasi-periodic covariance functions. We confirm the existence of two populations in the activity-period diagram. We find only one significant trend in the inactive population, namely that the cycle periods get shorter with increasing rotation. This is in contrast with earlier studies, that postulate the existence of trends in both of the populations. In terms of rotation to cycle period ratio, our data is consistent with only two activity branches such that the active branch merges together with the transitional one. The retrieved stellar cycles are uniformly distributed over the R'HK activity index, indicating that the operation of stellar large-scale dynamos carries smoothly over the Vaughan-Preston gap. At around the solar activity index, however, indications of a disruption in the cyclic dynamo action are seen. Our study shows that stellar cycle estimates depend significantly on the model applied. Such model-dependent aspects include the improper treatment of linear trends, while the assumption of strict harmonicity can result in the appearance of double cyclicities that seem more likely to be explained by the quasi-periodicity of the cycles.

Period estimation is one of the central topics in astronomical time series analysis, where data is often unevenly sampled. Especially challenging are studies of stellar magnetic cycles, as there the periods looked for are of the order of the same length than the datasets themselves. The datasets often contain trends, the origin of which is either a real long-term cycle or an instrumental effect, but these effects cannot be reliably separated, while they can lead to erroneous period determinations if not properly handled. In this study we aim at developing a method that can handle the trends properly, and by performing extensive set of testing, we show that this is the optimal procedure when contrasted with methods that do not include the trend directly to the model. The effect of the form of the noise (whether constant or heteroscedastic) on the results is also investigated. We introduce a Bayesian Generalised Lomb-Scargle Periodogram with Trend (BGLST), which is a probabilistic linear regression model using Gaussian priors for the coefficients and uniform prior for the frequency parameter. We show, using synthetic data, that when there is no prior information on whether and to what extent the true model of the data contains a linear trend, the introduced BGLST method is preferable to the methods which either detrend the data or leave the data untrended before fitting the periodic model. Whether to use noise with different than constant variance in the model depends on the density of the data sampling as well as on the true noise type of the process.

This paper contributes to the field of modeling and hindcasting of the total solar irradiance (TSI) based on different proxy data that extend further back in time than the TSI that is measured from satellites. We introduce a simple method to analyze persistent frequency-dependent correlations (FDCs) between the time series and use these correlations to hindcast missing historical TSI values. We try to avoid arbitrary choices of the free parameters of the model by computing them using an optimization procedure. The method can be regarded as a general tool for pairs of data sets, where correlating and anticorrelating components can be separated into non-overlapping regions in frequency domain. Our method is based on low-pass and band-pass filtering with a Gaussian transfer function combined with de-trending and computation of envelope curves. We find a major controversy between the historical proxies and satellite-measured targets: a large variance is detected between the low-frequency parts of targets, while the low-frequency proxy behavior of different measurement series is consistent with high precision. We also show that even though the rotational signal is not strongly manifested in the targets and proxies, it becomes clearly visible in FDC spectrum. The application of the new method to solar data allows us to obtain important insights into the different TSI modeling procedures and their capabilities for hindcasting based on the directly observed time intervals.

Many real world systems exhibit cyclic behavior that is, for example, due to the nearly harmonic oscillations being perturbed by the strong fluctuations present in the regime of significant non-linearities. For the investigation of such sys- tems special techniques relaxing the assumption to periodicity are required. In this paper, we present the generalization of one of such techniques, namely the D2 phase dispersion statistic, to multidimensional datasets, especially suited for the analysis of the outputs from three-dimensional numerical simulations of the full magnetohydrodynamic equations. We present the motivation and need for the usage of such a method with simple test cases, and present an application to a solar-like semi-global numerical dynamo simulation covering nearly 150 magnetic cycles.

We study the alignment of galaxies relative to their local environment in SDSS-DR8 and, using these data, we discuss evolution scenarios for different types of galaxies. We defined a vector field of the direction of anisotropy of the local environment of galaxies. We summed the unit direction vectors of all close neighbours of a given galaxy in a particular way to estimate this field. We found the alignment angles between the spin axes of disc galaxies, or the minor axes of elliptical galaxies, and the direction of anisotropy. The distributions of cosines of these angles are compared to the random distributions to analyse the alignment of galaxies. Sab galaxies show perpendicular alignment relative to the direction of anisotropy in a sparse environment, for single galaxies and galaxies of low luminosity. Most of the parallel alignment of Scd galaxies comes from dense regions, from 2...3 member groups and from galaxies with low luminosity. The perpendicular alignment of S0 galaxies does not depend strongly on environmental density nor luminosity; it is detected for single and 2...3 member group galaxies, and for main galaxies of 4...10 member groups. The perpendicular alignment of elliptical galaxies is clearly detected for single galaxies and for members of < 11 member groups; the alignment increases with environmental density and luminosity. We confirm the existence of fossil tidally induced alignment of Sab galaxies at low z. The alignment of Scd galaxies can be explained via the infall of matter to filaments. S0 galaxies may have encountered relatively massive mergers along the direction of anisotropy. Major mergers along this direction can explain the alignment of elliptical galaxies. Less massive, but repeated mergers are possibly responsible for the formation of elliptical galaxies in sparser areas and for less luminous elliptical galaxies.

Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mechanisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years, during which epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. The DNS employed is a semi-global (wedge) magnetoconvection model. For data analysis we use Ensemble Empirical Mode Decomposition (EEMD) and phase dispersion ($D^2$) methods. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like. This mode is accompanied by a higher frequency mode near the surface and a low-frequency mode in the bottom of the convection zone. The overall behavior of the dynamo solution is very complex exhibiting variable cycle lengths, epochs of disturbed and even ceased surface activity, and strong short-term hemispherical asymmetries. Surprisingly, the most prominent suppressed surface activity epoch is actually a global magnetic energy maximum. We interpret the overall irregular behavior to be due to the interplay of the different dynamo modes showing different equatorial symmetries, especially the smoother part of the irregular variations being related to the variations of the mode strengths, evolving with different and variable cycle lengths. The abrupt low activity epoch in the dominant dynamo mode near the surface is related to a strong maximum of the bottom toroidal field strength, which causes abrupt disturbances especially in the differential rotation profile via the suppression of the Reynolds stresses.

Late-type stars rotate differentially owing to anisotropic turbulence in their outer convection zones. The rotation is called solar-like (SL) when the equator rotates fastest and anti-solar (AS) otherwise. Hydrodynamic simulations show a transition from SL to AS rotation as the influence of rotation on convection is reduced, but the opposite transition occurs at a different point in the parameter space. The system is bistable, i.e., SL and AS rotation profiles can both be stable. We study the effect of a dynamo-generated magnetic field on the large-scale flows, particularly on the possibility of bistable behavior of differential rotation. We solve the hydromagnetic equations numerically in a rotating spherical shell for a set of different radiative conductivities controlling the relative importance of convection. In agreement with earlier findings, our models display SL rotation profiles when the rotational influence on convection is strong and a transition to AS when the rotational influence decreases. We find that dynamo-generated magnetic fields help to produce SL differential rotation compared to the hydrodynamic simulations. We do not observe any bistable states of differential rotation. In the AS cases we get coherent single-cell meridional circulation, whereas in SL cases we get multi-cellular patterns. In both cases, we obtain poleward circulation near the surface with a magnitude close to that observed in the Sun. Moreover, both differential rotation and meridional circulation have significant magnetic cycle-related variations that are similar in strength to those of the Sun. Purely hydrodynamic simulations of differential rotation and meridional circulation are shown to be of limited relevance as magnetic fields, self-consistently generated by dynamo action, significantly affect the flows.

We study LQ Hya photometry for 1982-2014 with the carrier fit (CF) -method and compare our results to earlier photometric analysis and recent Doppler imaging maps. We first utilize different types of statistical methods to estimate various candidates for the carrier period for the CF method. Secondly, a global fit to the whole data set and local fits to shorter segments are computed with the period that is found to be the optimal one. The harmonic least-squares analysis of all the available data reveals a short period close to 1.6 days as a limiting value for a set of significant frequencies. We interpret this as the rotation period of the spots near the equatorial region. In addition, the distribution of the significant periods is found to be bimodal, hinting of a longer-term modulating period, which we set out to study with a two-harmonic CF model. Weak modulation signal is, indeed retrieved, with a period of roughly 6.9 years. The phase dispersion analysis gives a clear symmetric minimum for coherence times lower than and around 100 days. We interpret this as the mean rotation period of the spots (1.60514 days), and this value is chosen to be used as the carrier period for the CF analysis. With the CF method we seek for any systematic trends in the spot distribution in the global time frame, and locally look for abrupt phase changes earlier reported in rapidly rotating objects. During 2005-2008 the global CF reveals a coherent structure rotating with a period of 1.6037 days, while during most other times the spot distribution appears rather random in phase. The evolution of the spot distribution of the object is found to be very chaotic, with no clear signs of an azimuthal dynamo wave that would persist over longer time scales, although the short-lived coherent structures observed occasionally do not rotate with the same speed as the mean spot distribution.

Context. Galaxies in the Universe form chains (filaments) that connect groups and clusters of galaxies. The filamentary network includes nearly half of the galaxies and is visually the most striking feature in cosmological maps. Aims. We study the distribution of galaxies along the filamentary network, trying to find specific patterns and regularities. Methods. Galaxy filaments are defined by the Bisous model, a marked point process with interactions. We use the two-point correlation function and the Rayleigh Z-squared statistic to study how galaxies and galaxy groups are distributed along the filaments. Results. We show that galaxies and groups are not uniformly distributed along filaments, but tend to form a regular pattern. The characteristic length of the pattern is around 7 Mpc/h. A slightly smaller characteristic length 4 Mpc/h can also be found, using the Z-squared statistic. Conclusions. We find that galaxy filaments in the Universe are like pearl necklaces, where the pearls are galaxy groups distributed more or less regularly along the filaments. We propose that this well defined characteristic scale could be used to test various cosmological models and to probe environmental effects on the formation and evolution of galaxies.

According to earlier Doppler images of the magnetically active primary giant component of the RS CVn binary II Peg, the surface of the star was dominated by one single active longitude that was clearly drifting in the rotational frame of the binary system during 1994-2002; later imaging for 2004-2010, however, showed decreased and chaotic spot activity, with no signs of the drift pattern. Here we set out to investigate from a more extensive photometric dataset whether such a drift is a persistent phenomenon, in which case it could be due to either an azimuthal dynamo wave or an indication of the binary system orbital synchronization still being incomplete. We analyse the datasets using the Carrier Fit method (hereafter CF), especially suitable for analyzing time series in which a fast clocking frequency (such as the rotation of the star) is modulated with a slower process (such as the stellar activity cycle). We combine all collected photometric data into one single data set, and analyze it with the CF method. As a result, we confirm the earlier results of the spot activity having been dominated by one primary spotted region almost through the entire data set, and the existence of a persistent, nearly linear drift. Disruptions of the linear trend and complicated phase behavior are also seen, but the period analysis reveals a rather stable periodicity with P(spot)=6.71054d plus/minus 0.00005d. After the linear trend is removed from the data, we identify several abrupt phase jumps, three of which are analyzed in more detail with the CF method. These phase jumps closely resemble what is called flip-flop event, but the new spot configurations do not, in most cases, persist for longer than a few months.

FK Comae is a rapidly rotating magnetically active star, the light curve of which is modulated by cool spots on its surface. It was the first star where the "flip-flop" phenomenon was discovered. Since then, flip-flops in the spot activity have been reported in many other stars. Therefore, it is of interest to perform a more thorough study of the evolution of the spot activity in FK Com. In this study, we analyse 15 years of photometric observations with two different time series analysis methods, with a special emphasis on detecting flip-flop type events from the data. We apply the continuous period search and carrier fit methods on long-term standard Johnson-Cousins V-observations from the years 1995--2010. The observations were carried out with two automated photometric telescopes, Phoenix-10 and Amadeus T7 located in Arizona. We identify complex phase behaviour in 6 of the 15 analysed data segments. We identify five flip-flop events and two cases of phase jumps, where the phase shift is \Delta \phi < 0.4. In addition we see two mergers of spot regions and two cases where the apparent phase shifts are caused by spot regions drifting with respect to each other. Furthermore we detect variations in the rotation period corresponding to a differential rotation coefficient of |k|>0.031. The flip-flop cannot be interpreted as a single phenomenon, where the main activity jumps from one active longitude to another. In some of our cases the phase shifts can be explained by differential rotation: Two spot regions move with different angular velocity and even pass each other. Comparison between the methods show that the carrier fit utility is better in retrieving slow evolution especially from a low amplitude light curve, while the continuous period search is more sensitive in case of rapid changes.

In this paper we first discuss observational evidence of longitudinal concentrations of magnetic activity in the Sun and rapidly rotating late-type stars with outer convective envelopes. Scenarios arising from the idea of rotationally influenced anisotropic convective turbulence being the key physical process generating these structures are then presented and discussed - such effects include the turbulent dynamo mechanism, negative effective magnetic pressure instability (NEMPI) and hydrodynamical vortex instability. Finally, we discuss non-axisymmetric stellar mean-field dynamo models, the results obtained with them, and compare those with the observational information gathered up so far. We also present results from a pure alpha-squared mean-field dynamo model, which show that time-dependent behavior of the dynamo solutions can occur both in the form of an azimuthal dynamo wave and/or oscillatory behavior related to the alternating energy levels of the active longitudes.

Knowing the variability of total solar irradiance (TSI) at the top of the atmosphere is crucial for specifying solar influence to climate variability. Satellite measured TSI data are available only for the last 32 years. But there is an opportunity to estimate approximate daily TSI values on the basis of the observed variability of solar activity. A common approach is based using data for Wolf sunspot numbers. Series of daily data from the DAVOS based TSI starting from 1978 have been compared with daily values of Wolf numbers and with the corresponding estimates for daily summary sunspot area. Our goal is to use established correlations to hindcast earlier TSI values. The available solar activity data enabled us to reproduce rough TSI estimates back to the 1870's.

The light curves of variable stars are commonly described using simple trigonometric models, that make use of the assumption that the model parameters are constant in time. This assumption, however, is often violated, and consequently, time series models with components that vary slowly in time are of great interest. In this paper we introduce a class of data analysis and visualization methods which can be applied in many different contexts of variable star research, for example spotted stars, variables showing the Blazhko effect, and the spin-down of rapid rotators. The methods proposed are of explorative type, and can be of significant aid when performing a more thorough data analysis and interpretation with a more conventional method.Our methods are based on a straightforward decomposition of the input time series into a fast "clocking" periodicity and smooth modulating curves. The fast frequency, referred to as the carrier frequency, can be obtained from earlier observations (for instance in the case of photometric data the period can be obtained from independently measured radial velocities), postulated using some simple physical principles (Keplerian rotation laws in accretion disks), or estimated from the data as a certain mean frequency. The smooth modulating curves are described by trigonometric polynomials or splines. The data approximation procedures are based on standard computational packages implementing simple or constrained least-squares fit-type algorithms.

Estimation of time delays from a noisy and gapped data is one of the simplest data analysis problems in astronomy by its formulation. But as history of real experiments show, the work with observed data sets can be quite complex and evolved. By analysing in detail previous attempts to build delay estimation algorithms we try to develop an automatic and robust procedure to perform the task. To evaluate and compare different variants of the algorithms we use real observed data sets which have been objects of past controversies. In this way we hope to select the methods and procedures which have highest probability to succeed in complex situations. As a result of our investigations we propose an estimation procedure which can be used as a method of choice in large photometric experiments. We can not claim that proposed methodology works with any reasonably well sampled input data set. But we hope that the steps taken are in correct direction and developed software is truly useful for practising astronomers.

The sunspots and other solar activity indicators tend to cluster on the surface of the Sun.These clusters very often occur at certain longitudes that persist in time.It is of general interest to find new and simple ways to characterize the observed distributions of different indicators and their behaviour in time. In the present work we use Greenwich sunspot data to evaluate statistical but not totally coherent stability of sunspot distribution along latitudes as well as longitudes. The aim was to obtain information on the longitudinal distribution of the underlying spot-generating mechanism rather than on the distribution and migration of sunspots or sunspot groups on the solar surface. Therefore only sunspot groups were included in the analysis, and only the time of their first appearance was used. We use simple nonparametric approach to reveal sunspot migration patterns and their persistency. Our analysis shows that regions where spots are generated tend to rotate differentially as the spots and spot groups themselves do. The activity areas, however, tend to break down relatively fast, during 7-15 solar rotations. This study provides a challenge for solar dynamo models, as our results are consistent with the presence of a non-axisymmetric spot-generating mechanism experiencing differential rotation (known as phase mixing in dynamo theory). The new nonparametric method introduced here, completely independent of the choice of the longitudinal distribution of sunspots, was found to be a very powerful tool for spatio-temporal analysis of surface features.

Notions and limits from standard time series analysis must be modified when treating series which are measured irregularly and contain long gaps. Classical Nyquist criterion to estimate frequency range which is potentially recoverable must be modified to handle this more complex situation. When basic exposition of the modified criterion is given in earlier papers, some minor problems and caveats are treated here. Using simple combinatorial arguments we show that for small sample sizes the modified Nyquist limit may overestimate the obtainable frequency range. On the other hand we will demonstrate that very high Nyquist limit values which are typical to irregularly sampled data can often be taken seriously and using proper observational techniques the frequency ranges for "time spectroscopy" can be significantly widened.

Long time photometric monitoring programs of gravitational lens systems are often carried on using modest equipment. The resolution of such observations is limited and some of the images may remain unresolved. It may be still possible to find a full set of time delays from such a blended data. We discuss here a particular but interesting case when we have two light curves that both are blends. A suitable computational algorithm is developed and tested to work with computer-generated model light curves. Our method combines both blended sequences using the hypothetical time delays between the initial components of the light curves as free input parameters. The combined curves are then compared using statistical distance estimation. It occurs that using an assumption of equal magnification ratios between the components of the blends, we can indeed recover the whole set of time delays.

Longtime monitoring of gravitational lens systems is often done using telescopes and recording equipment with modest resolution. Still, it would be interesting to get as much information as possible from the measured lightcurves. From high resolution images we know that the recorded quasar images are often blends and that the corresponding time series are not pure shifted replicas of the source variability. In this paper we will develop an algorithm to unscramble this kind of blended data. The proposed method is based on a simple idea. We use one of the photometric curves, which is supposedly a simple shifted replica of the source curve, to build different artificial combined curves. Then we compare these artificial curves with the blended curves. Proper solutions for a full set of time delays are then obtained by varying free input parameters and estimating statistical distances between the artificial and blended curves. We performed a check of feasibility and applicability of the new algorithm. For numerically generated data sets the time delay systems were recovered for a wide range of setups. Application of the new algorithm to the classical double quasar QSO 0957+561 A,B lightcurves shows a clear splitting of one of the images. This is an unexpected result and extremely interesting, especially in the context of the recent controversy about the exact time delay value for the system. The proposed method allows us to properly analyse the data from low resolution observations that have long time coverages. There are a number of gravitational lens monitoring programmes that can make use of the new algorithm.

It has recently been claimed that analysis of Greenwich sunspot data over 120 years reveals that sunspot activity clusters around two longitudes separated by 180 degrees (``active longitudes'') with clearly defined differential rotation during activity cycles.In the present work we extend this critical examination of methodology to the actual Greenwich sunspot data and also consider newly proposed methods of analysis claiming to confirm the original identification of active longitudes. Our analysis revealed that values obtained for the parameters of differential rotation are not stable across different methods of analysis proposed to track persistent active longitudes. Also, despite a very thorough search in parameter space, we were unable to reproduce results claiming to reveal the century-persistent active longitudes. We can therefore say that strong and well substantiated evidence for an essential and century-scale persistent nonaxisymmetry in the sunspot distribution does not exist.

Using Greenwich sunspot data for 120 years it was recently observed that activity regions on the Sun's surface tend to lie along smoothly changing longitude strips 180 degrees apart from each other. However, numerical experiments with random input data show that most, if not all, of the observed longitude discrimination can be looked upon as an artifact of the analysis method.

We present a short re-evaluation of a recently published time delay estimate for the gravitational lens system HE 1104-1805 with emphasis on important methodological aspects: bias of the statistics, inconsistency of the methods and use of the purposeful selection of data points(or so-called "cleaning") at the preprocessing stage. We show how the inadequate use of simple analysis methods can lead to too strong conclusions. Our analysis shows that there are indications for the time delay in HE 1104-1805 to be between -0.9 and -0.7 years, but still with a large uncertainty.

We present optical lightcurves of the gravitationally lensed components A (=A1+A2+A3) and B of the quadruple quasar RX J0911.4+0551 (z = 2.80). The observations were primarily obtained at the Nordic Optical Telescope between 1997 March and 2001 April and consist of 74 I-band data points for each component. The data allow the measurement of a time delay of 146 +- 8 days (2 sigma) between A and B, with B as the leading component. This value is significantly shorter than that predicted from simple models and indicates a very large external shear. Mass models including the main lens galaxy and the surrounding massive cluster of galaxies at z = 0.77, responsible for the external shear, yield H_0 = 71 +- 4 (random, 2 sigma) +- 8 (systematic) km/s/Mpc. The systematic model uncertainty is governed by the surface-mass density (convergence) at the location of the multiple images.

The classical T Tauri star RW Aur A is an irregular variable with a large amplitude in all photometric bands. In an extended series of photometric data we found small-amplitude periodic variations in the blue colours of the star, with a period of 2.64 days. The period was relatively stable over several years. The amplitude of the periodic signal is 0.21 mag in U-V, 0.07 mag in B-V, and about 0.02 mag in V-R and V-I. No periodicity was found in the V magnitude. The relevance of this photometric period to the recently discovered periodicity in spectral features of the star is discussed, and the hypothesis of a hot spot is critically considered.

We present optical I-band light curves of the gravitationally lensed double QSO B1600+434 from observations obtained at the Nordic Optical Telescope (NOT) between April 1998 and November 1999. The photometry has been performed by simultaneous deconvolution of all the data frames, involving a numerical lens galaxy model. Four methods have been applied to determine the time delay between the two QSO components, giving a mean estimate of \Delta_t = 51+/-4 days (95% confidence level). This is the fourth optical time delay ever measured. Adopting a Omega=0.3, Lambda=0 Universe and using the mass model of Maller et al. (2000), this time-delay estimate yields a Hubble parameter of H_0=52 (+14, -8) km s^-1 Mpc^-1 (95% confidence level) where the errors include time-delay as well as model uncertainties. There are time-dependent offsets between the two (appropriately shifted) light curves that indicate the presence of external variations due to microlensing.

From regular monitoring of the Double Quasar QSO 0957+561 A,B there is now general agreement on a time delay of about 416 days. This has made it possible to determine the microlensing residual in the light-curve, see Pelt et al. (1998). We have used two significant microlensing features: 1) A ``quiet'' period with a variability less than 0.05 mag lasting about 8 years, and 2) A change in the residual of 0.25 mag during a time interval of about 5 years. The first feature gives a lower limit for the lens mass, 'M', for a given normalized source radius, 'r', whereas the second feature gives an upper limit. We have considered the amount of mass in a continuum to be a free parameter with possible values between 0% (all mass in lenses) and 90%. At a significance level of 1% the mass can only be constrained within a rather wide range (1e-06 M_solar to 5 M_solar). For the radius 'R' of the source an upper limit of 1e16 cm is found, whereas the normalized source radius 'r' is restricted to be smaller than 30. At a level of 10%, however, the range of possible masses is much narrower (2e-03 M_solar to 0.5 M_solar), and the upper limit of R is about 6e15 cm, whereas the value of 'r' is restricted to be less than 2. We have used an effective transverse velocity 'V' equal to 600 km/s.

We present a new analysis of the presently available photometric data for the gravitationally lensed quasar 0957+561 A,B with the aim of determining the time delay between its two images. The basic method used is the dispersion estimation technique. Even by using the simplest non-parametric form of our method we can convincingly rule out a time delay near 536 days and show that the time delay is in the vicinity of 420 days. We then introduce refinements to our method in order to get a stable and reliable result for the time delay independent of high frequency noise in the data and sampling errors. Our best result for the time delay, checked by various statistical tests and using bootstrap error estimates, is 423 +/- 6 days. We furthermore confirm our earlier result that the radio data are compatible with this value. Using the best available model for the mass distribution in the lensing galaxy and cluster, our result for the time delay constrains the Hubble parameter to be smaller than 70 km/(s Mpc).

We present a new analysis of previously published optical and radio data sets of the gravitationally lensed quasar 0957+561 A,B with the aim of determining the time delay between its two images. We use a non-parametric estimate of the dispersion of the combined data set where, however, we only make use of alternating neighbours in order to avoid windowing effects. From the optical data a time delay of (415 $\pm$ 32) days and from the radio data a delay of (409 $\pm$ 23) days is suggested. We demonstrate a considerable sensitivity of different delay estimation procedures against the removal of only a few observational data points or against smoothing or detrending of the original data sets. The radio data give us formally a slightly more precise value for the time delay than the optical data. Also, due to the lack of windowing effects, the result obtained for the radio data can be considered as somewhat more reliable than the delay determined from the optical data.