Showing 1–35 of 35 results for author: Sudar, D

Atacama Large Aperture Submillimeter Telescope (AtLAST) Science: Solar and stellar observations

Authors: Sven Wedemeyer, Miroslav Barta, Roman Brajsa, Yi Chai, Joaquim Costa, Dale Gary, Guillermo Gimenez de Castro, Stanislav Gunar, Gregory Fleishman, Antonio Hales, Hugh Hudson, Mats Kirkaune, Atul Mohan, Galina Motorina, Alberto Pellizzoni, Maryam Saberi, Caius L. Selhorst, Paulo J. A. Simoes, Masumi Shimojo, Ivica Skokic, Davor Sudar, Fabian Menezes, Stephen White, Mark Booth, Pamela Klaassen , et al. (13 additional authors not shown)

Abstract: Observations at (sub-)millimeter wavelengths offer a complementary perspective on our Sun and other stars, offering significant insights into both the thermal and magnetic composition of their chromospheres. Despite the fundamental progress in (sub-)millimeter observations of the Sun, some important aspects require diagnostic capabilities that are not offered by existing observatories. In particul… ▽ More Observations at (sub-)millimeter wavelengths offer a complementary perspective on our Sun and other stars, offering significant insights into both the thermal and magnetic composition of their chromospheres. Despite the fundamental progress in (sub-)millimeter observations of the Sun, some important aspects require diagnostic capabilities that are not offered by existing observatories. In particular, simultaneously observations of the radiation continuum across an extended frequency range would facilitate the mapping of different layers and thus ultimately the 3D structure of the solar atmosphere. Mapping large regions on the Sun or even the whole solar disk at a very high temporal cadence would be crucial for systematically detecting and following the temporal evolution of flares, while synoptic observations, i.e., daily maps, over periods of years would provide an unprecedented view of the solar activity cycle in this wavelength regime. As our Sun is a fundamental reference for studying the atmospheres of active main sequence stars, observing the Sun and other stars with the same instrument would unlock the enormous diagnostic potential for understanding stellar activity and its impact on exoplanets. The Atacama Large Aperture Submillimeter Telescope (AtLAST), a single-dish telescope with 50\,m aperture proposed to be built in the Atacama desert in Chile, would be able to provide these observational capabilities. Equipped with a large number of detector elements for probing the radiation continuum across a wide frequency range, AtLAST would address a wide range of scientific topics including the thermal structure and heating of the solar chromosphere, flares and prominences, and the solar activity cycle. In this white paper, the key science cases and their technical requirements for AtLAST are discussed. △ Less

Submitted 6 March, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 14 pages, 4 figures, submitted to Open Research Europe as part of a collection on the Atacama Large Aperture Submillimeter Telescope (AtLAST)

Calculated brightness temperatures of solar structures compared with ALMA and Metsähovi measurements

Authors: F. Matković, R. Brajša, M. Kuhar, A. O. Benz, H. -G. Ludwig, C. L. Selhorst, I. Skokić, D. Sudar, A. Hanslmeier

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) allows for solar observations in the wavelength range of 0.3$-$10 mm, giving us a new view of the chromosphere. The measured brightness temperature at various frequencies can be fitted with theoretical models of density and temperature versus height. We use the available ALMA and Metsähovi measurements of selected solar structures (quiet sun… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) allows for solar observations in the wavelength range of 0.3$-$10 mm, giving us a new view of the chromosphere. The measured brightness temperature at various frequencies can be fitted with theoretical models of density and temperature versus height. We use the available ALMA and Metsähovi measurements of selected solar structures (quiet sun (QS), active regions (AR) devoid of sunspots, and coronal holes (CH)). The measured QS brightness temperature in the ALMA wavelength range agrees well with the predictions of the semiempirical Avrett$-$Tian$-$Landi$-$Curdt$-$Wülser (ATLCW) model, better than previous models such as the Avrett$-$Loeser (AL) or Fontenla$-$Avrett$-$Loeser model (FAL). We scaled the ATLCW model in density and temperature to fit the observations of the other structures. For ARs, the fitted models require 9%$-$13% higher electron densities and 9%$-$10% higher electron temperatures, consistent with expectations. The CH fitted models require electron densities 2%$-$40% lower than the QS level, while the predicted electron temperatures, although somewhat lower, do not deviate significantly from the QS model. Despite the limitations of the one-dimensional ATLCW model, we confirm that this model and its appropriate adaptations are sufficient for describing the basic physical properties of the solar structures. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: 11 pages, 3 figures, 2 tables, accepted and published in Astronomische Nachrichten/Astronomical Notes

Harmonizing the Generation and Pre-publication Stewardship of FAIR Image Data

Authors: Nikki Bialy, Frank Alber, Brenda Andrews, Michael Angelo, Brian Beliveau, Lacramioara Bintu, Alistair Boettiger, Ulrike Boehm, Claire M. Brown, Mahmoud Bukar Maina, James J. Chambers, Beth A. Cimini, Kevin Eliceiri, Rachel Errington, Orestis Faklaris, Nathalie Gaudreault, Ronald N. Germain, Wojtek Goscinski, David Grunwald, Michael Halter, Dorit Hanein, John W. Hickey, Judith Lacoste, Alex Laude, Emma Lundberg , et al. (22 additional authors not shown)

Abstract: Together with the molecular knowledge of genes and proteins, biological images promise to significantly enhance the scientific understanding of complex cellular systems and to advance predictive and personalized therapeutic products for human health. For this potential to be realized, quality-assured image data must be shared among labs at a global scale to be compared, pooled, and reanalyzed, thu… ▽ More Together with the molecular knowledge of genes and proteins, biological images promise to significantly enhance the scientific understanding of complex cellular systems and to advance predictive and personalized therapeutic products for human health. For this potential to be realized, quality-assured image data must be shared among labs at a global scale to be compared, pooled, and reanalyzed, thus unleashing untold potential beyond the original purpose for which the data was generated. There are two broad sets of requirements to enable image data sharing in the life sciences. One set of requirements is articulated in the companion White Paper entitled Enabling Global Image Data Sharing in the Life Sciences, which is published in parallel and addresses the need to build the cyberinfrastructure for sharing the digital array data. In this White Paper, we detail a broad set of requirements, which involves collecting, managing, presenting, and propagating contextual information essential to assess the quality, understand the content, interpret the scientific implications, and reuse image data in the context of the experimental details. We start by providing an overview of the main lessons learned to date through international community activities, which have recently made considerable progress toward generating community standard practices for imaging Quality Control (QC) and metadata. We then provide a clear set of recommendations for amplifying this work. The driving goal is to address remaining challenges and democratize access to everyday practices and tools for a spectrum of biomedical researchers, regardless of their expertise, access to resources, and geographical location. △ Less

Submitted 30 August, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: This manuscript is published with a closely related companion entitled, Enabling Global Image Data Sharing in the Life Sciences, which can be found at the following link, arXiv:2401.13023 [q-bio.OT]

Differences in physical properties of coronal bright points and their ALMA counterparts within and outside coronal holes

Authors: F. Matković, R. Brajša, M. Temmer, S. G. Heinemann, H. -G. Ludwig, S. H. Saar, C. L. Selhorst, I. Skokić, D. Sudar

Abstract: This study investigates and compares brightness and area of coronal bright points (CBPs) inside and outside of coronal holes (CHs) using the single-dish Band 6 observations by the Atacama Large Millimeter/submillimeter Array (ALMA), combined with extreme-ultraviolet (EUV) 193 $\overset{\circ}{\mathrm{A}}$ filtergrams obtained by the Atmospheric Imaging Assembly (AIA) and magnetograms obtained by t… ▽ More This study investigates and compares brightness and area of coronal bright points (CBPs) inside and outside of coronal holes (CHs) using the single-dish Band 6 observations by the Atacama Large Millimeter/submillimeter Array (ALMA), combined with extreme-ultraviolet (EUV) 193 $\overset{\circ}{\mathrm{A}}$ filtergrams obtained by the Atmospheric Imaging Assembly (AIA) and magnetograms obtained by the Helioseismic and Magnetic Imager (HMI), both on board Solar Dynamics Observatory (SDO). The CH boundaries were extracted from the SDO/AIA images using the Collection of Analysis Tools for Coronal Holes (CATCH) and CBPs were identified in the SDO/AIA, SDO/HMI, and ALMA data. Measurements of brightness and areas in both ALMA and SDO/AIA images were conducted for CBPs within CHs and quiet Sun regions outside CHs. A statistical analysis of the measured physical properties resulted in a lower average CBP brightness in both ALMA and SDO/AIA data for CBPs within the CHs. Depending on the CBP sample size, the difference in intensity for the SDO/AIA data, and brightness temperature for the ALMA data, between the CBPs inside and outside CHs ranged from between 2$σ$ and 4.5$σ$, showing a statistically significant difference between those two CBP groups. For CBP areas, CBPs within the CH boundaries showed smaller areas on average, with the observed difference between the two CBP groups between 1$σ$ and 2$σ$ for the SDO/AIA data, and up to 3.5$σ$ for the ALMA data, indicating that CBP areas are also significantly different. Given the measured properties, we conclude that the CBPs inside CHs tend to be less bright on average, but also smaller in comparison to those outside of CHs. This conclusion might point to the specific physical conditions and properties of the local CH region around a CBP limiting the maximum achievable intensity (temperature) and size of a CBP. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 27 pages, 28 figures, accepted in Astronomy and Astrophysics

Journal ref: A&A 670, A146 (2023)

Flares detected in ALMA single-dish images of the Sun

Authors: I. Skokić, A. O. Benz, R. Brajša, D. Sudar, F. Matković, M. Bárta

Abstract: The (sub)millimeter radiation of solar flares is poorly understood. Without spatial resolution, it cannot be compared easily to flare emissions in other wavelengths. The Atacama Large Millimeter-submillimeter Array (ALMA) offers sufficient resolution for the first time. However, used as an interferometer, its field of view is smaller than an active region and ALMA cannot observe on demand when a f… ▽ More The (sub)millimeter radiation of solar flares is poorly understood. Without spatial resolution, it cannot be compared easily to flare emissions in other wavelengths. The Atacama Large Millimeter-submillimeter Array (ALMA) offers sufficient resolution for the first time. However, used as an interferometer, its field of view is smaller than an active region and ALMA cannot observe on demand when a flare occurs. We use readily available large scale single-dish ALMA observations of solar millimeter flares and compare them to well-known features observed in other wavelengths. The properties of these other flare emissions, correlating in space and time, may then be used to interpret the millimeter brightenings and vice versa. The aim is to obtain reliable associations, limited by the time and space resolution of single-dish observations. We collected ALMA observations at 3 mm and 1 mm and searched for millimeter brightenings during times given in a flare catalog. We found five events with 9 or more images that can be used for comparison in time and space. The millimeter brightenings are associated with a variety of flare features in cool (H$α$, 30.4 nm), intermediate (17.1 nm), and hot (9.4 nm) lines. In several cases, the millimeter brightening peaked at the footpoint of a hot flare loop. In other cases the peak coincided with the top or footpoint of an active Hα filament. We found correlations also with post-flare loops and tops of a hot loop, and in some cases to no features at all. The wide field of view provided by the single-dish observations allowed for completely overviewing the flare activity in millimeter waves for the first time. The associated phenomena often changed during the flare in type and location, and may explain the sometimes bewildering behavior of millimeter flare emissions observed previously without spatial resolution. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 13 pages, 10 figures, accepted in Astronomy and Astrophysics

Journal ref: A&A 669, A156 (2023)

V1294 Aql = HD 184279: A bad boy among Be stars or an important clue to the Be phenomenon?

Authors: P. Harmanec, H. Božić, P. Koubský, S. Yang, D. Ruždjak, D. Sudar, M. Šlechta, M. Wolf, D. Korčáková, P. Zasche, A. Oplištilová, D. Vršnak, H. Ak, P. Eenens, H. Bakiş, V. Bakiş, S. Otero, R. Chini, T. Demsky, B. N. Barlow, P. Svoboda, J. Jonák, K. Vitovský, A. Harmanec

Abstract: A reliable determination of the basic physical properties and variability patterns of hot emission-line stars is important for understanding the Be phenomenon and ultimately, the evolutionary stage of Be stars. This study is devoted to one of the most remarkable Be stars, V1294 Aql = HD 184279. We collected and analysed spectroscopic and photometric observations covering a time interval of about 2… ▽ More A reliable determination of the basic physical properties and variability patterns of hot emission-line stars is important for understanding the Be phenomenon and ultimately, the evolutionary stage of Be stars. This study is devoted to one of the most remarkable Be stars, V1294 Aql = HD 184279. We collected and analysed spectroscopic and photometric observations covering a time interval of about 25000 d (68 yr). We present evidence that the object is a single-line 192.9 d spectroscopic binary and estimate that the secondary probably is a hot compact object with a mass of about 1.1-1.2 solar masses. We found and documented very complicated orbital and long-term spectral, light, and colour variations, which must arise from a combination of several distinct variability patterns. Attempts at modelling them are planned for a follow-up study. We place the time behaviour of V1294 Aql into context with variations known for some other systematically studied Be stars and discuss the current ideas about the nature of the Be phenomenon. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Comments: 15 pages, 15 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 666, A136 (2022)

Variation in solar differential rotation and activity in the period 1964-2016 determined by the Kanzelhöhe data set

Authors: I. Poljančić Beljan, R. Jurdana-Šepić, T. Jurkić, R. Brajša, I. Skokić, D. Sudar, D. Ruždjak, D. Hržina, W. Pötzi, A. Hanslmeier, A. M. Veronig

Abstract: We determined the differential rotation (DR) parameters $A$ and $B$ (corresponding to the equatorial rotation velocity and the gradient of the solar DR) by tracing sunspot groups in sunspot drawings of the Kanzelhöhe Observatory for Solar and Environmental Research (KSO; 1964-2008, for solar cycles (SC) 20-23) and KSO white-light images (2009-2016, for SC 24). We used different statistical methods… ▽ More We determined the differential rotation (DR) parameters $A$ and $B$ (corresponding to the equatorial rotation velocity and the gradient of the solar DR) by tracing sunspot groups in sunspot drawings of the Kanzelhöhe Observatory for Solar and Environmental Research (KSO; 1964-2008, for solar cycles (SC) 20-23) and KSO white-light images (2009-2016, for SC 24). We used different statistical methods and approaches to analyse cycle related variations, solar cycle phase-related variations and long-term variations of the DR. $A$ and $B$ show statistically significant periodic variability. The changes in $A$ related to solar cycle phase are in accordance with previously reported theoretical and experimental results (higher $A$ during solar minimum, lower $A$ during the maximum of activity), while changes in $B$ differ from the theoretical predictions as we observe more negative values of $B$, that is, a more pronounced DR during activity maximum. The main result of this paper for the long-term variations in $A$ is the detection of a phase shift between the activity flip (in the 1970s) and the equatorial rotation velocity flip (in the early 1990s). During this time period both $A$ and activity show a secular decreasing trend, indicating their correlation. Therefore, the theoretical model fails in the phase-shift time period that occurs after the modern Gleissberg maximum, while in the time period thereafter (after the 1990s), theoretical and experimental results are consistent. The long-term variations in $B$ in general yield an anticorrelation of $B$ and activity, as a rise of $B$ is observed during the entire time period (1964-2016) we analysed, during which activity decreased. We study for the first time the variation in solar DR and activity based on 53 years of KSO data. Our results agree well with the results related to the solar cycle phase from corona observations. △ Less

Submitted 15 April, 2022; originally announced April 2022.

Journal ref: A&A 663, A24 (2022)

Correlation between the solar magnetic field strength and the millimeter brightness temperature

Authors: I. Skokić, D. Sudar, R. Brajša

Abstract: Images of the Sun at millimeter wavelengths obtained by ALMA show a significant correspondence with the magnetograms. In this paper, we investigate this correspondence by comparing ALMA full-disk solar image taken at 1.2 mm with a SDO/HMI magnetogram and analyze their correlation. It is found that chromospheric network and active regions show a positive correlation where brightness temperature is… ▽ More Images of the Sun at millimeter wavelengths obtained by ALMA show a significant correspondence with the magnetograms. In this paper, we investigate this correspondence by comparing ALMA full-disk solar image taken at 1.2 mm with a SDO/HMI magnetogram and analyze their correlation. It is found that chromospheric network and active regions show a positive correlation where brightness temperature is increasing with the line-of-sight magnetic field strength, while sunspots have a negative correlation. Quiet Sun regions do not show any dependence of the brightness temperature with the magnetic field. Thermal bremsstrahlung is given as the best explanation for the observed correlations. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 11 pages, 4 figures, published in Central European Astrophysical Bulletin, 2020

Journal ref: Central European Astrophysical Bulletin, 44, (2020) 2, 1-11

A prediction for the 25th solar cycle maximum amplitude

Authors: R. Brajša, G. Verbanac, M. Bandić, A. Hanslmeier, I. Skokić, D. Sudar

Abstract: The minimum - maximum method, belonging to the precursor class of the solar activity forecasting methods, is based on a linear relationship between relative sunspot number in the minimum and maximum epochs of solar cycles. In the present analysis we apply a modified version of this method using data not only from the minimum year, but also from a couple of years before and after the minimum. The r… ▽ More The minimum - maximum method, belonging to the precursor class of the solar activity forecasting methods, is based on a linear relationship between relative sunspot number in the minimum and maximum epochs of solar cycles. In the present analysis we apply a modified version of this method using data not only from the minimum year, but also from a couple of years before and after the minimum. The revised 13-month smoothed monthly total sunspot number data set from SILSO/SIDC is used. Using data for solar cycle nos. 1-24 the largest correlation coefficient (CC) is obtained when correlating activity level 3 years before solar cycle minimum with the subsequent maximum (CC = 0.82), independent of inclusion or exclusion of the solar cycle no. 19. For the next solar maximum of the cycle no. 25 we predict: Rmax = 121 +- 33. Our results indicate that the next solar maximum (of the cycle no. 25) will be of the similar amplitude as the previous one, or even something lower. This is in accordance with the general middle-term lowering of the solar activity after the secular maximum in the 20th century and consistent with the Gleissberg period of the solar activity. The reliability of the 3 years before the minimum predictor is experimentally justified by the largest correlation coefficient and verified with the Student t-test. It is satisfactorily explained with the two empirical well-known findings: the extended solar cycle and the Waldmeier effect. Finally, we successfully reproduced the maxima of the last four solar cycles, nos. 21-25, using the 3 years before the minimum method. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 10 pages, 4 figures, accepted by Astronomische Nachrichten / Astronomical Notes

MITI Minimum Information guidelines for highly multiplexed tissue images

Authors: Denis Schapiro, Clarence Yapp, Artem Sokolov, Sheila M. Reynolds, Yu-An Chen, Damir Sudar, Yubin Xie, Jeremy L. Muhlich, Raquel Arias-Camison, Sarah Arena, Adam J. Taylor, Milen Nikolov, Madison Tyler, Jia-Ren Lin, Erik A. Burlingame, Human Tumor Atlas Network, Young H. Chang, Samouil L Farhi, Vésteinn Thorsson, Nithya Venkatamohan, Julia L. Drewes, Dana Pe'er, David A. Gutman, Markus D. Herrmann, Nils Gehlenborg , et al. (14 additional authors not shown)

Abstract: The imminent release of tissue atlases combining multi-channel microscopy with single cell sequencing and other omics data from normal and diseased specimens creates an urgent need for data and metadata standards that guide data deposition, curation and release. We describe a Minimum Information about highly multiplexed Tissue Imaging (MITI) standard that applies best practices developed for genom… ▽ More The imminent release of tissue atlases combining multi-channel microscopy with single cell sequencing and other omics data from normal and diseased specimens creates an urgent need for data and metadata standards that guide data deposition, curation and release. We describe a Minimum Information about highly multiplexed Tissue Imaging (MITI) standard that applies best practices developed for genomics and other microscopy data to highly multiplexed tissue images and traditional histology. △ Less

Submitted 23 February, 2022; v1 submitted 21 August, 2021; originally announced August 2021.

Probabilistic Drag-Based Ensemble Model (DBEM) Evaluation for Heliospheric Propagation of CMEs

Authors: Jaša Čalogović, Mateja Dumbović, Davor Sudar, Bojan Vršnak, Karmen Martinić, Manuela Temmer, Astrid Veronig

Abstract: The Drag-based Model (DBM) is a 2D analytical model for heliospheric propagation of Coronal Mass Ejections (CMEs) in ecliptic plane predicting the CME arrival time and speed at Earth or any other given target in the solar system. It is based on the equation of motion and depends on initial CME parameters, background solar wind speed, $w$ and the drag parameter $γ$. A very short computational time… ▽ More The Drag-based Model (DBM) is a 2D analytical model for heliospheric propagation of Coronal Mass Ejections (CMEs) in ecliptic plane predicting the CME arrival time and speed at Earth or any other given target in the solar system. It is based on the equation of motion and depends on initial CME parameters, background solar wind speed, $w$ and the drag parameter $γ$. A very short computational time of DBM ($<$ 0.01s) allowed us to develop the Drag-Based Ensemble Model (DBEM) that takes into account the variability of model input parameters by making an ensemble of n different input parameters to calculate the distribution and significance of the DBM results. Thus the DBEM is able to calculate the most likely CME arrival times and speeds, quantify the prediction uncertainties and determine the confidence intervals. A new DBEMv3 version is described in detail and evaluated for the first time determing the DBEMv3 performance and errors by using various CME-ICME lists as well as it is compared with previous DBEM versions. The analysis to find the optimal drag parameter $γ$ and ambient solar wind speed $w$ showed that somewhat higher values ($γ\approx 0.3 \times 10^{-7}$ km$^{-1}$, $w \approx$ 425 km\,s$^{-1}$) for both of these DBEM input parameters should be used for the evaluation compared to the previously employed ones. Based on the evaluation performed for 146 CME-ICME pairs, the DBEMv3 performance with mean error (ME) of -11.3 h, mean absolute error (MAE) of 17.3 h was obtained. There is a clear bias towards the negative prediction errors where the fast CMEs are predicted to arrive too early, probably due to the model physical limitations and input errors (e.g. CME launch speed). △ Less

Submitted 14 July, 2021; originally announced July 2021.

Journal ref: Solar Physics 2021

ALMA small-scale features in the quiet Sun and active regions

Authors: R. Brajsa, I. Skokic, D. Sudar, A. O. Benz, S. Krucker, H. -G. Ludwig, S. H. Saar, C. L. Selhorst

Abstract: Aims. The main aim of the present analysis is to decipher (i) the small-scale bright features in solar images of the quiet Sun and active regions obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) and (ii) the ALMA correspondence of various known chromospheric structures visible in the H-alpha images of the Sun. Methods. Small-scale ALMA bright features in the quiet Sun region w… ▽ More Aims. The main aim of the present analysis is to decipher (i) the small-scale bright features in solar images of the quiet Sun and active regions obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) and (ii) the ALMA correspondence of various known chromospheric structures visible in the H-alpha images of the Sun. Methods. Small-scale ALMA bright features in the quiet Sun region were analyzed using single-dish ALMA observations (1.21 mm, 248 GHz) and in an active region using interferometric ALMA measurements (3 mm, 100 GHz). With the single-dish observations, a full-disk solar image is produced, while interferometric measurements enable the high-resolution reconstruction of part of the solar disk, including the active region. The selected quiet Sun and active regions are compared with the H-alpha (core and wing sum), EUV, and soft X-ray images and with the magnetograms. Results. In the quiet Sun region, enhanced emission seen in the ALMA is almost always associated with a strong line-of-sight (LOS) magnetic field. Four coronal bright points were identified, while other small-scale ALMA bright features are most likely associated with magnetic network elements and plages. In the active region, in 14 small-scale ALMA bright features randomly selected and compared with other images, we found five good candidates for coronal bright points, two for plages, and five for fibrils. Two unclear cases remain: a fibril or a jet, and a coronal bright point or a plage. A comparison of the H-alpha core image and the 3 mm ALMA image of the analyzed active region showed that the sunspot appears dark in both images (with a local ALMA radiation enhancement in sunspot umbra), the four plage areas are bright in both images and dark small H-alpha filaments are clearly recognized as dark structures of the same shape also in ALMA. △ Less

Submitted 8 May, 2021; originally announced May 2021.

Comments: 12 pages, 4 figures, to be published in Astronomy and Astrophysics

Journal ref: A&A 651, A6 (2021)

Drag-based model (DBM) tools for forecast of coronal mass ejection arrival time and speed

Authors: Mateja Dumbovic, Jasa Calogovic, Karmen Martinic, Bojan Vrsnak, Davor Sudar, Manuela Temmer, Astrid Veronig

Abstract: Forecasting the arrival time of CMEs and their associated shocks is one of the key aspects of space weather research. One of the commonly used models is, due to its simplicity and calculation speed, the analytical drag-based model (DBM) for heliospheric propagation of CMEs. DBM relies on the observational fact that slow CMEs accelerate whereas fast CMEs decelerate, and is based on the concept of M… ▽ More Forecasting the arrival time of CMEs and their associated shocks is one of the key aspects of space weather research. One of the commonly used models is, due to its simplicity and calculation speed, the analytical drag-based model (DBM) for heliospheric propagation of CMEs. DBM relies on the observational fact that slow CMEs accelerate whereas fast CMEs decelerate, and is based on the concept of MHD drag, which acts to adjust the CME speed to the ambient solar wind. Although physically DBM is applicable only to the CME magnetic structure, it is often used as a proxy for the shock arrival. In recent years, the DBM equation has been used in many studies to describe the propagation of CMEs and shocks with different geometries and assumptions. Here we give an overview of the five DBM versions currently available and their respective tools, developed at Hvar Observatory and frequently used by researchers and forecasters. These include: 1) basic 1D DBM, a 1D model describing the propagation of a single point (i.e. the apex of the CME) or concentric arc (where all points propagate identically); 2) advanced 2D self-similar cone DBM, a 2D model which combines basic DBM and cone geometry describing the propagation of the CME leading edge which evolves self-similarly; 3) 2D flattening cone DBM, a 2D model which combines basic DBM and cone geometry describing the propagation of the CME leading edge which does not evolve self-similarly; 4) DBEM, an ensemble version of the 2D flattening cone DBM which uses CME ensembles as an input and 5) DBEMv3, an ensemble version of the 2D flattening cone DBM which creates CME ensembles based on the input uncertainties. All five versions have been tested and published in recent years and are available online or upon request. We provide an overview of these five tools, of their similarities and differences, as well as discuss and demonstrate their application. △ Less

Submitted 26 March, 2021; originally announced March 2021.

Comments: 21 pages, 4 figures, 2 tables, accepted in Frontiers March 24th 2021

QUAREP-LiMi: A community-driven initiative to establish guidelines for quality assessment and reproducibility for instruments and images in light microscopy

Authors: Glyn Nelson, Ulrike Boehm, Steve Bagley, Peter Bajcsy, Johanna Bischof, Claire M Brown, Aurelien Dauphin, Ian M Dobbie, John E Eriksson, Orestis Faklaris, Julia Fernandez-Rodriguez, Alexia Ferrand, Laurent Gelman, Ali Gheisari, Hella Hartmann, Christian Kukat, Alex Laude, Miso Mitkovski, Sebastian Munck, Alison J North, Tobias M Rasse, Ute Resch-Genger, Lucas C Schuetz, Arne Seitz, Caterina Strambio-De-Castillia , et al. (75 additional authors not shown)

Abstract: In April 2020, the QUality Assessment and REProducibility for Instruments and Images in Light Microscopy (QUAREP-LiMi) initiative was formed. This initiative comprises imaging scientists from academia and industry who share a common interest in achieving a better understanding of the performance and limitations of microscopes and improved quality control (QC) in light microscopy. The ultimate goal… ▽ More In April 2020, the QUality Assessment and REProducibility for Instruments and Images in Light Microscopy (QUAREP-LiMi) initiative was formed. This initiative comprises imaging scientists from academia and industry who share a common interest in achieving a better understanding of the performance and limitations of microscopes and improved quality control (QC) in light microscopy. The ultimate goal of the QUAREP-LiMi initiative is to establish a set of common QC standards, guidelines, metadata models, and tools, including detailed protocols, with the ultimate aim of improving reproducible advances in scientific research. This White Paper 1) summarizes the major obstacles identified in the field that motivated the launch of the QUAREP-LiMi initiative; 2) identifies the urgent need to address these obstacles in a grassroots manner, through a community of stakeholders including, researchers, imaging scientists, bioimage analysts, bioimage informatics developers, corporate partners, funding agencies, standards organizations, scientific publishers, and observers of such; 3) outlines the current actions of the QUAREP-LiMi initiative, and 4) proposes future steps that can be taken to improve the dissemination and acceptance of the proposed guidelines to manage QC. To summarize, the principal goal of the QUAREP-LiMi initiative is to improve the overall quality and reproducibility of light microscope image data by introducing broadly accepted standard practices and accurately captured image data metrics. △ Less

Submitted 27 January, 2021; v1 submitted 21 January, 2021; originally announced January 2021.

Comments: 17 pages, 3 figures, shortened abstract, Co-Lead Authors: Glyn Nelson and Ulrike Boehm, Corresponding author: Roland Nitschke

Journal ref: J. Microsc. 2021;1-18

A perspective on Microscopy Metadata: data provenance and quality control

Authors: Maximiliaan Huisman, Mathias Hammer, Alex Rigano, Ulrike Boehm, James J. Chambers, Nathalie Gaudreault, Alison J. North, Jaime A. Pimentel, Damir Sudar, Peter Bajcsy, Claire M. Brown, Alexander D. Corbett, Orestis Faklaris, Judith Lacoste, Alex Laude, Glyn Nelson, Roland Nitschke, David Grunwald, Caterina Strambio-De-Castillia

Abstract: The application of microscopy in biomedical research has come a long way since Antonie van Leeuwenhoek discovered unicellular organisms. Countless innovations have positioned light microscopy as a cornerstone of modern biology and a method of choice for connecting omics datasets to their biological and clinical correlates. Still, regardless of how convincing published imaging data looks, it does n… ▽ More The application of microscopy in biomedical research has come a long way since Antonie van Leeuwenhoek discovered unicellular organisms. Countless innovations have positioned light microscopy as a cornerstone of modern biology and a method of choice for connecting omics datasets to their biological and clinical correlates. Still, regardless of how convincing published imaging data looks, it does not always convey meaningful information about the conditions in which it was acquired, processed, and analyzed. Adequate record-keeping, reporting, and quality control are therefore essential to ensure experimental rigor and data fidelity, allow experiments to be reproducibly repeated, and promote the proper evaluation, interpretation, comparison, and re-use. To this end, microscopy images should be accompanied by complete descriptions detailing experimental procedures, biological samples, microscope hardware specifications, image acquisition parameters, and image analysis procedures, as well as metrics accounting for instrument performance and calibration. However, universal, community-accepted Microscopy Metadata standards and reporting specifications that would result in Findable Accessible Interoperable and Reproducible (FAIR) microscopy data have not yet been established. To understand this shortcoming and to propose a way forward, here we provide an overview of the nature of microscopy metadata and its importance for fostering data quality, reproducibility, scientific rigor, and sharing value in light microscopy. The proposal for tiered Microscopy Metadata Specifications that extend the OME Data Model put forth by the 4D Nucleome Initiative and by Bioimaging North America [1-3] as well as a suite of three complementary and interoperable tools are being developed to facilitate the process of image data documentation and are presented in related manuscripts [4-6]. △ Less

Submitted 31 May, 2021; v1 submitted 24 October, 2019; originally announced October 2019.

Centre to limb brightness variations from ALMA full disk solar images

Authors: Davor Sudar, Roman Brajša, Ivica Skokić, Arnold O. Benz

Abstract: Science Verification (SV) data of solar observations with Atacama Large Millimeter-submillimeter Array (ALMA) telescope were released to the scientific community. Understanding the centre to limb brightness function is necessary to compare features in full disk images. Our goals are to find the empirical centre to limb brightness functions in two available spectral bands and create flattened image… ▽ More Science Verification (SV) data of solar observations with Atacama Large Millimeter-submillimeter Array (ALMA) telescope were released to the scientific community. Understanding the centre to limb brightness function is necessary to compare features in full disk images. Our goals are to find the empirical centre to limb brightness functions in two available spectral bands and create flattened images with centre to limb brightness variations removed. We used second-order polynomial fit of the cosine of incidence angle to data points as a function of radial distance to the centre of the solar disk. The method also includes iterative removal of outliers based on the interquartile range. Fitting functions for all available images proved to adequately describe the data with comparatively small errors in the fitting coefficients. In both bands we found brightening towards the limb which is a consequence of increase in electron temperatures with radial distance in this region of the solar atmosphere. This study found that the $T_{b}$ of an active region has about 180 K difference between with and without the limb brightening at radial distance $\approx$0.75$R_{\odot}$ in Band 6. We also made flattened images with limb brightening removed. The limb brightening effect in ALMA images is significant enough (of the order of 10% for Band 3 and about 15% in Band 6) that it can not be neglected in further analyses. Since the effect of the side lobes was not included in this study, these values probably represent the lower limit of the limb brightening. The shape of the limb brightening function can also be used to constrain electron densities and temperatures in various layers of the solar atmosphere. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Properties and nature of Be stars 31. The binary nature, light variability, physical elements, and emission-line changes of HD~81357

Authors: P. Koubský, P. Harmanec, M. Brož, L. Kotková, S. Yang, H. Božić, D. Sudar, Y. Frémat, D. Korčákov'a, V. Votruba, P. Škoda, M. Šlechta, D. Ruždjak

Abstract: Reliable determination of the basic physical properties of hot emission-line binaries with Roche-lobe filling secondaries is important for developing the theory of mass exchange in binaries. It is not easy, however, due to the presence of circumstellar matter. Here, we report the first detailed investigation of a new representative of this class of binaries, HD~81357, based on the analysis of spec… ▽ More Reliable determination of the basic physical properties of hot emission-line binaries with Roche-lobe filling secondaries is important for developing the theory of mass exchange in binaries. It is not easy, however, due to the presence of circumstellar matter. Here, we report the first detailed investigation of a new representative of this class of binaries, HD~81357, based on the analysis of spectra and photometry from several observatories. HD~81357 was found to be a double-lined spectroscopic binary and an ellipsoidal variable seen under an intermediate orbital inclination of $\sim(63\pm5)^\circ$, having an orbital period of 33\fd77445(41) and a~circular orbit. From an automated comparison of the observed and synthetic spectra, we estimate the component's effective temperatures to be 12930(540)~K and 4260(24)~K. The combined light-curve and orbital solutions, also constrained by a very accurate Gaia Data Release 2 parallax, give the following values of the basic physical properties: masses $3.36\pm0.15$ and $0.34\pm0.04$~\Mnom, radii $3.9\pm0.2$ and 13.97\pm0.05$~\Rnom, and a~mass ratio $10.0\pm0.5$. Evolutionary modelling of the system including the phase of mass transfer between the components indicated that HD~81357 is a~system observed in the final slow phase of the mass exchange after the mass-ratio reversal. Contrary to what has been seen for similar binaries like AU~Mon, no cyclic light variations were found on a~time scale an~order of magnitude longer than the orbital period. 243,1 15% △ Less

Submitted 7 August, 2019; originally announced August 2019.

Comments: 16 pages, 9 figures; accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 629, A105 (2019)

Observations of the solar chromosphere with ALMA and comparison with theoretical models

Authors: R. Brajša, D. Sudar, I. Skokic, A. O. Benz, M. Kuhar, A. Kobelski, S. Wedemeyer, S. M. White, H. -G. Ludwig, M. Temmer, S. H. Saar, C. L. Selhorst

Abstract: In this work we use solar observations with the ALMA radio telescope at the wavelength of 1.21 mm. The aim of the analysis is to improve understanding of the solar chromosphere, a dynamic layer in the solar atmosphere between the photosphere and corona. The study has an observational and a modeling part. In the observational part full-disc solar images are analyzed. Based on a modified FAL atmosph… ▽ More In this work we use solar observations with the ALMA radio telescope at the wavelength of 1.21 mm. The aim of the analysis is to improve understanding of the solar chromosphere, a dynamic layer in the solar atmosphere between the photosphere and corona. The study has an observational and a modeling part. In the observational part full-disc solar images are analyzed. Based on a modified FAL atmospheric model, radiation models for various observed solar structures are developed. Finally, the observational and modeling results are compared and discussed. △ Less

Submitted 18 December, 2018; originally announced December 2018.

Comments: 4 pages, presented at The 20th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, July 29 - Aug 3 2018, Boston / Cambridge, USA

Meridional Motions and Reynolds Stress Determined by Using Kanzelhöhe Drawings and White Light Solar Images from 1964 to 2016

Authors: Domagoj Ruždjak, Davor Sudar, Roman Brajša, Ivica Skokić, Ivana Poljančić Beljan, Rajka Jurdana-Šepić, Arnold Hanslmeier, Astrid Veronig, Werner Pötzi

Abstract: Sunspot position data obtained from Kanzelhöhe Observatory for Solar and Environmental Research (KSO) sunspot drawings and white light images in the period 1964 to 2016 were used to calculate the rotational and meridional velocities of the solar plasma. Velocities were calculated from daily shifts of sunspot groups and an iterative process of calculation of the differential rotation profiles was u… ▽ More Sunspot position data obtained from Kanzelhöhe Observatory for Solar and Environmental Research (KSO) sunspot drawings and white light images in the period 1964 to 2016 were used to calculate the rotational and meridional velocities of the solar plasma. Velocities were calculated from daily shifts of sunspot groups and an iterative process of calculation of the differential rotation profiles was used to discard outliers. We found a differential rotation profile and meridional motions in agreement with previous studies using sunspots as tracers and conclude that the quality of the KSO data is appropriate for analysis of solar velocity patterns. By analysing the correlation and covariance of meridional velocities and rotation rate residuals we found that the angular momentum is transported towards the solar equator. The magnitude and latitudinal dependence of the horizontal component of the Reynolds stress tensor calculated is sufficient to maintain the observed solar differential rotation profile. Therefore, our results confirm that the Reynolds stress is the dominant mechanism responsible for transport of angular momentum towards the solar equator. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: 14 pages, 6 figures

First analysis of solar structures in 1.21 mm full-disc ALMA image of the Sun

Authors: R. Brajša, D. Sudar, A. O. Benz, I. Skokić, M. Bárta, B. De Pontieu, S. Kim, A. Kobelski, M. Kuhar, M. Shimojo, S. Wedemeyer, S. White, P. Yagoubov, Y. Yan

Abstract: Various solar features can be seen on maps of the Sun in the mm and sub-mm wavelength range. The recently installed Atacama Large Millimeter/submillimeter Array (ALMA) is capable of observing the Sun in that wavelength range with an unprecedented spatial, temporal and spectral resolution. To interpret solar observations with ALMA the first important step is to compare ALMA maps with simultaneous i… ▽ More Various solar features can be seen on maps of the Sun in the mm and sub-mm wavelength range. The recently installed Atacama Large Millimeter/submillimeter Array (ALMA) is capable of observing the Sun in that wavelength range with an unprecedented spatial, temporal and spectral resolution. To interpret solar observations with ALMA the first important step is to compare ALMA maps with simultaneous images of the Sun recorded in other spectral ranges. First we identify different structures in the solar atmosphere seen in the optical, IR and EUV parts of the spectrum (quiet Sun (QS), active regions (AR), prominences on the disc, magnetic inversion lines (IL), coronal holes (CH) and coronal bright points (CBPs)) in a full disc solar ALMA image. The second aim is to measure the intensities (brightness temperatures) of those structures and compare them with the corresponding QS level. A full disc solar image at 1.21 mm obtained on December 18, 2015 during a CSV-EOC campaign with ALMA is calibrated and compared with full disc solar images from the same day in Hα, in He I 1083 nm core, and with SDO images (AIA at 170 nm, 30.4 nm, 21.1 nm, 19.3 nm, and 17.1 nm and HMI magnetogram). The brightness temperatures of various structures are determined by averaging over corresponding regions of interest in the ALMA image. Positions of the QS, ARs, prominences on the disc, ILs, CHs and CBPs are identified in the ALMA image. At 1.21 mm ARs appear as bright areas (but sunspots are dark), while prominences on the disc and CHs are not discernible from the QS background, although having slightly less intensity than surrounding QS regions. ILs appear as large, elongated dark structures and CBPs correspond to ALMA bright points. These results are in general agreement with sparse earlier measurements at similar wavelengths. The identification of CBPs represents the most important new result. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 9 pages, 3 figures

Journal ref: A&A 613, A17 (2018)

A Relationship between the Solar Rotation and Activity Analysed by Tracing Sunspot Groups

Authors: Domagoj Ruždjak, Roman Brajša, Davor Sudar, Ivica Skokić, Ivana Poljančić Beljan

Abstract: The sunspot position from Greenwich Photoheliographic Results (GPR), US Air Force Solar Optical Observing Network and National Oceanic and Atmospheric Administration (USAF/NOAA), and Debrecen Photoheliographic Data (DPD) data bases in the period 1874 to 2016 were used to calculate yearly values of the solar differential-rotation parameters $A$ and $B$. The calculated differential-rotation paramete… ▽ More The sunspot position from Greenwich Photoheliographic Results (GPR), US Air Force Solar Optical Observing Network and National Oceanic and Atmospheric Administration (USAF/NOAA), and Debrecen Photoheliographic Data (DPD) data bases in the period 1874 to 2016 were used to calculate yearly values of the solar differential-rotation parameters $A$ and $B$. The calculated differential-rotation parameters were compared with the solar-activity level. We found that the Sun rotates more differentially at the minimum than at the maximum of activity during the 1977 - 2016 epoch. An inverse correlation between equatorial rotation and solar activity was found using the recently revised sunspot number. The secular decrease of equatorial rotation rate accompanying the increase of activity stopped in the last part of the 20th century. It was noted that when a significant peak of equatorial rotation velocity is observed during minimum of activity, the strength of the next maximum is smaller than the previous one. △ Less

Submitted 10 November, 2017; originally announced November 2017.

Comments: 13pages, 3 figures

Solar differential rotation in the period 1964 - 2016 determined by the Kanzelhöhe data set

Authors: I. Poljančić Beljan, R. Jurdana-Šepić, R. Brajša, D. Sudar, D. Ruždjak, D. Hržina, W. Pötzi, A. Hanslmeier, A. Veronig, I. Skokić, H. Wöhl

Abstract: The main aim of this work is to determine the solar differential rotation by tracing sunspot groups during the period 1964-2016, using the Kanzelhöhe Observatory for Solar and Environmental Research (KSO) sunspot drawings and white light images. Two procedures for the determination of the heliographic positions were applied: an interactive procedure on the KSO sunspot drawings (1964 - 2008, solar… ▽ More The main aim of this work is to determine the solar differential rotation by tracing sunspot groups during the period 1964-2016, using the Kanzelhöhe Observatory for Solar and Environmental Research (KSO) sunspot drawings and white light images. Two procedures for the determination of the heliographic positions were applied: an interactive procedure on the KSO sunspot drawings (1964 - 2008, solar cycles nos. 20 - 23) and an automatic procedure on the KSO white light images (2009 - 2016, solar cycle no. 24). For the determination of the synodic angular rotation velocities two different methods have been used: a daily shift (DS) method and a robust linear least-squares fit (rLSQ) method. Afterwards, the rotation velocities had to be converted from synodic to sidereal, which were then used in the least-squares fitting for the solar differential rotation law. For the test data from 2014, we found the rLSQ method for calculating rotational velocities to be more reliable than the DS method. The best fit solar differential rotation profile for the whole time period is $ω(b)$ = (14.47 $\pm$ 0.01) - (2.66 $\pm$ 0.10) $\sin^2b$ (deg/day) for the DS method and $ω(b)$ = (14.50 $\pm$ 0.01) - (2.87 $\pm$ 0.12) $\sin^2b$ (deg/day) for the rLSQ method. A barely noticeable north - south asymmetry is observed for the whole time period 1964 - 2016 in the present paper. Rotation profiles, using different data sets (e.g. Debrecen Photoheliographic Data, Greenwich Photoheliographic Results), presented by other authors for the same time periods and the same tracer types, are in good agreement with our results. Therefore, the KSO data set is suitable for the investigation of the long-term variabilities in the solar rotation profile. △ Less

Submitted 27 July, 2017; v1 submitted 25 July, 2017; originally announced July 2017.

Journal ref: A&A 606, A72 (2017)

Meridional Motion and Reynolds Stress from Debrecen Photoheliographic Data

Authors: Davor Sudar, Roman Brajša, Ivica Skokić, Ivana Poljančić Beljan, Hubertus Wöhl

Abstract: The Debrecen Photoheliographic Data catalogue is a continuation of the Greenwich Photoheliographic Results providing daily positions of sunspots and sunspot groups. We analyse the data for sunspot groups focusing on meridional motions and transfer of angular momentum towards the solar equator. Velocities are calculated with a daily shift method including an automatic iterative process of removing… ▽ More The Debrecen Photoheliographic Data catalogue is a continuation of the Greenwich Photoheliographic Results providing daily positions of sunspots and sunspot groups. We analyse the data for sunspot groups focusing on meridional motions and transfer of angular momentum towards the solar equator. Velocities are calculated with a daily shift method including an automatic iterative process of removing the outliers. Apart from the standard differential rotation profile, we find meridional motion directed towards the zone of solar activity. The difference in measured meridional flow in comparison to Doppler measurements and some other tracer measurements is interpreted as a consequence of different flow patterns inside and outside of active regions. We also find a statistically significant dependence of meridional motion on rotation velocity residuals confirming the transfer of angular momentum towards the equator. Analysis of horizontal Reynolds stress reveals that the transfer of angular momentum is stronger with increasing latitude up to about 40$^{\circ}$ where there is a possible maximum in absolute value. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: 16 pages, 5 figures

Meridional motions and Reynolds stress from SDO/AIA coronal bright points data

Authors: Davor Sudar, Steven H. Saar, Ivica Skokić, Ivana Poljančić Beljan, Roman Brajša

Abstract: Context. It is possible to detect and track coronal bright points (CBPs) in SDO/AIA images. Combination of high resolution and high cadence provides a wealth of data that can be used to determine velocity flows on the solar surface with very high accuracy. Aims. We derived a very accurate solar rotation profile and investigated meridional flows, torsional oscillations and horizontal Reynolds stres… ▽ More Context. It is possible to detect and track coronal bright points (CBPs) in SDO/AIA images. Combination of high resolution and high cadence provides a wealth of data that can be used to determine velocity flows on the solar surface with very high accuracy. Aims. We derived a very accurate solar rotation profile and investigated meridional flows, torsional oscillations and horizontal Reynolds stress based on $\approx$6 months of SDO/AIA data. Methods. We used a segmentation algorithm to detect CBPs in SDO/AIA images. We also used invariance of the solar rotation profile with central meridian distance (CMD) to determine the height of CBPs in 19.3 nm channel. Results. Best fit solar rotation profile is given by $ω(b)=(14.4060\pm0.0051 + (-1.662\pm0.050)\sin^{2}b + (-2.742\pm0.081)\sin^{4}b)${\degr} day$^{-1}$. Height of CBPs in SDO/AIA 19.3 nm channel was found to be $\approx$6500 km. Meridional motion is predominantly poleward for all latitudes, while solar velocity residuals show signs of torsional oscillations. Horizontal Reynolds stress was found to be small compared to similar works, but still showing transfer of angular momentum towards the solar equator. Conclusion. Most of the results are consistent with Doppler measurements rather than tracer measurements. Fairly small calculated value of horizontal Reynolds stress might be due to the particular phase of the solar cycle. Accuracy of the calculated rotation profile indicates that it is possible to measure changes in the profile as the solar cycle evolves. Analysis of further SDO/AIA CBP data will also help in better understanding of the temporal behaviour of the rotation velocity residuals, meridional motions and Reynolds stress. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Journal ref: A&A 587, A29 (2016)

Predicting Coronal Mass Ejections transit times to Earth with neural network

Authors: D. Sudar, B. Vršnak, M. Dumbović

Abstract: Predicting transit times of Coronal Mass Ejections (CMEs) from their initial parameters is a very important subject, not only from the scientific perspective, but also because CMEs represent a hazard for human technology. We used a neural network to analyse transit times for 153 events with only two input parameters: initial velocity of the CME, $v$, and Central Meridian Distance, CMD, of its asso… ▽ More Predicting transit times of Coronal Mass Ejections (CMEs) from their initial parameters is a very important subject, not only from the scientific perspective, but also because CMEs represent a hazard for human technology. We used a neural network to analyse transit times for 153 events with only two input parameters: initial velocity of the CME, $v$, and Central Meridian Distance, CMD, of its associated flare. We found that transit time dependence on $v$ is showing a typical drag-like pattern in the solar wind. The results show that the speed at which acceleration by drag changes to deceleration is $v\approx$500 km s$^{-1}$. Transit times are also found to be shorter for CMEs associated with flares on the western hemisphere than those originating on the eastern side of the Sun. We attribute this difference to the eastward deflection of CMEs on their path to 1 AU. The average error of the NN prediction in comparison to observations is $\approx$12 hours which is comparable to other studies on the same subject. △ Less

Submitted 24 November, 2015; originally announced November 2015.

Steps toward a high precision solar rotation profile: Results from SDO/AIA coronal bright point data

Authors: Davor Sudar, Ivica Skokić, Roman Brajša, Steven H. Saar

Abstract: Coronal bright points (CBP) are ubiquitous small brightenings in the solar corona associated with small magnetic bipoles. We derive the solar differential rotation profile by tracing the motions of CBPs detected by the Atmospheric Imaging Assembly (AIA) instrument aboard the Solar Dynamics Observatory (SDO). We also investigate problems related to detection of coronal bright points resulting from… ▽ More Coronal bright points (CBP) are ubiquitous small brightenings in the solar corona associated with small magnetic bipoles. We derive the solar differential rotation profile by tracing the motions of CBPs detected by the Atmospheric Imaging Assembly (AIA) instrument aboard the Solar Dynamics Observatory (SDO). We also investigate problems related to detection of coronal bright points resulting from instrument and detection algorithm limitations. To determine the positions and identification of coronal bright points we used a segmentation algorithm. A linear fit of their central meridian distance and latitude versus time was utilised to derive velocities. We obtained 906 velocity measurements in a time interval of only 2 days. The differential rotation profile can be expressed as $ω_{rot} = (14.47\pm 0.10 + (0.6\pm 1.0)\sin^{2}(b) + (-4.7\pm 1.7)\sin^{4}(b))$\degr day$^{-1}$. Our result is in agreement with other work and it comes with reasonable errors in spite of the very short time interval used. This was made possible by the higher sensitivity and resolution of the AIA instrument compared to similar equipment as well as high cadence. The segmentation algorithm also played a crucial role by detecting so many CBPs, which reduced the errors to a reasonable level. Data and methods presented in this paper show a great potential to obtain very accurate velocity profiles, both for rotation and meridional motion and, consequently, Reynolds stresses. The amount of coronal bright point data that could be obtained from this instrument should also provide a great opportunity to study changes of velocity patterns with a temporal resolution of only a few months. Other possibilities are studies of evolution of CBPs and proper motions of magnetic elements on the Sun. △ Less

Submitted 6 January, 2015; originally announced January 2015.

Journal ref: A&A 575, A63 (2015)

Geoeffectiveness of Coronal Mass Ejections in the SOHO era

Authors: Mateja Dumbovic, Andy Devos, Bojan Vrsnak, Davor Sudar, Luciano Rodriguez, Domagoj Ruzdjak, Kristoffer Leer, Susanne Vennerstrom, Astrid Veronig

Abstract: The main objective of the study is to determine the probability distributions of the geomagnetic Dst index as a function of the coronal mass ejection (CME) and solar flare parameters for the purpose of establishing a probabilistic forecast tool for the geomagnetic storm intensity. Several CME and flare parameters as well as the effect of successive-CME occurrence in changing the probability for a… ▽ More The main objective of the study is to determine the probability distributions of the geomagnetic Dst index as a function of the coronal mass ejection (CME) and solar flare parameters for the purpose of establishing a probabilistic forecast tool for the geomagnetic storm intensity. Several CME and flare parameters as well as the effect of successive-CME occurrence in changing the probability for a certain range of Dst index values, were examined. The results confirm some of already known relationships between remotely-observed properties of solar eruptive events and geomagnetic storms, namely the importance of initial CME speed, apparent width, source position, and the associated solar flare class. In this paper we quantify these relationships in a form to be used for space weather forecasting in future. The results of the statistical study are employed to construct an empirical statistical model for predicting the probability of the geomagnetic storm intensity based on remote solar observations of CMEs and flares. △ Less

Submitted 13 October, 2014; originally announced October 2014.

Revised physical elements of the astrophysically important O9.5+O9.5V eclipsing binary system Y Cyg

Authors: P. Harmanec, D. E. Holmgren, M. Wolf, H. Božić, E. F. Guinan, Y. W. Kang, P. Mayer, G. P. McCook, J. Nemravová, S. Yang, M. Šlechta, D. Ruždjak, D. Sudar, P. Svoboda

Abstract: Thanks to its long and rich observational history and rapid apsidal motion, the massive eclipsing binary Y Cyg represents one of the cornestones to critical tests of stellar evolution theory for massive stars. Yet, the determination of the basic physical properties is less accurate than it could be given the existing number of spectral and photometric observations. Our goal is to analyze all these… ▽ More Thanks to its long and rich observational history and rapid apsidal motion, the massive eclipsing binary Y Cyg represents one of the cornestones to critical tests of stellar evolution theory for massive stars. Yet, the determination of the basic physical properties is less accurate than it could be given the existing number of spectral and photometric observations. Our goal is to analyze all these data simultaneously with the new dedicated series of our own spectral and photometric observations from observatories widely separated in longitude. We obtained new series of UBV observations at three observatories separated in local time to obtain complete light curves of Y Cyg for its orbital period close to 3 days. This new photometry was reduced and carefully transformed to the standard UBV system using the HEC22 program. We also obtained new series of red spectra secured at two observatories and re-analyzed earlier obtained blue electronic spectra. Our analyses provide the most accurate so far published value of the apsidal period of 47.805 +/- 0.030 yrs and the following physical elements: M1=17.72+/-0.35$ Msun, M2=17.73+/-0.30 Msun, R1=5.785+/-0.091 Rsun, and R2=5.816+/-0.063 Rsun. The disentangling thus resulted in the masses, which are somewhat higher than all previous determinations and virtually the same for both stars, while the light curve implies a slighly higher radius and luminosity for star 2. The above empirical values imply the logarithm of the internal structure constant log k2 = -1.937. A comparison with Claret's stellar interior models implies an age close to 2 millions yrs for both stars. The claimed accuracy of modern element determination of 1-2 per cent seems still a bit too optimistic and obtaining new high-dispersion and high-resolution spectra is desirable. △ Less

Submitted 4 March, 2014; originally announced March 2014.

Comments: 13 pages; accepted for publication in Astronomy and Astrophysics

Tracing sunspot groups to determine angular momentum transfer on the Sun

Authors: D. Sudar, I. Skokić, D. Ruždjak, R. Brajša, H. Wöhl

Abstract: The goal of this paper is to investigate Reynolds stresses and to check if it is plausible that they are responsible for angular momentum transfer toward the solar equator. We also analysed meridional velocity, rotation velocity residuals and correlation between the velocities. We used sunspot groups position measurements from GPR (Greenwich Photographic Result) and SOON/USAF/NOAA (Solar Observing… ▽ More The goal of this paper is to investigate Reynolds stresses and to check if it is plausible that they are responsible for angular momentum transfer toward the solar equator. We also analysed meridional velocity, rotation velocity residuals and correlation between the velocities. We used sunspot groups position measurements from GPR (Greenwich Photographic Result) and SOON/USAF/NOAA (Solar Observing Optical Network/United States Air Force/National Oceanic and Atmospheric Administration) databases covering the period from 1878 until 2011. In order to calculate velocities we used daily motion of sunspot groups. The sample was also limited to $\pm$58\degr in Central Meridian Distance (CMD) in order to avoid solar limb effects. We mainly investigated velocity patterns depending on solar cycle phase and latitude. We found that meridional motion of sunspot groups is toward the centre of activity from all available latitudes and in all phases of the solar cycle. The range of meridional velocities is $\pm10$ m s$^{-1}$. Horizontal Reynolds stress is negative at all available latitudes and indicates that there is a minimum value ($q\approx$ - 3000 m$^2$ s$^{-2}$) located at $b\approx\pm$30$^{0}$. In our convention this means that angular momentum is transported toward the solar equator in agreement with the observed rotational profile of the Sun. △ Less

Submitted 23 January, 2014; v1 submitted 22 January, 2014; originally announced January 2014.

Exploring the Variable Sky with LINEAR. III. Classification of Periodic Light Curves

Authors: Lovro Palaversa, Željko Ivezić, Laurent Eyer, Domagoj Ruždjak, Davor Sudar, Mario Galin, Andrea Kroflin, Martina Mesarić, Petra Munk, Dijana Vrbanec, Hrvoje Božić, Sarah Loebman, Branimir Sesar, Lorenzo Rimoldini, Nicholas Hunt-Walker, Jacob VanderPlas, David Westman, J. Scott Stuart, Andrew C. Becker, Gregor Srdoč, Przemyslaw Wozniak, Hakeem Oluseyi

Abstract: We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0… ▽ More We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0.03 mag at $r=15$ to ~0.20 mag at r=18. Light curves include on average 250 data points, collected over about a decade. Using SDSS-based photometric recalibration of the LINEAR data for about 25 million objects, we selected ~200,000 most probable candidate variables and visually confirmed and classified approximately 7,000 periodic variables using phased light curves. The reliability and uniformity of visual classification across eight human classifiers was calibrated and tested using a SDSS Stripe 82 region variable star catalog, and verified using an unsupervised machine learning approach. The resulting sample of periodic LINEAR variables is dominated by 3,900 RR Lyrae stars and 2,700 eclipsing binary stars of all subtypes, and includes small fractions of relatively rare populations such as asymptotic giant branch stars and SX Phoenicis stars. We discuss the distribution of these mostly uncataloged variables in various diagrams constructed with optical-to-infrared SDSS, 2MASS and WISE photometry, and with LINEAR light curve features. An interesting side result is a robust and precise quantitative description of a strong correlation between the light-curve period and color/spectral type for close and contact eclipsing binary stars. These large samples of robustly classified variable stars will enable detailed statistical studies of the Galactic structure and physics of binary and other stars, and we make them publicly available. △ Less

Submitted 1 August, 2013; originally announced August 2013.

Comments: 71 pages, 21 figures. For additional machine readable tables and light curves, see http://www.astro.washington.edu/users/ivezic/linear/PaperIII/PLV.html . Accepted by AJ

The DWARF project: Eclipsing binaries - precise clocks to discover exoplanets

Authors: T. Pribulla, M. Vaňko, M. Ammler - von Eiff, M. Andreev, A. Aslantürk, N. Awadalla, D. Baluďanský, A. Bonanno, H. Božić, G. Catanzaro, L. Çelik, P. E. Christopoulou, E. Covino, F. Cusano, D. Dimitrov, P. Dubovský, E. M. Esmer, A. Frasca, Ľ. Hambálek, M. Hanna, A. Hanslmeier, B. Kalomeni, D. P. Kjurkchieva, V. Krushevska, I. Kudzej , et al. (31 additional authors not shown)

Abstract: We present a new observational campaign, DWARF, aimed at detection of circumbinary extrasolar planets using the timing of the minima of low-mass eclipsing binaries. The observations will be performed within an extensive network of relatively small to medium-size telescopes with apertures of ~20-200 cm. The starting sample of the objects to be monitored contains (i) low-mass eclipsing binaries with… ▽ More We present a new observational campaign, DWARF, aimed at detection of circumbinary extrasolar planets using the timing of the minima of low-mass eclipsing binaries. The observations will be performed within an extensive network of relatively small to medium-size telescopes with apertures of ~20-200 cm. The starting sample of the objects to be monitored contains (i) low-mass eclipsing binaries with M and K components, (ii) short-period binaries with sdB or sdO component, and (iii) post-common-envelope systems containing a WD, which enable to determine minima with high precision. Since the amplitude of the timing signal increases with the orbital period of an invisible third component, the timescale of project is long, at least 5-10 years. The paper gives simple formulas to estimate suitability of individual eclipsing binaries for the circumbinary planet detection. Intrinsic variability of the binaries (photospheric spots, flares, pulsation etc.) limiting the accuracy of the minima timing is also discussed. The manuscript also describes the best observing strategy and methods to detect cyclic timing variability in the minima times indicating presence of circumbinary planets. First test observation of the selected targets are presented. △ Less

Submitted 28 June, 2012; originally announced June 2012.

Comments: 12 pages, 2 figures, submitted to Astron. Nachrichten

V2368 Oph: An eclipsing and double-lined spectroscopic binary used as a photometric comparison star for U Oph

Authors: P. Harmanec, H. Božić, P. Mayer, P. Eenens, M. Brož, M. Wolf, S. Yang, M. Šlechta, D. Ruždjak, D. Sudar, H. Ak

Abstract: The A-type star HR 6412 = V2368 Oph was used by several investigators as a photometric comparison star for the known eclipsing binary U Oph but was found to be variable by three independent groups, including us. By analysing series of new spectral and photometric observations and a critical compilation of available radial velocities, we were able to find the correct period of light and radial-velo… ▽ More The A-type star HR 6412 = V2368 Oph was used by several investigators as a photometric comparison star for the known eclipsing binary U Oph but was found to be variable by three independent groups, including us. By analysing series of new spectral and photometric observations and a critical compilation of available radial velocities, we were able to find the correct period of light and radial-velocity variations and demonstrate that the object is an eclipsing and double-lined spectroscopic binary moving in a highly eccentric orbit. We derived a linear ephemeris T min.I = HJD (2454294.67 +/- 0.01) + (38.32712 +/- 0.00004)d x E and estimated preliminary basic physical properties of the binary. The dereddened UBV magnitudes and effective temperatures of the primary and secondary, based on our light- and velocity-curve solutions, led to distance estimates that agree with the Hipparcos distance within the errors. We find that the mass ratio must be close to one, but the limited number and wavelength range of our current spectra does not allow a truly precise determination of the binary masses. Nevertheless, our results show convincingly that both binary components are evolved away from the main sequence, which makes this system astrophysically very important. There are only a few similarly evolved A-type stars among known eclipsing binaries. Future systematic observations and careful analyses can provide very stringent tests for the stellar evolutionary theory. △ Less

Submitted 11 May, 2011; originally announced May 2011.

Comments: 10 pages, 7 figs, in press 2011 A&A

Spectral and photometric analysis of the eclipsing binary epsilon Aurigae prior to and during the 2009-2011 eclipse

Authors: Pavel Chadima, Petr Harmanec, Philip D. Bennett, Brian Kloppenborg, Robert Stencel, Stevenson Yang, Hrvoje Bozic, Miroslav Slechta, Lenka Kotkova, Marek Wolf, Petr Skoda, Viktor Votruba, Jeff L. Hopkins, Christian Buil, Davor Sudar

Abstract: A series of 353 red electronic spectra obtained between 1994 and 2010, and of 171 UBV photometric observations of the 2010 eclipse, were analyzed in an effort to better understand the eclipsing binary eps Aur. The main results follow. (1) We attempted to recover a spectrum of the companion by disentangling the observed spectra of the eps Aur binary failed, but we were able to disentangle the spect… ▽ More A series of 353 red electronic spectra obtained between 1994 and 2010, and of 171 UBV photometric observations of the 2010 eclipse, were analyzed in an effort to better understand the eclipsing binary eps Aur. The main results follow. (1) We attempted to recover a spectrum of the companion by disentangling the observed spectra of the eps Aur binary failed, but we were able to disentangle the spectrum of telluric lines and obtain a mean spectrum of the F-type primary star. The latter was then compared to a grid of synthetic spectra for a number of plausible values of T(eff) and log(g), but a reasonably good match was not found. However, we conclude that the observed spectrum is that of a low-gravity star. (2) We examined changes in the complex H-alpha line profiles over the past 16 years, with particular emphasis on the 2009-2011 eclipse period, by subtracting a mean out-of-eclipse H-alpha profile (appropriately shifted in radial velocity) from the observed spectra. We find that the dark disk around the unseen companion has an extended 'atmosphere' that manifests itself via blueshifted and redshifted H-alpha 'shell' absorptions seen projected against the F star. Significantly, the H-alpha shell line first appeared three years before first contact of the optical eclipse when the system was not far past maximum separation. (3) Analyses of radial velocities and central intensities of several strong, unblended spectral lines, as well as UBV photometry, demonstrated that these observables showed apparent multiperiodic variability during eclipse. The dominant period of 66.21 was common to all the observables, but with different phase shifts between these variables. This result strongly supports our earlier suggestion that the photometric variability seen during eclipse is intrinsic to the F star, and therefore, the idea of a central brightening due to a hole in the disk should be abandoned. △ Less

Submitted 30 April, 2011; originally announced May 2011.

Comments: accepted for publication in A&A (April 5, 2011)

Journal ref: 2011A&A...530A.146C

UX Monocerotis as a W Serpentis binary

Authors: D. Sudar, P. Harmanec, H. Lehmann, S. Yang, H. Božić, D. Ruždjak

Abstract: Using our new photometric and spectroscopic observations as well as all available published data, we present a new interpretation of the properties of the peculiar emission-line binary UX Mon. We conclude that this binary is in a rare phase of fast mass transfer between the binary components prior to the mass ratio reversal. We firmly establish that the orbital period is secularly decreasing at a… ▽ More Using our new photometric and spectroscopic observations as well as all available published data, we present a new interpretation of the properties of the peculiar emission-line binary UX Mon. We conclude that this binary is in a rare phase of fast mass transfer between the binary components prior to the mass ratio reversal. We firmly establish that the orbital period is secularly decreasing at a rate of $\dot{P}=-0.260$ seconds per year. From several lines of reasoning, we show that the mass ratio of the component losing mass to the mass-gaining component $q$ must be larger than 1 and find our most probable value to be $q=1.15\pm0.1$. The BINSYN suite of programs and the steepest descent method were used to perform the final modeling. We modeled the star as a W Ser star with a thick disk around its primary. Although the remaining uncertainties in some of the basic physical elements describing the system in our model are not negligible, the model is in fair agreement with available observations. Only the nature of the light variations outside the primary eclipse remains unexplained. △ Less

Submitted 28 June, 2013; v1 submitted 9 March, 2011; originally announced March 2011.

A new ephemeris and an orbital solution of ε Aurigae

Authors: Pavel Chadima, Petr Harmanec, Stephenson Yang, Phillip D. Bennett, Hrvoje Božić, Domagoj Ruždjak, Davor Sudar, Petr Škoda, Miroslav Šlechta, Marek Wolf, Martin Lehký, Pavol Dubovský

Abstract: We collected rich series of RV measurements covering last 110 years and photometric observations from the past 6 primary eclipses, complemented them by our new observations and derived a new precise ephemeris and an orbital solution of epsilon Aur. We collected rich series of RV measurements covering last 110 years and photometric observations from the past 6 primary eclipses, complemented them by our new observations and derived a new precise ephemeris and an orbital solution of epsilon Aur. △ Less

Submitted 3 May, 2010; originally announced May 2010.

Comments: accepted and published in IBVS (paper No. 5937)

Journal ref: 2010IBVS.5937....1C