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.

An online software tool for the easy preparation of ephemerides of the solar surface features is presented. It was developed as a helper tool for the preparation of observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), but it can be used at other observatories as well. The tool features an easy to use point-and-click graphical user interface with the possibility to enter or adjust input parameters, while the result is a table of predicted positions in the celestial equatorial coordinate system, suitable for import into the ALMA Observing Tool software. The tool has been successfully used for the preparation and execution of solar observations with ALMA.

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.

Polar brightening of the Sun at radio frequencies has been studied for almost fifty years and yet a disagreement persists between solar atmospheric models and observations. Some observations reported brightening values much smaller than the expected values obtained from the models, with discrepancies being particularly large at millimeter wavelengths. New clues to calibrate the atmospheric models can be obtained with the advent of the Atacama Large Millimeter/submillimeter Array (ALMA) radio interferometer. In this work, we analyzed the lower limit of the polar brightening observed at 100 and 230 GHz by ALMA, during its Science Verification period, 2015 December 16-20. We find that the average polar intensity is higher than the disk intensity at 100 and 230 GHz, with larger brightness intensities at the South pole in eight of the nine maps analyzed. The observational results were compared with calculations of the millimetric limb brightnening emission for two semi-empirical atmospheric models, FAL- C (Fontenla et al. 1993) and SSC (Selhorst et al. 2005a). Both models presented larger limb intensities than the averaged observed values. The intensities obtained with the SSC model were closer to the observations, with polar brightenings of 10.5% and 17.8% at 100 and 230 GHz, respectively. This discrepancy may be due to the presence of chromospheric features (like spicules) at regions close to the limb.

Sunspot position data obtained from Kanzelh\"{o}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.

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\alpha, 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.

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.

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\"ohe 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 $\omega(b)$ = (14.47 $\pm$ 0.01) - (2.66 $\pm$ 0.10) $\sin^2b$ (deg/day) for the DS method and $\omega(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.

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.

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 $\omega(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.

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 $\omega_{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.

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.

Existing methods for conversion between synodic and sidereal rotation velocities of the Sun are tested for validity using state of the art ephemeris data. It is found that some of them are in good agreement with ephemeris calculations, while the other ones show a discrepancy of almost 0.01 deg per day. The discrepancy is attributed to a missing factor and a new corrected relation is given.