330 likes | 424 Views
Temperature effect of muon component and practical questions of its account in real time for Global muon detector Network. 1 7- 1 9 October 2011 Ottawa. Observations of the CR secondary components Anisotropy - one of the most important Space Weather parameters.
E N D
Temperature effect of muoncomponent and practicalquestions of its account in realtime for Global muon detector Network 17-19 October 2011 Ottawa
Observations of the CR secondary components Anisotropy - one of the most important Space Weather parameters. Muon telescopes are adjusted to the CR anisotropy observation much better than NM because of their narrow direction and specific construction (location at one point several detectors of different directions). In some cases such a telescope may substitute a group of detectors (for example, NMS) located in different points at the globe. But use of muon detector data very oftenrestrains presence at this data of the big temperature effect inherent muonic component of secondary CR. What is the nature of temperature effect? Layer of the muon generation in the atmosphere
Main stages. Meteorological effects were discovered just with the beginning of systematic study of the variations of secondary CR. The nature of barometric effect was determined relatively quickly (in 1926!), but an estimation of temperature effect took several tens. In 1932 P.Blacket, basing on the hypotetic muon (which was experimentally discovered in 1936), explained successfully the negative temperature effect. M.Forro found also positive temperature effect in 1947, which he was happened to explain with two-meson model (in 1947 the existence of peon was experimentally confirmed). Temperature effect was also investigated by V.Hess, 100years Anniversary of whom is coming, but in 1940 he could not explain the experimental results obtained. V. Hess, "On the Seasonal and the Atmospheric Temperature Effect in Cosmic Radiation", Phys. Rev., V. 57, May 1, 1940.
Method of crossed telescopes [Elliot, 1950] . It allows to get rid of variations of an atmospheric origin, having kept almost completely anisotropic variations. However, we cannot define the CR density in this case, or estimate spectrum of the CR variations. The seasonal temperature effect (~5%) is comparable with the CR variations of non terrestrial origin (for example, 11-year variations), and it is impossible to use muon data in a full volume if it uncorrected for temperature effect.
Temperature effect and long term variations. Ionisation chamber in Yakutsk(wide directed detector with lead screen to cut off the soft component of the CR). Ionization chamber In Yakutsk Upper panel – uncorrected and corrected for temperature counting rates. Comparison with the data from NMs Huancayo и Haleakala. From 1995-drief.
Temperature effect and long term variations. Muon telescope NAGOYA, vertical Upper panel – uncorrected and corrected for temperature counting rates. Comparison with the data from NMs Huancayo и Haleakala.
Methods of the Temperature Effect exclusion. 1) Method Duperier (A. Duperier,1949) (empirical) Where αH decay coefficient (%/km) – negative effect, αTpositive temperature coefficient (%/C)(empirical) 2) The integral method (L. Dorman, 1954; Maeda & Wada, 1954; Olbert, 1953) Where δT(h)=TB-T and WT(h) – density of temperature coefficient. 3) The method of the effective temperature (P. Barrett et. al., 1952) 4) The method of mean-mass temperature (Yu. Krestyannikov, 1976)
Density of temperature coefficient We see here Mean-mass temperature For extra high energies P.H.Barrett et. al.,. Rev. Mod. Phys., 24: 133, 1952.L.V.Volkova, Nucl. Phys., 12 (2): 347-359, 1970.L.Dorman, Meteorological Effects of Cosmic Rays, Nauka, 1972.; Dorman & Yanke, 1971.A.Dmitrieva, Astroparticle Physics 34, 401–411, 2011.Gaisser, T., Cosmic Rays and Particle Physics, Cambridge University Press, Cambridge, Chapter 6, 1990.Maeda& Wada, 1954.
Where do we take Temperature data ? We use the result of the Global Forecast System (GFS) temperature model representing by the National Centers for Environmental Prediction — NCEP (USA). On the basis of GFS the system GEFS (Global Ensemble Forecast System) submits temperature data and prognosis at 28 vertical levels every 6 hours (at 00, 06,12, and 18 UT). Data are interpolated on the grid with resolution 1°x1° with 6-hour interval). Getting data of GEFS model is going in real time by means of distributed system of access to the geophysical data bases (http://www.ngdc.noaa.gov/wdc/wdcmain.html). Weather serverhttp://esse.wdcb.ru; Atmosphere temperature profile in real timehttp://phoenix.wdcb.ru, (mirror http://dimm.wdcb.ru)
Temperature Data in Real Time Updating of database for reanalysis is carried up to 5 days. Simultaneously the prognosos is performed for the next 5 days. Its data on the current day T(h) is used for calculation and exclusion of temperature effect in real time.
Temperature Data 2011-08-03 00:00:00+00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 101314.95 12.06 9.50 7.25 5.31 3.682.35 1.32 0.62 0.26 0.26 0.651.45 2.66 4.18 5.93 7.77 9.62 11.35 12.87 14.10 14.97 15.42 15.37 100017.15 16.07 15.06 14.13 13.28 12.52 11.85 11.27 10.80 10.41 10.13 9.949.85 9.86 9.99 10.27 10.70 11.32 12.15 13.17 14.27 15.30 16.11 16.54 925 12.15 11.24 10.46 9.79 9.24 8.80 8.45 8.19 8.00 7.85 7.73 7.60 7.45 7.27 7.09 6.98 6.98 7.15 7.55 8.20 9.00 9.82 10.54 11.03 850 5.75 5.41 5.14 4.92 4.76 4.64 4.55 4.49 4.44 4.38 4.27 4.10 3.85 3.50 3.11 2.72 2.41 2.23 2.25 2.49 2.90 3.38 3.84 4.20 700 4.25 4.59 4.84 5.01 5.10 5.115.05 4.91 4.71 4.45 4.13 3.76 3.35 2.90 2.40 1.84 1.21 0.48 -0.35 -1.29 -2.28 -3.29 -4.26 -5.13 600-3.75 -3.91 -4.09 -4.29 -4.49 -4.71-4.95 -5.20 -5.46 -5.73 -6.00 -6.28-6.55 -6.82 -7.10 -7.40 -7.73 -8.11-8.55 -9.06 -9.62 -10.23 -10.84 -11.46 500-15.25 -14.99 -14.88 -14.90 -15.03 -15.26-15.55-15.89 -16.26 -16.62 -16.94 -17.20-17.35-17.39 -17.32 -17.16 -16.95 -16.71-16.45 -16.20 -15.99 -15.85 -15.79 -15.85 400-27.25 -26.95 -26.82 -26.83 -26.98 -27.22-27.55 -27.93 -28.34 -28.73 -29.07 -29.32-29.45 -29.43 -29.27 -29.04 -28.74 -28.44-28.15 -27.91 -27.73 -27.59 -27.48 -27.41 300-38.45 -38.55 -38.59 -38.58 -38.51 -38.40-38.25 -38.06 -37.84 -37.57 -37.27 -36.93-36.55 -36.13 -35.70 -35.29 -34.92 -34.63-34.45 -34.39 -34.45 -34.58 -34.78 -35.01 250-42.85 -42.92 -42.93 -42.89 -42.81 -42.69 -42.55 -42.39 -42.23 -42.09 -41.98 -41.93-41.95 -42.05 -42.22 -42.41 -42.60 -42.76-42.85 -42.86 -42.81 -42.73 -42.66 -42.62 200-48.65 -48.98 -49.22 -49.36 -49.42 -49.41-49.35 -49.25 -49.13 -49.01 -48.93 -48.90-48.95 -49.09 -49.31 -49.55 -49.77 -49.96-50.05 -50.03 -49.94 -49.81 -49.68 -49.62 150-57.25 -57.69 -58.04 -58.33 -58.56 -58.73-58.85 -58.93 -58.99 -59.02 -59.03 -59.04 -59.05 -59.07 -59.10 -59.16 -59.25 -59.38-59.55 -59.77 -60.01 -60.24 -60.45 -60.59 100-64.35 -64.12 -64.00 -63.98 -64.06 -64.22-64.45 -64.74 -65.06 -65.39 -65.72 -66.01-66.25 -66.42 -66.55 -66.67 -66.80 -66.98-67.25 -67.61 -68.03 -68.45 -68.79 -69.01 70-68.35 -68.08 -67.82 -67.54 -67.27 -67.01-66.75 -66.50 -66.27 -66.06 -65.88 -65.74-65.65 -65.60 -65.60 -65.61 -65.64 -65.65-65.65 -65.61 -65.49 -65.24 -64.84 -64.22 50-60.25 -60.73 -61.06 -61.27 -61.36 -61.35-61.25 -61.08 -60.87 -60.64 -60.41 -60.20-60.05 -59.97 -59.93 -59.92 -59.91 -59.86-59.75 -59.56 -59.32 -59.04 -58.75 -58.48 30-52.85 -53.11 -53.34 -53.54 -53.70 -53.84-53.95 -54.03 -54.10 -54.16 -54.21 -54.27-54.35 -54.44 -54.53 -54.59 -54.60 -54.53-54.35 -54.06 -53.69 -53.29 -52.94 -52.67 20-50.35 -50.53 -50.67 -50.78 -50.86 -50.92-50.95 -50.96 -50.96 -50.94 -50.91 -50.88-50.85 -50.82 -50.80 -50.76 -50.71 -50.64-50.55 -50.42 -50.27 -50.08 -49.86 -49.62 10-47.45 -47.91 -48.28 -48.54 -48.72 -48.82-48.85 -48.81 -48.72 -48.58 -48.39 -48.18-47.95 -47.70 -47.44 -47.18 -46.91 -46.63-46.35 -46.08 -45.83 -45.64 -45.52 -45.52 Tm-19.23 -19.35 -19.48 -19.61 -19.74 -19.88-20.02 -20.17 -20.31 -20.45 -20.59 -20.72-20.85-20.97 -21.09 -21.19 -21.27 -21.33-21.37 -21.38 -21.39 -21.40 -21.43 -21.51 To obtain temperature values between 0, 6, 12, 18 hrs the spline approximation is used.
Temperature Data Comparison of the temperature profiles in the atmosphere with the model GFS for Greifswald point in2009 (data accuracy).
Temperature data January March June September DecemberComparison of the temperature profiles in the atmosphere with the model GFS for Moscow in 1999. Black circles-measurements, red triangles- GFS model).
Temperature Data in Real Time Deflection of forecasted distributionTfor on the nearest day from the model value T calculated in 5 days, is about 0.5 ºС and doesn’t exceed 1 ºС.
Points where on the basis of GFS(Global Forecast System) model the hourly vertical profiles of temperature at 17 standard isobaric levels of the atmosphere are obtained.
Map of points where on the basis of GFS(Global Forecast System) model the hourly vertical profiles of temperature at 17 standard isobaric levels of the atmosphere are obtained.
μddb— muon detector data basehttp://cr20.izmiran.ru/phpMyAdmin Data are kept in three tables: muon data, atmospheric pressure and temperature distribution. Hourly Data Pressure Data Temperаture (P) Data
Example for Greifswald, 2009 Grey-hourly data, red-daily) 2009! Gaussian distribution Uncorrected data – Two peaks!
Example for Greifswald, 2009 2009! Gaussian distribution Uncorrected data – Two peaks!
Example for YangBaJing, 2009 For YangBaJing station (altitude 4300 m) the applied method did not give a chance to correct completely data for temperature effect. Experimentally found temperature coefficient is twice more than calculated and has no physical background. This may be caused by a quality of current data at this station.
Example for Baksan, 2009 For demonstration of the method ability the pictures for the underground detector Baksan are presented (energy threshold is 220 GeV). The usual CR variations should be negligible small at this detector, so, the observed variations reflect only local fluctuations and temperature effect. Temperature effect can be estimated on this station for any other year, not only for the quiet period. In this case (underground detector) the temperature effect is positive in contrast with above demonstrated stations.
Conclusions • Altitudinal distribution of temperature in the atmosphere, byGFS model obtained, is sufficient for correction for temperature effect the observable muon data in real time. • For all the points with muon detectors the collection of hourly data on the altitudinal distribution of temperature in the atmosphere is carried out. • Calculations of the density of temperature coefficient are necessary to perform for all the directions of existed detectors accounting its real geometry. • It is more preferable to find the temperature coefficients experimentally, using method of effective temperature.
Many thanks to:- the Workshop Organizers for their kindly interest.- all of your for attentions.
Для возможного канадского проекта с точки зрения методики исключения метеорологических эффектов необходимо: 1) Минутные данные. Прецизионный ДД (точность несколько десятых mb и такая же долговременная стабильность). Необходимость минутных данных обсудить, может достаточно часовых. 2) Часовые данные. Формирование вертикального температурного распределения по данным модели атмосферы для 17 стандартных изобарических уровней. Измерение локальной приземной температуры T2 (точность несколько десятых градусов и такая же долговременная стабильность).
RESULTS on correction for temperature effect На рисунке приведены: Uncorrected data (blue), Corrected for temperature effect| by integral method (red); Mean-mass T method-black; Temperature variations Computed by integral method-brown, Median method (orange); Duperje method- green