<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "https://jats.nlm.nih.gov/nlm-dtd/publishing/3.0/journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0" article-type="research-article"><?xmltex \bartext{Research article}?>
  <front>
    <journal-meta><journal-id journal-id-type="publisher">ACP</journal-id><journal-title-group>
    <journal-title>Atmospheric Chemistry and Physics</journal-title>
    <abbrev-journal-title abbrev-type="publisher">ACP</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Atmos. Chem. Phys.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">1680-7324</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/acp-23-15365-2023</article-id><title-group><article-title>The role of a low-level jet for stirring the stable atmospheric surface layer in the Arctic</article-title><alt-title>The role of a low-level jet for stirring the stable atmospheric surface layer in the Arctic</alt-title>
      </title-group><?xmltex \runningtitle{The role of a low-level jet for stirring the stable atmospheric surface layer in the Arctic}?><?xmltex \runningauthor{U. Egerer et al.}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="yes" rid="aff1 aff4">
          <name><surname>Egerer</surname><given-names>Ulrike</given-names></name>
          <email>ulrike.egerer@nrel.gov</email>
        <ext-link>https://orcid.org/0000-0001-6107-612X</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Siebert</surname><given-names>Holger</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>Hellmuth</surname><given-names>Olaf</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff3">
          <name><surname>Sørensen</surname><given-names>Lise Lotte</given-names></name>
          
        <ext-link>https://orcid.org/0000-0002-9823-589X</ext-link></contrib>
        <aff id="aff1"><label>1</label><institution>Atmospheric Microphysics, Leibniz Institute for Tropospheric Research, <?xmltex \hack{\break}?>Permoserstr. 15,  04318 Leipzig, Germany</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>Modeling of Atmospheric Processes, Leibniz Institute for Tropospheric Research, <?xmltex \hack{\break}?>Permoserstr. 15, 04318 Leipzig, Germany</institution>
        </aff>
        <aff id="aff3"><label>3</label><institution>Department of Environmental Science, iClimate, Arctic Research Centre ARC, Aarhus University, Frederiksborgvej 399, 4000 Roskilde, Denmark</institution>
        </aff>
        <aff id="aff4"><label>a</label><institution>now at: National Renewable Energy Laboratory (NREL), 15013 Denver W Parkway,<?xmltex \hack{\break}?> Golden, CO 80401, USA</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Ulrike Egerer (ulrike.egerer@nrel.gov)</corresp></author-notes><pub-date><day>15</day><month>December</month><year>2023</year></pub-date>
      
      <volume>23</volume>
      <issue>24</issue>
      <fpage>15365</fpage><lpage>15373</lpage>
      <history>
        <date date-type="received"><day>24</day><month>March</month><year>2023</year></date>
           <date date-type="rev-request"><day>13</day><month>July</month><year>2023</year></date>
           <date date-type="rev-recd"><day>29</day><month>October</month><year>2023</year></date>
           <date date-type="accepted"><day>1</day><month>November</month><year>2023</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2023 </copyright-statement>
        <copyright-year>2023</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://acp.copernicus.org/articles/.html">This article is available from https://acp.copernicus.org/articles/.html</self-uri><self-uri xlink:href="https://acp.copernicus.org/articles/.pdf">The full text article is available as a PDF file from https://acp.copernicus.org/articles/.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d1e135">In this study, we analyze the transition of a stable atmospheric boundary layer (ABL) with a low-level jet (LLJ) to a traditional stable ABL with a classic Ekman helix in the late-winter central Arctic. Vertical profiles in the ABL were measured with a hot-wire anemometer on a tethered balloon during a 15 h period in March 2018 in northeast Greenland. The tethered balloon allows high-resolution turbulence observations from the ground to the top of the ABL. The core of the LLJ was observed at about 150 m altitude, and its height and strength were associated with the temperature inversion. Increased turbulence was observed in the vicinity of the LLJ, but most of the turbulence does not reach down to the surface, thus decoupling the LLJ from the surface. Only when the LLJ collapses and the ABL again exhibits a more classical Ekman spiral is a coupling to the surface re-established. The LLJ might enhance both advective and turbulent vertical transport of passive tracers such as aerosol particles or moisture in the often stably stratified Arctic ABL.</p>
  </abstract>
    
<funding-group>
<award-group id="gs1">
<funding-source>Deutsche Forschungsgemeinschaft</funding-source>
<award-id>268020496 - TRR 172</award-id>
</award-group>
</funding-group>
</article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d1e147">Low-level jets (LLJs) are vertically more or less bounded wind fields with local maximum wind velocities exceeding the geostrophic wind. An LLJ, especially in conjunction with stably stratified boundary layers, is usually found in the upper region of the atmospheric boundary layer (ABL) or even partly just inside or above the inversion. LLJs are often related to inertia oscillations and a decoupling of the flow dynamics from the surface-layer friction <xref ref-type="bibr" rid="bib1.bibx5 bib1.bibx17" id="paren.1"/>. There is a lively debate about the origin of LLJs <xref ref-type="bibr" rid="bib1.bibx20 bib1.bibx11 bib1.bibx19 bib1.bibx9 bib1.bibx6" id="paren.2"/>, and details of the mechanism behind this are still not completely understood. Polar regions are preferable locations for the occurrence of LLJs <xref ref-type="bibr" rid="bib1.bibx19 bib1.bibx12" id="paren.3"/> due to frequently observed (extremely) stably stratified and shallow ABLs, in particular during cloudless (wintertime) conditions with strongly negative thermal-infrared net irradiance. In northeast coastal Greenland, the area of interest for this study, LLJs occur 60 %–70 % of the time, are observed at very low altitudes and are attributed to katabatic flows <xref ref-type="bibr" rid="bib1.bibx19" id="paren.4"/>.</p>
      <?pagebreak page15366?><p id="d1e162">The stably stratified boundary layer, as defined by an increase in potential temperature with height, usually exhibits comparably low turbulence. When neglecting cloud-related effects such as radiative cooling at the cloud top, the only significant source of turbulent kinetic energy (TKE) is the surface roughness. However, due to the strong local wind shear below and above the wind maximum, an LLJ can introduce a significant amount of TKE <xref ref-type="bibr" rid="bib1.bibx17" id="paren.5"/> and modulate the vertical distribution of TKE <xref ref-type="bibr" rid="bib1.bibx11 bib1.bibx4" id="paren.6"/>. These wind shear zones are usually considered a source of TKE production as long as the damping effect of the inversion does not dominate the TKE production. Boundary-layer flows with a source of turbulence at the top of the inversion generated by vertical wind shear have also been referred to as “upside down” <xref ref-type="bibr" rid="bib1.bibx14" id="paren.7"/>. The relative importance of turbulence-generating wind shear and damping inversion is quantified by the dimensionless Richardson number (Ri) defined later on. In addition to TKE production, LLJs may play an important role in the advection of momentum, turbulence, and aerosol particles and precursor gases <xref ref-type="bibr" rid="bib1.bibx18 bib1.bibx2" id="paren.8"/>.</p>
      <p id="d1e177">Here, we study the ABL evolution in terms of an LLJ, as observed by a set of subsequent vertical wind and temperature profiles measured with a tethered balloon in the framework of the Polar Airborne Measurements and Arctic Regional Climate Model SImulation Project (PAMARCMiP) at Station Nord in northeast Greenland. The central question of this work is how the LLJ might affect the vertical distribution of turbulence and whether increased turbulence is observed at the surface when the LLJ occurs or collapses. This could indicate that properties advected with the LLJ, such as increased aerosol concentration or precursor gases, could also be mixed down to the surface after being advected over a certain distance inside the LLJ.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Observations</title>
<sec id="Ch1.S2.SS1">
  <label>2.1</label><title>PAMARCMiP campaign</title>
      <p id="d1e195">The PAMARCMiP field campaign was conducted at the Villum Research Station (VRS) on the military outpost “Station Nord” in the northeast corner of Greenland on the small peninsula of Princess Ingeborg (81<inline-formula><mml:math id="M1" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>36<inline-formula><mml:math id="M2" display="inline"><mml:msup><mml:mi/><mml:mo>′</mml:mo></mml:msup></mml:math></inline-formula> N, 16<inline-formula><mml:math id="M3" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>40<inline-formula><mml:math id="M4" display="inline"><mml:msup><mml:mi/><mml:mo>′</mml:mo></mml:msup></mml:math></inline-formula> W). Balloon observations were made from 10 March to 7 April 2018 about 2 <inline-formula><mml:math id="M5" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">km</mml:mi></mml:mrow></mml:math></inline-formula> south of VRS at a small observation hut (“Flygers Hut”) to minimize the station's influence. The site is located on the coastline between the Greenland Ice Sheet to the south and the Arctic Ocean, which is covered with sea ice for most of the year at this location. A glacier is located about 20 <inline-formula><mml:math id="M6" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">km</mml:mi></mml:mrow></mml:math></inline-formula> to the south.</p>
      <p id="d1e251">The LLJ and ABL structure observations are mainly based on the tethered-balloon system BELUGA <xref ref-type="bibr" rid="bib1.bibx7" id="paren.9"/>. Continuously running meteorological measurements at 9 <inline-formula><mml:math id="M7" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> altitude at the VRS and three-dimensional wind data from an ultrasonic anemometer (USA-1, manufactured by METEK GmbH, Germany) mounted on a tower at 65 <inline-formula><mml:math id="M8" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> height support the analysis of the profile measurements.</p>
</sec>
<sec id="Ch1.S2.SS2">
  <label>2.2</label><title>The BELUGA setup</title>
      <p id="d1e281">BELUGA is a modular system consisting of tethered balloons of various sizes and a variety of sensor packages. During PAMARCMiP, a 9 <inline-formula><mml:math id="M9" display="inline"><mml:mrow class="unit"><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> balloon with a maximum payload of 3 <inline-formula><mml:math id="M10" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">kg</mml:mi></mml:mrow></mml:math></inline-formula> was used with multiple sensor packages. The instrument packages were deployed in a rotating sequence aboard BELUGA to ensure adherence to the payload limit; a flight overview table with deployed instruments is provided in <xref ref-type="bibr" rid="bib1.bibx8" id="text.10"/>. In this work, data from two sensor packages are used. (i) A turbulence probe based essentially on a single-component hot-wire anemometer measuring wind speed with a sampling frequency of 500 <inline-formula><mml:math id="M11" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">Hz</mml:mi></mml:mrow></mml:math></inline-formula>. A Pitot static tube served as a reference for the hot-wire sensors. Fast temperature measurements were made with a cold-wire sensor.  In addition to a calibrated thermometer as a reference for the cold-wire sensor, the data logger provided additional pressure measurements for barometric altitude. (ii) A standard meteorological probe (“StdMeteo”)  based on a Graw DFM-09 radiosonde with an additional Pitot static probe provided wind speed and direction using a compass. This package measured wind, temperature and relative humidity at a frequency of 1 <inline-formula><mml:math id="M12" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">Hz</mml:mi></mml:mrow></mml:math></inline-formula> and was intended for regular use on all flights. The absolute accuracy of the wind speed measurement with the Pitot static tube is determined by the accuracy of the differential pressure sensor and is additionally inversely proportional to the wind speed itself. For typical wind speeds of 5 <inline-formula><mml:math id="M13" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, the absolute accuracy is only in the range of 0.5 <inline-formula><mml:math id="M14" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> which is essentially noticeable by a zero offset and can therefore be partially corrected; however, the relative resolution is higher by a factor of 10 for the Pitot static tube and the spectrally determined resolution of the hot-wire anemometer is in the range of 0.2 <inline-formula><mml:math id="M15" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">cm</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> <xref ref-type="bibr" rid="bib1.bibx7" id="paren.11"/>.</p>
</sec>
<sec id="Ch1.S2.SS3">
  <label>2.3</label><title>Measurements and synoptical situation</title>
      <p id="d1e385">This study focuses on nine subsequent flights between 10:00 UTC on 29 March and 01:00 UTC on 30 March 2018, during the transition between polar night and day.  The period was mainly characterized by the persistent occurrence of an LLJ, followed by a transition to a more classic stable ABL which is characterized by the wind speed increasing mostly logarithmically from the ground with height.  The turbulence is generated by the wind shear and is therefore strongest near the ground and decreases continuously with height, depending on the damping effect of the temperature inversion. Synoptic conditions during this period were influenced by high air pressure over central Greenland yielding calm weather conditions at ground level. All flights occurred under cloudless conditions with only thin clouds well above 1000 <inline-formula><mml:math id="M16" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>. Therefore, possible icing was not an issue.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1"><?xmltex \currentcnt{1}?><?xmltex \def\figurename{Figure}?><label>Figure 1</label><caption><p id="d1e398">Time series for 29–30 March 2018 of <bold>(a)</bold> BELUGA flight altitude from the two instrument packages “StdMeteo” and “Hotwire” and <bold>(b)</bold> VRS near-surface measurements of temperature <inline-formula><mml:math id="M17" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> and specific humidity <inline-formula><mml:math id="M18" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> as well as <bold>(c)</bold> wind velocity <inline-formula><mml:math id="M19" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> and wind direction dd.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/15365/2023/acp-23-15365-2023-f01.png"/>

        </fig>

      <?pagebreak page15367?><p id="d1e438">Figure <xref ref-type="fig" rid="Ch1.F1"/> provides an overview of the nine BELUGA flight profiles up to altitudes of 600 to 1400 <inline-formula><mml:math id="M20" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> in combination with near-surface measurements. The near-surface conditions are quite variable in terms of temperature <inline-formula><mml:math id="M21" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> and specific humidity <inline-formula><mml:math id="M22" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula>, with <inline-formula><mml:math id="M23" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> ranging between <inline-formula><mml:math id="M24" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>35 and <inline-formula><mml:math id="M25" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>26 <inline-formula><mml:math id="M26" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C and <inline-formula><mml:math id="M27" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> ranging between 0.15 to 0.3 <inline-formula><mml:math id="M28" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">kg</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>. Qualitatively, there is an obvious high correlation between <inline-formula><mml:math id="M29" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M30" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula>. Wind speeds <inline-formula><mml:math id="M31" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> were below 2 <inline-formula><mml:math id="M32" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and decreased towards the end of the observation period. The wind direction dd is mainly from the west but turns to being south and then to north from 16:00 UTC and back to west after 19:00 UTC. Between 13:00 and 14:00 UTC, <inline-formula><mml:math id="M33" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M34" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> drop abruptly by <inline-formula><mml:math id="M35" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>8 <inline-formula><mml:math id="M36" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">K</mml:mi></mml:mrow></mml:math></inline-formula> and 0.2 <inline-formula><mml:math id="M37" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">kg</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, respectively, along with a change in wind direction (southeast to north)  and reach their previous values about 1 h later.</p>
</sec>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Data analysis methods</title>
      <p id="d1e615">In this study, we use the balloon-borne vertical profiles to relate the properties of an LLJ to the vertical structure of stability and turbulence. The literature reveals various definitions of an LLJ, mostly based on a local low-altitude wind velocity maximum greater than around 2 <inline-formula><mml:math id="M38" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> <xref ref-type="bibr" rid="bib1.bibx19 bib1.bibx3 bib1.bibx5" id="paren.12"/>. We adopt and slightly modify this definition and define the following criteria for an LLJ: (i) the wind velocity maximum occurs below 250 <inline-formula><mml:math id="M39" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> altitude, and (ii) the difference between the maximum wind velocity in the LLJ core <inline-formula><mml:math id="M40" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mtext>LLJ</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> and the wind minimum <inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mtext>min</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> above  and below (commonly the near-surface wind velocity) exceeds 2.5 <inline-formula><mml:math id="M42" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>. The LLJ strength is then defined as <inline-formula><mml:math id="M43" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>u</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mtext>LLJ</mml:mtext></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mtext>min</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> with <inline-formula><mml:math id="M44" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mtext>min</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> being the higher value either above or below the jet core. The height of the LLJ core <inline-formula><mml:math id="M45" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mtext>LLJ</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> is typically located at the maximum height of the (strong) temperature inversion <inline-formula><mml:math id="M46" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. We tailor some of the criteria used in the literature to our specific observations and the unique characteristics of this particular LLJ.</p>
      <p id="d1e743">The dimensionless gradient Richardson number <inline-formula><mml:math id="M47" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the ratio of buoyancy to shear:
          <disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M48" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>g</mml:mi><mml:mover accent="true"><mml:mi mathvariant="italic">θ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>∂</mml:mo><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>/</mml:mo><mml:mo>∂</mml:mo><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mo>∂</mml:mo><mml:mi>u</mml:mi><mml:mo>/</mml:mo><mml:mo>∂</mml:mo><mml:mi>z</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mstyle><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        with the vertical gradient of the mean flow speed <inline-formula><mml:math id="M49" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mo>∂</mml:mo><mml:mi>u</mml:mi><mml:mo>/</mml:mo><mml:mo>∂</mml:mo><mml:mi>z</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>≈</mml:mo><mml:mo>(</mml:mo><mml:mo>∂</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mi>x</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mo>∂</mml:mo><mml:mi>z</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>(</mml:mo><mml:mo>∂</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mi>y</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mo>∂</mml:mo><mml:mi>z</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>, the vertical potential temperature gradient <inline-formula><mml:math id="M50" display="inline"><mml:mrow><mml:mo>∂</mml:mo><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>/</mml:mo><mml:mo>∂</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:math></inline-formula>, and the acceleration of gravity <inline-formula><mml:math id="M51" display="inline"><mml:mi>g</mml:mi></mml:math></inline-formula>. When stratification dominates over wind shear and a critical Richardson number <inline-formula><mml:math id="M52" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mtext>c</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula>  – estimated to be between 0.25 and 1 <xref ref-type="bibr" rid="bib1.bibx15 bib1.bibx1" id="paren.13"/> – is reached, turbulence decreases. Thus, the contribution of buoyancy and shear to the turbulence profile throughout the ABL can be analyzed by the shear and buoyancy terms of the <inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. For estimating <inline-formula><mml:math id="M54" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M55" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M56" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> profiles are smoothed with a 10 <inline-formula><mml:math id="M57" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">s</mml:mi></mml:mrow></mml:math></inline-formula> rolling mean before calculating the local gradients. The final <inline-formula><mml:math id="M58" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> profile is then smoothed again with a rolling mean with a 10 <inline-formula><mml:math id="M59" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">s</mml:mi></mml:mrow></mml:math></inline-formula> window. As an approximation for the range between the surface and the core of the LLJ, the bulk Richardson number <inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> <xref ref-type="bibr" rid="bib1.bibx13" id="paren.14"/> is calculated for the height <inline-formula><mml:math id="M61" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula>:
          <disp-formula id="Ch1.E2" content-type="numbered"><label>2</label><mml:math id="M62" display="block"><mml:mrow><mml:msub><mml:mtext>Ri</mml:mtext><mml:mi mathvariant="normal">b</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>g</mml:mi><mml:mover accent="true"><mml:mi mathvariant="italic">θ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo><mml:mo>⋅</mml:mo><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>u</mml:mi><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
        Here, <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, or <inline-formula><mml:math id="M64" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>u</mml:mi><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, is the difference between <inline-formula><mml:math id="M65" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula>, or <inline-formula><mml:math id="M66" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula>, at <inline-formula><mml:math id="M67" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> and close to the surface.</p>
      <p id="d1e1128">The local turbulence is characterized by means of the local energy dissipation rate <inline-formula><mml:math id="M68" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula>, which is derived from the second-order structure function of <inline-formula><mml:math id="M69" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> applying inertial sub-range scaling:
          <disp-formula id="Ch1.E3" content-type="numbered"><label>3</label><mml:math id="M70" display="block"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi>u</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mo>*</mml:mo></mml:msup><mml:mo>)</mml:mo><mml:mo>-</mml:mo><mml:mi>u</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mi>C</mml:mi><mml:mo>⋅</mml:mo><mml:msup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup><mml:mo>⋅</mml:mo><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>*</mml:mo></mml:msup><mml:mo>⋅</mml:mo><mml:mover accent="true"><mml:mi>u</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:mrow></mml:mfenced><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        with <inline-formula><mml:math id="M71" display="inline"><mml:mrow><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:math></inline-formula> and a time lag <inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> for the longitudinal wind velocity <inline-formula><mml:math id="M73" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> <xref ref-type="bibr" rid="bib1.bibx16" id="paren.15"><named-content content-type="post">Eqs. 10 and 11;</named-content></xref>. Here, the overline in Eq. (<xref ref-type="disp-formula" rid="Ch1.E3"/>) denotes time averaging; an averaging interval of 2 <inline-formula><mml:math id="M74" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">s</mml:mi></mml:mrow></mml:math></inline-formula> was selected which yields robust estimates <xref ref-type="bibr" rid="bib1.bibx7" id="paren.16"/>. The minimum resolvable <inline-formula><mml:math id="M75" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> is estimated to be about <inline-formula><mml:math id="M76" display="inline"><mml:mrow><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> <inline-formula><mml:math id="M77" display="inline"><mml:mrow class="unit"><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>.</p>
      <?pagebreak page15368?><p id="d1e1318">To better classify both the integration and averaging times in the determination of parameters as discussed before but also to better interpret the vertical profiles in the next sections, knowledge of the length scales involved is of certain interest. Integral length scales of a turbulent flow, describing the size of the largest energy-containing eddies, can be defined according to
          <disp-formula id="Ch1.E4" content-type="numbered"><label>4</label><mml:math id="M78" display="block"><mml:mrow><mml:mi mathvariant="script">L</mml:mi><mml:mo>=</mml:mo><mml:mi mathvariant="italic">β</mml:mi><mml:mo>⋅</mml:mo><mml:msubsup><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msubsup><mml:mo>/</mml:mo><mml:mi mathvariant="italic">ε</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
        <xref ref-type="bibr" rid="bib1.bibx22" id="paren.17"><named-content content-type="pre">p. 16;</named-content></xref> with a constant <inline-formula><mml:math id="M79" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula> of  <inline-formula><mml:math id="M80" display="inline"><mml:mrow><mml:mi mathvariant="script">O</mml:mi><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>. We call the ratio of individual, locally fluctuating <inline-formula><mml:math id="M81" display="inline"><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:mrow><mml:mi mathvariant="italic">ε</mml:mi></mml:mfrac></mml:mstyle></mml:math></inline-formula> values “local length scales” <inline-formula><mml:math id="M82" display="inline"><mml:mi mathvariant="script">L</mml:mi></mml:math></inline-formula>. Here, <inline-formula><mml:math id="M83" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> is interpreted as a local parameter derived from 30 s long sub-records to be consistent with <inline-formula><mml:math id="M84" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> (with <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> being the standard deviation of <inline-formula><mml:math id="M86" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula>) averaged over 30 <inline-formula><mml:math id="M87" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">s</mml:mi></mml:mrow></mml:math></inline-formula>, since <inline-formula><mml:math id="M88" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is determined by the largest contributing scales. Earlier studies <xref ref-type="bibr" rid="bib1.bibx7 bib1.bibx16" id="paren.18"/> have shown that the estimation of dissipation rates is insensitive to the choice of the averaging time. Figure <xref ref-type="fig" rid="Ch1.F2"/> shows a scatterplot of <inline-formula><mml:math id="M89" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> versus <inline-formula><mml:math id="M90" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> based on all BELUGA data observed on 29 March 2018. For each measuring point, the corresponding <inline-formula><mml:math id="M91" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> value is represented by a color code to better classify the values. The relation of observed <inline-formula><mml:math id="M92" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M93" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> is almost linear as predicted by Eq. (<xref ref-type="disp-formula" rid="Ch1.E4"/>) but shows some variations. A majority of resulting local length scales (98 <inline-formula><mml:math id="M94" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">%</mml:mi></mml:mrow></mml:math></inline-formula> of the data points, assuming <inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>) are in the range between 0.1 to 10 <inline-formula><mml:math id="M96" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>  with somewhat larger <inline-formula><mml:math id="M97" display="inline"><mml:mi mathvariant="script">L</mml:mi></mml:math></inline-formula> in regions of comparably high turbulence (blue–purple colors). A few data points with <inline-formula><mml:math id="M98" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mo>&gt;</mml:mo><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> (yellow color) indicate small length scales below around 1 <inline-formula><mml:math id="M99" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>. The ratio of the assumed mean wind velocity of <inline-formula><mml:math id="M100" display="inline"><mml:mrow><mml:mover accent="true"><mml:mi>u</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mo>≈</mml:mo></mml:mrow></mml:math></inline-formula> 5 <inline-formula><mml:math id="M101" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to the averaging time of 30 <inline-formula><mml:math id="M102" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">s</mml:mi></mml:mrow></mml:math></inline-formula> for estimating <inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> yields an averaging length of 150 <inline-formula><mml:math id="M104" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>, which is about 10 to 100 times larger than the <inline-formula><mml:math id="M105" display="inline"><mml:mi mathvariant="script">L</mml:mi></mml:math></inline-formula> estimates. Therefore, we can safely conclude that the averaging period covers enough eddies and the results are statistically robust.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2"><?xmltex \currentcnt{2}?><?xmltex \def\figurename{Figure}?><label>Figure 2</label><caption><p id="d1e1643">Relation of turbulence parameters for all BELUGA flights on 29 March 2018 with the hot-wire: gradient Richardson number <inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>u</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> and local dissipation rate <inline-formula><mml:math id="M108" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula>. Each dot represents a 30 <inline-formula><mml:math id="M109" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">s</mml:mi></mml:mrow></mml:math></inline-formula> time interval of averaged parameters. The dashed line shows a linear regression fitted to the data points. </p></caption>
        <?xmltex \igopts{width=227.622047pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/15365/2023/acp-23-15365-2023-f02.png"/>

      </fig>

      <p id="d1e1691">For the 65 <inline-formula><mml:math id="M110" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> mast wind data, <inline-formula><mml:math id="M111" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> is estimated analogously to the vertical profiles but averaged over 5 <inline-formula><mml:math id="M112" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">min</mml:mi></mml:mrow></mml:math></inline-formula> segments and using the vertical wind velocity component. The vertical wind velocity of the 65 <inline-formula><mml:math id="M113" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> mast sonic anemometer for the flux calculation is tilt-corrected using the double-rotation algorithm described by <xref ref-type="bibr" rid="bib1.bibx21" id="text.19"/>. For the determination of the turbulent fluxes of virtual sensible heat, <inline-formula><mml:math id="M114" display="inline"><mml:mrow><mml:mi>H</mml:mi><mml:mo>=</mml:mo><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo>⋅</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mover accent="true"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">θ</mml:mi><mml:mi>v</mml:mi><mml:mo>′</mml:mo></mml:msubsup><mml:msup><mml:mi>w</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>, and momentum, <inline-formula><mml:math id="M115" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi>u</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo>⋅</mml:mo><mml:mover accent="true"><mml:mrow><mml:msup><mml:mi>u</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:msup><mml:mi>w</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>, the fluctuations are calculated by applying linear detrending and averaging over 30 <inline-formula><mml:math id="M116" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">min</mml:mi></mml:mrow></mml:math></inline-formula> periods at 5 <inline-formula><mml:math id="M117" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">min</mml:mi></mml:mrow></mml:math></inline-formula> time steps.</p>
</sec>
<sec id="Ch1.S4">
  <label>4</label><title>Observed vertical structure of the ABL and LLJ</title>
<sec id="Ch1.S4.SS1">
  <label>4.1</label><title>Evolution of the LLJ and ABL structure</title>
      <p id="d1e1828">We first analyze the evolution of the mean ABL structure and the LLJ throughout the observation period. Figure <xref ref-type="fig" rid="Ch1.F3"/> provides an overview of the period based on BELUGA profiles (time–height contour plots), VRS observations in terms of meteorological measurements at 9 <inline-formula><mml:math id="M118" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> and the sonic data at 65 <inline-formula><mml:math id="M119" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> height. The observation period can be divided into three sub-periods: (i) an LLJ period, (ii) a transition period and (iii) a standard stable ABL. By a “standard stable ABL” we refer to a cloud-free, shallow stable ABL in which terrestrial radiation causes a surface-based temperature inversion, comparable to a nocturnal boundary layer at midlatitudes.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3"><?xmltex \currentcnt{3}?><?xmltex \def\figurename{Figure}?><label>Figure 3</label><caption><p id="d1e1851">Time–height contour plots and time series observed at 9 <inline-formula><mml:math id="M120" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> (green) and 65 <inline-formula><mml:math id="M121" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> (blue) altitude for <bold>(a)</bold> wind velocity <inline-formula><mml:math id="M122" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula>, <bold>(b)</bold> wind direction dd, and <bold>(c)</bold> temperature <inline-formula><mml:math id="M123" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> or potential temperature <inline-formula><mml:math id="M124" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula>. The data are from three different sources: (i) the time–height profiles are from BELUGA, (ii) the data shown at 9 <inline-formula><mml:math id="M125" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> altitude are from VRS and (iii) the data at 65 <inline-formula><mml:math id="M126" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> altitude are from the mast. The dashed vertical lines represent the time of the example flights in Fig. <xref ref-type="fig" rid="Ch1.F4"/>.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/15365/2023/acp-23-15365-2023-f03.png"/>

        </fig>

      <p id="d1e1925">In the first period (10:00 to 16:00 UTC), a clear LLJ emerges with a wind velocity maximum of about 10 <inline-formula><mml:math id="M127" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> at a height of around 100 to 150 <inline-formula><mml:math id="M128" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>, while the wind velocity in the near-surface layer remains below 2 <inline-formula><mml:math id="M129" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>. The wind direction in the lower 400 <inline-formula><mml:math id="M130" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> is west to northwest with the highest variability at surface level. A very stable surface layer develops with near-surface temperatures around <inline-formula><mml:math id="M131" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>30 <inline-formula><mml:math id="M132" display="inline"><mml:mrow class="unit"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mi mathvariant="normal">C</mml:mi></mml:mrow></mml:math></inline-formula>, which increase by 10 <inline-formula><mml:math id="M133" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">K</mml:mi></mml:mrow></mml:math></inline-formula> to 100 <inline-formula><mml:math id="M134" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> height. The sharp surface temperature drop at 13:15 is not obvious at 65 <inline-formula><mml:math id="M135" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> altitude and qualitatively correlates with a wind rotation to the north.</p>
      <p id="d1e2024">Between about 16:30 and 21:00 UTC, the wind velocity becomes more constant with height but generally decreases throughout the profile from about 5 <inline-formula><mml:math id="M136" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to less than 2 <inline-formula><mml:math id="M137" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> with a highly variable wind direction. The LLJ disappears almost completely, and the strong surface temperature inversion now extends only to the lowermost tens of meters. At 17:00, however, another smaller LLJ occurs with a maximum at a lower altitude compared to the previous LLJs of the first period. The temperature throughout the entire profile decreases by about 5 to 10 <inline-formula><mml:math id="M138" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">K</mml:mi></mml:mrow></mml:math></inline-formula>. This period is labeled as a<?pagebreak page15369?> transition between the LLJ period and a standard stable ABL structure.</p>
      <p id="d1e2069">After 21:00 UTC, the wind velocity increases with a local maximum at around 100 <inline-formula><mml:math id="M139" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> but is much less pronounced compared to the LLJ structure observed in the first period. The wind direction is almost constant throughout the profile. The temperature profile is similar to the transition phase with a temperature inversion situated at a lower altitude.</p>
      <p id="d1e2080">For a more detailed analysis of the stratification during the three periods, four selected individual profiles of <inline-formula><mml:math id="M140" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M141" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M142" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> are plotted in Fig. <xref ref-type="fig" rid="Ch1.F4"/>. A well-developed LLJ is observed at 14:31 UTC. The wind velocity maximum of 9 <inline-formula><mml:math id="M143" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> in 100 <inline-formula><mml:math id="M144" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> coincides with the top of the strong surface-based temperature inversion. Above this inversion, the ABL is almost adiabatically stratified up to 150 <inline-formula><mml:math id="M145" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> – the region with decreasing wind velocity. This region between 100  and 150 <inline-formula><mml:math id="M146" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> shows the highest local variability in wind velocity, although the wind shear is lower compared to the height range between 40  and 100 <inline-formula><mml:math id="M147" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>. This observation is visible in the local energy dissipation, which is highest – apart from the lowermost surface layer – in the upper part of the LLJ. Here, with almost neutral stratification, the wind shear term in Eq. (<xref ref-type="disp-formula" rid="Ch1.E1"/>) dominates the buoyancy term by 1 order of magnitude. Below the LLJ core at 40 to 100 <inline-formula><mml:math id="M148" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>, the turbulence generation due to wind shear is reduced (compared to above the LLJ) by the influence of the temperature inversion. A local minimum of <inline-formula><mml:math id="M149" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> is located between the surface and the base of the LLJ, where the damping effect of the strong temperature inversion coincides with a height-constant wind velocity. Above 160 <inline-formula><mml:math id="M150" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M151" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> again increases slightly with height and <inline-formula><mml:math id="M152" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> drops by 2 orders of magnitude.</p>
      <p id="d1e2196">The profile observed at 19:12 UTC in the transition phase shows the lowest wind velocity (<inline-formula><mml:math id="M153" display="inline"><mml:mo lspace="0mm">≈</mml:mo></mml:math></inline-formula> 1 <inline-formula><mml:math id="M154" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) between 50 <inline-formula><mml:math id="M155" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> and 160 <inline-formula><mml:math id="M156" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>. The maximum wind velocity is above 3 <inline-formula><mml:math id="M157" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> in the near-surface layer, coinciding with a temperature inversion up to 40 <inline-formula><mml:math id="M158" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>. However, just above this inversion, there is a shallow, 20 <inline-formula><mml:math id="M159" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> thick layer where the wind shear is strong enough to develop some turbulence indicated by a local maximum of <inline-formula><mml:math id="M160" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> up to <inline-formula><mml:math id="M161" display="inline"><mml:mrow><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mrow class="unit"><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mrow></mml:math></inline-formula>. Above and below this maximum, stable stratification in combination with weak wind shear results in reduced turbulence. In the transition phase, one single profile includes an LLJ again, with a slightly different structure than the first period observations (Fig. <xref ref-type="fig" rid="Ch1.F4"/>, upper right). Here, the wind velocity increases with height directly above  the surface with <inline-formula><mml:math id="M162" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mtext>LLJ</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:math></inline-formula> m s<inline-formula><mml:math id="M163" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> at 60 <inline-formula><mml:math id="M164" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> height. The lowermost 20 <inline-formula><mml:math id="M165" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> are neutrally stratified followed by a  temperature inversion at <inline-formula><mml:math id="M166" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">100</mml:mn></mml:mrow></mml:math></inline-formula> m, which coincides with the wind minimum above the LLJ. Here, the <inline-formula><mml:math id="M167" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> profile is different from the LLJ cases, with maximum turbulence at the surface, gradually decreasing towards the upper bound of the LLJ, allowing turbulent mixing between the surface and the LLJ core. The upper part of the LLJ is still located inside the temperature inversion, leading to less turbulence in that region than for the LLJ phase. Above the LLJ, the turbulence intensity remains low and slightly increases at higher altitudes, with the wind velocity increasing again.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4"><?xmltex \currentcnt{4}?><?xmltex \def\figurename{Figure}?><label>Figure 4</label><caption><p id="d1e2380">Example vertical profiles for each phase:  wind velocity <inline-formula><mml:math id="M168" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> with the definition of LLJ strength <inline-formula><mml:math id="M169" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>u</mml:mi></mml:mrow></mml:math></inline-formula> and LLJ height <inline-formula><mml:math id="M170" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mtext>LLJ</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula>, potential temperature <inline-formula><mml:math id="M171" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> with the temperature inversion strength <inline-formula><mml:math id="M172" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:mrow></mml:math></inline-formula> and inversion height <inline-formula><mml:math id="M173" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mtext>i</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> and dissipation rate <inline-formula><mml:math id="M174" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula>. For the transition phase, one profile with and one profile without LLJ is shown.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/15365/2023/acp-23-15365-2023-f04.png"/>

        </fig>

      <p id="d1e2454">In the period with a standard stable ABL structure, the surface layer up to 50 <inline-formula><mml:math id="M175" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> shows the largest increase in <inline-formula><mml:math id="M176" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M177" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula>, thus resulting in generally low values of <inline-formula><mml:math id="M178" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> compared to the previous periods. Above 50 <inline-formula><mml:math id="M179" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>, the ABL is slightly stably stratified with an almost height-constant wind velocity, weakest turbulence at the top of the surface inversion and some variable turbulence at higher altitudes. Throughout the profile, the values for buoyancy and shear increase or decrease to a similar degree, resulting in turbulence predominantly induced by surface roughness.</p>
      <p id="d1e2494">As a next step, we examine the temporal evolution of LLJ and ABL parameters. Figure <xref ref-type="fig" rid="Ch1.F5"/> shows the time series of LLJ and temperature inversion strength and height, as well as <inline-formula><mml:math id="M180" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> between surface and <inline-formula><mml:math id="M181" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> derived from the profiles, combined with continuously measured turbulence parameters at 65 <inline-formula><mml:math id="M182" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>. Except for a few cases in the transition phase, the LLJ strength correlates with the temperature inversion strength (Pearson correlation coefficient <inline-formula><mml:math id="M183" display="inline"><mml:mrow><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.52</mml:mn></mml:mrow></mml:math></inline-formula>). The LLJ height correlates even more clearly with the inversion height (<inline-formula><mml:math id="M184" display="inline"><mml:mrow><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.64</mml:mn></mml:mrow></mml:math></inline-formula>). <inline-formula><mml:math id="M185" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is increased in the transition period and around Ri=1 (weak turbulence) for the other two periods. Energy dissipation <inline-formula><mml:math id="M186" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> in the LLJ period is 1 order of magnitude higher than in phases II and III and shows much more variability. The turbulent fluxes of heat <inline-formula><mml:math id="M187" display="inline"><mml:mi>H</mml:mi></mml:math></inline-formula> and momentum <inline-formula><mml:math id="M188" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi>u</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> show increased values and variability and individual events with high flux magnitude in the LLJ period, compared to low fluxes in the transition and standard ABL phase. The LLJ, as observed at 17:04 UTC, coincides with a short period<?pagebreak page15370?> of weak upward-oriented momentum fluxes and downward heat fluxes.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5"><?xmltex \currentcnt{5}?><?xmltex \def\figurename{Figure}?><label>Figure 5</label><caption><p id="d1e2592">Temporal development of <bold>(a)</bold> the LLJ parameters strength <inline-formula><mml:math id="M189" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>u</mml:mi></mml:mrow></mml:math></inline-formula> and height <inline-formula><mml:math id="M190" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mtext>LLJ</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula>, <bold>(b)</bold> the temperature inversion strength <inline-formula><mml:math id="M191" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:mrow></mml:math></inline-formula> and height <inline-formula><mml:math id="M192" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mtext>i</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula>, <bold>(c)</bold> bulk Richardson number <inline-formula><mml:math id="M193" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (between surface and <inline-formula><mml:math id="M194" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>), <bold>(d)</bold> dissipation rate <inline-formula><mml:math id="M195" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> at 65 <inline-formula><mml:math id="M196" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> (5 <inline-formula><mml:math id="M197" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">min</mml:mi></mml:mrow></mml:math></inline-formula> averages), and <bold>(e)</bold> turbulent fluxes <inline-formula><mml:math id="M198" display="inline"><mml:mi>H</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M199" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mtext>u</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> at 65 <inline-formula><mml:math id="M200" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> (calculated every 5 <inline-formula><mml:math id="M201" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">min</mml:mi></mml:mrow></mml:math></inline-formula> for a 30 <inline-formula><mml:math id="M202" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">min</mml:mi></mml:mrow></mml:math></inline-formula> period).</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/15365/2023/acp-23-15365-2023-f05.png"/>

        </fig>

</sec>
<sec id="Ch1.S4.SS2">
  <label>4.2</label><title>Normalized vertical profiles</title>
      <p id="d1e2756">Each of the three periods introduced in Sect. <xref ref-type="sec" rid="Ch1.S4.SS1"/> features a distinct vertical structure of thermodynamic and turbulence parameters. For each of the three periods, Fig. <xref ref-type="fig" rid="Ch1.F6"/> shows normalized vertical profiles. Whereas the height is normalized with <inline-formula><mml:math id="M203" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M204" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> is normalized with <inline-formula><mml:math id="M205" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mtext>max</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> (below 300 <inline-formula><mml:math id="M206" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>), and <inline-formula><mml:math id="M207" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M208" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> are normalized with their values at <inline-formula><mml:math id="M209" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. The box plots are assigned to height intervals and include median values for each profile within the respective period. The wind velocity profile shows a characteristic shape for the LLJ period and the standard stable ABL period, with low variability within the profiles. The transition phase exhibits much higher variability, and wind velocity increases with height. The <inline-formula><mml:math id="M210" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M211" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> profiles are similar for all phases, with the lowest variability below <inline-formula><mml:math id="M212" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> for the LLJ period. In each period, there is a turbulence maximum in <inline-formula><mml:math id="M213" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> close to the surface, which indicates surface-driven turbulence for all cases. In relation to the wind profile, the vertical turbulence profile is different for each period and will be discussed in more detail.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F6" specific-use="star"><?xmltex \currentcnt{6}?><?xmltex \def\figurename{Figure}?><label>Figure 6</label><caption><p id="d1e2861">Normalized vertical profiles of <inline-formula><mml:math id="M214" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M215" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M216" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M217" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M218" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> for phase I (with LLJ, top), phase II (transition, center) and phase III (profiles without LLJ, bottom). The box plots show variations within the individual <inline-formula><mml:math id="M219" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula> profiles in each phase. The height <inline-formula><mml:math id="M220" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> is normalized with the temperature inversion base height <inline-formula><mml:math id="M221" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. The wind velocity is normalized with <inline-formula><mml:math id="M222" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mtext>max</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> (<inline-formula><mml:math id="M223" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mtext>LLJ</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> or maximum <inline-formula><mml:math id="M224" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> below 300 <inline-formula><mml:math id="M225" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>). The quantities <inline-formula><mml:math id="M226" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M227" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> are normalized with their value at <inline-formula><mml:math id="M228" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. The boxes include the lower and upper quartile values of the data, with the orange line at the median. </p></caption>
          <?xmltex \igopts{width=398.338583pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/15365/2023/acp-23-15365-2023-f06.png"/>

        </fig>

      <p id="d1e2998">In the LLJ period, the average <inline-formula><mml:math id="M229" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> structure has a characteristic shape with two local maxima: near the surface and at <inline-formula><mml:math id="M230" display="inline"><mml:mrow><mml:mi>z</mml:mi><mml:mo>/</mml:mo><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>≈</mml:mo><mml:mn mathvariant="normal">1.5</mml:mn></mml:mrow></mml:math></inline-formula> with reduced values around <inline-formula><mml:math id="M231" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> itself. A local minimum of <inline-formula><mml:math id="M232" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> at <inline-formula><mml:math id="M233" display="inline"><mml:mrow><mml:mi>z</mml:mi><mml:mo>/</mml:mo><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>≈</mml:mo><mml:mn mathvariant="normal">0.3</mml:mn></mml:mrow></mml:math></inline-formula> suggests decoupling of the LLJ from the surface. In the transition period, the average  <inline-formula><mml:math id="M234" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> profile is almost height-constant but shows high internal variability. The <inline-formula><mml:math id="M235" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> profile is highly variable with height without any characteristic structure. The only common features are minimum values around <inline-formula><mml:math id="M236" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and maximum values close to the surface.</p>
      <p id="d1e3091">The last period is characterized by a gradual increase in the normalized velocity from the surface up to <inline-formula><mml:math id="M237" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> followed by a slight decrease above. Turbulence shows comparably high values only close to the surface with a clear minimum around <inline-formula><mml:math id="M238" display="inline"><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, where values close to the resolution limit are observed. The low variability within the profiles results partly from only two profiles with turbulence measurements in this phase. The profiles of <inline-formula><mml:math id="M239" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> generally match the <inline-formula><mml:math id="M240" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> profiles with low <inline-formula><mml:math id="M241" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">Ri</mml:mi><mml:mi mathvariant="normal">g</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, correlating with high <inline-formula><mml:math id="M242" display="inline"><mml:mi mathvariant="italic">ε</mml:mi></mml:math></inline-formula> but generally showing a high variability.</p>
      <p id="d1e3153">These normalized profiles show that the presence of an LLJ enhances turbulent mixing directly above and below the jet maximum compared to a stable ABL without an LLJ. The turbulence increase is more pronounced above the jet core. However, a stably stratified  region close to the surface with height-constant wind speeds decouples the LLJ from the surface. The enhanced turbulent mixing by the LLJ might be important for the vertical mixing of advected long-range-transported tracers.</p>
</sec>
</sec>
<sec id="Ch1.S5" sec-type="conclusions">
  <label>5</label><title>Summary and discussion</title>
      <p id="d1e3166">This study presents the observation of an LLJ based on tethered-balloon measurements in the late-winter central Arctic in March 2018 in northeast Greenland at the Villum Research Station/Station Nord. The measurements span a 15 h period with a transition from a stable ABL with a prominent LLJ to a classically stable nocturnal ABL without an LLJ. The observations include measurements of mean standard meteorological measurements as well as turbulence measurements from the ground to well above the ABL (typically 300 <inline-formula><mml:math id="M243" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula>). The   balloon observations were supported by continuous turbulence measurements on a mast at 65 <inline-formula><mml:math id="M244" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:math></inline-formula> height – corresponding to the lower part of the LLJ.</p>
      <p id="d1e3185">During the LLJ phase, observations indicate increased turbulence (increased local dissipation rates) within the LLJ, but the increased turbulence does not reach the surface. Only during the transition phase to a more classical boundary layer structure do these increased turbulence intensities also reach the surface, allowing thorough mixing to occur. These observations lead to the hypothesis that within an  LLJ, a passive tracer can be “trapped” and transported over a long distance, without vertical mixing beyond the LLJ boundaries greatly reducing the tracer concentration. Only the dissolution of the LLJ can finally lead to an increased concentration also at the ground.</p>
      <?pagebreak page15372?><p id="d1e3188">An accompanying analytical model  study <xref ref-type="bibr" rid="bib1.bibx10" id="paren.20"/> shows that LLJs may impact  aerosol formation and evolution under Arctic conditions. However, the one-dimensional analytical LLJ model revealed difficulties in the reproduction of the observed LLJ wind peak (underestimation of the sharpness of the observed vertical wind velocity), which leads to a misprediction of the shear-induced contribution to the TKE budget. As a consequence, the TKE production at the LLJ level is underestimated. The removal of this shortcoming would require the consideration of momentum advection, baroclinicity and terrain-induced cold-air flow. Nevertheless, within the framework of this conceptual study, we demonstrate that an LLJ can effectively promote the horizontal transport of a passive tracer released close to the ground. This result, in turn, suggests that an LLJ may have a significant impact on both in situ and ex situ aerosol formation and evolution under Arctic conditions. In future research, the simulation of the intrinsically three-dimensional nature of this phenomenon with a more sophisticated model requires further investigation through appropriate observations of the parameters involved, including vertical profiling.</p>
</sec>

      
      </body>
    <back><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d1e3198">The observational data for this paper are publicly available (<xref ref-type="bibr" rid="bib1.bibx8" id="altparen.21"/>, <ext-link xlink:href="https://doi.org/10.1594/PANGAEA.900240" ext-link-type="DOI">10.1594/PANGAEA.900240</ext-link>).</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d1e3210">UE and HS performed the balloon measurements and processed and analyzed the data. LLS contributed the meteorological mast measurements. OH contributed to the data evaluation and discussion. All co-authors contributed to drafting the paper.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d1e3216">The contact author has declared that none of the authors has any competing interests.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d1e3222">Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.</p>
  </notes><ack><title>Acknowledgements</title><p id="d1e3228">We gratefully acknowledge the funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – project number 268020496 – TRR 172, within the Transregional Collaborative Research Center “ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes, and Feedback Mechanisms (AC)<inline-formula><mml:math id="M245" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:math></inline-formula>” in sub-project A02. We acknowledge the essential help of Andreas Herber, Jens Voigtländer, Frank Stratmann, PAMARCMiP scientists, and the crew of Station Nord and Villum Research Station during the balloon operations. Sven-Erik Gryning provided and advised us on ceilometer and wind lidar data at VRS. We wish to thank Andreas Massling, Daniel Charles Thomas and Jakob Boyd Pernov for discussing VRS data.</p></ack><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d1e3242">This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. 268020496 – TRR 172).<?xmltex \hack{\newline}?><?xmltex \hack{\newline}?>The publication of this article was funded by the <?xmltex \notforhtml{\newline}?> Open Access Fund of the Leibniz Association.</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d1e3253">This paper was edited by Geraint Vaughan and reviewed by Geraint Vaughan and one anonymous referee.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bibx1"><?xmltex \def\ref@label{{Abarbanel et~al.(1984)Abarbanel, Holm, Marsden, and
Ratiu}}?><label>Abarbanel et al.(1984)Abarbanel, Holm, Marsden, and Ratiu</label><?label Abarbanel1984?><mixed-citation>Abarbanel, H. D. I., Holm, D. D., Marsden, J. E., and Ratiu, T.: Richardson number criterion for the nonlinear stability of three-dimensional stratified flow, Phys. Rev. Lett., 52, 2352–2355, <ext-link xlink:href="https://doi.org/10.1103/PhysRevLett.52.2352" ext-link-type="DOI">10.1103/PhysRevLett.52.2352</ext-link>, 1984.</mixed-citation></ref>
      <ref id="bib1.bibx2"><?xmltex \def\ref@label{{Algarra et~al.(2019)Algarra, Eiras-Barca, Miguez-Macho, Nieto, and
Gimeno}}?><label>Algarra et al.(2019)Algarra, Eiras-Barca, Miguez-Macho, Nieto, and Gimeno</label><?label algarra:2019?><mixed-citation>Algarra, I., Eiras-Barca, J., Miguez-Macho, G., Nieto, R., and Gimeno, L.: On the assessment of the moisture transport by the Great Plains low-level jet, Earth Syst. Dynam., 10, 107–119, <ext-link xlink:href="https://doi.org/10.5194/esd-10-107-2019" ext-link-type="DOI">10.5194/esd-10-107-2019</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx3"><?xmltex \def\ref@label{{Andreas et~al.(2000)Andreas, Claffey, and Makshtas}}?><label>Andreas et al.(2000)Andreas, Claffey, and Makshtas</label><?label Andreas2000?><mixed-citation>Andreas, E. L., Claffey, K. J., and Makshtas, A. P.: Low-level atmospheric jets and inversions over the western Weddell Sea, Bound.-Lay. Meteorol., 97, 459–486, <ext-link xlink:href="https://doi.org/10.1023/A:1002793831076" ext-link-type="DOI">10.1023/A:1002793831076</ext-link>, 2000.</mixed-citation></ref>
      <ref id="bib1.bibx4"><?xmltex \def\ref@label{{Banta et~al.(2006)Banta, Pichugina, and Brewer}}?><label>Banta et al.(2006)Banta, Pichugina, and Brewer</label><?label Banta2006?><mixed-citation>Banta, R. M., Pichugina, Y. L., and Brewer, W. A.: Turbulent Velocity-Variance Profiles in the Stable Boundary Layer Generated by a Nocturnal Low-Level Jet, J. Atmos. Sci., 63, 2700–2719, <ext-link xlink:href="https://doi.org/10.1175/JAS3776.1" ext-link-type="DOI">10.1175/JAS3776.1</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bibx5"><?xmltex \def\ref@label{{Blackadar(1957)}}?><label>Blackadar(1957)</label><?label Blackadar1957?><mixed-citation>Blackadar, A. K.: Boundary Layer Wind Maxima and Their Significance for the Growth of Nocturnal Inversions, B. Am. Meteorol. Soc., 38, 283–290, <ext-link xlink:href="https://doi.org/10.1175/1520-0477-38.5.283" ext-link-type="DOI">10.1175/1520-0477-38.5.283</ext-link>, 1957.</mixed-citation></ref>
      <ref id="bib1.bibx6"><?xmltex \def\ref@label{{Chechin and Lüpkes(2019)}}?><label>Chechin and Lüpkes(2019)</label><?label Chechin_2019?><mixed-citation>Chechin, D. G. and Lüpkes, C.: Baroclinic low-level jets in Arctic marine cold-air outbreaks, IOP Conf. Ser. Earth Env., 231, 012011, <ext-link xlink:href="https://doi.org/10.1088/1755-1315/231/1/012011" ext-link-type="DOI">10.1088/1755-1315/231/1/012011</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx7"><?xmltex \def\ref@label{{Egerer et~al.(2019a)Egerer, Gottschalk, Siebert, Ehrlich, and
Wendisch}}?><label>Egerer et al.(2019a)Egerer, Gottschalk, Siebert, Ehrlich, and Wendisch</label><?label Egerer2019?><mixed-citation>Egerer, U., Gottschalk, M., Siebert, H., Ehrlich, A., and Wendisch, M.: The new BELUGA setup for collocated turbulence and radiation measurements using a tethered balloon: first applications in the cloudy Arctic boundary layer, Atmos. Meas. Tech., 12, 4019–4038, <ext-link xlink:href="https://doi.org/10.5194/amt-12-4019-2019" ext-link-type="DOI">10.5194/amt-12-4019-2019</ext-link>, 2019a.</mixed-citation></ref>
      <ref id="bib1.bibx8"><?xmltex \def\ref@label{{{Egerer} et~al.(2019b){Egerer}, {Siebert}, {Voigtl\"{a}nder}, and
{Gottschalk}}}?><label>Egerer et al.(2019b)Egerer, Siebert, Voigtländer, and Gottschalk</label><?label Egerer2019_pangaea_villum?><mixed-citation>Egerer, U., Siebert, H., Voigtländer, J., and Gottschalk, M.: Tethered balloon-borne measurements of turbulence and radiation during the Arctic field campaign PAMARCMiP in March/April 2018, PANGAEA [data set], <ext-link xlink:href="https://doi.org/10.1594/PANGAEA.900240" ext-link-type="DOI">10.1594/PANGAEA.900240</ext-link>, 2019b.</mixed-citation></ref>
      <ref id="bib1.bibx9"><?xmltex \def\ref@label{{Guest et~al.(2018)Guest, Persson, Wang, Jordan, Jin, Blomquist, and
Fairall}}?><label>Guest et al.(2018)Guest, Persson, Wang, Jordan, Jin, Blomquist, and Fairall</label><?label Guest2018?><mixed-citation>Guest, P., Persson, P. O. G., Wang, S., Jordan, M., Jin, Y., Blomquist, B., and Fairall, C.: Low-Level Baroclinic Jets Over the New Arctic Ocean, J. Geophys. Res.-Oceans, 123, 4074–4091, <ext-link xlink:href="https://doi.org/10.1002/2018JC013778" ext-link-type="DOI">10.1002/2018JC013778</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx10"><?xmltex \def\ref@label{{Hellmuth et~al.(2023)Hellmuth, Egerer, Siebert, and
Sørensen}}?><label>Hellmuth et al.(2023)Hellmuth, Egerer, Siebert, and Sørensen</label><?label hellmuth:2023?><mixed-citation>Hellmuth, O., Egerer, U., Siebert, H., and Sørensen, L. L.: PAMARCMiP Contribution: An analytic model companion based on observations: The role of low-level jets in the advection of passive tracers in the high Arctic, Tech. rep., Zenodo, <ext-link xlink:href="https://doi.org/10.5281/zenodo.7689308" ext-link-type="DOI">10.5281/zenodo.7689308</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bibx11"><?xmltex \def\ref@label{{Jakobson et~al.(2013)Jakobson, Vihma, Jakobson, Palo, M\"{a}nnik, and
Jaagus}}?><label>Jakobson et al.(2013)Jakobson, Vihma, Jakobson, Palo, Männik, and Jaagus</label><?label Jakobson2013?><mixed-citation>Jakobson, L., Vihma, T<?pagebreak page15373?>., Jakobson, E., Palo, T., Männik, A., and Jaagus, J.: Low-level jet characteristics over the Arctic Ocean in spring and summer, Atmos. Chem. Phys., 13, 11089–11099, <ext-link xlink:href="https://doi.org/10.5194/acp-13-11089-2013" ext-link-type="DOI">10.5194/acp-13-11089-2013</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx12"><?xmltex \def\ref@label{{L{\'{o}}pez-Garc{\'{\i}}a et~al.(2022)L{\'{o}}pez-Garc{\'{\i}}a, Neely,
Dahlke, and Brooks}}?><label>López-García et al.(2022)López-García, Neely, Dahlke, and Brooks</label><?label lopez2022?><mixed-citation>López-García, V., Neely, R. R., Dahlke, S., and Brooks, I. M.: Low-level jets over the Arctic Ocean during MOSAiC, Elementa: Science of the Anthropocene, 10, 00063, <ext-link xlink:href="https://doi.org/10.1525/elementa.2022.00063" ext-link-type="DOI">10.1525/elementa.2022.00063</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bibx13"><?xmltex \def\ref@label{{Mahrt(1985)}}?><label>Mahrt(1985)</label><?label mahrt:85?><mixed-citation> Mahrt, L.: Vertical structure and turbulence in the very stable boundary layer, J. Atmos. Sci., 42, 2333–2349, 1985.</mixed-citation></ref>
      <ref id="bib1.bibx14"><?xmltex \def\ref@label{{Mahrt(1999)}}?><label>Mahrt(1999)</label><?label mahrt:1999?><mixed-citation> Mahrt, L.: Stratified Atmospheric Boundary Layers, Bound.-Lay. Meteorol., 90, 375–396, 1999.</mixed-citation></ref>
      <ref id="bib1.bibx15"><?xmltex \def\ref@label{{Miles(1961)}}?><label>Miles(1961)</label><?label Miles1961?><mixed-citation>Miles, J. W.: On the stability of heterogeneous shear flows, J. Fluid Mech., 10, 496–508, <ext-link xlink:href="https://doi.org/10.1017/S0022112061000305" ext-link-type="DOI">10.1017/S0022112061000305</ext-link>, 1961.</mixed-citation></ref>
      <ref id="bib1.bibx16"><?xmltex \def\ref@label{{Siebert et~al.(2006)Siebert, Lehmann, and Wendisch}}?><label>Siebert et al.(2006)Siebert, Lehmann, and Wendisch</label><?label siebert:06?><mixed-citation> Siebert, H., Lehmann, K., and Wendisch, M.: Observations of small scale turbulence and energy dissipation rates in the cloudy boundary layer, J. Atmos. Sci., 63, 1451–1466, 2006.</mixed-citation></ref>
      <ref id="bib1.bibx17"><?xmltex \def\ref@label{{Smedman et~al.(1993)Smedman, Tjernstr{\"{o}}m, and
H{\"{o}}grstr{\"{o}}m}}?><label>Smedman et al.(1993)Smedman, Tjernström, and Högrström</label><?label smedman:1993?><mixed-citation>Smedman, A.-S., Tjernström, M., and Högrström, U.: Analysis of the turbulence structure of a marine Low-Level Jet, Bound.-Lay. Meteorol., 66, 195–126, 1993.  </mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bibx18"><?xmltex \def\ref@label{{Stensrud(1996)}}?><label>Stensrud(1996)</label><?label stensrud:1996?><mixed-citation> Stensrud, D. J.: Importance of Low-Level Jets to Climate: A Review, J. Climate, 9, 1698–1711, 1996.</mixed-citation></ref>
      <ref id="bib1.bibx19"><?xmltex \def\ref@label{{Tuononen et~al.(2015)Tuononen, Sinclair, and Vihma}}?><label>Tuononen et al.(2015)Tuononen, Sinclair, and Vihma</label><?label Tuononen2015?><mixed-citation>Tuononen, M., Sinclair, V. A., and Vihma, T.: A climatology of low-level jets in the mid-latitudes and polar regions of the Northern Hemisphere, Atmos. Sci. Lett., 16, 492–499, <ext-link xlink:href="https://doi.org/10.1002/asl.587" ext-link-type="DOI">10.1002/asl.587</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx20"><?xmltex \def\ref@label{{Vihma et~al.(2011)Vihma, Kilpeläinen, Manninen, Sjöblom, Jakobson,
Palo, Jaagus, and Maturilli}}?><label>Vihma et al.(2011)Vihma, Kilpeläinen, Manninen, Sjöblom, Jakobson, Palo, Jaagus, and Maturilli</label><?label Vihma2011?><mixed-citation>Vihma, T., Kilpeläinen, T., Manninen, M., Sjöblom, A., Jakobson, E., Palo, T., Jaagus, J., and Maturilli, M.: Characteristics of Temperature and Humidity Inversions and Low-Level Jets over Svalbard Fjords in Spring, Adv. Meteorol., 2011, 1–14, <ext-link xlink:href="https://doi.org/10.1155/2011/486807" ext-link-type="DOI">10.1155/2011/486807</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx21"><?xmltex \def\ref@label{{Wilczak et~al.(2001)Wilczak, Oncley, and Stage}}?><label>Wilczak et al.(2001)Wilczak, Oncley, and Stage</label><?label Wilczak2001?><mixed-citation>Wilczak, J. M., Oncley, S. P., and Stage, S. A.: Sonic anemometer tilt correction algorithms, Bound.-Lay. Meteorol., 99, 127–150, <ext-link xlink:href="https://doi.org/10.1023/A:1018966204465" ext-link-type="DOI">10.1023/A:1018966204465</ext-link>, 2001.</mixed-citation></ref>
      <ref id="bib1.bibx22"><?xmltex \def\ref@label{{Wyngaard(2010)}}?><label>Wyngaard(2010)</label><?label Wyngaard2010?><mixed-citation>Wyngaard, J. C.: Turbulence in the atmosphere, Cambridge University Press, <ext-link xlink:href="https://doi.org/10.1017/CBO9780511840524" ext-link-type="DOI">10.1017/CBO9780511840524</ext-link>, 2010.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>The role of a low-level jet for stirring the stable atmospheric surface layer in the Arctic</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>Abarbanel et al.(1984)Abarbanel, Holm, Marsden, and
Ratiu</label><mixed-citation>
      
Abarbanel, H. D. I., Holm, D. D., Marsden, J. E., and Ratiu, T.: Richardson
number criterion for the nonlinear stability of three-dimensional stratified
flow, Phys. Rev. Lett., 52, 2352–2355, <a href="https://doi.org/10.1103/PhysRevLett.52.2352" target="_blank">https://doi.org/10.1103/PhysRevLett.52.2352</a>,
1984.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>Algarra et al.(2019)Algarra, Eiras-Barca, Miguez-Macho, Nieto, and
Gimeno</label><mixed-citation>
      
Algarra, I., Eiras-Barca, J., Miguez-Macho, G., Nieto, R., and Gimeno, L.: On the assessment of the moisture transport by the Great Plains low-level jet, Earth Syst. Dynam., 10, 107–119, <a href="https://doi.org/10.5194/esd-10-107-2019" target="_blank">https://doi.org/10.5194/esd-10-107-2019</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>Andreas et al.(2000)Andreas, Claffey, and Makshtas</label><mixed-citation>
      
Andreas, E. L., Claffey, K. J., and Makshtas, A. P.: Low-level atmospheric jets
and inversions over the western Weddell Sea, Bound.-Lay. Meteorol.,
97, 459–486, <a href="https://doi.org/10.1023/A:1002793831076" target="_blank">https://doi.org/10.1023/A:1002793831076</a>, 2000.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>Banta et al.(2006)Banta, Pichugina, and Brewer</label><mixed-citation>
      
Banta, R. M., Pichugina, Y. L., and Brewer, W. A.: Turbulent Velocity-Variance
Profiles in the Stable Boundary Layer Generated by a Nocturnal Low-Level Jet,
J. Atmos. Sci., 63, 2700–2719, <a href="https://doi.org/10.1175/JAS3776.1" target="_blank">https://doi.org/10.1175/JAS3776.1</a>,
2006.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>Blackadar(1957)</label><mixed-citation>
      
Blackadar, A. K.: Boundary Layer Wind Maxima and Their Significance for the
Growth of Nocturnal Inversions, B. Am. Meteorol.
Soc., 38, 283–290, <a href="https://doi.org/10.1175/1520-0477-38.5.283" target="_blank">https://doi.org/10.1175/1520-0477-38.5.283</a>, 1957.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>Chechin and Lüpkes(2019)</label><mixed-citation>
      
Chechin, D. G. and Lüpkes, C.: Baroclinic low-level jets in Arctic marine
cold-air outbreaks, IOP Conf. Ser. Earth Env.,
231, 012011, <a href="https://doi.org/10.1088/1755-1315/231/1/012011" target="_blank">https://doi.org/10.1088/1755-1315/231/1/012011</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>Egerer et al.(2019a)Egerer, Gottschalk, Siebert, Ehrlich, and
Wendisch</label><mixed-citation>
      
Egerer, U., Gottschalk, M., Siebert, H., Ehrlich, A., and Wendisch, M.: The new BELUGA setup for collocated turbulence and radiation measurements using a tethered balloon: first applications in the cloudy Arctic boundary layer, Atmos. Meas. Tech., 12, 4019–4038, <a href="https://doi.org/10.5194/amt-12-4019-2019" target="_blank">https://doi.org/10.5194/amt-12-4019-2019</a>, 2019a.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>Egerer et al.(2019b)Egerer, Siebert, Voigtländer, and
Gottschalk</label><mixed-citation>
      
Egerer, U., Siebert, H., Voigtländer, J., and Gottschalk, M.:
Tethered balloon-borne measurements of turbulence and radiation during the
Arctic field campaign PAMARCMiP in March/April 2018, PANGAEA [data set],
<a href="https://doi.org/10.1594/PANGAEA.900240" target="_blank">https://doi.org/10.1594/PANGAEA.900240</a>, 2019b.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>Guest et al.(2018)Guest, Persson, Wang, Jordan, Jin, Blomquist, and
Fairall</label><mixed-citation>
      
Guest, P., Persson, P. O. G., Wang, S., Jordan, M., Jin, Y., Blomquist, B., and
Fairall, C.: Low-Level Baroclinic Jets Over the New Arctic Ocean, J.
Geophys. Res.-Oceans, 123, 4074–4091,
<a href="https://doi.org/10.1002/2018JC013778" target="_blank">https://doi.org/10.1002/2018JC013778</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>Hellmuth et al.(2023)Hellmuth, Egerer, Siebert, and
Sørensen</label><mixed-citation>
      
Hellmuth, O., Egerer, U., Siebert, H., and Sørensen, L. L.: PAMARCMiP
Contribution: An analytic model companion based on observations: The role of
low-level jets in the advection of passive tracers in the high Arctic, Tech.
rep., Zenodo, <a href="https://doi.org/10.5281/zenodo.7689308" target="_blank">https://doi.org/10.5281/zenodo.7689308</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>Jakobson et al.(2013)Jakobson, Vihma, Jakobson, Palo, Männik, and
Jaagus</label><mixed-citation>
      
Jakobson, L., Vihma, T., Jakobson, E., Palo, T., Männik, A., and Jaagus, J.: Low-level jet characteristics over the Arctic Ocean in spring and summer, Atmos. Chem. Phys., 13, 11089–11099, <a href="https://doi.org/10.5194/acp-13-11089-2013" target="_blank">https://doi.org/10.5194/acp-13-11089-2013</a>, 2013.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>López-García et al.(2022)López-García, Neely,
Dahlke, and Brooks</label><mixed-citation>
      
López-García, V., Neely, R. R., Dahlke, S., and Brooks, I. M.:
Low-level jets over the Arctic Ocean during MOSAiC, Elementa: Science
of the Anthropocene, 10, 00063, <a href="https://doi.org/10.1525/elementa.2022.00063" target="_blank">https://doi.org/10.1525/elementa.2022.00063</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>Mahrt(1985)</label><mixed-citation>
      
Mahrt, L.: Vertical structure and turbulence in the very stable boundary layer,
J. Atmos. Sci., 42, 2333–2349, 1985.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>Mahrt(1999)</label><mixed-citation>
      
Mahrt, L.: Stratified Atmospheric Boundary Layers, Bound.-Lay. Meteorol., 90,
375–396, 1999.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>Miles(1961)</label><mixed-citation>
      
Miles, J. W.: On the stability of heterogeneous shear flows, J. Fluid
Mech., 10, 496–508, <a href="https://doi.org/10.1017/S0022112061000305" target="_blank">https://doi.org/10.1017/S0022112061000305</a>, 1961.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>Siebert et al.(2006)Siebert, Lehmann, and Wendisch</label><mixed-citation>
      
Siebert, H., Lehmann, K., and Wendisch, M.: Observations of small scale
turbulence and energy dissipation rates in the cloudy boundary layer, J.
Atmos. Sci., 63, 1451–1466, 2006.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>Smedman et al.(1993)Smedman, Tjernström, and
Högrström</label><mixed-citation>
      
Smedman, A.-S., Tjernström, M., and Högrström, U.: Analysis of the
turbulence structure of a marine Low-Level Jet, Bound.-Lay. Meteorol., 66,
195–126, 1993.


    </mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>Stensrud(1996)</label><mixed-citation>
      
Stensrud, D. J.: Importance of Low-Level Jets to Climate: A Review, J.
Climate, 9, 1698–1711, 1996.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>Tuononen et al.(2015)Tuononen, Sinclair, and Vihma</label><mixed-citation>
      
Tuononen, M., Sinclair, V. A., and Vihma, T.: A climatology of low-level jets
in the mid-latitudes and polar regions of the Northern Hemisphere,
Atmos. Sci. Lett., 16, 492–499,
<a href="https://doi.org/10.1002/asl.587" target="_blank">https://doi.org/10.1002/asl.587</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>Vihma et al.(2011)Vihma, Kilpeläinen, Manninen, Sjöblom, Jakobson,
Palo, Jaagus, and Maturilli</label><mixed-citation>
      
Vihma, T., Kilpeläinen, T., Manninen, M., Sjöblom, A., Jakobson, E., Palo,
T., Jaagus, J., and Maturilli, M.: Characteristics of Temperature and
Humidity Inversions and Low-Level Jets over Svalbard Fjords in Spring,
Adv. Meteorol., 2011, 1–14, <a href="https://doi.org/10.1155/2011/486807" target="_blank">https://doi.org/10.1155/2011/486807</a>, 2011.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>Wilczak et al.(2001)Wilczak, Oncley, and Stage</label><mixed-citation>
      
Wilczak, J. M., Oncley, S. P., and Stage, S. A.: Sonic anemometer tilt
correction algorithms, Bound.-Lay. Meteorol., 99, 127–150,
<a href="https://doi.org/10.1023/A:1018966204465" target="_blank">https://doi.org/10.1023/A:1018966204465</a>, 2001.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>Wyngaard(2010)</label><mixed-citation>
      
Wyngaard, J. C.: Turbulence in the atmosphere, Cambridge University Press,
<a href="https://doi.org/10.1017/CBO9780511840524" target="_blank">https://doi.org/10.1017/CBO9780511840524</a>, 2010.

    </mixed-citation></ref-html>--></article>
