SciELO - Scientific Electronic Library Online

 
vol.39 issue1Study of coagulating/flocculating characteristics of organic polymers extracted from biowaste for water treatmentSupport for incident management in optical networks through critical points identification author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Ingeniería e Investigación

Print version ISSN 0120-5609

Ing. Investig. vol.39 no.1 Bogotá Jan./Apr. 2019  Epub Nov 05, 2019

https://doi.org/10.15446/ing.investig.v39n1.71670 

Original articles

Assessment of groundwater level variations using multivariate statistical methods

Evaluación de cambios en el nivel freático mediante métodos estadísticos multivariados

Fausto Molina-Gómez1 

Lenin A. Bulla-Cruz2 

Luis Á. Moreno-Anselmi3 

Juan C. Ruge4 

Carol Arévalo-Daza5 

1 Civil Engineer, Universidad Militar Nueva Granada, Colombia. M.Sc., Universidad de los Andes, Colombia. Affiliation: Universidad Militar Nueva Granada, Colombia. E-mail: fausto.molina@unimilitar.edu.co

2Civil Engineer, Universidad Militar Nueva Granada, Colombia. M.Sc., Universidad Nacional de Colombia, Colombia. Affiliation: Universidad Nacional de Colombia, Colombia. E-mail: labullac@unal.edu.co

3 Civil Engineer, Universidad Francisco de Paula Santander, Colombia. M.Sc., Universidad de los Andes, Colombia. Affiliation: Universidad Militar Nueva Granada, Colombia. E-mail: luis.moreno@unimilitar.edu.co

4 Civil Engineer, Universidad Francisco de Paula Santander, Colombia. Ph.D., Universidade de Brasilia, Brasil. Affiliation: Universidad Militar Nueva Granada, Colombia. E-mail: juan.ruge@unimilitar.edu.co

5 Civil Engineer, Universidad Francisco de Paula Santander, Colombia. M.Sc., Universidad de los Andes, Colombia. Affiliation: Universidad Militar Nueva Granada, Colombia. E-mail: carol.arevalo@unimilitar.edu.co.


ABSTRACT

Fluctuation of groundwater level induces changes in pore-water pressure of soil. However, this variation is not considered for underground constructions. This article explores the application of a statistical method to evaluate the groundwater level variation in geotechnical designs. The methodology included: (i) data collection, (ii) statistic formulation, and (iii) statistic data analysis. We collected information from the technical studies of the project "Metro de Bogotá", and selected four boreholes spanning 160 m, approximately, where the 1o de Mayo metro station will be built, in the south of the city. We used groundwater level readings reported by different piezometers for 30 days and data variance was assessed using a multivariate statistical method: analysis of repeated measures profiles. Results present a procedure to estimate the groundwater level fluctuation during a short monitoring period. We concluded that the analysis of repeated measures profiles allows estimating the groundwater level variation under a significance level 1-α.

Keywords: Analysis of repeated measures profiles; boreholes; infrastructure projects; Metro of Bogotá city

RESUMEN

La fluctuación del nivel freático induce cambios en la presión de poros del suelo. Sin embargo, esta variación no se contempla en construcciones subterráneas. Este documento explora la aplicación de un método estadístico para evaluar la variación del nivel freático en diseños geotécnicos. La metodología incluyó: (i) recolección de datos, (ii) formulación estadística y (iii) análisis estadístico de datos. Se recopiló información de los estudios técnicos del proyecto "Metro de Bogotá". Se seleccionaron cuatro sondeos que abarcan 160 m, aproximadamente, donde se construirá la estación 1 o de Mayo, al sur de la ciudad. Se utilizaron lecturas de nivel freático reportadas por varios piezómetros, durante 30 días y la variabilidad de los datos se evaluó utilizando el método estadístico multivariado: análisis de perfiles de medidas repetidas. Los resultados presentan un procedimiento para estimar la fluctuación del nivel freático durante un período corto de monitoreo. Se concluyó que el análisis de perfiles de medidas repetidas permite estimar la variación del nivel freático bajo un nivel de significancia 1- α.

Palabras clave: Análisis de perfiles de medidas repetidas; Metro de la ciudad de Bogotá; proyectos de infraestructura; sondeos

Introduction

Geotechnical explorations often require to register the position of the water along the soil profile. This registration is a parameter known as groundwater level. The identification of groundwater level allows calculating the effective stress. Likewise, such reference point is used to establish the possible drainage conditions of the soil in the structure design. Moreover, groundwater level depends on the hydro-geological conditions of the ground (Gonzalez de Vallejo and Ferrer, 2011) and its variation is conditioned by the weather and the hydraulic properties of soil in the hydraulic parameters of the soil layers (Ruge, Da Cunha, Colmenares, and Mendoza, 2017).

However, in many cases, the variability of this state is not contemplated as part of the design process. Osterberg (2004) states that better exploration and sample practices are necessary to improve the building quality. Hence, in order to avoid extra laboratory tests, the groundwater level variation must be included in the designs of underground constructions (Tristá, Sotolongo, Cristía, and Fernández, 2016).

The variation of water content affects the stiffness of the material and increases its strains during reload stages (Molina-Gómez, Camacho-Tauta, and Reyes-Ortiz, 2016). In addition, fluctuations in the groundwater level induce variations in the lateral pressure, and those changes can be considered in the design of retaining wall construction (Ruge, 2014). Ausilio and Conte (2005) affirm that the groundwater level position is an issue when computing the bearing capacity of shallow foundations. Therefore, the groundwater level may affect the stability of geotechnical structures, such as tunnels.

Nevertheless, the study of the effect of the groundwater level variation has acquired relevance, especially, in the slope stability analysis. Reddi and Wu (1991) and Cascini, Calvello, and Grimaldi (2010) proposed a model to derive the time-dependent shear strength along the main slip surfaces. Conte and Troncone (2011) developed a method, based on a simple sliding-block model, to estimate the probability of failure in slopes, induced by the increment of pore water pressure during the rising of groundwater level. During seismic events, water under the surface controls the saturation degree, which can affect the soil strength due to the liquefaction phenomena (Soares and Viana da Fonseca, 2016).

The fluctuation of the groundwater level can be assessed, mainly, by statistical methods. Zhao, Li, Zhang, and Wang (2016) used a regression model to calculate the position of the water table. They validated the results using field measurements and found that the equation can predict the water level variations, with good precision. Han et al. (2016) implemented a groundwater level modelling framework through the coupling of two spatial and temporal clustering techniques. In addition, their procedure used self-organizing map technique to identify spatially homogeneous clusters of groundwater level piezometers. Yoon et al. (2016) predicted the long-term groundwater level fluctuation using a time series model and artificial neural network to evaluate the effect of rainfall on the soil.

In this study, we applied the analysis of repeated measures profiles to the estimation of the groundwater fluctuation. The remainder of this paper has four sections. The first section corresponds to the compilation of theoretical background and the description of the statistical method. The second section shows the data source and describes the soil composition of the study case. The third section refers to the validation of the null hypothesis of the technique. Finally, the fourth section presents the analysis and conclusions of this research.

Repeated measures profiles

Multivariate analysis allows to solve problems based on analytical criteria, which include all the variables involved (Gatingon, 2013). Hence, those methods can help to interpret better any type of information (Johnson and Wichern, 2007). Statistics methods do not only analyze numerical information, but they can also interpret graphical data with a quantitative approach. One of those methods is the repeated measures profiles, which analyzes graphical results.

Through a repeated measures profiles design, it is possible to estimate the variation of a response variable subjected to different treatments. As stated by Tabachnick and Fidell (2013), the data set may come from a dependent variable measured several times under the same pattern, i.e. the same independent variable is considered. This technique is an application of multivariate analysis of variance (MANOVA), where all samples (n) are measured within a fixed or constant variable. It focuses mainly on the comparison of variances, hence, the mean vectors of an specific treatment measured at the same level are compared (Friendly, 2010).

We evaluate three null hypotheses, according to Davis (2002). Those hypotheses are parallelism, flatness and coincidence. Harrar and Kong (2016) affirm that the technique seeks to respond the following questions: (i) whether there is interaction effect between-subjects and within subject factors, (ii) whether there is a between-subject factor effect, and (iii) whether there is a within-subject factor effect. Figure 1 presents a graphical representation of the null hypothesis.

Source: Authors

Figure 1 Null hypotheses: (a) parallelism; (b) flatness; (c) 

Mathematically, the profiles of repeated measures analyze the variance or covariance of the data (Johnson and Wichern, 2007). Therefore, the method compares the matrix of the slope parameters. Timm (2004) suggested Equations (1-3) to describe the null hypotheses (H0).

where μ represents the media of the measurement in each repetition; ρ the number of the measurement and i the number of the profile.

Moreover, all the hypotheses validation has been programmed in several statistical software, due to the amount of information collected. Through those routines, it is possible to calculate the variance and significance level between data. Those criteria can be estimated via p-value, which is a parameter used to reject or not reject the null hypothesis in any statistical model (Habiger, 2015). In addition, the ρ-value provides results with a confidence level 1 - α. Therefore, if the p-value is lower than a, the null hypothesis is rejected (Wackerly, Mendenhall, and Scheaffer, 2008). However, it is important to establish previously if the hypothesis must be rejected or not. Bulut and Desjardins (2017) stated that the repeated measures profiles are parallel, flatness and coincident when the p-value < α.

Data collecting

Information comes from the geotechnical exploration of the site designated as 1o de Mayo metro station. Such station is part of the most important infrastructure project in the city: the metro of Bogotá. For this project, an exhaustive ground identification, which included 15 months of geological-geotechnical exploration was performed. The Institute for Urban Development of the city affirms that they hired the drilling of 563 boreholes in the 27 km of the metro line (IDU, 2015a). Each perforation had 50 m depth and were approximately 100 meters from each other. Likewise, the studies identified the physical, mechanical and dynamical properties of the subsoil, through more than 2000 laboratory tests. The data used in this article comes from the records reported in the geotechnical study for the metro of Bogotá (IDU, 2015b) for the boreholes SE1-25, SE1-26, SE1-27 and SE1-28. Figure 2 presents the localization of the boreholes.

According to the information provided in the Decree 523 of 2010, "Microzonificación Sísmica de Bogotá", the soil of the area in Figure 2 corresponds to an alluvial material (Secretaría General de la Alcaldía Mayor de Bogotá D. C., 2010). Molina-Gómez, Moreno-Anselmi, and Arévalo-Daza (2016) explain that this type of soil has medium to high load-carrying capacity, low compressibility, medium liquefaction susceptibility and could be unstable in open excavations. In addition, Caicedo, Mendoza, López, and Lizcano, (2018) indicate that these deposits are located in plain areas composed of loose to compacted clayey sands. Table 1 shows the soil lithology of the zone based on the classification proposed by the Unified System of Soils Classification (USCS).

Source: Authors adapted from Google Earth®

Figure 2 Location of the boreholes of the 1o de Mayo Station. 

Table 1 Soil classification and its physical properties 

Source: Authors adapted from IDU (2015b)

The zone for the 1 o de Mayo metro station is located in the area of Tunjuelo River, which is a basin in the locality of Kennedy. The underground water flow of such river moves in southwest to northeast direction along 73 km distance. In addition, the slope of this aquifer changes 15 to 3, with an average slope of 5 at the south of Bogotá. According to geological-geotechnical reports for the first line of metro (Oteo-Mazo, 2015), the piezometric level of the aquifer is between 2560 and 2540 masl. Figure 3 presents the hydrogeological profile of the 1o de Mayo Station.

Source: Authors

Figure 3 Hydrogeological profile of the 1o de Mayo Station. 

In the area for the 1o de Mayo Station construction, piezometers were positioned after drilling at the four research points for the future metro station (Figure 4). Piezometric readings were registered daily. For this research, readings at the same specific time were selected, in order to ensure repeated measures in the entire exploration site.

Source: IDU (2015b)

Figure 4 Piezometric space for groundwater level reading into borehole SE1-26. 

Results and analysis

We used RStudio to process the data, which is a free statistical software based on an object-oriented algorithm. This software allows data plotting and information modelling by several techniques. In addition, RStudio has different packages to estimate the data variance using different multivariate techniques.

In this research, we used profileR package, which was proposed by Bulut, Davison, and Rodriguez (2017). Computation procedure covers an experiment design, including calculation of variance through sums squares matrix and vector products of the Equations (1-3). This tool provides a set of multivariate methods and data visualization options to implement profile analysis and cross-validation techniques described by Bulut (2013) and Davison and Davenport (2002). Likewise, it includes routines to perform criterion-related profile analysis, profile analysis via multidimensional scaling, moderated profile analysis, profile analysis by group, and a within-person factor model to derive score profiles. In addition, it allows to compare simultaneously the effect of treatments by univariate techniques as the Hotelling's T2 test.

Furthermore, we used four different techniques to assess the parallelism of groundwater level along the boreholes. Those statistics were Pillai trace, Wilks' Lambda, Hotelling-Lawley trace and Roy's largest root. The aforementioned statistical models were described by Molina-Gómez et al., (2016). We validated the null hypothesis using the application of the procedure proposed by Bulut and Desjardins (2017), which evaluates the profiles of repeated measures with the profileR package.

Based on the exploration records, we selected thirty different measures of groundwater level during thirty different days. The measures started on March, 2014 and finished on April, 2014. During the monitoring period, there are no values of groundwater level at the surface. Topographic records indicate that all the research points and piezometers are at the elevation position, 2558 masl. Figure 5 shows the descriptive exploration of the groundwater level readings in the research points.

Source: Authors

Figure 5 Box plot ofthe measures. 

Results revealed an interaction of information between the boreholes SE1-25 with SE1-26 and SE1-27 with SE1-28. We found similar inter-quartile distribution in all measures. However, we identified an outlier in SE1-26 and SE1-28, which is the first measure in SE1-26, and a change during monitoring in SE1-28.

In order to satisfy the homogeneity of variances assumption, using the median, we performed the Levene's test before executing the analysis of repeated measures profiles. The null hypothesis of this test assesses if all the populations' variances are equal and the alternative hypothesis considers that at least two of them differ. Results showed that the groundwater variances of the profiles are homogenous, under a confidence level of 99 (α = 0,01). If data satisfy the homogeneity condition, it is valid to apply the analysis of repeated measures profiles technique. Table 2 presents the Levene's model outcomes.

Table 2 Results of the homogeneity of variance evaluation 

Source: Authors

Boreholes 3 3,6153 1,53 X 10-2

Residuals 116

Figure 6 presents the repeated measures profiles of the groundwater measurements. From a qualitative viewpoint, we observed that all the profiles remain constant along the readings. Nevertheless, we confirmed the outlier of Figure 6 and a decreasing groundwater level during day23, as well as a lower fluctuation in SE1-28. Possible reasons for the variation are the soil permeability and the material heterogeneity. According to the soil classification presented in Table 1 and the contrast with the typical values of permeability coefficient presented in the literature by Warrick (2003), it is possible to affirm that the soil has intermediate hydraulic conductivity and could experience fluctuations during short periods of time.

Source: Authors

Figure 6 Groundwater level profiles. 

Table 3 presents the analysis of the repeated measures results. We assessed the three null hypotheses of the method, under a confidence level of 99 (α= 0,01). The value of a selected is based on the experiments design in the literature by Kuehl (2000) and Ramachandran and Tsokos (2009). Since the p-value obtained is less than a, there is statistical evidence to reject all the hypotheses (Habiger, 2015; Wackerly et al., 2008) under a confidence level of 99. Therefore, results suggest that the profiles of groundwater level are flat (horizontal), but not parallel and neither coincident.

Statistical results showed that the groundwater level of the aquifer remains constant during the period of study. In addition, under the area for the station construction, the water flows through stratified zones where the piezometric readings change from a research point to the other due to the effect of the equivalent hydraulic conductivity of the layered soil (Dulcey-Leal, Molina-Gómez, and Bulla-Cruz, 2018). However, for underground constructions, like the metro, it is necessary to ensure a constant groundwater level using a set of monitoring/control procedures, in order to avoid additional pore-water pressures that can increase soil lateral stresses. In this context, flatness hypothesis (horizontality) is the most important null hypothesis, H02, since it establishes the uniformity of groundwater level during the monitoring period.

Conclusions

This article addressed a statistical assessment of ground-water variations for underground constructions. Results present a quantitative procedure to estimate the uniformity of the groundwater level position. We analyzed four different research points using data collected for the design of the infrastructure project "Metro de Bogotá". The points are continuous and cover a distance of 160 m. From results, the following conclusions were drawn:

Table 3 Validation results ofnull hypotheses and variance 

Source: Authors

  1. We obtained a profile of repeated measures. We analyzed data that came from the geotechnical studies of the "Metro de Bogotá" project, including thirty different measures of groundwater level position. By visual inspection, not all the boreholes presented variations in the groundwater level position, except in day 23. However, we identified that the plot does not indicate the degree or quantitative value of the variation of such level.

  2. We implemented a multivariate graphical statistical method, the analysis of repeated measures profiles, and assessed the variation of groundwater level for underground constructions. We evaluated the null hypotheses of the repeated measures analysis. Outcomes showed that the profiles are flat (horizontal), but are not parallel and not coincident under a confidence level of 99 (a = 0,01). Thus, we found no variations of groundwater level in a period of 30 days, which indicates no probable pore-water pressure build-up and no increments in the lateral stress of the soil during the station construction.

  3. Analysis of repeated measures profiles showed that, when assessing groundwater variations, the most important null hypothesis is H02, which determines the horizontality of the profiles. Hence, the measurements of the four boreholes have a between-subject factor effect. In this way, we found that there are no variations of the groundwater level. Therefore, we can affirm that the groundwater level is statistically uniform, even if we observed some variations graphically.

Acknowledgements

This work is part of the Research Project INV-DIS-2068 supported by the Vicerrectoría de Investigaciones of the Universidad Militar Nueva Granada.

References

Ausilio, E., and Conte, E. (2005). Influence of groundwater on the bearing capacity of shallow foundations. Canadian Geotechnical Journal, 42(2), 663-672. DOI: 10.1139/T04-084 [ Links ]

Bulut, O. (2013). Between-person and Within-person Subscore Reliability: Comparison of Unidimensional and Multidimensional IRT Models. (Ph.D. thesis, University of Minnesota). Retrieved from: http://citeseerx.ist.psu.edu/viewdoc/download?.doi=10.1.1.816.4794&rep=rep1&type=pdf. [ Links ]

Bulut, O., Davison, M. L., and Rodriguez, M. C. (2017). Estimating Between-Person and Within-Person Subscore Reliability with Profile Analysis. Multivariate Behavioral Research, 52(1), 86-104. DOI: 10.1080/00273171.2016.1253452. [ Links ]

Bulut, O., and Desjardins, C. D. (2017). Profile Analysis of Multivariate Data in R. R Library [R Package version 0.3-5]. Retrieved from https://cran.r-project.org/web/packages/profileR/profileR.pdf. [ Links ]

Caicedo, B., Mendoza, C., López, F., and Lizcano, A. (2018). Behavior of diatomaceous soil in lacustrine deposits of Bogotá, Colombia. Journal of Rock Mechanics and Geotechnical Engineering, 10(2), 367-379. DOI: 10.1016/J.JRMGE.2017.10.005. [ Links ]

Cascini, L., Calvello, M., and Grimaldi, G. M. (2010). Groundwater Modeling for the Analysis of Active Slow-Moving Landslides. Journal of Geotechnical and Geoenvironmental Engineering, 136(9), 1220-1230. DOI: 10.1061/(ASCE)GT.1943-5606.0000323. [ Links ]

Conte, E., and Troncone, A. (2011). Analytical Method for Predicting the Mobility of Slow-Moving Landslides owing to Groundwater Fluctuations. Journal of Geotechnical and Geoenvironmental Engineering , 137(8), 777-784. DOI: 10.1061/(ASCE)GT.1943-5606.0000486. [ Links ]

Davis, C. S. (2002). Statistical methods for the analysis of repeated measurements (2nd ed.). New York: Springer. [ Links ]

Davison, M. L., and Davenport, E. C. (2002). Identifying criterion-related patterns of predictor scores using multiple regression. Psychological Methods, 7(4), 468-484. Retrieved from http://www.ncbi.nlm.nih.gov/pubmed/12530704. [ Links ]

Dulcey-Leal, E., Molina-Gómez, F. A., and Bulla-Cruz, L. A. (2018). Hydraulic conductivity in layered saturated soils assessed through a novel physical model. DYNA, 85(205), 119-124. DOI: 10.15446/dyna.v85n205.64473. [ Links ]

Friendly, M. (2010). HE Plots for Repeated Measures Designs. Journal of Statistical Software, 37(4), 1-40. DOI: 10.18637/jss.v037.i04. [ Links ]

Gatingon, H. (2013). Statistical Analysis of Management Data (3rd ed.). New York: Springer Science & Business Media. [ Links ]

González de Vallejo, L. I., and Ferrer, M. (2011). Geological Engineering. London: CRC Press/Balkema. DOI: 10.1201/b11745. [ Links ]

Habiger, J. D. (2015). Multiple test functions and adjusted p-values for test statistics with discrete distributions. Journal of Statistical Planning and Inference, 167, 1-13. DOI: 10.1016/j.jspi.2015.06.003. [ Links ]

Han, J.-C., Huang, Y., Li, Z., Zhao, C., Cheng, G., and Huang, P. (2016). Groundwater level prediction using a SOM-aided stepwise cluster inference model. Journal of Environmental Management, 182, 308-321. DOI: 10.1016/j.jenvman.2016.07.069. [ Links ]

Harrar, S. W., and Kong, X. (2016). High-dimensional multivariate repeated measures analysis with unequal covariance matrices. Journal of Multivariate Analysis, 145, 1-21. DOI: 10.1016/j.jmva.2015.11.012. [ Links ]

Instituto de Desarrollo Urbano - IDU. (2015a). Estudios geológicos y geotécnicos. Unpublished raw data. [ Links ]

Instituto de Desarrollo Urbano - IDU. (2015b). Información Geotécnica PLMB Tramo I. Unpublished raw data. [ Links ]

Johnson, R. A., and Wichern, D. W. (2007). Applied Multivariate Statistical Analysis (6th ed.). New Yersey: Pearson Education International. [ Links ]

Kuehl, R. (2000). Design of experiments: statistical principles of research design and analysis (2nd ed.). Pacific Grove (Calif.): Duxbury/Thomson Learning. [ Links ]

Molina-Gómez, F. A., Moreno-Anselmi, L. A., and Arévalo-Daza, C. E. (2016). Aplicación del análisis de perfiles de medidas repetidas en exploraciones geotécnicas, Paper presented at the XV Congreso Colombiano de Geotecnia (p. 7). Cartagena, Sociedad Colombiana de Geotecnia. [ Links ]

Molina-Gómez, F. A., Camacho-Tauta, J. F., and Reyes-Ortiz, O. J. (2016). Stiffness of a granular base under optimum and saturated water contents. Revista Tecnura, 20(49), 75-85. DOI: 10.14483/udistrital.jour.tecnura.2016.3.a05. [ Links ]

Osterberg, J. O. (2004). Geotechnical Engineers, Wake-up - The Soil Exploration Process Needs Drastic Change. In J.P. Turner, and P.W. Mayne (Eds.), GeoSupport 2004: Drilled Shafts, Micropiling, Deep Mixing, Remedial Methods, and Specialty Foundation Systems (pp. 450-459). Reston, VA: American Society of Civil Engineers. DOI: 10.1061/40713(2004)38. [ Links ]

Oteo-Mazo, C. (2015). Condiciones geotécnicas del Proyecto de la Línea 1. Bogotá. Unpublished raw data. [ Links ]

Ramachandran, K. M., and Tsokos, C. P. (2009). Mathematical statistics with applications. San Diego, CA: Else-vier/Academic Press. [ Links ]

Reddi, L. N., and Wu, T. H. (1991). Probabilistic Analysis of Ground Water Levels in Hillside Slopes. Journal of Geotechnical Engineering, 117(6), 872-890. DOI: 10.1061/(ASCE)0733-9410(1991)117:6(872). [ Links ]

Ruge, J. C. (2014). Análise do comportamento de cortina de estacas executada em solo poroso metaestável mediante o uso de um modelo constitutivo hipoplástico considerando a resposta não saturada. (Ph.D. thesis, Universidade de Brasília). Retrieved from: http://repositorio.unb.br/handle/10482/17558?mode=full. [ Links ]

Ruge, J. C., Da Cunha, R., Colmenares, J. E., and Mendoza, C. C. (2017). Class A prediction of a retaining structure made by a pile curtain wall executed on a tropical soil. DYNA, 84(202), 278-288. DOI: 10.15446/dyna.v84n202.63965. [ Links ]

Secretaría General de la Alcaldía Mayor de Bogotá D. C. (September 16, 2010). Por el cual se adopta la Micro-zonificación Sísmica de Bogotá D. C. [Decreto 523 de 2010]. Retrieved from http://www.alcaldiabogota.gov.co/sisjur/normas/Norma1.jsp?i=40984. [ Links ]

Soares, M., and Viana da Fonseca, A. (2016). Factors Affecting Steady State Locus in Triaxial Tests. Geotechnical Testing Journal, 39(6), 20150228. DOI: 10.1520/GTJ20150228. [ Links ]

Tabachnick, B. G., and Fidell, L. S. (2013). Using multivariate statistics (6th ed.). Boston: Pearson Education. [ Links ]

Timm, N. H. (2004). Applied Multivariate Analysis. New York, NY: Springer New York. DOI: 10.1007/b98963. [ Links ]

Tristá, J. G., Sotolongo, G. Q., Cristía, W. C., and Fernández, C. G. (2016). Simulation of the soil-water characteristic curve using laboratory tests and empirical models for Cubans unsaturated soils. Revista de la Construcción, 15(1), 22-31.Retrieved from: http://revistadelaconstruccion.uc.cl/index.php/rdlc/article/view/494/127. [ Links ]

Wackerly, D. D., Mendenhall, W., and Scheaffer, R. L. (2008). Mathematical statistics with applications. Belmont, CA: Thomson Brooks/Cole. [ Links ]

Warrick, A. W. (2003). Soil Water Dynamics (1st ed.). New York: Oxford University Press. [ Links ]

Yoon, H., Hyun, Y., Ha, K., Lee, K.-K., and Kim, G.-B. (2016). A method to improve the stability and accuracy of ANN- and SVM-based time series models for long-term groundwater level predictions. Computers & Geosciences, 90, 144-155. DOI: 10.1016/j.cageo.2016.03.002. [ Links ]

Zhao, Y., Li, Y., Zhang, L., and Wang, Q. (2016). Groundwater level prediction of landslide based on classification and regression tree. Geodesy and Geodynamics, 7(5), 348-355. DOI: 10.1016/j.geog.2016.07.005. [ Links ]

How to cite: Molina-Gómez, F., Bulla-Cruz, L. A., Moreno-Anselmi, L. A., Ruge, J. C., and Arévalo-Daza, C. (2019). Assessment of groundwater level variations using multivariate statistical methods. Ingeniería e Investigación, 39(1), 36-42. DOI: 10.15446/ing.investig.v39n1.71670

Attribution 4.0 International (CC BY 4.0) Share - Adapt

Received: April 12, 2018; Accepted: February 28, 2019

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License