Monitoring a Landslide with High Accuracy by Total Station: A DTM-Based Model to Correct for the Atmospheric Effects

: For the monitoring of large landslides, total stations equipped with an Electronic Distance Meter (EDM) are widely used. To obtain the atmospheric parameters, required along the line of sight of every measure, the data collected by a weather station close to the instrument are usually adopted. Even after these corrections, the results obtained in the monitoring of areas with complex topography don’t reach the accuracies theoretically attainable by the high-end instruments. The article proposes a method for removing the errors due to the inﬂuence of microclimate on the measurements obtained by a high-end EDM, in order to get the maximum accuracy obtainable from such instruments. The method is based on an atmospheric model, set up by using the climatic data and a digital terrain model (DTM) of the landslide area. The methodology has been applied to a landslide in southern Italy. Over 38,000 distances, acquired for each monitored point, were used. The results demonstrate the effectiveness of the method: the standard deviations of the distances after their correction, show a reduction, ranging from 20% to 50%, with respect to the most diffused procedures; furthermore, the obtained accuracy equals the one declared by the manufacturer of the instrument for measurements in optimal conditions.


Introduction
The techniques of geomatics are widely used for monitoring landslides. The monitoring involves the points inside the whole body and along the crown and the flanks of the landslide. It is in general extended to the structures situated near the crown and the ridges, and it is performed by using areal and punctual sensors. Among the areal techniques, the most commonly used is the terrestrial laser scanner [1][2][3]; Synthetic Aperture Radar (SAR) interferometry is used both from satellites [4] and from the ground [5]. Among the punctual sensors, Global Navigation Satellite Systems (GNSS) receivers have been used for several years and recent applications also use low-end instruments [6]. Low cost sensors, based on the measurements of position and tilt, are often used for activating early warning procedures [7].
In spite of the increasing use of new technologies, total station-the integration of a theodolite and an Electronic Distance Meter (EDM)-remains a fundamental instrument for monitoring [8][9][10][11][12][13]. Accurate angle and range measurements are useful when describing the surface of the ground, measuring the displacements of selected points and evaluating morphological evolutions. The accuracy of EDMs, above all if a continuously operating total station can be used, effectively contributes to early detect the trigger of a landslide and to activate warning procedures.
An accurate knowledge of the refractive index of air is essential for length interferometry or geodetic surveying; for this purposes in recent decades, several authors proposed formulae [14,15]. altitude is 250 m above sea level. The climate is characterized by hot, dry summers and rainy winters. The average annual temperature in 2012 was about 14.9 • C. The prevailing winds are from the West and North-West. On 15 February 2010 at 14:30 hours local time, the failure stage of the landslide started after a long, slow deformation phase, involving mainly the middle sector of the slope, close to the west side of the town, where the main access road was bisected. This landslide, widely described in several articles [23][24][25][26][27][28], can be considered as a first detachment phenomenon [24].
During the night between 20 February and 21 February 2010, a second slide occurred, resulting in a retrogression of the central part of the head scarp of about 80 m, with a width of about 200 m and a depth of 25 m. The landslide has an area of about 0.3 km 2 ; the width of the crown is about 500 m and the distance from the top of the crown to the tip is about 1400 m. The volume of the displaced material has been estimated at 8.0 million cubic meters [25,26]. Presently, the landslide is moving at a very low rate (mm/month) [23,27].
The landslide zone is characterized by a concave shape. The main scarp and the right flank have a constant height of approximately 50 m; the height of the left flank, beyond which the built zone is located, decreases to the south-east. The prevailing winds are sheltered by the coastal range and by the hills surrounding the main scarp. These characteristics cause the presence of calm air layers at the bottom of the landslide; in this zone the daily climatic extremes are wider, especially in the summer, i.e., temperatures are lower at night and much higher during the day. In fact, since the circulation is limited by the slope, the air becomes stagnant; it is heated by insolation during the day and cooled by terrestrial radiation during the clear nights, in sharp contrast to the wind-exposed, well-mixed air layers of the upper slopes, where heat transfer due to convection is preeminent. The limited presence of vegetation, in addition, contributes to further increase the heat exchange by radiation between air and soil and reduces the so called thermal flywheel effect, i.e., the smoothing of the average air temperature due to the long term fluctuation of the ground temperature.
By observing the head scarp of the landslide, it is possible to recognize the youngest part of the stratigraphic succession, also detectable on the geological map of Italy [28] and confirmed by the geotechnical investigations performed in the area [29].
According with [25], below the vegetal soil, we can observe (from top to bottom) the following rock strata ( Figure 2): (1) Pleistocene eluvio-colluvial deposits of a typical red color with a thicknesses of 0-10 m; (2) Middle Pliocene deposits of a light grey or brown color with a thicknesses of 5-25 m; (3) highly porous Upper Miocene evaporitic limestone, with a thicknesses varying between 25 and 45 m. The rock strata reduce in thickness to the west. Taking into account the above described composition, and the interpretations of lithology and stratigraphy from borehole data [29], we can evaluate that the soil in the landslide zone has the following characteristics [30,31]: thermal conductivity k = 1.0 [W/m K], density ρ = 1600 [kg/m 3 ], heat capacity c = 1100 [J/kg K].

The Total Station Monitoring: Layout and Instruments
Landslide monitoring using total station is generally performed by using three components: (1) a reference network materialized with prisms fixed in a stable position that can be observed from the instrument, (2) one or more stations established on stable ground at locations from which the landslide surface is visible, (3) a group of monitoring prisms at the likely unstable slope zone or area of interest.
By following the suggestions of the geologists, a permanent station was installed on a balcony of a building, located on the opposite side of the landslide, from which it is possible to observe almost the entire crown of the landslide and the town (Figure 3). The permanent station is positioned on the right side of the Scotrapiti ditch, near the border but outside the DSGSD. Two points, outside the landslide area and visible from the total station, are used for referencing aims and to control global movements. These reference points are periodically surveyed by GNSS receivers. The stability of the station was confirmed by the data collected by a Topcon GPT-8003 total station, which had operated in the same location for eight months after the failure stage.
The topographic monitoring of the landslide involves 20 points, 12 of which are located along the crown and eight on some surrounding buildings. On each point, a cube corner prism was installed for the execution of automated measurements. Along the crown of the landslide, the prisms were positioned and installed on iron bars, about one meter in height, grounded with concrete and braced; the prisms on the walls of the surrounding buildings were fixed by raw plugs. Figure 4 shows the layout of the total station monitoring; point S is the station, black points are located along the landslide crown, blue points are prisms fixed on buildings; yellow squares are points periodically surveyed by GNSS receivers, I1 and I2 are two boreholes which house inclinometers. In accordance with the geologists' indications, a prism was positioned on a building, one km away, to have a stable reference (point E1). A second stable point E2, facing the same side of the monitored buildings, is used as control point, to evaluate the stability and effectiveness of the total station.   Figure 5 shows the skyline observed from the total station: we can see that the height declines to the Scotrapiti River. This is confirmed by the present cross-sections shown in Figure 6, where the distances and the heights on the left side are referred to the total station, while the dashed line is horizontal.  For the continuous monitoring, a total station Leica TS30 was used, connected to a panel for the remote control and for data exchange. The instrument was positioned on a support consisting of an iron bracket specially designed, which was affixed on a wall of the building. For the protection against wind and rain, the bracket is equipped with a Plexiglas cover, which does not interfere with the measurements. The main characteristics of the total station are summarized in Table 1. The panel for remote control and data exchange houses: (1) a processing unit; (2) a data transmission module; (3) a device for electrical protection; (4) a switch-box controllable remotely for the restart of the total station in case of malfunction; and (5) an uninterruptible power supply to ensure proper functioning in the absence of electrical network for a time varying from 3 to 6 h.
A weather station was positioned close to the total station; it allows to obtain atmospheric parameters (pressure, relative humidity, and temperature). All measurements are performed every 20 min and data is sent to the Geomatics Lab of the University of Calabria for processing.

Analysis of the First Results
In order to fully utilize the performances of EDM, the atmospheric parameters obtained by our weather station are used. The refractive index of air can be obtained by using known models [15,32]. In our case, to correct the EDM measurements for the atmospheric effects, the following formula [18] based on the Barrel and Sears model [33] and obtained for a wavelength of 658 nm was applied: where ∆D is the atmospheric correction in parts per million (ppm), p is the atmospheric pressure in mbar, h r is the relative humidity, T is the temperature in • C, α = 1/273.15 and x = (7.5T/(237.3 + T)) + 0.7857. Given a raw measure D raw , the weather corrected distance D c is therefore: Taking into account the climatic data provided by the Regional Agency for the Environment of Calabria Region [34], it is reasonable to adopt the same values for atmospheric pressure and a variation of the relative humidity less than 15% along the path followed by the EDM laser beam in the whole landslide area; in this case, we obtain a correction of almost 1 ppm per degree of temperature increment. A commonly used procedure takes advantage of the use of some stable points well distributed, but in our case we cannot consider stable points in the town and on the left flank of the landslide; in this way, the only reference points could be E1 and E2, decentered with respect to the others. Figure 7A shows the raw values of the slope distance from station to E1. The red lines represent the results after the application of a 300-point (four days) moving average filter. Figure 7B,D shows the slope distances from station to E1 and E2 after applying Formula (2). Figure 7C is an enlargement of Figure 7A. We can observe the same oscillations of the results: actually, the terrain cross sections are very similar. Thus E1 was chosen as reference point and E2 as control point.
To calibrate the monitoring system, the average weather corrected value obtained for the first year of acquisition for point E1 was computed and was chosen as reference value. Following the procedure, for every measure, the residual is transformed in part per million (ppm) of the distance: where D 1raw is the raw distance to E1 and D 1fix is the reference value. The correction is then applied to the other measured distances; therefore, the corrected distances D j are: where D jraw is the raw measure from total station to target j. A similar method can be applied to obtain the correction of the vertical angles, used for height measurement [26]. Figure 8 shows the slope distance from station to point E4 (a) and E9 (b) after applying the correction in ppm obtained by assuming fixed the distance to E1. Despite the correction, we can observe both short and long term fluctuations. Short oscillations have an identical behavior: the average amplitude is equal to ±0.7 mm for E4 (1.4 ppm) and to ±1 mm for E9 (1.4 ppm). Long term fluctuations are quite different. If we consider the smoothed values (red lines) we can observe that the amplitude is equal to ±1.3 mm for E4 (1.8 ppm) and to ±0.6 mm for E9 (0.8 ppm).   Figure 9 shows the daily temperature variations and the slope distance from station to point E9. We can observe that the amplitude of the short oscillations is correlated to the daily temperature ranges.   The nearly linear relation between amplitude of short oscillations and daily temperature range is shown in Figure 11 for points E4, E5 and E9. The regression coefficients are 0.137 for E4, 0.140 for E5 and 0.147 for E9. For the whole data, the regression coefficient is equal to 0.141 ppm/ • C.
To obtain quantitative values for the period of short and long term fluctuations, we applied Fast Fourier Transform to the values obtained after applying the correction in ppm derived by assuming fixed the distance to E1. In this way, the searched periods have been detected, corresponding to the annual and daily oscillations: a first peak is obtained for the frequency of 3.307 × 10 −8 Hz, corresponding to a period of about 350 days; a second peak is detected for the frequency 1.159 × 10 −5 Hz, which corresponding period is 24 h. The analysis of the results obtained after applying the correction derived by fixing the distance of target E1, suggests the following remarks: (1) After the correction, the distances show both annual and daily oscillations, almost sinusoidal; this is confirmed by the results of the Fast Fourier Transform; (2) The amplitude of the annual oscillations is maximum for target E4, it decreases in the east, almost following the skyline as seen from the station; (3) An offset of about three months can be detected between maximum values of annual oscillations and equinoxes; (4) It is possible to observe a correlation between the oscillation of the distances and the daily temperature ranges; their ratio is almost constant and equal to 0.14 ppm/ • C. The beginning of the daily oscillation is at 20 o'clock.
The methods based on distance ratios and/or local scale parameters work best when all trajectories connecting the total station to the monitored point's lines have similar longitudinal profiles. This is not true in general. In our case, the use of the distance ratios introduces a time varying error because of the final 25-30% of the measurement line S-E1 which runs close and almost parallel to the ground which is not the case for the monitoring sites E4 to E9. Furthermore, the daily temperature variations are probably connected to wind and total irradiance variations during the day, which drives the amount of error introduced by the final portion of S-E1. The highly time-varying component of the atmospheric effects is taken into account by the only local scale factor (and for this aim the current meteo measurements are not needed), the spatial variability mainly depends on height above ground because the soil/surface is very homogeneous, humidity variations are negligible because of relatively dry air and/or low temperatures, and the difference between temperature close to the ground and up to 80 m above ground is well reflected by the assumed soil radiation model and by the logarithmic temperature variation with height.
The previous considerations suggest a methodology for improving the accuracy of measurements, by using a temperature model derived by taking into account the geometric characteristics of the landslide, with a peculiar attention to the DTM. Figure 12 shows the flow chart of the proposed algorithm. Starting with the climatic data, the solar radiation values on a horizontal surface and by using a 3D TIN model, the average daily solar energy incident on the ground in the landslide area is computed. A computer code allows to model the ground surface temperature in the landslide zone. The maximum and minimum temperatures are then obtained for each month.
Once computed the maximum and minimum annual temperatures on the soil surface, and by taking into account a convective heat transfer coefficient depending on the laser beam azimuth, the layering of the air temperature in the landslide zone is obtained. The average temperatures along the trajectories connecting the total station to the monitored points are then computed, along with the maximum and minimum ppm corrections. For all trajectories, day by day and hour by hour, the ppm corrections are computed, by using a sinusoidal variation and the daily temperature ranges. These corrections are applied to the distances previously modified using the correction in ppm obtained by assuming fixed the distance to E1.
To apply the proposed methodology, the ground surface temperature and the layering of the air temperature must be obtained.

The Solar Radiation in the Landslide
Solar radiation and air temperature, as well as wind and rain, are the most important meteorological factors that influence the surface and subsurface temperature [19]. The heat exchange between the atmosphere and the ground depends also on vegetation and on geometric characteristics of the terrain. Slope orientation and inclination determine the angle between solar rays and normal to the surface and, consequently, the amount of solar power incident on the ground and the net radiation balance [35].
Below the ground surface, almost all thermal energy is transferred by conduction; thus the most important properties to use are heat capacity, thermal conductivity and soil latent heat [36,37]. The wind speed is the most important information for modeling the convection heat transfer between air and ground and between air layers. The effect of evaporation may not be considered in our case, given that there is no vegetation in the landslide.
Electromagnetic radiation is the most important factor in the total energy budget. During the day, the shortwave radiation, due to the sun, is dominant, and the incoming energy, for a lightened surface, is higher than the outgoing one; at night the energy balance is dominated by long-wave radiation [19] and the outgoing energy, most of all in clear nights and in absence of dense vegetation, is greater than the incoming one.
The values of the main atmospheric parameters for the zone of Maierato are reported in Table 2. Data are the average values of the last 40 years and was provided by the Regional Agency for the Environment of Calabria Region's historical database [34]. It should be underlined that the values of solar energy are referred to a horizontal surface with a flat skyline; slope and the orientation of the terrain along with the shadowing due to the surrounding elements (buildings, hills, vegetation, etc.) can be taken into account by using a DTM of the landslide zone.
To obtain the post-failure DTM, a terrestrial laser scanner survey has been performed; the instrument used is the LMS-Z420i manufactured by RIEGL GmbH, Horn, Austria, with a measurement range of about 1000 m and an accuracy of 10 mm. Figure 13 shows the 3D model of the present situation superimposed to an aerial photo: due to the landslide, a zone in the form of a basin was created. Among the several formulas available in literature to obtain the daily solar energy H incident on a tilted surface, we used the one suggested by UNI, the Italian Organization for Standardization [38]. The known elements are: (a) the global and the diffuse daily solar energy on a horizontal surface (H h , H d ) attainable by Table 2, (b) the latitude φ of Maierato (38 • 42 16" N), (c) the angle β the surface makes with the horizontal, (d) the surface azimuth γ, defined as the deviation of the projection on a horizontal plane of the normal to the surface from the local meridian, with zero due south, east negative, west positive, (e) the solar constant G o = 1,353 W/m 2 , (f) the sun's declination angle δ, (g) the albedo ρ (reflection coefficient of the soil surface) that, in our case, can be considered equal to 0.2, (h) the sun hour angle at sunrise for a horizontal surface s , (i) the sun hour angles at sunrise and sunset for the considered surface and ". The daily solar energy H incident on a tilted surface is given by: To obtain the coefficient R, the geometric characteristics of the area are taken into account. Topographic shading of direct radiation is accounted for and diffuse radiation is modified by the proportion of sky in view: the irradiation deficit due to shadowing is taken into account through the angles and ". Since in the landslide area there are not artificial shadings like in the towns, the actual sunrise and sunset times give sufficient information. Furthermore, for our aims we do not need a very detailed microclimate model, so a simplified method can be used.
For β = 0, the coefficients T h , U h and V h can be obtained. The direct solar energy on the tilted surface H b and the direct solar energy on a horizontal surface H bh can be computed, along with the ratio R b : Finally, the coefficient R is given by: (12) Figure 14 shows the instantaneous solar lighting in the landslide and in the surroundings and the average daily solar energy incident on the ground for different dates and times, obtained by using the data of Table 2 and the TIN model. Red lines are the paths of the EDM laser beams, starting from the station on the hill. In the first row, from left to right, we can observe the radiation incident upon the ground on 21 December, winter solstice, for different times (7:30, noon and 16:30); the light areas are lit by the sun, the dark ones are in shadow; the fourth image represents the daily solar energy collected by the ground in a color scale, ranging from 5 to 26 MJ/m 2 day. In the second and third rows, the radiation and daily energy values are shown on 21 March, the vernal equinox, and on 21 June, summer solstice.
We can make the following remarks: (1) during winter, when the solar angle is minimum, the zone where the targets E4-E9 are located, gets solar energy only a few hours, due to the low sun angle and to the shading caused by the hill where the total station is positioned. The prevailing winds (SW, SE) are shielded by the same hill. Thus, the temperature on the soil and near the soil is lower than outside the landslide, during the day, due to the low solar radiation and at night (above all in clear nights) due to the long-wave radiation and the absence of wind. (2) During summer, the mentioned zone gets almost the same solar energy as the external one, due to the high sun angle. The prevailing winds (SW, NW) are shielded by the same hill and by the main scarp of the landslide. Thus, the temperature on and near the soil is higher than outside the landslide due to the greater solar energy gain in absence of vegetation and wind. (3) Furthermore, by taking into account the skyline (Figure 4) and the cross sections ( Figure 5) we can observe that the shielding effect due to the hill decreases from E4 to E9 and an increasing convective heat exchange must be foreseen.

Modeling the Ground Surface Temperature in the Landslide Zone
A software program was developed by using Matlab ® software and the Simulink ® libraries. The calculations are based on the following parameters' input: thermal conductivity, density, heat capacity, air temperature, relative humidity, solar radiation, and convective heat transfer coefficient. For the first three parameters, the values reported in Section 2.1 were used. The daily air temperature and the relative humidity were obtained by interpolating the values of Table 2, while the solar radiation was computed as previously described, and it was considered constant every month. During the night, the radiative cooling was taken into account. The net outgoing radiation varies from about 15 W/m 2 for an overcast night, to 100 W/m 2 for a clear sky night [39], so we adopted values ranging from −45 W/m 2 in January (40% of clear sky nights) to −87 W/m 2 in July (90% of clear sky nights). A convective heat transfer coefficient equal to 6 W/m 2 · K was chosen for the calm zone; the coefficient increases with the azimuth up to 24 W/m 2 · K, going from the direction to E4, to the direction to the bottom of the Scotrapiti River.
The procedure described in detail in [40] was applied. The soil model make use of two main equations for heat and moisture transfer within soil. For heat conduction the one-dimensional diffusion equation is: ∂T ∂t where ν, c s , and ρ s , are the thermal conductivity, specific heat capacity, and soil density, respectively (obtained in the section geological settings). The expression ν/ρ s c s = k s is called the thermal diffusivity.
In several cases, k s can be considered constant along z-axis [41]; the equation then becomes: The vertical one-dimensional equation of diffusive and gravitational motions for soil moisture transfer is: ∂n ∂t = ∂ ∂z D n ∂n ∂z + ∂K n ∂z (15) where D n and K n are diffusional and hydraulic conductivity, respectively. The heat budget equation for the interface-spanning layer with average temperature T g is written as: where R n is net radiation, H is sensible heat flux, L ν E is latent heat flux, and G is soil heat flux, G = ν∂T/∂z. (Radiation fluxes are defined positive toward the surface, other fluxes positive away from the surface.) The software program uses the finite difference method. The maximum and minimum temperatures have been obtained, respectively, in July and January.

Average Temperature Along the Paths of the Laser Beams and ppm Corrections
According with the literature [19], it is reasonable to hypothesize that the air temperature varies with the height in a binary logarithm way. Thus, for our simulation, we can consider air layers starting from each target, characterized by an average difference of temperature of one • C and by doubled thicknesses. Figure 15 shows, in a range of colors, the difference between air temperatures in the landslide zone and external air temperature calculated for July. The values range from 0 • C (blue) to 5 • C (red). The heights on the left side refer to the total station S, the dashed line is the horizontal, the yellow line is the EDM laser beam path. For target E4, the layers are shown, separated by white lines. To obtain the height of the layers, we use the following procedure. We define: ∆T j = temperature difference between the target E j and the total station; ∆H j = height difference between the target E j and the total station; n = the integer part of ∆T j ; We divide the atmosphere between each target and the total station into n + 1 layers. For the thickness of the ith layer of the jth trajectory (starting from the soil and for i = 1 to n), we can state the proportionality: For the thickness of the last layer we can write: The total thickness, equal to the height difference between the target E j and the total station, is: The thickness of the ith layer is given by: The thickness of the last layer is: If we neglect the curvature of the laser beam of EDM, the path from the total station to the target E j can be divided in n + 1 segments, whose lengths are directly proportional to the thickness of the relevant layer. The average temperature of the atmosphere along a laser beam path is then given by: where T Ej is the temperature of the air at the point E j . The difference of the average temperatures between the path to Ej and the path to E 1 is given by: Actually, the air layers close to the ground follow the topography. Thus, for each layer, the length of the beam path depends not only on the thickness, but also on its angle of incidence, that has been neglected in our model. Nevertheless, this simplification is justified if we consider that: (a) even for the lower layers, the separation surfaces are roughly horizontal in most cases; (b) the lower layers have a very limited thickness compared to the higher ones. Therefore, a more accurate modeling of the laser beam path could be requested just for critical conditions, e.g., for a signal path almost parallel to a very steep slope.
As stated in Section 2.3, it is reasonable to consider a correction of 1 ppm per degree of temperature difference. The absolute values of maximum correction A j (in ppm), obtained with the procedure described previously, are reported in Table 3. The values range from 2.10 for the target E4, to 0.65 for the target E9. If we do not consider the daily oscillation of air temperature, the correction for the slope distance is sinusoidal and it ranges between the values A j and −A j , given by Table 3 for each target. Given that the annual oscillations begin around 15 September (close to the fall equinox), the correction for the Nth day of the year, expressed in ppm, is given by: where: A j = maximum absolute value of ppm correction of the slope distance from total station to E j ; N = day of the year; Day 0 = day of the beginning of the yearly oscillation (258 = 15 September). Taking into account the daily ranges of air temperature and the correlation described in Section 2.3, the correction for the slope distance, for the H th hour of the N th day of the year, expressed in ppm, is given by: where: H r = hour of the day; H r0 = 20, hour of the beginning of the daily oscillation; ∆T N = temperature range of Nth day. The distances corrected are:

The Atmospheric Correction of the EDM Measurements
To obtain the atmospheric correction for the EDM measurements, the following steps are followed: (a) We compute the maximum and minimum annual temperatures on the soil surface, by taking into account low-medium values of the ground thermal diffusivity and a convective heat transfer coefficient, which depends on the laser beam azimuth; (b) We consider calm air and the stratification described in Section 2.6; (c) Taking into account the maximum and minimum temperatures of the soil surface, we compute the average temperatures along the trajectories connecting the total station to the target E1 and to targets E4-E9 using Equation (22); (d) We calculate the maximum and minimum temperature differences between the trajectory to E1 and the trajectories to E4-E9 using Equation (23) and, consequently, the relevant maximum and minimum ppm corrections; (e) For all trajectories we compute, day by day and hour by hour, the ppm corrections using Equations (24) and (25). We compute distances corrected using Equation (26). Figure 16 shows the corrections, expressed in ppm, of the slope distances, computed for target E4. In Figure 16a the corrections for the whole monitoring period are represented; in Figure 16b the corrections for a period of 30 days, circled in red in Figure 13a, are shown, along with the daily range of external air temperature.

Results and Discussion
In the following, we can observe a comparison between the results obtained by applying the corrections in ppm deduced by fixing the distance to E1, and the results obtained by adopting the methodology previously described. The procedure has been applied to the measurements carried out over a 15 month period, from June 2012 to September 2013, during which over 38,000 distances have been acquired for each target. Figure 17 shows the results obtained for target E4 (a,b); the red lines represent the results after the application of a 300-point moving average filter. The values of slope distances for the whole monitoring period are shown; some power failures caused the loss of measurements for short periods, easily detectable in the graphs. On the left we can see the results of applying the corrections in ppm obtained by fixing the distance E1: the residual seasonal oscillations are noticeable. On the right the distances obtained by adopting the methodology described in Section 2.6 are shown; two results have been obtained: (1) the long term oscillations are almost eliminated; (2) the range between maximum and minimum values is dramatically reduced.   The overall results obtained during the whole monitoring period for each target are summarized in Table 4, where: Di is the average slope distance in m; σ raw is the standard deviation of the raw acquisitions in mm; σ temp is the standard deviation of the distances obtained by using Equation (1) and the atmospheric data collected by the weather station; σ E1 is the standard deviation of the distances with the corrections in ppm obtained by fixing the distance to the target E1; σ tot is the standard deviation of the distances obtained by adopting the described methodology; Declared EDM accuracy is the value of accuracy, at 95% confidence level, which can be expected in ideal conditions, according to the characteristics declared by the manufacturer and reported in Table 1 (0.6 mm + 1 ppm). We can underline the progressive reduction of the standard deviation. By using the atmospheric data acquired by the weather station, the reduction is about 60%. A further reduction is obtained by fixing the distance to target E1; the improvement is less for the targets E4 and E5, which are more influenced by the layering of the air in the landslide area. By using the method described in the paper, we obtained the best results and the influence of the microclimate is almost corrected: the reduction of the standard deviation is about 90%, while the final standard deviation is about proportional to the distances and does not depend on the directions. In the last column we report the accuracy, at 95% confidence level, which can be expected in ideal conditions, according to the characteristics declared by the manufacturer (Table 1): we can observe that these values are almost equal to 2σ tot ; this confirms the effectiveness of the described methodology.
The fluctuation of the results that still remain even after the application of the proposed methodology is probably due to changes in the parameters not considered in the model. In any case, the effects of these variations are less than the stated accuracy for the total station used, so a more sophisticated model might not be very useful.
With reference to the usability and feasibility, the method is more effective when the longitudinal profiles along the trajectories connecting the total station to the monitored points are different, which can produce different meteo-climatic parameters (wind, temperature range). We observe that the methodology presented could be largely adopted with very low costs. In fact, the requested information is a DTM, along with historical climatic data. Given that the same robotic stations used for monitoring can perform automatic scanning and the recent EDMs can measure ranges up to 2000 m without targets, a DTM can be obtained easily and inexpensively. The historical climatic data can be easily found in the websites of public and private institutions.

Conclusions
A method for the removal of the errors due to the influence of microclimate in the monitoring of a large landslide by a high-end total station has been described. The method is based on the realization of an atmospheric model by using climatic data and the DTM of the landslide area. The novel aspect of the study lies in the possibility to obtain high EDM accuracy without the necessity to install several weather stations distributed in the whole landslide area.
The methodology has been applied to the large landslide at Maierato, southern Italy. The results of the test carried out demonstrate the effectiveness of the method. The standard deviations show a reduction, ranging from 20% to 50%, with respect to the results obtained by using the most diffused procedures, still affected by both long and short term oscillations. The more evident improvements are obtained for the points situated inside the landslide zone, characterized by a basin shape. The obtained standard deviations are is in line with the values declared by the manufacturer for measurements in optimal conditions. The methodology presented could be largely adopted, above all when high accuracies are requested (monitoring of structures, cliffs, etc.), with a low cost.