On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width
Abstract The standard model (e.g., Hocking in Earth Planets Space 51:525–541, 1999), $$\varepsilon = c_{0} \sigma^{2} N$$ ε=c0σ2N , (where $$\sigma$$ σ is the radar spectral width assumed to be equal to vertical turbulence velocity fluctuation $$\sqrt {\overline{{w^{2} }} }$$ w2¯ , N is the buoyancy...
| Main Authors: | , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
SpringerOpen
2018-12-01
|
| Series: | Earth, Planets and Space |
| Subjects: | |
| Online Access: | http://link.springer.com/article/10.1186/s40623-018-0957-7 |
| _version_ | 1818363370044129280 |
|---|---|
| author | Lakshmi Kantha Hubert Luce Hiroyuki Hashiguchi |
| author_facet | Lakshmi Kantha Hubert Luce Hiroyuki Hashiguchi |
| author_sort | Lakshmi Kantha |
| collection | DOAJ |
| description | Abstract The standard model (e.g., Hocking in Earth Planets Space 51:525–541, 1999), $$\varepsilon = c_{0} \sigma^{2} N$$ ε=c0σ2N , (where $$\sigma$$ σ is the radar spectral width assumed to be equal to vertical turbulence velocity fluctuation $$\sqrt {\overline{{w^{2} }} }$$ w2¯ , N is the buoyancy frequency, and c 0 is a constant), derived from Weinstock (J Atmos Sci 35:1022–1027, 1978; J Atmos Sci 38:880–883, 1981) formulation, has been used extensively for estimating the turbulence kinetic energy (TKE) dissipation rate $$\varepsilon$$ ε under stable stratification from VHF radar Doppler spectral width $$\sigma$$ σ . The Weinstock model can be derived by simply integrating the TKE spectrum in the wavenumber space from the buoyancy wavenumber $$k_{\text{B}} = \frac{N}{\sigma }$$ kB=Nσ to $$\infty$$ ∞ . However, it ignores the radar volume dimensions and hence its spatial weighting characteristics. Labitt (Some basic relations concerning the radar measurements of air turbulence, MIT Lincoln Laboratory, ATC Working Paper NO 46WP-5001, 1979) and White et al. (J Atmos Ocean Technol 16:1967–1972, 1999) formulations do take into account the radar spatial weighting characteristics, but assume that the wavenumber range in the integration of TKE spectrum extends from 0 to $$\infty$$ ∞ . The White et al. model accounts for wind speed effects, whereas the other two do not. More importantly, all three formulations make the assumption that k −5/3 spectral shape of TKE spectrum extends across the entire wavenumber range of integration. It is traditional to use Weinstock formulation for $$k_{\text{B}}^{ - 1} < 2a,2b$$ kB-1<2a,2b (where a and b are radar volume dimensions in the horizontal and vertical directions) and White et al. formulation (without wind advection) for $$k_{\text{B}}^{ - 1} > 2a,2b$$ kB-1>2a,2b . However, there is no need to invoke these asymptotic limits. We present here a numerical model, which is valid for all values of buoyancy wavenumber $$k_{\text{B}}$$ kB and transitions from $$\varepsilon \sim\sigma^{2}$$ ε∼σ2 behavior at lower values of $$\sigma$$ σ in accordance with Weinstock’s model, to $$\varepsilon \sim\sigma^{3}$$ ε∼σ3 at higher values of $$\sigma$$ σ , in agreement with Chen (J Atmos Sci 31:2222–2225, 1974) and Bertin et al. (Radio Sci 32:791–804, 1997). It can also account for the effects of wind speed, as well the beam width and altitude. Following Hocking (J Atmos Terr Phys 48:655–670, 1986, Earth Planets Space 51:525–541, 1999), the model also takes into account contributions of velocity fluctuations beyond the inertial subrange. The model has universal applicability and can also be applied to convective turbulence in the atmospheric column. It can also be used to explore the parameter space and hence the influence of various parameters and assumptions on the extracted $$\varepsilon$$ ε values. In this note, we demonstrate the utility of the numerical model and make available a MATLAB code of the model for potential use by the radar community. The model results are also compared against in situ turbulence measurements using an unmanned aerial vehicle (UAV) flown in the vicinity of the MU radar in Shigaraki, Japan, during the ShUREX 2016 campaign. |
| first_indexed | 2024-12-13T21:47:23Z |
| format | Article |
| id | doaj.art-2860913b8afa494cae92a5ea9a0a059c |
| institution | Directory Open Access Journal |
| issn | 1880-5981 |
| language | English |
| last_indexed | 2024-12-13T21:47:23Z |
| publishDate | 2018-12-01 |
| publisher | SpringerOpen |
| record_format | Article |
| series | Earth, Planets and Space |
| spelling | doaj.art-2860913b8afa494cae92a5ea9a0a059c2022-12-21T23:30:22ZengSpringerOpenEarth, Planets and Space1880-59812018-12-0170111810.1186/s40623-018-0957-7On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral widthLakshmi Kantha0Hubert Luce1Hiroyuki Hashiguchi2Department of Aerospace Engineering Sciences, University of ColoradoMediterranean Institute of Oceanography (MIO), UM 110, UMR 7294, Université de ToulonResearch Institute for Sustainable Humanosphere, Kyoto UniversityAbstract The standard model (e.g., Hocking in Earth Planets Space 51:525–541, 1999), $$\varepsilon = c_{0} \sigma^{2} N$$ ε=c0σ2N , (where $$\sigma$$ σ is the radar spectral width assumed to be equal to vertical turbulence velocity fluctuation $$\sqrt {\overline{{w^{2} }} }$$ w2¯ , N is the buoyancy frequency, and c 0 is a constant), derived from Weinstock (J Atmos Sci 35:1022–1027, 1978; J Atmos Sci 38:880–883, 1981) formulation, has been used extensively for estimating the turbulence kinetic energy (TKE) dissipation rate $$\varepsilon$$ ε under stable stratification from VHF radar Doppler spectral width $$\sigma$$ σ . The Weinstock model can be derived by simply integrating the TKE spectrum in the wavenumber space from the buoyancy wavenumber $$k_{\text{B}} = \frac{N}{\sigma }$$ kB=Nσ to $$\infty$$ ∞ . However, it ignores the radar volume dimensions and hence its spatial weighting characteristics. Labitt (Some basic relations concerning the radar measurements of air turbulence, MIT Lincoln Laboratory, ATC Working Paper NO 46WP-5001, 1979) and White et al. (J Atmos Ocean Technol 16:1967–1972, 1999) formulations do take into account the radar spatial weighting characteristics, but assume that the wavenumber range in the integration of TKE spectrum extends from 0 to $$\infty$$ ∞ . The White et al. model accounts for wind speed effects, whereas the other two do not. More importantly, all three formulations make the assumption that k −5/3 spectral shape of TKE spectrum extends across the entire wavenumber range of integration. It is traditional to use Weinstock formulation for $$k_{\text{B}}^{ - 1} < 2a,2b$$ kB-1<2a,2b (where a and b are radar volume dimensions in the horizontal and vertical directions) and White et al. formulation (without wind advection) for $$k_{\text{B}}^{ - 1} > 2a,2b$$ kB-1>2a,2b . However, there is no need to invoke these asymptotic limits. We present here a numerical model, which is valid for all values of buoyancy wavenumber $$k_{\text{B}}$$ kB and transitions from $$\varepsilon \sim\sigma^{2}$$ ε∼σ2 behavior at lower values of $$\sigma$$ σ in accordance with Weinstock’s model, to $$\varepsilon \sim\sigma^{3}$$ ε∼σ3 at higher values of $$\sigma$$ σ , in agreement with Chen (J Atmos Sci 31:2222–2225, 1974) and Bertin et al. (Radio Sci 32:791–804, 1997). It can also account for the effects of wind speed, as well the beam width and altitude. Following Hocking (J Atmos Terr Phys 48:655–670, 1986, Earth Planets Space 51:525–541, 1999), the model also takes into account contributions of velocity fluctuations beyond the inertial subrange. The model has universal applicability and can also be applied to convective turbulence in the atmospheric column. It can also be used to explore the parameter space and hence the influence of various parameters and assumptions on the extracted $$\varepsilon$$ ε values. In this note, we demonstrate the utility of the numerical model and make available a MATLAB code of the model for potential use by the radar community. The model results are also compared against in situ turbulence measurements using an unmanned aerial vehicle (UAV) flown in the vicinity of the MU radar in Shigaraki, Japan, during the ShUREX 2016 campaign.http://link.springer.com/article/10.1186/s40623-018-0957-7VHF radarMU radarUnmanned aerial vehicles (UAV)Turbulence kinetic energy (TKE)TKE dissipation rateRadar model |
| spellingShingle | Lakshmi Kantha Hubert Luce Hiroyuki Hashiguchi On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width Earth, Planets and Space VHF radar MU radar Unmanned aerial vehicles (UAV) Turbulence kinetic energy (TKE) TKE dissipation rate Radar model |
| title | On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width |
| title_full | On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width |
| title_fullStr | On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width |
| title_full_unstemmed | On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width |
| title_short | On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width |
| title_sort | on a numerical model for extracting tke dissipation rate from very high frequency vhf radar spectral width |
| topic | VHF radar MU radar Unmanned aerial vehicles (UAV) Turbulence kinetic energy (TKE) TKE dissipation rate Radar model |
| url | http://link.springer.com/article/10.1186/s40623-018-0957-7 |
| work_keys_str_mv | AT lakshmikantha onanumericalmodelforextractingtkedissipationratefromveryhighfrequencyvhfradarspectralwidth AT hubertluce onanumericalmodelforextractingtkedissipationratefromveryhighfrequencyvhfradarspectralwidth AT hiroyukihashiguchi onanumericalmodelforextractingtkedissipationratefromveryhighfrequencyvhfradarspectralwidth |