Uncertainty in the measurement of indoor temperature and humidity in naturally ventilated dairy buildings as influenced by measurement technique and data variability

[EN] The microclimatic conditions in dairy buildings affect animal welfare and gaseous emissions. Measurements are highly variable due to the inhomogeneous distribution of heat and humidity sources (related to farm management) and the turbulent inflow (associated with meteorologic boundary conditions). The selection of the measurement strategy (number and position of the sensors) and the analysis methodology adds to the uncertainty of the applied measurement technique. To assess the suitability of different sensor positions, in situations where monitoring in the direct vicinity of the animals is not possible, we collected long-term data in two naturally ventilated dairy barns in Germany between March 2015 and April 2016 (horizontal and vertical profiles with 10 to 5 min temporal resolution). Uncertainties related to the measurement setup were assessed by comparing the device outputs under lab conditions after the on-farm experiments. We found out that the uncertainty in measurements of relative humidity is of particular importance when assessing heat stress risk and resulting economic losses in terms of temperature-humidity index. Measurements at a height of approximately 3 m-3.5 m turned out to be a good approximation for the microclimatic conditions in the animal occupied zone (including the air volume close to the emission active zone). However, further investigation along this cross-section is required to reduce uncertainties related to the inhomogeneous distribution of humidity. In addition, a regular sound cleaning (and if possible recalibration after few months) of the measurement devices is crucial to reduce the instrumentation uncertainty in long-term monitoring of relative humidity in dairy barns (C) 2017 The Authors. Published by Elsevier Ltd on behalf of IAgrE.The work was financially supported by the German Federal Ministry of Food and Agriculture (BMEL) through the Federal Office for Agriculture and Food (BLE), grant number 2814ERA02C.Hempel, S.; König, M.; Menz, C.; Janke, D.; Amon, B.; Banhazi, T.; Estellés, F.... (2018). Uncertainty in the measurement of indoor temperature and humidity in naturally ventilated dairy buildings as influenced by measurement technique and data variability. Biosystems Engineering. 166:58-75. https://doi.org/10.1016/j.biosystemseng.2017.11.004S587516

Retrieval of characteristic parameters for water vapour transmittance in the development of ground based Sun-Sky radiometric measurements of columnar water vapour

Sunsky radiometers are instruments created for aerosol study, but they can measure in the water vapour absorption band allowing the estimation of columnar water vapour in clear sky simultaneously with aerosol characteristics, with high temporal resolution. A new methodology is presented for estimating calibration parameters (i.e. characteristic parameters of the atmospheric transmittance and solar calibration constant) directly from the sunsky radiometer measurements. The methodology is based on the hypothesis that characteristic parameters of the atmospheric transmittance are dependent on vertical profiles of pressure, temperature and moisture occurring at each site of measurement. To obtain the parameters from the proposed methodology some seasonal independent measurements of columnar water vapour taken over a large range of solar zenith angle simultaneously with the sunsky radiometer measurements, are needed. In this work high time resolution columnar water vapour measurements by GPS were used as independent data set, but also the case when such measurements are not available was considered by developing the surface humidity method (SHM). This methodology makes it possible to retrieve the needed independent data set of columnar water vapour using the standard surface meteorological observations (temperature, pressure and relative humidity) more readily available. The time pattern of columnar water vapour from sunsky radiometer retrieved using both the methodologies was compared with simultaneous measurements from microwave radiometer, radiosondings and GPS. Water vapour from sunsky radiometer, obtained using GPS independent measurements, was characterized by an error varying from 1% up to 5%, whereas water vapour from SHM showed an error from 1% up to 11%, depending on the local columnar water occurring at the site during the year. These errors were estimated by comparing water vapour series from sunsky radiometer against measurements taken by GPS at a nearby station. The accordance between retrievals from sunsky radiometer and simultaneous measurements from the other instruments was found always within the error both in the case of SHM and of the GPS independent data set. Water vapour obtained using characteristic parameters of the atmospheric transmittance dependent on water vapour was also compared against GPS retrievals, showing a clear improvement with respect to the case when these parameters are kept fixed

Evaluation of the new ESR network software for the retrieval of direct sun products from CIMEL CE318 and PREDE POM01 sun-sky radiometers

The European Skynet Radiometers network (EuroSkyRad or ESR) has been recently established as a research network of European PREDE sun-sky radiometers. Moreover, ESR is federated with SKYNET, an international network of PREDE sun-sky radiometers mostly present in East Asia. In contrast to SKYNET, the European network also integrates users of the CIMEL CE318 sky–sun photometer. Keeping instrumental duality in mind, a set of open source algorithms has been developed consisting of two modules for (1) the retrieval of direct sun products (aerosol optical depth, wavelength exponent and water vapor) from the sun extinction measurements; and (2) the inversion of the sky radiance to derive other aerosol optical properties such as size distribution, single scattering albedo or refractive index. In this study we evaluate the ESR direct sun products in comparison with the AERosol RObotic NETwork (AERONET) products. Specifically, we have applied the ESR algorithm to a CIMEL CE318 and PREDE POM simultaneously for a 4-yr database measured at the Burjassot site (Valencia, Spain), and compared the resultant products with the AERONET direct sun measurements obtained with the same CIMEL CE318 sky–sun photometer. The comparison shows that aerosol optical depth differences are mostly within the nominal uncertainty of 0.003 for a standard calibration instrument, and fall within the nominal AERONET uncertainty of 0.01–0.02 for a field instrument in the spectral range 340 to 1020 nm. In the cases of the Ångström exponent and the columnar water vapor, the differences are lower than 0.02 and 0.15 cm, respectively. Therefore, we present an open source code program that can be used with both CIMEL and PREDE sky radiometers and whose results are equivalent to AERONET and SKYNET retrievals

How far may rabbit cage’s space recommendations reach: the gap between science and regulations

Animal welfare is a major concern in livestock production. The impact of housing conditions on animal well-being is under study in several species such as chickens, laying hens, or pigs and, some European Directives have been developed in order to protect animal welfare. In these Directives aspects related to housing conditions have been regulated. Several attempts have been done to propose a European Directive for rabbit farming in terms of welfare protection, but any agreement have been reached for now. One of the main hot spots is related to cages design. Strong efforts are being made nowadays in order to reach a compromise to define cages characteristics which really benefit rabbits’ welfare. In this context, the main aim of this work is to develop a critical review about rabbit does space needs in relation to cages design. It becomes necessary to determine the behavioural needs of rabbits in terms of welfare, so that the space requirements can be known. This might lead to define appropriate housing conditions, although the needs of the animal can change according to age, learning, diurnal rhythm, season and genetic relations. In general, regarding the height of the cages (as well as width), few papers have been found and a European Food Safety Authorities’ report recognises, that in the absence of scientific evidence concerning these needs of rabbits, it may be important for growing rabbits to be able to sit and stand with ears erect, as well as rear up occasionally. Regarding the use of platforms, enriching the cage with raised platforms aims at satisfying the doe’s need for isolation from the litter, rather than stimulating exercise. Nevertheless, the use of platforms may cause hygiene problems which have to be solved. Therefore, elevated platforms might be considered as environmental enrichment elements and not as structural needs. The patent general lack of information in the topic must be considered when new housing conditions are being proposed, and only those aspects which could lead to real improvement of rabbit’s welfare have to be taken into account. The authors propose a cage size for lactating does of 45×75×38 cm (height, length and width) with a minimum surface of 3,500 cm2 (nest-box not included)

Intercomparison of spectroradiometers and Sun photometers for the determination of the aerosol optical depth during the VELETA-2002 field campaign

[ 1] In July 2002 the VELETA-2002 field campaign was held in Sierra Nevada ( Granada) in the south of Spain. The main objectives of this field campaign were the study of the influence of elevation and atmospheric aerosols on measured UV radiation. In the first stage of the field campaign, a common calibration and intercomparison between Licor-1800 spectroradiometers and Cimel-318 Sun photometers was performed in order to assess the quality of the measurements from the whole campaign. The intercomparison of the Licor spectroradiometers showed, for both direct and global irradiances, that when the comparisons were restricted to the visible part of the spectrum the deviations were within the instruments' nominal accuracies which allows us to rely on these instruments for measuring physical properties of aerosols at the different measurement stations. A simultaneous calibration on AOD data was performed for the Cimel-318 Sun photometers. When a common calibration and methodology was applied, the deviation was lowered to much less than 0.01 for AOD. At the same time an intercomparison has been made between the AOD values given by the spectroradiometers and the Sun photometers, with deviations obtained from 0.01 to 0.03 for the AOD in the visible range, depending on the channel. In the UVA range, the AOD uncertainty was estimated to be around 0.02 and 0.05 for Cimel and Licor respectively. In general the experimental differences were in agreement with this uncertainty estimation. In the UVB range the AOD measurements should not be used due to maximum instrumental uncertainties

Towards an integrated crowdsourcing definition

Crowdsourcing is a relatively recent concept that encompasses many practices. This diversity leads to the blurring of the limits of crowdsourcing that may be identified virtually with any type of internet-based collaborative activity, such as co-creation or user innovation. Varying definitions of crowdsourcing exist, and therefore some authors present certain specific examples of crowdsourcing as paradigmatic, while others present the same examples as the opposite. In this article, existing definitions of crowdsourcing are analysed to extract common elements and to establish the basic characteristics of any crowdsourcing initiative. Based on these existing definitions, an exhaustive and consistent definition for crowdsourcing is presented and contrasted in 11 cases.Estelles Arolas, E.; González-Ladrón-De-Guevara, F. (2012). Towards an integrated crowdsourcing definition. Journal of Information Science. 32(2):189-200. doi:10.1177/0165551512437638S189200322Vukovic, M., & Bartolini, C. (2010). Towards a Research Agenda for Enterprise Crowdsourcing. Leveraging Applications of Formal Methods, Verification, and Validation, 425-434. doi:10.1007/978-3-642-16558-0_36Brabham, D. C. (2008). Crowdsourcing as a Model for Problem Solving. Convergence: The International Journal of Research into New Media Technologies, 14(1), 75-90. doi:10.1177/1354856507084420Vukovic, M. (2009). Crowdsourcing for Enterprises. 2009 Congress on Services - I. doi:10.1109/services-i.2009.56Doan, A., Ramakrishnan, R., & Halevy, A. Y. (2011). Crowdsourcing systems on the World-Wide Web. Communications of the ACM, 54(4), 86. doi:10.1145/1924421.1924442Brabham, D. C. (2008). Moving the crowd at iStockphoto: The composition of the crowd and motivations for participation in a crowdsourcing application. First Monday, 13(6). doi:10.5210/fm.v13i6.2159Huberman, B. A., Romero, D. M., & Wu, F. (2009). Crowdsourcing, attention and productivity. Journal of Information Science, 35(6), 758-765. doi:10.1177/0165551509346786Andriole, S. J. (2010). Business impact of Web 2.0 technologies. Communications of the ACM, 53(12), 67. doi:10.1145/1859204.1859225Denyer, D., Tranfield, D., & van Aken, J. E. (2008). Developing Design Propositions through Research Synthesis. Organization Studies, 29(3), 393-413. doi:10.1177/0170840607088020Egger, M., Smith, G. D., & Altman, D. G. (Eds.). (2001). Systematic Reviews in Health Care. doi:10.1002/9780470693926Tatarkiewicz, W. (1980). A History of Six Ideas. doi:10.1007/978-94-009-8805-7Cosma, G., & Joy, M. (2008). Towards a Definition of Source-Code Plagiarism. IEEE Transactions on Education, 51(2), 195-200. doi:10.1109/te.2007.906776Brabham, D. C. (2009). Crowdsourcing the Public Participation Process for Planning Projects. Planning Theory, 8(3), 242-262. doi:10.1177/1473095209104824Alonso, O., & Lease, M. (2011). Crowdsourcing 101. Proceedings of the fourth ACM international conference on Web search and data mining - WSDM ’11. doi:10.1145/1935826.1935831Bederson, B. B., & Quinn, A. J. (2011). Web workers unite! addressing challenges of online laborers. Proceedings of the 2011 annual conference extended abstracts on Human factors in computing systems - CHI EA ’11. doi:10.1145/1979742.1979606Grier, D. A. (2011). Not for All Markets. Computer, 44(5), 6-8. doi:10.1109/mc.2011.155Heer, J., & Bostock, M. (2010). Crowdsourcing graphical perception. Proceedings of the 28th international conference on Human factors in computing systems - CHI ’10. doi:10.1145/1753326.1753357Heymann, P., & Garcia-Molina, H. (2011). Turkalytics. Proceedings of the 20th international conference on World wide web - WWW ’11. doi:10.1145/1963405.1963473Kazai, G. (2011). In Search of Quality in Crowdsourcing for Search Engine Evaluation. Advances in Information Retrieval, 165-176. doi:10.1007/978-3-642-20161-5_17La Vecchia, G., & Cisternino, A. (2010). Collaborative Workforce, Business Process Crowdsourcing as an Alternative of BPO. Lecture Notes in Computer Science, 425-430. doi:10.1007/978-3-642-16985-4_40Liu, E., & Porter, T. (2010). Culture and KM in China. VINE, 40(3/4), 326-333. doi:10.1108/03055721011071449Oliveira, F., Ramos, I., & Santos, L. (2010). Definition of a Crowdsourcing Innovation Service for the European SMEs. Lecture Notes in Computer Science, 412-416. doi:10.1007/978-3-642-16985-4_37Porta, M., House, B., Buckley, L., & Blitz, A. (2008). Value 2.0: eight new rules for creating and capturing value from innovative technologies. Strategy & Leadership, 36(4), 10-18. doi:10.1108/10878570810888713Ribiere, V. M., & Tuggle, F. D. (Doug). (2010). Fostering innovation with KM 2.0. VINE, 40(1), 90-101. doi:10.1108/03055721011024955Sloane, P. (2011). The brave new world of open innovation. Strategic Direction, 27(5), 3-4. doi:10.1108/02580541111125725Wexler, M. N. (2011). Reconfiguring the sociology of the crowd: exploring crowdsourcing. International Journal of Sociology and Social Policy, 31(1/2), 6-20. doi:10.1108/01443331111104779Whitla, P. (2009). Crowdsourcing and Its Application in Marketing Activities. Contemporary Management Research, 5(1). doi:10.7903/cmr.1145Yang, J., Adamic, L. A., & Ackerman, M. S. (2008). Crowdsourcing and knowledge sharing. Proceedings of the 9th ACM conference on Electronic commerce - EC ’08. doi:10.1145/1386790.1386829Brabham, D. C. (2010). MOVING THE CROWD AT THREADLESS. Information, Communication & Society, 13(8), 1122-1145. doi:10.1080/13691181003624090Giudice, K. D. (2010). Crowdsourcing credibility: The impact of audience feedback on Web page credibility. Proceedings of the American Society for Information Science and Technology, 47(1), 1-9. doi:10.1002/meet.14504701099Stewart, O., Huerta, J. M., & Sader, M. (2009). Designing crowdsourcing community for the enterprise. Proceedings of the ACM SIGKDD Workshop on Human Computation - HCOMP ’09. doi:10.1145/1600150.1600168Maslow, A. H. (1943). A theory of human motivation. Psychological Review, 50(4), 370-396. doi:10.1037/h0054346Veal, A. J. (Ed.). (2002). Leisure and tourism policy and planning. doi:10.1079/9780851995465.0000Dahlander, L., & Gann, D. M. (2010). How open is innovation? Research Policy, 39(6), 699-709. doi:10.1016/j.respol.2010.01.01

Results from the Fourth WMO Filter Radiometer Comparison for aerosol optical depth measurements

This study presents the results of the Fourth Filter Radiometer Comparison that was held in Davos, Switzerland, between 28 September and 16 October 2015. Thirty filter radiometers and spectroradiometers from 12 countries participated including reference instruments from global aerosol networks. The absolute differences of all instruments compared to the reference have been based on the World Meteorological Organization (WMO) criterion defined as follows: 95% of the measured data has to be within 0.005±0.001∕m (where m is the air mass). At least 24 out of 29 instruments achieved this goal at both 500 and 865nm, while 12 out of 17 and 13 out of 21 achieved this at 368 and 412nm, respectively. While searching for sources of differences among different instruments, it was found that all individual differences linked to Rayleigh, NO2, ozone, water vapor calculations and related optical depths and air mass calculations were smaller than 0.01 in aerosol optical depth (AOD) at 500 and 865nm. Different cloud-detecting algorithms used have been compared. Ångström exponent calculations showed relatively large differences among different instruments, partly because of the high calculation uncertainty of this parameter in low AOD conditions. The overall low deviations of these AOD results and the high accuracy of reference aerosol network instruments demonstrated a promising framework to achieve homogeneity, compatibility and harmonization among the different spectral AOD networks in the near future

Exploring Barriers and Opportunities in Adopting Crowdsourcing Based New Product Development in Manufacturing SMEs

Crowdsourcing is an innovative business practice of obtaining needed services, ideas, or content or even funds by soliciting contributions from a large group of people (the ‘Crowd’). The potential benefits of utilizing crowdsourcing in product design are well-documented, but little research exists on what are the barriers and opportunities in adopting crowdsourcing in new product development (NPD) of manufacturing SMEs. In order to answer the above questions, a Proof of Market study is carried out on crowdsourcing-based product design under an Innovate UK funded Smart project, which aims at identifying the needs, challenges and future development opportunities associated with adopting crowdsourcing strategies for NPD. The research findings from this study are reported here and can be used to guide future development of crowdsourcing-based collaborative design methods and tools and provide some practical references for industry to adopt this new and emerging collaborative design method in their business

Search of the early O3 LIGO data for continuous gravitational waves from the Cassiopeia A and Vela Jr. supernova remnants

R. Abbott et al.We present directed searches for continuous gravitational waves from the neutron stars in the Cassiopeia A (Cas A) and Vela Jr. supernova remnants. We carry out the searches in the LIGO detector data from the first six months of the third Advanced LIGO and Virgo observing run using the weave semicoherent method, which sums matched-filter detection-statistic values over many time segments spanning the observation period. No gravitational wave signal is detected in the search band of 20–976 Hz for assumed source ages greater than 300 years for Cas A and greater than 700 years for Vela Jr. Estimates from simulated continuous wave signals indicate we achieve the most sensitive results to date across the explored parameter space volume, probing to strain magnitudes as low as ∼6.3×10−26 for Cas A and ∼5.6×10−26 for Vela Jr. at frequencies near 166 Hz at 95% efficiency.The authors also gratefully acknowledge the support of the Science and Technology Facilities Council (STFC) of the United Kingdom,the Max-Planck-Society (MPS), and the State of Niedersachsen/Germany for support of the construction o Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. The authors gratefully acknowledge the Italian Istituto Nazionale di Fisica Nucleare (INFN), the French Centre National de la Recherche Scientifique (CNRS) and the Netherlands Organization for Scientific Research (NWO), for the construction and operation of the Virgo detector and the creation and support of the EGO consortium. The authors also gratefully acknowledge research support from these agencies as well as by the Council of Scientific and Industrial Research of India, the Department of Science and Technology, India, the Science and Engineering Research Board (SERB), India, the Ministry of Human Resource Development,India, the Spanish Agencia Estatal de Investigación (AEI), the Spanish Ministerio de Ciencia e Innovación and Ministerio de Universidades, the Conselleria de Fons Europeus, Universitati Cultura and the Direcció General de Política Universitaria i Recerca del Govern de les Illes Balears, the Conselleria d’Innovació, Universitats, Ciencia i Societat Digital de la Generalitat Valenciana and the CERCA Programme Generalitatd e Catalunya,Spain,the National Science Centre of Poland and the European Union—European Regional Development Fund; Foundation for Polish Science (FNP),the Swiss National Science Foundation(SNSF), the Russian Foundation for Basic Research, the Russian Science Foundation, the European Commission, the European Social Funds(ESF),the European Regional Development Funds(ERDF),the Royal Society,the Scottish Funding Council,the Scottish Universities Physics Alliance,the Hungarian ScientificResearch Fund (OTKA),the French Lyon Institute of Origins (LIO),the Belgian Fonds de la Recherche Scientifique (FRS-FNRS),Actions de Recherche Concertees(ARC) and Fonds Wetenschappelijk Onderzoek—Vlaanderen (FWO),Belgium,the ParisÎle-de-France Region,the National Research, Development and Innovation Office Hungary (NKFIH),the National Research Foundation of Korea,the Natural Science and Engineering Research Council Canada,Canadian Foundation for Innovation (CFI),the Brazilian Ministr of Science,Technology, and Innovations,the International Center for Theoretical Physics South American Institute for Fundamental Research(ICTP-SAIFR), the Research Grants Councilof Hong Kong,the National Natural Science Foundation of China(NSFC),the LeverhulmeTrust,the Research Corporation,the Ministry of Science and Technology (MOST),Taiwan, the United States Departmentof Energy, and the Kavli Foundation.The authors gratefully acknowledge thes upport of theNSF,STFC,INFN and CNRSf or provision of computational resources.Peer reviewe

Constraints on dark photon dark matter using data from LIGO's and Virgo's third observing run

LIGO Scientific Collaboration - Virgo Collaboration - KAGRA Collaboration: Abbott, R. et al.We present a search for dark photon dark matter that could couple to gravitational-wave interferometers using data from Advanced LIGO and Virgo’s third observing run. To perform this analysis, we use two methods, one based on cross-correlation of the strain channels in the two nearly aligned LIGO detectors, and one that looks for excess power in the strain channels of the LIGO and Virgo detectors. The excess power method optimizes the Fourier transform coherence time as a function of frequency, to account for the expected signal width due to Doppler modulations. We do not find any evidence of dark photon dark matter with a mass between mA∼10−14–10−11 eV/c2, which corresponds to frequencies between 10–2000 Hz, and therefore provide upper limits on the square of the minimum coupling of dark photons to baryons, i.e., U(1)B dark matter. For the cross-correlation method, the best median constraint on the squared coupling is ∼1.31×10−47 at mA∼4.2×10−13eV/c2; for the other analysis, the best constraint is ∼2.4×10−47 at mA∼5.7×10−13eV/c2. These limits improve upon those obtained in direct dark matter detection experiments by a factor of ∼100 for mA∼[2–4]×10−13eV/c2, and are, in absolute terms, the most stringent constraint so far in a large mass range mA∼2×10−13–8×10−12eV/c2.This material is based upon work supported by NSF’s LIGO Laboratory which is a major facility fully funded by the National Science Foundation. The authors also gratefully acknowledge the support of the Science and Technology Facilities Council (STFC) of the United Kingdom, the Max-Planck-Society (MPS), and the State of Niedersachsen/ Germany for support of the construction of Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. The authors gratefully acknowledge the Italian Istituto Nazionale di Fisica Nucleare (INFN), the French Centre National de la Recherche Scientifique (CNRS) and the Netherlands Organization for Scientific Research, for the construction and operation of the Virgo detector and the creation and support of the EGO consortium. The authors also gratefully acknowledge research support fromthese agencies as well as by the Council of Scientific and Industrial Research of India, the Department of Science and Technology, India, the Science & Engineering Research Board (SERB), India, the Ministry of Human Resource Development, India, the Spanish Agencia Estatal de Investigación, the Vicepresidencia i Conselleria d’Innovació, Recerca i Turisme and the Conselleria d’Educació i Universitat del Govern de les Illes Balears, the Conselleria d’Innovació, Universitats, Ciencia i Societat Digital de la Generalitat Valenciana and the CERCA Programme Generalitat de Catalunya, Spain, the National Science Centre of Poland and the Foundation for Polish Science (FNP), the Swiss National Science Foundation (SNSF), the Russian Foundation for Basic Research, the Russian Science Foundation, the European Commission, the European Regional Development Funds (ERDF), the Royal Society, the Scottish Funding Council, the Scottish Universities Physics Alliance, the Hungarian Scientific Research Fund (OTKA), the French Lyon Institute of Origins (LIO), the Belgian Fonds de la Recherche Scientifique (FRS-FNRS), Actions de Recherche Concert´ees (ARC) and Fonds Wetenschappelijk Onderzoek—Vlaanderen (FWO), Belgium, the Paris Île-de-France Region, the National Research, Development and Innovation Office Hungary (NKFIH), the National Research Foundation of Korea, the Natural Science and Engineering Research Council Canada, Canadian Foundation for Innovation (CFI), the Brazilian Ministry of Science, Technology, and Innovations, the International Center for Theoretical Physics South American Institute for Fundamental Research (ICTPSAIFR), the Research Grants Council of Hong Kong, the National Natural Science Foundation of China (NSFC), the Leverhulme Trust, the Research Corporation, the Ministry of Science and Technology (MOST), Taiwan, the United States Department of Energy, and the Kavli Foundation. The authors gratefully acknowledge the support of the NSF, STFC, INFN and CNRS for provision of computational resources. This work was supported by MEXT, JSPS Leading-edge Research Infrastructure Program, JSPS Grant-in-Aid for Specially Promoted Research 26000005, JSPS Grant-in-Aid for Scientific Research on Innovative Areas 2905: No. JP17H06358, No. JP17H06361 and No. JP17H06364, JSPS Core-to-Core Program A. Advanced Research Networks, JSPS Grant-in-Aid for Scientific Research (S) No. 17H06133, the joint research program of the Institute forCosmic Ray Research, University of Tokyo, National Research Foundation (NRF) and Computing Infrastructure Project of KISTI-GSDC in Korea, Academia Sinica (AS), AS Grid Center (ASGC) and the Ministry of Science and Technology (MoST) in Taiwan under grants including AS-CDA-105-M06, Advanced Technology Center (ATC) of NAOJ, and Mechanical Engineering Center of KEK.Peer reviewe