• Title/Summary/Keyword: relative error

Search Result 1,544, Processing Time 0.03 seconds

Calibration of Hargreaves Equation Coefficient for Estimating Reference Evapotranspiration in Korea (우리나라 기준증발산량 추정을 위한 Hargreaves 공식의 계수 보정)

  • Hwang, Seon-ah;Han, Kyung-hwa;Zhang, Yong-seon;Cho, Hee-rae;Ok, Jung-hun;Kim, Dong-Jin;Kim, Gi-sun;Jung, Kang-ho
    • Korean Journal of Agricultural and Forest Meteorology
    • /
    • v.21 no.4
    • /
    • pp.238-249
    • /
    • 2019
  • The evapotranspiration is estimated based on weather factors such as temperature, wind speed and humidity, and the Hargreaves equation is a simple equation for calculating evapotranspiration using temperature data. However, the Hargreaves equation tends to be underestimated in areas with wind speeds above 3 m s-1 and overestimated in areas with high relative humidity. The study was conducted to determine Hargreaves equation coefficient in 82 regions in Korea by comparing evapotranspiration determined by modified Hargreaves equation and the Penman-Monteith equation for the time period of 2008~2018. The modified Hargreaves coefficients for 50 inland areas were estimated to be 0.00173~0.00232(average 0.00196), which is similar to or lower than the default value 0.0023. On the other hand, there are 32 coastal areas, and the modified coefficients ranged from 0.00185 to 0.00303(average 0.00234). The east coastal area was estimated to be similar to or higher than the default value, while the west and south coastal areas showed large deviations by area. As results of estimating the evapotranspiration by the modified Hargreaves coefficient, root mean square error(RMSE) is reduced from 0.634~1.394(average 0.857) to 0.466~1.328(average 0.701), and Nash-Sutcliffe Coefficient(NSC) increased from -0.159~0.837(average 0.647) to -0.053~0.910(average 0.755) compared with original Hargreaves equation. Therefore, we confirmed that the Hargreaves equation can be overestimated or underestimated compared to the Penman-Monteith equation, and expected that it will be able to calculate the high accuracy evapotranspiration using the modified Hargreaves equation. This study will contribute to water resources planning, irrigation schedule, and environmental management.

The Three Dimensional Analysis on Nasal Airway Morphology in Class III Malocclusion (골격성 III급 부정교합자의 Nasal Airway 형태에 관한 3차원적 분석 연구)

  • Kim, Moon-Hwan;Lee, Jin-Woo;Cha, Kyung-Suk;Chung, Dong-Hwa
    • Journal of Dental Rehabilitation and Applied Science
    • /
    • v.24 no.4
    • /
    • pp.389-403
    • /
    • 2008
  • In Angle's Class III malocclusion, which has higher incidence in Korean than Western, depressed midfacial profile with protruded lower lips and mandible may give rise to many functional, esthetic, psychological, social problems. Due to the different malocclusion incidence according to racial differences, many previous studies focused on the relationship between Class II malocclusion and nasal airway obstruction. Previous studies used lateral cephalography which has limitations of 2 dimensional image with projection error and identification error. Therefore, the purpose of this study was to analyze morphologic differences in the nasal airway between normal occlusion and Angle's Class III malocclusion patients using 3-dimensional facial computed tomography. Thirteen normal occlusion(7 men and 6 women) and sixteen skeletal Class III(7 men and 9 women) patients were selected and 3-dimensional facial computed tomography taking was performed. Comparison between two group in volume and sectional area of nasal airway were carried out. The results were followed. 1. In the comparison of absolute nasal airway volume, oropharyngeal space of experimental group were larger than control group but there are no significant difference in other. 2. In the comparison of relative nasal airway volume, oropharyngeal space of experimental group were larger than control group but there are no significant difference in other. 3. In the oropharyngeal space width on frontal and lateral view, the similar tendency was revealed between two groups. 4. In the lateral curvature of nasal airway, the similar tendency was revealed between two groups.

The Effect of the Optical Points Difference between Finished-Reading Glasses and Dispensing Reading Glasses (완성품 돋보기와 조제가공된 돋보기가 광학적 요소에 미치는 영향)

  • Shim, Young-Cheol;Yoo, Gun-Chang;Kim, In-Suk
    • Journal of Korean Ophthalmic Optics Society
    • /
    • v.13 no.3
    • /
    • pp.65-71
    • /
    • 2008
  • Purpose: This paper studied the effect of eyes on the comparison between the distance optical centers problem of dispensing reading glasses made by optician and finished reading glasses in the current market. Methods: The method of this study has been measured by eleven different categories from +1.00D to +4.00D. This study also separated into three groups by their optical frame size and measured optical center point (O.C) and optical center height (O.H) with 200 peoples of man and females over 40 years old without ocular disease living in Gwang-san gu, Gwang-Ju city. Results: As a result, optical center point ranged from 57 mm to 80 mm and it turned to be most common range is from 61 mm to 65 mm (64.6%). Moreover, the optical center height ranged from 1 mm to 8mm and most common ranged (23%) were 4 mm. In other words, finished reading glasses have irregular optical ranges. After observing 200 people who are over 40 years old men and women, result shows that more than 75.5% (151 people) currently use finished reading glasses. Survey of 151 people, most common error between the finished reading glasses's O.C and the wearers P.D were 4 mm (45%). Furthermore, the most common error between the finished reading glassses's O.H and the wearers O.H ranged from 3 mm to 4 mm. Astonishingly, the entire 151 people who wear finished reading glasses appeal that they feel tiredness on their eyes when they wear finished reading glasses. 53 people (35%) claimed that they feel tiredness on their eyes after 10 to 20 minutes wearing finished reading glasses. Base on the research, We conducted more experiment to find the value of prism of optical centers err because it will tell us whether the finished reading glasses are good enough to wear or not. We multiplied diopter by the difference between finished reading glasses's O.C. and wearer's P.D. Consequently, We found out that the finished reading glasses counter to the German RAL-RG 915 policy. And We also found that it is relative to the diopter of lenses. In conclusion, based on the researched that wearing finished reading glasses have a dangerous factor for our vision. Therefore optician must need to recommend correctly made dispensing reading glasses based on the optical center point.

  • PDF

Development of Greenhouse Cooling and Heating Load Calculation Program Based on Mobile (모바일 기반 온실 냉난방 부하 산정 프로그램 개발)

  • Moon, Jong Pil;Bang, Ji Woong;Hwang, Jeongsu;Jang, Jae Kyung;Yun, Sung Wook
    • Journal of Bio-Environment Control
    • /
    • v.30 no.4
    • /
    • pp.419-428
    • /
    • 2021
  • In order to develope a mobile-based greenhouse energy calculation program, firstly, the overall thermal transmittance of 10 types of major covers and 16 types of insulation materials were measured. In addition, to estimate the overall thermal transmittance when the cover and insulation materials were installed in double or triple layers, 24 combinations of double installations and 59 combinations of triple installations were measured using the hotbox. Also, the overall thermal transmittance value for a single material and the thermal resistance value were used to calculate the overall thermal transmittance value at the time of multi-layer installation of covering and insulating materials, and the linear regression equation was derived to correct the error with the measured values. As a result of developing the model for estimating thermal transmittance when installing multiple layers of coverings and insulating materials based on the value of overall thermal transmittance of a single-material, the model evaluation index was 0.90 (good when it is 0.5 or more), indicating that the estimated value was very close to the actual value. In addition, as a result of the on-site test, it was evaluated that the estimated heat saving rate was smaller than the actual value with a relative error of 2%. Based on these results, a mobile-based greenhouse energy calculation program was developed that was implemented as an HTML5 standard web-based mobile web application and was designed to work with various mobile device and PC browsers with N-Screen support. It had functions to provides the overall thermal transmittance(heating load coefficient) for each combination of greenhouse coverings and thermal insulation materials and to evaluate the energy consumption during a specific period of the target greenhouse. It was estimated that an energy-saving greenhouse design would be possible with the optimal selection of coverings and insulation materials according to the region and shape of the greenhouse.

Observation of Ice Gradient in Cheonji, Baekdu Mountain Using Modified U-Net from Landsat -5/-7/-8 Images (Landsat 위성 영상으로부터 Modified U-Net을 이용한 백두산 천지 얼음변화도 관측)

  • Lee, Eu-Ru;Lee, Ha-Seong;Park, Sun-Cheon;Jung, Hyung-Sup
    • Korean Journal of Remote Sensing
    • /
    • v.38 no.6_2
    • /
    • pp.1691-1707
    • /
    • 2022
  • Cheonji Lake, the caldera of Baekdu Mountain, located on the border of the Korean Peninsula and China, alternates between melting and freezing seasonally. There is a magma chamber beneath Cheonji, and variations in the magma chamber cause volcanic antecedents such as changes in the temperature and water pressure of hot spring water. Consequently, there is an abnormal region in Cheonji where ice melts quicker than in other areas, freezes late even during the freezing period, and has a high-temperature water surface. The abnormal area is a discharge region for hot spring water, and its ice gradient may be used to monitor volcanic activity. However, due to geographical, political and spatial issues, periodic observation of abnormal regions of Cheonji is limited. In this study, the degree of ice change in the optimal region was quantified using a Landsat -5/-7/-8 optical satellite image and a Modified U-Net regression model. From January 22, 1985 to December 8, 2020, the Visible and Near Infrared (VNIR) band of 83 Landsat images including anomalous regions was utilized. Using the relative spectral reflectance of water and ice in the VNIR band, unique data were generated for quantitative ice variability monitoring. To preserve as much information as possible from the visible and near-infrared bands, ice gradient was noticed by applying it to U-Net with two encoders, achieving good prediction accuracy with a Root Mean Square Error (RMSE) of 140 and a correlation value of 0.9968. Since the ice change value can be seen with high precision from Landsat images using Modified U-Net in the future may be utilized as one of the methods to monitor Baekdu Mountain's volcanic activity, and a more specific volcano monitoring system can be built.

Pipetting Stability and Improvement Test of the Robotic Liquid Handling System Depending on Types of Liquid (용액에 따른 자동분주기의 분주능력 평가와 분주력 향상 실험)

  • Back, Hyangmi;Kim, Youngsan;Yun, Sunhee;Heo, Uisung;Kim, Hosin;Ryu, Hyeonggi;Lee, Guiwon
    • The Korean Journal of Nuclear Medicine Technology
    • /
    • v.20 no.2
    • /
    • pp.62-68
    • /
    • 2016
  • Purpose In a cyclosporine experiment using a robotic liquid handing system has found a deviation of its standard curve and low reproducibility of patients's results. The difference of the test is that methanol is mixed with samples and the extractions are used for the test. Therefore, we assumed that the abnormal test results came from using methanol and conducted this test. In a manual of a robotic liquid handling system mentions that we can choose several setting parameters depending on the viscosity of the liquids being used, the size of the sampling tips and the motor speeds that you elect to use but there's no exact order. This study was undertaken to confirm pipetting ability depending on types of liquids and investigate proper setting parameters for the optimum dispensing ability. Materials and Methods 4types of liquids(water, serum, methanol, PEG 6000(25%)) and $TSH^{125}I$ tracer(515 kBq) are used to confirm pipetting ability. 29 specimens for Cyclosporine test are used to compare results. Prepare 8 plastic tubes for each of the liquids and with multi pipette $400{\mu}l$ of each liquid is dispensed to 8 tubes and $100{\mu}l$ of $TSH^{125}I$ tracer are dispensed to all of the tubes. From the prepared samples, $100{\mu}l$ of liquids are dispensed using a robotic liquid handing system, counted and calculated its CV(%) depending on types of liquids. And then by adjusting several setting parameters(air gap, dispense time, delay time) the change of the CV(%)are calcutated and finds optimum setting parameters. 29 specimens are tested with 3 methods. The first(A) is manual method and the second(B) is used robotic liquid handling system with existing parameters. The third(C) is used robotic liquid handling system with adjusted parameters. Pipetting ability depending on types of liquids is assessed with CV(%). On the basis of (A), patients's test results are compared (A)and(B), (A)and(C) and they are assessed with %RE(%Relative error) and %Diff(%Difference). Results The CV(%) of the CPM depending on liquid types were water 0.88, serum 0.95, methanol 10.22 and PEG 0.68. As expected dispensing of methanol using a liquid handling system was the problem and others were good. The methanol's dispensing were conducted by adjusting several setting parameters. When transport air gap 0 was adjusted to 2 and 5, CV(%) were 20.16, 12.54 and when system air gap 0 was adjusted to 2 and 5, CV(%) were 8.94, 1.36. When adjusted to system air gap 2, transport air gap 2 was 12.96 and adjusted to system air gap 5, Transport air gap 5 was 1.33. When dispense speed was adjusted 300 to 100, CV(%) was 13.32 and when dispense delay was adjusted 200 to 100 was 13.55. When compared (B) to (A), the result increased 99.44% and %RE was 93.59%. When compared (C-system air gap was adjusted 0 to 5) to (A), the result increased 6.75% and %RE was 5.10%. Conclusion Adjusting speed and delay time of aspiration and dispense was meaningless but changing system air gap was effective. By adjusting several parameters proper value was found and it affected the practical result of the experiment. To optimize the system active efforts are needed through the test and in case of dispensing new types of liquids proper test is required to check the liquid is suitable for using the equipment.

  • PDF

An Intelligent Decision Support System for Selecting Promising Technologies for R&D based on Time-series Patent Analysis (R&D 기술 선정을 위한 시계열 특허 분석 기반 지능형 의사결정지원시스템)

  • Lee, Choongseok;Lee, Suk Joo;Choi, Byounggu
    • Journal of Intelligence and Information Systems
    • /
    • v.18 no.3
    • /
    • pp.79-96
    • /
    • 2012
  • As the pace of competition dramatically accelerates and the complexity of change grows, a variety of research have been conducted to improve firms' short-term performance and to enhance firms' long-term survival. In particular, researchers and practitioners have paid their attention to identify promising technologies that lead competitive advantage to a firm. Discovery of promising technology depends on how a firm evaluates the value of technologies, thus many evaluating methods have been proposed. Experts' opinion based approaches have been widely accepted to predict the value of technologies. Whereas this approach provides in-depth analysis and ensures validity of analysis results, it is usually cost-and time-ineffective and is limited to qualitative evaluation. Considerable studies attempt to forecast the value of technology by using patent information to overcome the limitation of experts' opinion based approach. Patent based technology evaluation has served as a valuable assessment approach of the technological forecasting because it contains a full and practical description of technology with uniform structure. Furthermore, it provides information that is not divulged in any other sources. Although patent information based approach has contributed to our understanding of prediction of promising technologies, it has some limitations because prediction has been made based on the past patent information, and the interpretations of patent analyses are not consistent. In order to fill this gap, this study proposes a technology forecasting methodology by integrating patent information approach and artificial intelligence method. The methodology consists of three modules : evaluation of technologies promising, implementation of technologies value prediction model, and recommendation of promising technologies. In the first module, technologies promising is evaluated from three different and complementary dimensions; impact, fusion, and diffusion perspectives. The impact of technologies refers to their influence on future technologies development and improvement, and is also clearly associated with their monetary value. The fusion of technologies denotes the extent to which a technology fuses different technologies, and represents the breadth of search underlying the technology. The fusion of technologies can be calculated based on technology or patent, thus this study measures two types of fusion index; fusion index per technology and fusion index per patent. Finally, the diffusion of technologies denotes their degree of applicability across scientific and technological fields. In the same vein, diffusion index per technology and diffusion index per patent are considered respectively. In the second module, technologies value prediction model is implemented using artificial intelligence method. This studies use the values of five indexes (i.e., impact index, fusion index per technology, fusion index per patent, diffusion index per technology and diffusion index per patent) at different time (e.g., t-n, t-n-1, t-n-2, ${\cdots}$) as input variables. The out variables are values of five indexes at time t, which is used for learning. The learning method adopted in this study is backpropagation algorithm. In the third module, this study recommends final promising technologies based on analytic hierarchy process. AHP provides relative importance of each index, leading to final promising index for technology. Applicability of the proposed methodology is tested by using U.S. patents in international patent class G06F (i.e., electronic digital data processing) from 2000 to 2008. The results show that mean absolute error value for prediction produced by the proposed methodology is lower than the value produced by multiple regression analysis in cases of fusion indexes. However, mean absolute error value of the proposed methodology is slightly higher than the value of multiple regression analysis. These unexpected results may be explained, in part, by small number of patents. Since this study only uses patent data in class G06F, number of sample patent data is relatively small, leading to incomplete learning to satisfy complex artificial intelligence structure. In addition, fusion index per technology and impact index are found to be important criteria to predict promising technology. This study attempts to extend the existing knowledge by proposing a new methodology for prediction technology value by integrating patent information analysis and artificial intelligence network. It helps managers who want to technology develop planning and policy maker who want to implement technology policy by providing quantitative prediction methodology. In addition, this study could help other researchers by proving a deeper understanding of the complex technological forecasting field.

Evaluation of the Usefulness of MapPHAN for the Verification of Volumetric Modulated Arc Therapy Planning (용적세기조절회전치료 치료계획 확인에 사용되는 MapPHAN의 유용성 평가)

  • Woo, Heon;Park, Jang Pil;Min, Jae Soon;Lee, Jae Hee;Yoo, Suk Hyun
    • The Journal of Korean Society for Radiation Therapy
    • /
    • v.25 no.2
    • /
    • pp.115-121
    • /
    • 2013
  • Purpose: Latest linear accelerator and the introduction of new measurement equipment to the agency that the introduction of this equipment in the future, by analyzing the process of confirming the usefulness of the preparation process for applying it in the clinical causes some problems, should be helpful. Materials and Methods: All measurements TrueBEAM STX (Varian, USA) was used, and a file specific to each energy, irradiation conditions, the dose distribution was calculated using a computerized treatment planning equipment (Eclipse ver 10.0.39, Varian, USA). Measuring performance and cause errors in MapCHECK 2 were analyzed and measured against. In order to verify the performance of the MapCHECK 2, 6X, 6X-FFF, 10X, 10X-FFF, 15X field size $10{\times}10$ cm, gantry $0^{\circ}$, $180^{\circ}$ direction was measured by the energy. IGRT couch of the CT values affect the measurements in order to confirm, CT number values : -800 (Carbon) & -950 (COUCH in the air), -100 & 6X-950 in the state for FFF, 15X of the energy field sizes $10{\times}10$, gantry $180^{\circ}$, $135^{\circ}$, $275^{\circ}$ directionwas measured at, MapPHAN allocated to confirm the value of HU were compared, using the treatment planning computer for, Measurement error problem by the sharp edges MapPHAN Learn gantry direction MapPHAN of dependence was measured in three ways. GANTRY $90^{\circ}$, $270^{\circ}$ in the direction of the vertically erected settings 6X-FFF, 15X respectively, and Setting the state established as a horizontal field sizes $10{\times}10$, $90^{\circ}$, $45^{\circ}$, $315^{\circ}$, $270^{\circ}$ of in the direction of the energy-6X-FFF, 15X, respectively, were measured. Without intensity modulated beam of the third open arc were investigated. Results: Of basic performance MapCHECK confirm the attenuation measured by Couch, measured from the measured HU values that are assigned to the MAP-PHAN, check for calculation accuracy for the angled edge of the MapPHAN all come in a range of valid measurement errors do not affect the could see. three ways for the Gantry direction dependence, the first of the meter built into the value of the Gantry $270^{\circ}$ (relative $0^{\circ}$), $90^{\circ}$ (relative $180^{\circ}$), 6X-FFF, 15X from each -1.51, 0.83% and -0.63, -0.22% was not affected by the AP/PA direction represented. Setting the meter horizontally Gantry $90^{\circ}$, $270^{\circ}$ from the couch, Energy 6X-FFF 4.37, 2.84%, 15X, -9.63, -13.32% the difference. By-side direction measurements MapPHAN in value is not within the valid range can not, because that could be confirmed as gamma pass rate 3% of the value is greater than the value shown. You can check the Open Arc 6X-FFF, 15X energy, field size $10{\times}10$ cm $360^{\circ}$ rotation of the dose distribution in the state to look at nearly 90% pass rate to emerge. Conclusion: Based on the above results, the MapPHAN gantry direction dependence by side in the direction of the beam relative dose distribution suitable for measuring the gamma value, but accurate measurement of the absolute dose can not be considered is. this paper, a more accurate treatment plan in order to confirm, Reduce the tolerance for VMAT, such as lateral rotation investigation in order to measure accurate absolute isodose using a combination of IMF (Isocentric Mounting Fixture) MapCHEK 2, will be able to minimize the impact due to the angular dependence.

  • PDF

Thermal Properties of Granite from the Central Part of Korea (한국 중부 지역의 화강암 열물성)

  • Kim, Jongchan;Lee, Youngmin;Koo, Min-Ho
    • Economic and Environmental Geology
    • /
    • v.47 no.4
    • /
    • pp.441-453
    • /
    • 2014
  • Thermal and physical properties were measured on 206 Jurassic granite samples obtained from three boreholes in the central part of Korea. Thermal conductivity(${\lambda}$), thermal diffusivity(${\alpha}$), and specific heat(Cp) were measured in a laboratory; the average values are ${\lambda}$=2.813 W/mK, ${\alpha}=1.296mm^2/sec$, and Cp=0.816 J/gK, respectively. In addition, porosity(${\phi}$), and dry and saturated density(${\rho}$) were measured in the laboratory; the average values are ${\phi}$=0.01, ${\rho}(dry)=2.662g/cm^3$ and ${\rho}(saturated)=2.67g/cm^3$, respectively. Thermal diffusivity of 10 granite samples were measured with increasing temperature from $25^{\circ}C$ to $200^{\circ}C$. In this study, we found that thermal diffusivity at $200^{\circ}C$ is about 30% lower than thermal diffusivity at $25^{\circ}C$. In correlation analysis, thermal conductivity increases with increasing thermal diffusivity. However, thermal conductivity does not show good correlation with porosity and density. Consequently, we know that thermal conductivity of granite would be more influenced by mineral composition than by porosity. We also derived ${\rho}=-2.393{\times}{\phi}+2.705$ from density and porosity data. XRD and XRF analysis were performed to investigate effects of mineral and chemical composition on thermal conductivity. From those results, we found that thermal conductivity increases with increasing quartz and $SiO_2$, and decreases with increasing albite and $Al_2O_3$. Regression analysis using those mineral and chemical composition were carried out ; we found $K=0.0294V_{Quartz}+1.93$ for quartz, $K=0.237W_{SiO_2}-14.09$ for $SiO_2$, and $K=0.053W_{SiO_2}-0.476W_{Al_2O_3}+6.52$ for $SiO_2$ and $Al_2O_3$. Specific gravities were measured on 10 granite samples in the laboratory. The measured specific gravity depends on chemical compositions of granite. Therefore, specific gravity can be estimated by the felsic-mafic index(F) that is calculated from chemical composition. The estimated specific gravity ranges from 2.643 to 2.658. The average relative error between measured and estimated specific gravities is 0.677%.

Error Analysis of Three Types of Satellite-observed Surface Skin Temperatures in the Sea Ice Region of the Northern Hemisphere (북반구 해빙 지역에서 세 종류 위성관측 표면온도에 대한 오차분석)

  • Kang, Hee-Jung;Yoo, Jung-Moon
    • Journal of the Korean earth science society
    • /
    • v.36 no.2
    • /
    • pp.139-157
    • /
    • 2015
  • We investigated the relative errors of satellite-observed Surface Skin Temperature (SST) data caused by sea ice in the northern hemispheric ocean ($30-90^{\circ}N$) during April 16-24, 2003-2014 by intercomparing MODerate Resolution Imaging Spectroradiometer (MODIS) Ice Surface Temperature (IST) data with two types of Atmospheric Infrared Sounder (AIRS) SST data including one with the AIRS/Advanced Microwave Sounding Unit-A (AMSU) and the other with 'AIRS only'. The MODIS temperatures, compared to the AIRS/AMSU, were systematically up to ~1.6 K high near the sea ice boundaries but up to ~2 K low in the sea ice regions. The main reason of the difference of skin temperatures is that the MODIS algorithm used infrared channels for the sea ice detection (i.e., surface classification), while microwave channels were additionally utilized in the AIRS/AMSU. The 'AIRS only' algorithm has been developed from NASA's Goddard Space Flight Center (NASA/GSFC) to prepare for the degradation of AMSU-A by revising part of the AIRS/AMSU algorithm. The SST of 'AIRS only' compared to AIRS/AMSU showed a bias of 0.13 K with RMSE of 0.55 K over the $30-90^{\circ}N$ region. The difference between AIRS/AMSU and 'AIRS only' was larger over the sea ice boundary than in other regions because the 'AIRS only' algorithm utilized the GCM temperature product (NOAA Global Forecast System) over seasonally-varying frozen oceans instead of the AMSU microwave data. Three kinds of the skin temperatures consistently showed significant warming trends ($0.23-0.28Kyr^{-1}$) in the latitude band of $70-80^{\circ}N$. The systematic disagreement among the skin temperatures could affect the discrepancies of their trends in the same direction of either warming or cooling.