1. Introduction
The atmospheric mixing layer (AML) is the turbulent layer of the atmosphere adjacent to the earth’s surface and determines the vertical mixing of air via thermal and/or mechanical turbulence [
1,
2]. The AML responds to variations in evapotranspiration and sensible heat fluxes within the timescales of an hour or less. In the context of the planetary boundary layer (PBL), the mixing layer height (MLH) represents a salient parameter because it has a considerable impact on the processes involved in the transportation and dispersion of pollutants from various sources. [
3,
4]. The development of the AML in a day is governed by a variety of parameters, such as cloud cover, water vapor content, the concentration of air pollutants, strength of synoptic wind patterns, soil moisture, nighttime cloud cover, and stratification in the free troposphere [
5]. During the daytime, the AML top represents the entrainment zone; hence, an understanding of the diurnal and seasonal variations in the atmospheric mixing layer height (MLH) is essential for discerning the mechanisms controlling the air quality, chemical processes, and numerical modeling of the lower atmosphere [
6,
7].
Automatic lidar ceilometers (ALC) with their compact design, and a high range of resolution (∼10 m), making them advantageous to many of the alternative systems for the MLH estimation. AML height estimation from ceilometers can reveal intricate features about the vertical structure of the atmosphere [
8,
9,
10].
2. Data and Methods
For the scope of this study, the first datasets of the Campbell Ceilometer C135, which is continuously operating in the National Kapodistrian University of Athens (Greece), the Department of Aerospace Science and Technology in the facilities of Euboea Island, were used. The ceilometer in this study (
Figure 1) used the method of Haij et al., 2007, to estimate the MLH [
11,
12]. At this point, it is noteworthy to point out that the area around the location of the ceilometer is rural with a significant agricultural sector and a small town with a population of 8000. Also, the location has a straight-line distance of 75 km from the international airport of Athens, which is the largest airport in Greece. This means that the information derived from the ceilometer can be very useful for understanding air transportation safety in the greater area.
All the measurements from the instrument (
Figure 1a) for the period from 1 June 2023 to 31 August 2023 were collected. During this period, a large number of valid measurements were examined. The initial dataset was filtered to keep only the clear sky measurements from days that were considered to have clear sky conditions. More specifically, as clear sky measurements were considered, in each instantaneous measurement in which there was a “cloud base”, as it was recorded by the instrument, a null value was produced. At this point, it must be stated that the ceilometer provides three different estimations of the MLH for each measurement, labeled according to their accuracy [
12]. Considering this information, a second filter was applied based on the difference between the MLH values, retaining only those measurements with a difference of less than 500 m.
In order to identify the mean daily variation in MLH, the final datasets were split by month (June, July, and August), and the mean hourly MLH was calculated (
Figure 2). Using the provided estimations of the MLH from the ceilometer, the weighted hourly mean MLH was calculated.
3. Results and Discussion
The results of the mean daily variation in MLH during the summer months (June, July, and August) can be seen in
Figure 2. This Figure depicts the fact that the lowest daily MLH values are recorded from 05:00 to 09:00, which was expected since the land stopped receiving heat from incoming solar radiation many hours before 05:00 (the previous night) and had not yet started heating from solar radiation of the current day (daily heat). The MLH started to rapidly increase after 12:00, with the whole summer period reaching a plateau of 1250–1500 m a little after dusk (around 20:00–22:00) when the impact of solar radiation stopped again due to the Earth’s daily heating cycle. The use of a third-order polynomial function to fit the mean hourly MLH values provided a high correlation coefficient (adjusted R
2) above 0.75. This type of function can efficiently represent the general shape of MLH variation. The maximum range of values of the MLH in all the examined months indicated a well-mixed atmosphere with significant heating of the land’s surface. The role of wind shear is also important, highlighting the need for such a dataset to enrich this finding. An interesting result of
Figure 2 is that the mean MLH remains at high values (above 100 m) in a significant part of the night (00:00–04:00), especially during June and July. This result indicates low surface cooling rates after the transition from day to night. This can be explained not only by the vertical stratification of the atmospheric boundary layer (ABL) itself but also by the differentiation of the land cover due to its cultivation in the surrounding areas [
13]. Nevertheless, this result needs more data (
Figure 2d) and a parallel study of land use/the land cover classification of the study area during the summer period to be evaluated.
Regarding the number of data used (
Figure 2d) to calculate the mean hourly MLH (
Figure 2a–c) during the summer months, it can be considered quite satisfactory from a statistical point of view to provide preliminary results, but a larger number of data must be used in a future study to record more efficiently the trends and the variations depicted in
Figure 2a–c. Also, alternative filters (and thresholds) to characterize clear-sky (and/or cloudy) conditions must be examined in a future study to conclude robust criteria that can represent different atmospheric conditions. Finally, the larger values of the standard deviation of the mean MLH are concentrated mainly (but not exclusively) during the day, which is consistent with the usual turbulence in the ABL afternoon, which, in turn, causes a large variation in the MLH. Finally, the role of land–sea breezes must also be examined in the affection of the MLH daily variation because five kilometers west of the ceilometer’s location is extended to the Gulf Euboea (
Figure 1b).
4. Conclusions
The central scope of this study is to examine the daily variations in the mixing layer height during the summer period in a rural environment (Euboea, Greece) during clear-sky conditions with the use of a ceilometer (Campbell ceilometer CS135). Among the main results of this study are the maximum height, which is located during the late afternoon, as well as the high values of the MLH during the night (00:00–04:00). This may be caused by the low cooling rate of the land’s surface and the lowest part of the troposphere after dusk. The daily profile for August differs from June and July, which may be caused by the land use/land cover profile due to the cultivation activities in the agricultural sector.
Moreover, the lowest MLH values are located in the transition period between day and night (05:00 to 09:00), with the August profile having significantly lower MLH values than the other two summer months.
The preliminary results of this study will be enriched in the near future with more data coming from the ceilometer, which is intended to be correlated with wind measurements and land cover data to examine in-depth the MLH variation in the study area.
Author Contributions
Conceptualization, S.K.; methodology, S.K.; software, S.K. and N.P.; validation, S.K. and N.P.; formal analysis, S.K. and N.P.; investigation, S.K..; resources, N.P.; data curation, S.K. and N.P.; writing—original draft preparation, S.K.; writing—review and editing, S.K. and N.P.; visualization, N.P.; supervision, S.K.; project administration, S.K. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to research restrictions and relative terms and conditions.
Acknowledgments
The authors acknowledge the National Kapodistrian University of Athens (Greece), Department of Aerospace Science and Technology, which supported the supply and the operation of two ceilometers (the second one will be operational by the end of the year 2023) to implement such valuable measurements.
Conflicts of Interest
The authors declare no conflicts of interest.
References
- Stull, R.B. An Introduction to Boundary Layer Meteorology; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1988. [Google Scholar]
- Seibert, P.; Beyrich, F.; Gryning, S.E.; Joffre, S.; Rasmussen, A.; Tercier, P. Review and intercomparison of operational methods for the determination of the mixing height. Atmos. Environ. 2000, 34, 1001–1027. [Google Scholar] [CrossRef]
- Yassin, M.F.; Al-Shatti, L.A.; Al Rashidi, M.S. Assessment of the atmospheric mixing layer height and its effects on pollutant dispersion. Environ Monit Assess. 2018, 190, 372. [Google Scholar] [CrossRef] [PubMed]
- Danchovski, V. Summertime Urban Mixing Layer Height over Sofia, Bulgaria. Atmosphere 2019, 10, 36. [Google Scholar] [CrossRef]
- Schween, J.H.; Hirsikko, A.; Löhnert, U.; Crewell, S. Mixing-layer height retrieval with ceilometer and Doppler lidar: From case studies to long-term assessment. Atmos. Meas. Tech. 2014, 7, 3685–3704. [Google Scholar] [CrossRef]
- Monks, P.; Granier, C.; Fuzzi, S.; Stohl, A.; Williams, M.; Akimoto, H.; Amann, M.; Baklanov, A.; Baltensperger, U.; Bey, I.; et al. Atmospheric composition change—Global and regional air quality. Atmos. Environ. 2009, 43, 5268–5350. [Google Scholar] [CrossRef]
- Bonin, T.A.; Choukulkar, A.; Brewer, W.A.; Sandberg, S.P.; Weickmann, A.M.; Pichugina, Y.; Banta, R.M.; Oncley, S.P.; Wolfe, D.E. Evaluation of Turbulence Measurement Techniques from a Single Doppler Lidar. Atmos. Meas. Tech. Discuss. 2017, 10, 3021–3039. [Google Scholar] [CrossRef]
- Yang, D.W.; Li, C.; Lau, A.K.-H.; Li, Y. Long-term measurement of daytime atmospheric mixing layer height over Hong Kong. J. Geophys. Res. Atmos. 2013, 118, 2422–2433. [Google Scholar] [CrossRef]
- Emeis, S.; Münkel, C.; Vogt, S.; Müller, W.J.; Schäfer, K. Atmospheric boundary-layer structure from simultaneous SODAR, RASS, and ceilometer measurements. Atmos. Environ. 2004, 38, 273–286. [Google Scholar] [CrossRef]
- Eresmaa, N.; Karppinen, A.; Joffre, S.M.; Räsänen, J.; Talvitie, H. Mixing height determination by ceilometer. Atmos. Chem. Phys. 2006, 6, 1485–1493. [Google Scholar] [CrossRef]
- de Haij, M.J.; Klein Baltink, H.; Wauben, W.M.F. Continuous Mixing Layer Height Determination Using the LD-40 Ceilometer: A Feasibility Study; Scientific Report WR 2007-01; KNMI: De Bilt, The Netherlands, 2007. [Google Scholar]
- Campbell Scientific. CS135 LIDAR Ceilometer Manual. 2018. Available online: https://s.campbellsci.com/documents/eu/manuals/cs135.pdf (accessed on 10 September 2023).
- Singh, N.; Solanki, R.; Ojha, N.; Janssen, R.H.H.; Pozzer, A.; Dhaka, S.K. Boundary layer evolution over the central Himalayas from radio wind profiler and model simulations. Atmos. Chem. Phys. 2016, 16, 10559–10572. [Google Scholar] [CrossRef]
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).