Radiative aerosols are known to influence the surface energy budget and hence the evolution of the planetary boundary layer. In this study, we develop a method to estimate the aerosol-induced reduction in the planetar...Radiative aerosols are known to influence the surface energy budget and hence the evolution of the planetary boundary layer. In this study, we develop a method to estimate the aerosol-induced reduction in the planetary boundary layer height (PBLH) based on two years of ground-based measurements at a site, the Station for Observing Regional Processes of the Earth System (SORPES), at Nanjing University, China, and radiosonde data from the meteorological station of Nanjing. The observations show that increased aerosol loads lead to a mean decrease of 67.1 W m-2 for downward shortwave radiation (DSR) and a mean increase of 19.2 W m-2 for downward longwave radiation (DLR), as well as a mean decrease of 9.6 W m-2 for the surface sensible heat flux (SHF) in the daytime. The relative variations of DSR, DLR and SHF are shown as a function of the increment of column mass concentration of particulate matter (PM2.5). High aerosol loading can significantly increase the atmospheric stability in the planetary boundary layer during both daytime and nighttime. Based on the statistical relationship between SHF and PM2.5 column mass concentrations, the SHF under clean atmospheric conditions (same as the background days) is derived. In this case, the derived SHF, together with observed SHF, are then used to estimate changes in the PBLH related to aerosols. Our results suggest that the PBLH decreases more rapidly with increasing aerosol loading at high aerosol loading. When the daytime mean column mass concentration of PM2.5 reaches 200 mg m-2, the decrease in the PBLH at 1600 LST (local standard time) is about 450 m.展开更多
The boundary layer structure and related heavy rainfall of Typhoon Fitow(2013), which made landfall in Zhejiang Province, China, are studied using the Advanced Research version of the Weather Research and Forecasting ...The boundary layer structure and related heavy rainfall of Typhoon Fitow(2013), which made landfall in Zhejiang Province, China, are studied using the Advanced Research version of the Weather Research and Forecasting model, with a focus on the sensitivity of the simulation to the planetary boundary layer parameterization. Two groups of experiments—one with the same surface layer scheme and including the Yonsei University(YSU), Mellor–Yamada–Nakanishi–Niino Level 2.5,and Bougeault and Lacarrere schemes; and the other with different surface layer schemes and including the Mellor–Yamada–Janjic′ and Quasi-Normal Scale Elimination schemes—are investigated. For the convenience of comparative analysis, the simulation with the YSU scheme is chosen as the control run because this scheme successfully reproduces the track, intensity and rainfall as a whole. The maximum deviations in the peak tangential and peak radial winds may account for 11% and 33%of those produced in the control run, respectively. Further diagnosis indicates that the vertical diffusivity is much larger in the first group, resulting in weaker vertical shear of the tangential and radial winds in the boundary layer and a deeper inflow layer therein. The precipitation discrepancies are related to the simulated track deflection and the differences in the simulated low-level convergent flow among all tests. Furthermore, the first group more efficiently transfers moisture and energy and produces a stronger ascending motion than the second, contributing to a deeper moist layer, stronger convection and greater precipitation.展开更多
The planetary boundary layer height(PBLH) was calculated using the radiosonde sounding data, including120 L-band operational sites and 8 GPS sites in China. The diurnal and seasonal variations of PBLH were analyzed us...The planetary boundary layer height(PBLH) was calculated using the radiosonde sounding data, including120 L-band operational sites and 8 GPS sites in China. The diurnal and seasonal variations of PBLH were analyzed using radiosonde sounding(OBS-PBLH) and ERA data(ERA-PBLH). Based on comparison and error analyses, we discussed the main error sources in these data. The frequency distributions of PBLH variations under different regimes(the convective boundary layer, the neutral residual layer, and the stable boundary layer) can be well fitted by a Gamma distribution and the shape parameter k and scale parameter s values were obtained for different regions of China. The variation characteristics of PBLH were found in summer under these three regimes for different regions. The relationships between PBLH and PM_(2.5) concentration generally follow a power law under very low or no precipitation conditions in the region of Beijing, Tianjin and Hebei in summer. The results usually deviated from this power distribution only under strong precipitation or high relative humidity conditions because of the effects of hygroscopic growth of aerosols or wet deposition. The OBS-PBLH provided a reasonable spatial distribution relative to ERA-PBLH.This indicates that OBS-PBLH has the potential for identifying the variation of PM_(2.5) concentration.展开更多
Few studies have specifically focused on the validation and spatiotemporal distribution of planetary boundary layer height (PBLH) and relative humidity (RH) data in China. In this analysis, continuous PBLH and sur...Few studies have specifically focused on the validation and spatiotemporal distribution of planetary boundary layer height (PBLH) and relative humidity (RH) data in China. In this analysis, continuous PBLH and surface-level RH data simulated from GEOS-5 between 2004 and 2012, were validated against ground-based observations. Overall, the simulated RH was consistent with the statistical data from meteorological stations, with a correlation coefficient of 0.78 and a slope of 0.9. However, the simulated PBLH was underestimated compared to LIDAR data by a factor of approximately two, which was primarily because of poor simulation in late summer and early autumn. We further examined the spatiotemporal distribution characteristics of two factors in four regions--North China, South China, Northwest China, and the Tibetan Plateau. The results showed that the annual PBLH trends in all regions were fairly moderate but sensitive to solar radiation and precipitation, which explains why the PBLH values were ranked in order from largest to smallest as follows: Tibetan Plateau, Northwest China, North China, and South China. Strong seasonal variation of the PBLH exhibited high values in summer and low values in winter, which was also consistent with the turbulent vertical exchange. Not surprisingly, the highest RH in South China and the lowest RH in desert areas of Northwest China (less than 30%). Seasonally, South China exhibited little variation, whereas Northwest China exhibited its highest humidity in winter and lowest humidity in spring, the maximum values in the other regions were obtained from July to September.展开更多
Based on different parameterization schemes of planetary boundary layer (PBL), the uncertainty of intensity and structure of the Super-strong Typhoon Rammasun (1409) is investigated using the WRF model (v3.4) with six...Based on different parameterization schemes of planetary boundary layer (PBL), the uncertainty of intensity and structure of the Super-strong Typhoon Rammasun (1409) is investigated using the WRF model (v3.4) with six PBL parameterization schemes. Results indicate that PBL uncertainty leads to the uncertainty in tropical cyclone (TC)prediction, which increases with forecast time. The uncertainty in TC prediction is mainly reflected in the uncertainty in TC intensity, with significant differences in the TC intensity forecasts using various PBL schemes. The uncertainty in TC prediction is also reflected in the uncertainty in TC structures. Greater intensity is accompanied by smaller vortex width,tighter vortex structure, stronger wind in the near-surface layer and middle and lower troposphere, stronger inflow(outflow) wind at the lower (upper) levels, stronger vertical upward wind, smaller thickness of the eye wall, smaller outward extension of the eye wall, and warmer warm core at the upper levels of eye. PBL height, surface upward heat flux and water vapor flux are important factors that cause the uncertainty in TC intensity and structure. The more surface upward heat flux and water vapor flux and the lower PBL height, the faster TC development and the stronger TC intensity.展开更多
Diurnal variations in the planetary boundary layer height(PBLH)at different latitudes over different surface characteristics are described,based on 45 years(1973−2017)of radiosonde observations.The PBLH is determined ...Diurnal variations in the planetary boundary layer height(PBLH)at different latitudes over different surface characteristics are described,based on 45 years(1973−2017)of radiosonde observations.The PBLH is determined from the radiosonde data by the bulk Richardson number(BRN)method and verified by the parcel method and the potential temperature gradient method.In general,the BRN method is able to represent the height of the convective boundary layer(BL)and neutral residual layer cases but has relatively large uncertainty in the stable BL cases.The diurnal cycle of the PBLH over land is quite different from the cycle over ocean,as are their seasonal variations.For stations over land,the PBLH shows an apparent diurnal cycle,with a distinct maximum around 15:00 LT,and seasonal variation,with higher values in summer.Compared with the PBLH over land,over oceans the PBLH diurnal cycles are quite mild,the PBLHs are much lower,and the seasonal changes are less pronounced.The seasonal variations in the median PBLH diurnal cycle are positively correlated with the near-surface temperature and negatively correlated with the near-surface relative humidity.Finally,although at most latitudes the daytime PBLH exhibits,over these 45 years,a statistically significant increasing trend at most hours between 12:00 LT and 18:00 LT over both land and ocean,there is no significant trend over either land or ocean in the nighttime PBLH for almost all the studied latitudes.展开更多
The global planetary boundary layer height(PBLH)estimated from 11 years(2007–17)of Integrated Global Radiosonde Archive(IGRA)data,Constellation Observing System for Meteorology,Ionosphere and Climate(COSMIC)soundings...The global planetary boundary layer height(PBLH)estimated from 11 years(2007–17)of Integrated Global Radiosonde Archive(IGRA)data,Constellation Observing System for Meteorology,Ionosphere and Climate(COSMIC)soundings,and European Center for Medium-Range Weather Forecasts(ECMWF)interim reanalysis(ERAInterim)data,are compared in this study.In general,the spatial distribution of global PBLH derived from ERAInterim is consistent with the one from IGRA,both at 1200 UTC and 0000 UTC.High PBLH occurs at noon local time,because of strong radiation energy and convective activity.There are larger differences between the results of COSMIC and the other two datasets.PBLHs derived from COSMIC are much higher than those from radiosonde and reanalysis data.However,PBLHs derived from the three datasets all exhibit higher values in the low latitudes and lower ones in the high latitudes.The latitudinal difference between IGRA and COSMIC ranges from−1700 m to−500 m,while it ranges from−500 m to 250 m for IGRA and ERA-Interim.It is found that the differences among the three datasets are larger in winter and smaller in summer for most studied latitudes.展开更多
New Reynolds' mean momentum equations including both turbulent viscosity and dispersion are used to analyze atmospheric balance motions of the planetary boundary layer. It is pointed out that turbulent dispersion ...New Reynolds' mean momentum equations including both turbulent viscosity and dispersion are used to analyze atmospheric balance motions of the planetary boundary layer. It is pointed out that turbulent dispersion with r 0 will increase depth of Ekman layer, reduce wind velocity in Ekman layer and produce a more satisfactory Ekman spiral lines fit the observed wind hodograph. The wind profile in the surface layer including tur-bulent dispersion is still logarithmic but the von Karman constant k is replaced by k1 = 1 -2/k, the wind increasesa little more rapidly with height.展开更多
The investigations on the dynamies of the PBL have been developed in recent years. Some authors emphasized macro-dynamics and others emphasized micro-structure of the PBL. In this paper, we study and review some main ...The investigations on the dynamies of the PBL have been developed in recent years. Some authors emphasized macro-dynamics and others emphasized micro-structure of the PBL. In this paper, we study and review some main characteristics of the wind field in the PBL from the view point connecting the macro-dynamics and micro-stucture of the PBL, thus providing the physical basis for the further research of the dynamics and the parameterization of the PBL.展开更多
In the Universidad Nacional de Colombia, Sede Medellin, the Lasers and Spectroscopy Group (GLEO) has been designed and manufactured a tropospheric lidar station based on a pulsed Nd:YAG laser operating at 532-nm wa...In the Universidad Nacional de Colombia, Sede Medellin, the Lasers and Spectroscopy Group (GLEO) has been designed and manufactured a tropospheric lidar station based on a pulsed Nd:YAG laser operating at 532-nm wavelength. The main scientific objective has been to evaluate the vertical structure of the Planetary Boundary Layer (PBL) in urban sites of Medellin-Colombia, South America, (Longitude 75°34′05″ West, Latitude 6°13′55″ North), using suspended aerosols as tracers of the atmospheric motion. This paper reports the design, manufactures and performance of an elastic lidar system, which was operated under varying air pollution and meteorological conditions and the observations presented here were performed in early 2010, over the metropolitan area of Medellin, city included in the Andean region of Colombia.展开更多
在20世纪90年代,有学者认为峨眉山大火成岩省(Emeishan Large Igneous Province,ELIP)大规模火山活动与二叠-三叠系之交(Permian-Triassic Boundary,P-TB)的生物大灭绝事件在时间上有耦合关系,随后的40Ar/39Ar同位素测年结果也显示峨眉...在20世纪90年代,有学者认为峨眉山大火成岩省(Emeishan Large Igneous Province,ELIP)大规模火山活动与二叠-三叠系之交(Permian-Triassic Boundary,P-TB)的生物大灭绝事件在时间上有耦合关系,随后的40Ar/39Ar同位素测年结果也显示峨眉山大火成岩省是晚二叠世形成的。但是,近些年大量的SHRIMP U-Pb测年结果表明,ELIP大规模火山喷发约在~260Ma;因此有研究认为,ELIP火山活动与中二叠世瓜德卢普期末(end-Guadalupian)的生物灭绝事件在时间上联系更加紧密。至于P-T界线生物大灭绝,现在多数学者认为是,由于西伯利亚大火成岩省火山强烈活动释放大量气体和火山灰所造成环境变化引起的。最近,我们在ELIP东部的贵州盘县峨眉山玄武岩系剖面中发现顶部发育厚度达近百米的凝灰岩层,其LA-ICP-MSU-Pb法测年结果为251.0±1.0Ma,与浙江煤山剖面中二叠系-三叠系边界处黏土层或火山灰层的锆石U-Pb年龄接近。因此,峨眉山玄武岩喷发结束的时间应该在P-T边界,与西伯利亚大火成岩省的主体喷发时间一致。新的测年结果暗示了ELIP火山活动与地球历史上最大的一次生物灭绝事件(P-T边界)可能存在着成因联系。展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 91544231)the State Key Research and Development Program of China (Grant No. 2016YFC0200500)+1 种基金Jiangsu Provincial Collaborative Innovation Center of Climate ChangeJun ZOU was also supported by the Program for Outstanding Ph D Candidates of Nanjing University
文摘Radiative aerosols are known to influence the surface energy budget and hence the evolution of the planetary boundary layer. In this study, we develop a method to estimate the aerosol-induced reduction in the planetary boundary layer height (PBLH) based on two years of ground-based measurements at a site, the Station for Observing Regional Processes of the Earth System (SORPES), at Nanjing University, China, and radiosonde data from the meteorological station of Nanjing. The observations show that increased aerosol loads lead to a mean decrease of 67.1 W m-2 for downward shortwave radiation (DSR) and a mean increase of 19.2 W m-2 for downward longwave radiation (DLR), as well as a mean decrease of 9.6 W m-2 for the surface sensible heat flux (SHF) in the daytime. The relative variations of DSR, DLR and SHF are shown as a function of the increment of column mass concentration of particulate matter (PM2.5). High aerosol loading can significantly increase the atmospheric stability in the planetary boundary layer during both daytime and nighttime. Based on the statistical relationship between SHF and PM2.5 column mass concentrations, the SHF under clean atmospheric conditions (same as the background days) is derived. In this case, the derived SHF, together with observed SHF, are then used to estimate changes in the PBLH related to aerosols. Our results suggest that the PBLH decreases more rapidly with increasing aerosol loading at high aerosol loading. When the daytime mean column mass concentration of PM2.5 reaches 200 mg m-2, the decrease in the PBLH at 1600 LST (local standard time) is about 450 m.
基金supported by the National Natural Science Foundation of China (Grant No. 41375056)the National Basic Research and Development Project (973 program) of China under contract no. 2015CB452805+2 种基金the National Key Technology R&D Program (Grant No. 2012BAC03)the Social Welfare Technology Development Projects of the Science and Technology Department of Zhejiang Province (Grant No. 2014C33056)the Key Project of Science and Technology Plan of Zhejiang Meteorological Provincial Bureau (2017ZD04)
文摘The boundary layer structure and related heavy rainfall of Typhoon Fitow(2013), which made landfall in Zhejiang Province, China, are studied using the Advanced Research version of the Weather Research and Forecasting model, with a focus on the sensitivity of the simulation to the planetary boundary layer parameterization. Two groups of experiments—one with the same surface layer scheme and including the Yonsei University(YSU), Mellor–Yamada–Nakanishi–Niino Level 2.5,and Bougeault and Lacarrere schemes; and the other with different surface layer schemes and including the Mellor–Yamada–Janjic′ and Quasi-Normal Scale Elimination schemes—are investigated. For the convenience of comparative analysis, the simulation with the YSU scheme is chosen as the control run because this scheme successfully reproduces the track, intensity and rainfall as a whole. The maximum deviations in the peak tangential and peak radial winds may account for 11% and 33%of those produced in the control run, respectively. Further diagnosis indicates that the vertical diffusivity is much larger in the first group, resulting in weaker vertical shear of the tangential and radial winds in the boundary layer and a deeper inflow layer therein. The precipitation discrepancies are related to the simulated track deflection and the differences in the simulated low-level convergent flow among all tests. Furthermore, the first group more efficiently transfers moisture and energy and produces a stronger ascending motion than the second, contributing to a deeper moist layer, stronger convection and greater precipitation.
基金National Key R&D Program Pilot Projects of China(2016YFC203300)Major Program of National Natural Science Foundation of China(91644223)+2 种基金Special Funding Project for Public Industry Research and Development of Ministry of Environmental Protection(201509001)National Natural Science Foundation of China(9133700041575008)
文摘The planetary boundary layer height(PBLH) was calculated using the radiosonde sounding data, including120 L-band operational sites and 8 GPS sites in China. The diurnal and seasonal variations of PBLH were analyzed using radiosonde sounding(OBS-PBLH) and ERA data(ERA-PBLH). Based on comparison and error analyses, we discussed the main error sources in these data. The frequency distributions of PBLH variations under different regimes(the convective boundary layer, the neutral residual layer, and the stable boundary layer) can be well fitted by a Gamma distribution and the shape parameter k and scale parameter s values were obtained for different regions of China. The variation characteristics of PBLH were found in summer under these three regimes for different regions. The relationships between PBLH and PM_(2.5) concentration generally follow a power law under very low or no precipitation conditions in the region of Beijing, Tianjin and Hebei in summer. The results usually deviated from this power distribution only under strong precipitation or high relative humidity conditions because of the effects of hygroscopic growth of aerosols or wet deposition. The OBS-PBLH provided a reasonable spatial distribution relative to ERA-PBLH.This indicates that OBS-PBLH has the potential for identifying the variation of PM_(2.5) concentration.
基金supported by the National Key R&D Program of China (2016YFC0201507)the National Natural Science Foundation of China (Grant Nos. 41471367, 91543128 and 41571417)
文摘Few studies have specifically focused on the validation and spatiotemporal distribution of planetary boundary layer height (PBLH) and relative humidity (RH) data in China. In this analysis, continuous PBLH and surface-level RH data simulated from GEOS-5 between 2004 and 2012, were validated against ground-based observations. Overall, the simulated RH was consistent with the statistical data from meteorological stations, with a correlation coefficient of 0.78 and a slope of 0.9. However, the simulated PBLH was underestimated compared to LIDAR data by a factor of approximately two, which was primarily because of poor simulation in late summer and early autumn. We further examined the spatiotemporal distribution characteristics of two factors in four regions--North China, South China, Northwest China, and the Tibetan Plateau. The results showed that the annual PBLH trends in all regions were fairly moderate but sensitive to solar radiation and precipitation, which explains why the PBLH values were ranked in order from largest to smallest as follows: Tibetan Plateau, Northwest China, North China, and South China. Strong seasonal variation of the PBLH exhibited high values in summer and low values in winter, which was also consistent with the turbulent vertical exchange. Not surprisingly, the highest RH in South China and the lowest RH in desert areas of Northwest China (less than 30%). Seasonally, South China exhibited little variation, whereas Northwest China exhibited its highest humidity in winter and lowest humidity in spring, the maximum values in the other regions were obtained from July to September.
基金National Natural Science Foundation of China(41575108,41475082)Scientific Research Program of Shanghai Science and Technology Commission(19dz1200101)。
文摘Based on different parameterization schemes of planetary boundary layer (PBL), the uncertainty of intensity and structure of the Super-strong Typhoon Rammasun (1409) is investigated using the WRF model (v3.4) with six PBL parameterization schemes. Results indicate that PBL uncertainty leads to the uncertainty in tropical cyclone (TC)prediction, which increases with forecast time. The uncertainty in TC prediction is mainly reflected in the uncertainty in TC intensity, with significant differences in the TC intensity forecasts using various PBL schemes. The uncertainty in TC prediction is also reflected in the uncertainty in TC structures. Greater intensity is accompanied by smaller vortex width,tighter vortex structure, stronger wind in the near-surface layer and middle and lower troposphere, stronger inflow(outflow) wind at the lower (upper) levels, stronger vertical upward wind, smaller thickness of the eye wall, smaller outward extension of the eye wall, and warmer warm core at the upper levels of eye. PBL height, surface upward heat flux and water vapor flux are important factors that cause the uncertainty in TC intensity and structure. The more surface upward heat flux and water vapor flux and the lower PBL height, the faster TC development and the stronger TC intensity.
基金This work was supported by the Meteorological Research Open Foundation of Huaihe Basin(HRM201604).
文摘Diurnal variations in the planetary boundary layer height(PBLH)at different latitudes over different surface characteristics are described,based on 45 years(1973−2017)of radiosonde observations.The PBLH is determined from the radiosonde data by the bulk Richardson number(BRN)method and verified by the parcel method and the potential temperature gradient method.In general,the BRN method is able to represent the height of the convective boundary layer(BL)and neutral residual layer cases but has relatively large uncertainty in the stable BL cases.The diurnal cycle of the PBLH over land is quite different from the cycle over ocean,as are their seasonal variations.For stations over land,the PBLH shows an apparent diurnal cycle,with a distinct maximum around 15:00 LT,and seasonal variation,with higher values in summer.Compared with the PBLH over land,over oceans the PBLH diurnal cycles are quite mild,the PBLHs are much lower,and the seasonal changes are less pronounced.The seasonal variations in the median PBLH diurnal cycle are positively correlated with the near-surface temperature and negatively correlated with the near-surface relative humidity.Finally,although at most latitudes the daytime PBLH exhibits,over these 45 years,a statistically significant increasing trend at most hours between 12:00 LT and 18:00 LT over both land and ocean,there is no significant trend over either land or ocean in the nighttime PBLH for almost all the studied latitudes.
基金supported by the Meteorological Research Open Foundation of Huaihe Basin grant number HRM201604。
文摘The global planetary boundary layer height(PBLH)estimated from 11 years(2007–17)of Integrated Global Radiosonde Archive(IGRA)data,Constellation Observing System for Meteorology,Ionosphere and Climate(COSMIC)soundings,and European Center for Medium-Range Weather Forecasts(ECMWF)interim reanalysis(ERAInterim)data,are compared in this study.In general,the spatial distribution of global PBLH derived from ERAInterim is consistent with the one from IGRA,both at 1200 UTC and 0000 UTC.High PBLH occurs at noon local time,because of strong radiation energy and convective activity.There are larger differences between the results of COSMIC and the other two datasets.PBLHs derived from COSMIC are much higher than those from radiosonde and reanalysis data.However,PBLHs derived from the three datasets all exhibit higher values in the low latitudes and lower ones in the high latitudes.The latitudinal difference between IGRA and COSMIC ranges from−1700 m to−500 m,while it ranges from−500 m to 250 m for IGRA and ERA-Interim.It is found that the differences among the three datasets are larger in winter and smaller in summer for most studied latitudes.
文摘New Reynolds' mean momentum equations including both turbulent viscosity and dispersion are used to analyze atmospheric balance motions of the planetary boundary layer. It is pointed out that turbulent dispersion with r 0 will increase depth of Ekman layer, reduce wind velocity in Ekman layer and produce a more satisfactory Ekman spiral lines fit the observed wind hodograph. The wind profile in the surface layer including tur-bulent dispersion is still logarithmic but the von Karman constant k is replaced by k1 = 1 -2/k, the wind increasesa little more rapidly with height.
文摘The investigations on the dynamies of the PBL have been developed in recent years. Some authors emphasized macro-dynamics and others emphasized micro-structure of the PBL. In this paper, we study and review some main characteristics of the wind field in the PBL from the view point connecting the macro-dynamics and micro-stucture of the PBL, thus providing the physical basis for the further research of the dynamics and the parameterization of the PBL.
文摘In the Universidad Nacional de Colombia, Sede Medellin, the Lasers and Spectroscopy Group (GLEO) has been designed and manufactured a tropospheric lidar station based on a pulsed Nd:YAG laser operating at 532-nm wavelength. The main scientific objective has been to evaluate the vertical structure of the Planetary Boundary Layer (PBL) in urban sites of Medellin-Colombia, South America, (Longitude 75°34′05″ West, Latitude 6°13′55″ North), using suspended aerosols as tracers of the atmospheric motion. This paper reports the design, manufactures and performance of an elastic lidar system, which was operated under varying air pollution and meteorological conditions and the observations presented here were performed in early 2010, over the metropolitan area of Medellin, city included in the Andean region of Colombia.