Solar modulation of galactic cosmic rays during 19–23 solar cycles

We make a detailed analysis of cross-correlation and time-lag between monthly data of galactic cosmic rays (GCRs) intensity and different solar activity indices (e.g., sunspot number, sunspot area, green coronal Fe line and 10.7 cm solar radio flux) during 19–23 solar cycles. GCRs time-series data from Kiel neutron monitor station and solar data from the last 50 years period, covering five solar cycles (19–23), and alternating solar polarity states (i.e., five A < 0 and four A > 0) have been investigated. We find a clear asymmetry in the cross-correlation between GCRs and solar activity indicators for both odd and even-numbered solar cycles. The time-lags between GCRs and solar parameters are found different in different solar cycles as well as in the opposite polarity states (A < 0 and A > 0) within the same solar cycle. Possible explanations of the observed results are discussed in light of modulation models, including drift effects.     

The orbital and superhump periods of the dwarf nova HS 0417+7445 in Camelopardalis

We present the 2005–2010 outburst history of the SU UMa-type dwarf HS 0417+7445, along with a detailed analysis of extensive time-series photometry obtained in March 2008 during the second recorded superoutburst of the system. The mean outburst interval is 197 ± 59 d, with a median of 193 d. The March 2008 superoutburst was preceded by a precursor outburst, had an amplitude of 4.2 magnitudes, and the whole event lasted about 16 days. No superhumps were detected during the decline from the precursor outburst, and our data suggests instead that orbital humps were present during that phase. Early superhumps detected during the rise to the superoutburst maximum exhibited an unusually large fractional period excess of ? = 0.137 (P_sh = 0.0856(88) d). Following the maximum, a linear decline in brightness followed, lasting at least 6 days. During this decline, a stable superhump period of P_sh = 0.07824(2) d was measured. Superimposed on the superhumps were orbital humps, which allowed us to accurately measure the orbital period of HS 0417+7445, P_orb = 0.07531(8) d, which was previously only poorly estimated. The fractional superhump period excess during the main phase of the outburst was ? = 0.037, which is typical for SU UMa dwarf novae with similar orbital period. Our observations are consistent with the predictions of the thermal-tidal instability model for the onset of superoutbursts, but a larger number of superoutbursts with extensive time-series photometry during the early phases of the outburst would be needed to reach a definite conclusion on the cause of superoutbursts.     

Hemispheric asymmetry of sunspot area in solar cycle 23 and rising phase of solar cycle 24: Comparison of three data sets

Analysis of long-term solar data from different observatories is required to compare and confirm the various level of solar activity in depth. In this paper, we study the north–south asymmetry of monthly mean sunspot area distribution during the cycle-23 and rising phase of cycle-24 using the data from Kodaikanal Observatory (KO), Michelson Doppler Imager (MDI) and Solar Optical Observing Network (SOON). Our analysis confirmed the double peak behavior of solar cycle-23 and the dominance of southern hemisphere in all the sunspot area data obtained from three different resources. The analysis also showed that there is a 5–6 months time delay in the activity levels of two hemispheres. Furthermore, the wavelet analysis carried on the same data sets showed several known periodicities (e.g., 170–180 days, 2.1 year) in the north–south difference of sunspot area data. The temporal occurrence of these periods is also the same in all the three data sets. These results could help in understanding the underlying mechanism of north–south asymmetry of solar activity.     

Scaling relation between Sunyaev–Zel’dovich effect and X-ray luminosity and scale-free evolution of cosmic baryon field

It has been revealed recently that, in the scale free range, i.e. from the scale of the onset of nonlinear evolution to the scale of dissipation, the velocity and mass density fields of cosmic baryon fluid are extremely well described by the self-similar log-Poisson hierarchy. As a consequence of this evolution, the relations among various physical quantities of cosmic baryon fluid should be scale invariant, if the physical quantities are measured in cells on scales larger than the dissipation scale, regardless the baryon fluid is in virialized dark halo, or in pre-virialized state. We examine this property with the relation between the Compton parameter of the thermal Sunyaev–Zel’dovich effect, y(r), and X-ray luminosity, L_x(r), where r being the scale of regions in which y and L_x are measured. According to the self-similar hierarchical scenario of nonlinear evolution, one should expect that (1) in the y(r) ? L_x(r) relation, y(r) = 10^A(r)[L_x(r)]^α(r), the coefficients A(r) and α(r) are scale-invariant; (2) The relation y(r) = 10^A(r)[L_x(r)]^α(r) given by cells containing collapsed objects is also available for cells without collapsed objects, only if r is larger than the dissipation scale. These two predictions are well established with a scale decomposition analysis of observed data, and a comparison of observed y(r) ? L_x(r) relation with hydrodynamic simulation samples. The implication of this result on the characteristic scales of non-gravitational heating is also addressed.     

Gleissberg cycle of solar activity over the last 7000 years

Long-term variations of solar activity significantly affect terrestrial phenomena. Studies have shown cyclic components in solar activity and geophysical phenomena (e.g., the Schwabe, Hale, Gleissberg, and Suess cycles, and a cycle of about 2300 years). In this paper, the wavelet technique is employed to investigate the Gleissberg cycle in solar variations during 5000 BC–1995 AD. Analysis shows time-variable characteristics in the Gleissberg solar cycle over the period; no obvious correlation between the Gleissberg and Suess cycles has been found.     

The near-contact binary star RZ Dra revisited

This paper presents the absolute parameters of RZ Dra. New CCD observations were made at the Mt. Suhora Observatory in 2007. Two photometric data sets (1990 BV and 2007 BVRI) were analysed using modern light-curve synthesis methods. Large asymmetries in the light curves may be explained in terms of a dark starspot on the primary component, an A6 type star. Due to this magnetic activity, the primary component would appear to belong to the class of Ap-stars and would show small amplitude with δ Scuti-type pulsations. With this in mind, a time-series analysis of the residual light curves was made. However, we found no evidence of pulsation behaviour in RZ Dra. Combining the solutions of our light curves and Rucinski et al. (2000)’s radial velocity curves, the following absolute parameters of the components were determined: M₁ = 1.63 ± 0.03 M_⊙, M₂ = 0.70 ± 0.02 M_⊙, R₁ = 1.65 ± 0.02R_⊙, R₂ = 1.15 ± 0.02 R_⊙, L₁ = 9.72 ± 0.30 L_⊙ and L₂ = 0.74 ± 0.10 L_⊙. The distance to RZ Dra was calculated as 400 ± 25 pc, taking into account interstellar extinction. The orbital period of the system was studied using updated O–C information. It was found that the orbital period varied in its long-period sinusoidal form, superimposed on a downward parabola. The parabolic term shows a secular period decrease at a slow rate of 0.06 ± 0.02 s per century and is explained by the mass loss via magnetized wind of the Ap-star primary. The tilted sinusoidal form of the period variation may be considered as an apparent change and may be interpreted in terms of the light-time effect due to the presence of a third body.     

Identification of Zirconium Oxide molecular lines in the sunspot umbral spectra

The ZrO molecule has been detected in sunspot umbrae through the identification of following laboratory molecular transitions: ¹Σ⁺ ? X¹Σ⁺ (0, 0), A³Φ₂ ? X²Δ₁ (0, 0), A³Φ₃ ? X²Δ₂ (0, 0), A³Φ₄ ? X²Δ₃ (0, 0), B³Π₂ ? X³Δ₃ (0, 0), B³Π₁ ? X³Δ₂ (0, 0) and B³Π₀ ? X³Δ₁ (0, 0) in red – infrared region using high resolution, visible range Fourier Transform Spectrum of sunspot umbra observed at the National Solar Observatory in Kitt Peak (NSO/KP). Much new identification has been made in the searched spectral wavenumber region from 16650 cm^?1 to 18007 cm^?1 of sunspot spectrum. Equivalent widths of well resolved lines, versus rotational quantum number J have been used to determine the effective rotational temperature for seven bands of the ZrO molecule. This result agrees well with the temperatures derived for other molecules’ presence in sunspot umbrae. It is evident that ZrO molecular lines are formed in higher layers of the atmosphere of relatively “cold” sunspots.     

Spectroscopic and photometric study of the contact binary BO CVn

We present the results of the study of the contact binary system BO CVn. We have obtained physical parameters of the components based on combined analysis of new, multi-color light curves and spectroscopic mass ratio. This is the first time the latter has been determined for this object. We derived the contact configuration for the system with a very high filling factor of about 88%. We were able to reproduce the observed light curve, namely the flat bottom of the secondary minimum, only if a third light has been added into the list of free parameters. The resulting third light contribution is significant, about 20–24%, while the absolute parameters of components are: M₁ = 1.16, M₂ = 0.39, R₁ = 1.62 and R₂ = 1.00 (in solar units).The O-C diagram shows an upward parabola which, under the conservative mass transfer assumption, would correspond to a mass transfer rate of dM/dt = 6.3 × 10^?8M_⊙/yr, matter being transferred from the less massive component to the more massive one. No cyclic, short-period variations have been found in the O-C diagram (but longer-term variations remain a possibility).     

Galactic longitude dependent galactic model parameters

We present the Galactic model parameters for thin disc estimated by Sloan Digital Sky Survey (SDSS) data of 14 940 stars with apparent magnitudes 16 < g₀ ⩽ 21 in six intermediate latitude fields in the first Galactic quadrant. Star/galaxy separation was performed by using the SDSS photometric pipeline and the isodensity contours in the (g − r)₀ − (r − i)₀ two colour diagram. The separation of thin disc stars is carried out by the bimodal distribution of stars in the (g − r)₀ histogram, and the absolute magnitudes were evaluated by a procedure presented in the literature (Bilir, S., Karaali, S., Tunçel, S. 2005. AN 326, 321). Exponential density law fits better to the derived density functions for the absolute magnitude intervals 8 < M(g) ⩽ 9 and 11 < M(g) ⩽ 12, whereas sech/sech² laws are more appropriate for absolute magnitude intervals 9 < M(g) ⩽ 10 and 10 < M(g) ⩽ 11. We showed that the scaleheight and scalelength are Galactic longitude dependent. The average values and ranges of the scaleheight and the scalelength are 〈H〉 = 220 pc (196 ⩽ H ⩽ 234 pc) and 〈H〉 = 1900 pc (1561 ⩽ h ⩽ 2280 pc) respectively. This result would be useful to explain different numerical values claimed for those parameters obtained by different authors for the fields in different directions of the Galaxy.     

An Estimate for the Size of Sunspot Cycle 24

For the sunspot cycles in the modern era (cycle?10 to the present), the ratio of R _Z(max)/R _Z(36th month) equals 1.26±0.22, where R _Z(max) is the maximum amplitude of the sunspot cycle?using smoothed monthly mean sunspot number and R _Z(36th month) is the smoothed monthly mean sunspot number 36 months after cycle?minimum. For the current sunspot cycle?24, the 36th month following the cycle?minimum occurred in November 2011, measuring?61.1. Hence, cycle?24 likely will have a maximum amplitude of about 77.0±13.4 (the one-sigma prediction interval), a value well below the average R _Z(max) for the modern era sunspot cycles (about 119.7±39.5).     

The detached eclipsing binary TX Her revisited

This paper presents new CCD Bessell BVRI light curves and photometric analysis of the Algol-type binary star TX Her. The CCD observations were carried out at Çanakkale Onsekiz Mart University Observatory in 2010. New BVRI light curves from this study and radial velocity curves from Popper (1970) were solved simultaneously using modern light and radial velocity curves synthesis methods. The general results show that TX Her is a well-detached eclipsing binary, however, both component stars fill at least half of their Roche lobes. A significant third light contribution to the total light of the system could not be determined. Using O–C residuals formed by the updated minima times, an orbital period study of the system was performed. It was confirmed that the tilted sinusoidal O–C variation corresponds to an apparent period variation caused by the light travel time effect due to an unseen third body. The following absolute parameters of the components were derived: M₁ = 1.62 ± 0.04 M_⊙, M₂ = 1.45 ± 0.03 M_⊙, R₁ = 1.69 ± 0.03 R_⊙, R₂ = 1.43 ± 0.03 R_⊙, L₁ = 8.21 ± 0.90 L_⊙ and L₂ = 3.64 ± 0.60 L_⊙. The distance to TX Her was calculated as 155 ± 10 pc, taking into account interstellar extinction. The position of the components of TX Her in the HR diagram are also discussed. The components are young stars with an age of ～500 Myr.     

What the Sunspot Record Tells Us About Space Climate

The records concerning the number, sizes, and positions of sunspots provide a direct means of characterizing solar activity over nearly 400 years. Sunspot numbers are strongly correlated with modern measures of solar activity including: 10.7-cm radio flux, total irradiance, X-ray flares, sunspot area, the baseline level of geomagnetic activity, and the flux of galactic cosmic rays. The Group Sunspot Number provides information on 27 sunspot cycles, far more than any of the modern measures of solar activity, and enough to provide important details about long-term variations in solar activity or “Space Climate.” The sunspot record shows: 1) sunspot cycles have periods of 131± 14 months with a normal distribution; 2) sunspot cycles are asymmetric with a fast rise and slow decline; 3) the rise time from minimum to maximum decreases with cycle amplitude; 4) large amplitude cycles are preceded by short period cycles; 5) large amplitude cycles are preceded by high minima; 6) although the two hemispheres remain linked in phase, there are significant asymmetries in the activity in each hemisphere; 7) the rate at which the active latitudes drift toward the equator is anti-correlated with the cycle period; 8) the rate at which the active latitudes drift toward the equator is positively correlated with the amplitude of the cycle after the next; 9) there has been a significant secular increase in the amplitudes of the sunspot cycles since the end of the Maunder Minimum (1715); and 10) there is weak evidence for a quasi-periodic variation in the sunspot cycle amplitudes with a period of about 90 years. These characteristics indicate that the next solar cycle should have a maximum smoothed sunspot number of about 145 ± 30 in 2010 while the following cycle should have a maximum of about 70 ± 30 in 2023.     

Measurement of the radial density gradient of cosmic ray in the heliosphere by the GRAPES-3 experiment

A radial anisotropy in the flux of cosmic rays in heliosphere was theoretically predicted by Parker and others within the framework of the diffusion–convection mechanism. The solar wind is responsible for sweeping out the galactic cosmic rays, creating a radial density gradient within the heliosphere. This gradient coupled with the interplanetary magnetic field induces a flow of charged particles perpendicular to the ecliptic plane which was measured and correctly explained by Swinson, and is hereafter referred as ‘Swinson flow’. The large area GRAPES-3 tracking muon telescope offers a powerful probe to measure the Swinson flow and the underlying radial density gradient of the galactic cosmic rays at a relatively high rigidity of ∼100 GV. The GRAPES-3 data collected over a period of six years (2000–2005) were analyzed and the amplitude of the Swinson flow was estimated to be (0.0644 ± 0.0008)% of cosmic ray flux which was an ∼80σ effect. The phase of the maximum flow was at a sidereal time of (17.70 ± 0.05) h which was 18 min earlier than the expected value of 18 h. This small 18 min phase difference had a significance of ∼6σ indicating the inherent precision of the GRAPES-3 measurement. The radial density gradient of the galactic cosmic rays at a median rigidity of 77 GV was found to be 0.65% AU⁻¹.     

Orbital period investigations of two short-period early-type overcontact binaries BH Cen and V701 Sco in two extremely young galactic clusters IC 2944 and NGC 6383

Both V701 Sco and BH Cen are two early-type short-period overcontact systems (P = 0.^d762 and P = 0.^d792, respectively). V701 Sco is a member of the young galactic cluster NGC 6383, while BH Cen is a component of a younger galactic cluster IC 2944 where star formation is in process. They provide good opportunity to understand the formation and evolution of binary stars. In the present paper, orbital period changes of the two binaries are investigated. It is discovered that the orbital period of BH Cen shows a long-term increase with a rate of dP/dt = +1.70(±0.39) × 10⁻⁷ days/year while it undergoes a cyclic oscillation with a period of 44.6 years and an amplitude of A₃ = 0.^d0216. For V701 Sco, its O-C curve reveals a periodic change with a period of 41.2 years and amplitude of A₃ = 0.^d0158. The mass ratio of BH Cen is 0.84, but V701 Sco contains twin B1-1.5V type stars with a mass ratio of unit. The continuous period increase of BH Cen is caused by the mass transfer from the less massive component to the more massive one at a rate of dM₂/dt = 3.5 × 10⁻⁶ days/year.The cyclic period changes of both systems can be plausibly explained as the results of light-travel time effects suggesting that they are triple systems. The astrophysical parameters of the unseen tertiary components in the two systems have been determined. We think that the invisible tertiary components in both binaries played an important role in the formations and evolutions of the overcontact configurations by bringing angular momentum out from the central systems. For BH Cen, this process created the initial short period and will support its evolution into an overcontact configuration via a Case A mass transfer within the life time of the extremely young cluster IC 2944. For V701 Sco, two identical zero-age main-sequence components in an overcontact configuration suggest that it may have been formed by fission, possibly by the fission of the third body. The fact that no long-term continuous period variations were found for V701 Sco may suggest that an overcontact binary with the mass ratio of unity can be in an equilibrium revealing that the original configuration of the binary was overcontact as is its present state. It has been reported that faint stars in the two extremely young clusters are relatively scare. From the present study, it is shown that faint stars in young clusters are usually formed as companions of OB stars (including binaries). It is very difficult to detect them because of their low luminosity when compared with the more luminous OB stars.     

Periodicity of flare index revisited using the Hilbert–Huang transform method

Using the Hilbert–Huang Transform (HHT) method, we investigate the periodicity in the monthly mean flare indices from 1966 to 2007 (corresponding to almost four complete solar cycles), calculated by T. Atac and A. Ozguc. The results show as following. (1) The periods of 9.37 ± 2.53, 11.8 ± 0.172 and 23.6 ± 0.316 years are found to be statistically significant in the flare index. The most eminent period is 9.37 ± 2.53 years. (2) Other periods of 0.237 ± 0.196 years (86.6 ± 71.6 days), 0.525 ± 0.0508 years (191 ± 18.5 days), 1.05 ± 0.478 years (383 ± 174 days) and 2.37 ± 0.395 years are below the 99% confidence level line, suggesting they are due to stochastic random noise.     

Short-Term Periodicities in Sunspot Activity and Flare Index Data during Solar Cycle 23

The short-term periodicities in sunspot numbers, sunspot areas, and flare index data are investigated in detail using the Date Compensated Discrete Fourier Transform (DCDFT) for the full disk of the Sun separately over the rising, the maximum, and the declining portions of solar cycle 23 (1996 – 2006). While sunspot numbers and areas show several significant periodicities in a wide range between 23.1 and 36.4 days, the flare index data do not exhibit any significant periodicity. The earlier conclusion of Pap, Tobiska, and Bouwer (1990, Solar Phys. 129, 165) and Kane (2003, J. Atmos. Solar-Terr. Phys. 65, 1169), that the 27-day periodicity is more pronounced in the declining portion of a solar cycle than in the rising and maximum ones, seems to be true for sunspot numbers and sunspot area data analyzed here during solar cycle 23.     

Latest results on Jovian disk X-rays from XMM-Newton

We present the results of a spectral study of the soft X-ray emission (0.2–2.5 keV) from low-latitude (‘disk’) regions of Jupiter. The data were obtained during two observing campaigns with XMM-Newton in April and November 2003. While the level of the emission remained approximately the same between April and the first half of the November observation, the second part of the latter shows an enhancement by about 40% in the 0.2–2.5 keV flux. A very similar, and apparently correlated increase, in time and scale, was observed in the solar X-ray and EUV flux.The months of October and November 2003 saw a period of particularly intense solar activity, which appears reflected in the behavior of the soft X-rays from Jupiter's disk. The X-ray spectra, from the XMM-Newton EPIC CCD cameras, are all well fitted by a coronal model with temperatures in the range 0.4–0.5 keV, with additional line emission from Mg XI (1.35 keV) and Si XIII (1.86 keV): these are characteristic lines of solar X-ray spectra at maximum activity and during flares.The XMM-Newton observations lend further support to the theory that Jupiter's disk X-ray emission is controlled by the Sun, and may be produced in large part by scattering, elastic and fluorescent, of solar X-rays in the upper atmosphere of the planet.     

The growth of wind-waves in Titan’s hydrocarbon seas

Ocean wave growth on Titan is considered. The classic Sverdrup–Munk theory for terrestrial wave growth is applied to Titan, and is compared with a simple energy balance model that exposes the effect of Titan’s environmental parameters (air density, gravity, and fluid density). These approaches are compared with the only previously-published (semi-empirical) model (Ghafoor, N.A.-L., Zarnecki, J.C., Challenor, P., Srokosz, M.A. [2000] J. Geophys. Res. 105, 12,077–12,091, hereafter G2k), and allow the impact of various parameters such as atmospheric density to be transparently explored.Our model, like G2k, suggests fully-developed significant wave heights on Titan H_s = 0.2 U², where U is the windspeed (SI units): in dimensionless terms this is rather close to H_s = 0.2 U²/g, a rule of thumb previously noted for terrestrial waves (we find various datasets where the prefactor varies by ～2). It is noted that liquid and air densities affect the growth rate of waves, but not their fully-developed height: for 1 m/s winds wave amplitude reaches 0.15 m (75% of fully-developed) with a fetch of only 1 km, rather faster than predicted by G2k. Liquid viscosity has no major effect on gravity wave growth, but does influence the threshold windspeed at which gravity–capillary waves form in the first place.The model is used to develop predicted ranges for wave height to guide the design of the Titan Mare Explorer (TiME), a proposed Discovery-class mission to float a capsule on Ligeia Mare in 2023. For the expected maximum 1 m/s winds, a significant wave height of 0.2 m and wavelength of ～4 m can be expected. Assuming that wave heights follow Rayleigh statistics as they do on Earth, then given the wave period of ～4 s, individual waves of ～0.6 m might be encountered over a 3 month period.For predicted Titan winds at Kraken Mare, significant wave heights may reach ～0.6 m in the peak of summer but do not exceed the tidal amplitude at its northern end, consistent with the area around Mayda Insula being a tidal flat, while elsewhere on Kraken and Ligeia and at Ontario Lacus, shorelines may be wave- or tidally-dominated, depending on the specific location.     

Absolute magnitudes of asteroids and a revision of asteroid albedo estimates from WISE thermal observations

We obtained estimates of the Johnson V absolute magnitudes (H) and slope parameters (G) for 583 main-belt and near-Earth asteroids observed at Ond?ejov and Table Mountain Observatory from 1978 to 2011. Uncertainties of the absolute magnitudes in our sample are <0.21 mag, with a median value of 0.10 mag. We compared the H data with absolute magnitude values given in the MPCORB, Pisa AstDyS and JPL Horizons orbit catalogs. We found that while the catalog absolute magnitudes for large asteroids are relatively good on average, showing only little biases smaller than 0.1 mag, there is a systematic offset of the catalog values for smaller asteroids that becomes prominent in a range of H greater than ～10 and is particularly big above H ～ 12. The mean (H_catalog ? H) value is negative, i.e., the catalog H values are systematically too bright. This systematic negative offset of the catalog values reaches a maximum around H = 14 where the mean (H_catalog ? H) is ?0.4 to ?0.5. We found also smaller correlations of the offset of the catalog H values with taxonomic types and with lightcurve amplitude, up to ～0.1 mag or less. We discuss a few possible observational causes for the observed correlations, but the reason for the large bias of the catalog absolute magnitudes peaking around H = 14 is unknown; we suspect that the problem lies in the magnitude estimates reported by asteroid surveys. With our photometric H and G data, we revised the preliminary WISE albedo estimates made by Masiero et al. (Masired, J.R. et al. [2011]. Astrophys. J. 741, 68–89) and Mainzer et al. (Mainzer, A. et al. [2011b]. Astrophys. J. 743, 156–172) for asteroids in our sample. We found that the mean geometric albedo of Tholen/Bus/DeMeo C/G/B/F/P/D types with sizes of 25–300 km is p_V = 0.057 with the standard deviation (dispersion) of the sample of 0.013 and the mean albedo of S/A/L types with sizes 0.6–200 km is 0.197 with the standard deviation of the sample of 0.051. The standard errors of the mean albedos are 0.002 and 0.006, respectively; systematic observational or modeling errors can predominate over the quoted formal errors. There is apparent only a small, marginally significant difference of 0.031 ± 0.011 between the mean albedos of sub-samples of large and small (divided at diameter 25 km) S/A/L asteroids, with the smaller ones having a higher albedo. The difference will have to be confirmed and explained; we speculate that it may be either a real size dependence of surface properties of S type asteroids or a small size-dependent bias in the data (e.g., a bias towards higher albedos in the optically-selected sample of asteroids). A trend of the mean of the preliminary WISE albedo estimates increasing with asteroid size decreasing from D ～ 30 down to ～5 km (for S types) showed in Mainzer et al. (Mainzer, A. et al. [2011a]. Astrophys. J. 741, 90–114) appears to be mainly due to the systematic bias in the MPCORB absolute magnitudes that progressively increases with H in the corresponding range H = 10–14.