Uncertainties of drag coefficient estimates above sea ice from field data

Surface turbulent exchanges play a key role on sea ice dynamics, on ocean and sea ice heat budgets and on the polar atmosphere. Uncertainties in parameterizations of surface turbulent fluxes are mostly held by the transfer coefficients and estimates of those transfer coefficients from field data are required for parameterization development. Measurement errors propagate through the computation of transfer coefficients and contribute to its total error together with the uncertainties in the empirical stability functions used to correct for stability effects. Here we propose a methodology to assess their contributions individually to each coefficient estimate as well as the total drag coefficient uncertainty and we apply this methodology on the example of the SHEBA campaign. We conclude that for most common drag coefficient values (between 1.0×10 -3 and 2.5×10 -3), the relative total uncertainty ranges from 25 and 50%. For stable or unstable conditions with a stability parameter |ζ|>1 on average, the total uncertainty in the neutral drag coefficient exceeds the neutral drag coefficient value itself, while for |ζ|<1 the total uncertainty is around 25% of the drag coefficient. For closer-to-neutral conditions, this uncertainty is dominated by measurement uncertainties in surface turbulent momentum fluxes which should therefore be the target of efforts in uncertainty reduction. We also propose an objective data-screening procedure for field data, which consists of retaining data for which the relative error on neutral drag coefficient does not exceed a given threshold. This method, in addition to the commonly used flux quality control procedure, allows for a reduction of the drag coefficient dispersion compared to other data-screening methods, which we take as an indication of better dataset quality

Experiment and modelling of stable stmopheric soundary layer on complex terrain : the turbulent katabatic wind process

La couche-limite atmosphérique turbulente stable, particulièrement en zone de relief, n'est pas totalement comprise. Elle est, donc, mal représentée par les modèles atmosphériques. En présence de pente et d'un refroidissement du sol, l'augmentation locale de masse volumique génère un écoulement catabatique. En région de montagne, le maximum de vent est généralement enregistré à une hauteur (z_j) de 1-10 m. Le jet de paroi engendre un changement de signe du flux de qdm ainsi qu'une variabilité du flux de chaleur sensible proche du sol. Ces variabilités de flux turbulents contredisent l'applicabilité de la théorie des similitudes de Monin-Obukhov (TSMO), pourtant utilisée de manière universelle dans les modèles atmosphériques. Si la TSMO est discutable pour les cas très stables, c'est en présence de pente qu'elle n'est naturellement plus valide puisqu'elle néglige le couplage entre les équations de vitesse et de température. Il est donc impossible de représenter correctement un écoulement catabatique (z_j O(1m)) par un modèle utilisant la TSMO et avec une résolution verticale de l'ordre de la hauteur du maximum de vent. L'objectif du travail de thèse est d'apporter une contribution dans la compréhension et la modélisation de ces écoulements.Afin de compléter les observations peu nombreuses, une campagne de mesure a été réalisée sur une pente raide (20-40 deg) : la pente ouest du Grand Colon (chaîne de Belledonne, Alpes). Les analyses spectrales témoignent de la sensibilité de l'écoulement local aux perturbations externes, même faibles. Les caractéristiques turbulentes classiques sont observées à haute-fréquence alors que des comportements moins standards sont observés aux fréquences intermédiaires ou basses et expliqués par la présence de perturbations turbulentes d'énergie du même ordre de grandeur que l'injection locale. Les cospectres montrent un comportement propre aux écoulements catabatiques: recouvrement progressif selon z des corrélations croisées 0. La TSMO est mise en défaut pour l'écoulement observé et une solution alternative est utilisée pour estimer les flux en surface, permettant une bonne description de la vitesse de frottement.Le modèle 1D de surface de ISBA (Météo-France) est modifié pour répondre à la modélisation des écoulements catabatiques. Dans un premier temps, le modèle est validé sur un cas standard: en comparaison avec un modèle de Prandtl adapté. Dans un second temps, les données in-situ sont modélisées, d'abord en fournissant des profils de diffusivités effectives puis en utilisant un modèle modifié de turbulence d'ordre 1.5. Les modélisations 1D représentent correctement les champs moyens de vitesse et température mais montrent cependant des comportements trop diffusifs. Le modèle de longueur de mélange est principalement remis en cause, y compris en utilisant des paramétrisations adaptées.Des simulations LES 3D réalistes (Meso-NH, Météo-France) sont effectuées à haute résolution pour représenter le cas d'étude. Ces modélisations représentent finement les variabilités spatiales de l'écoulement catabatique. Cependant, des biais sont engendrés principalement par l'utilisation de la TSMO en condition aux limites de surface. Malgré la forte résolution spatiale, l'utilisation de la TSMO repousse à seulement z=2 m la perception des termes sources de l'écoulement catabatique par le modèle, alors que la source de l'écoulement atteint son maximum précisément en surface. Les modèles analytiques d'écoulement catabatique (de type Prandtl, qui pourraient aisément être intégrés en conditions aux limites) nécessitent de connaître "a priori" les profils de diffusivité. Ceci implique l'utilisation d'un modèle de turbulence. Le couplage du modèle 1D de surface (précédemment modifié et validé "off-line") est donc proposé pour répondre au manque de description de la physique par les CaL classiques de surface. Le travail préliminaire du couplage est présenté et des solutions sont proposées en perspective.The stable atmospheric boundary layer, particularly in complex terrain, is not yet fully understood and it is thus still inadequately modelled. A surface cooling of a sloping terrain generates katabatic wind due to local density increase. This flow behaves as a wall-bounded turbulent jet, often simply modelled by a local balance between the buoyancy force and the turbulent friction. In mountainous regions, the wind maximum is typically observed at a height (z_j) of 1-10 m above the ground. The wall-bounded jet is responsible for a momentum-flux sign change and a heat-flux variability close to the ground. Those turbulent-flux variabilities are fully conflicting with the aplicability of the Monin-Obukhov similarity theory (MOST), which is nevertheless universally used in the atmospheric models to provide the surface boundary condition. If the MOST is already questionable for the very stable cases, it is obviously not valide over sloping surfaces because it neglectes the coupling of the wind and temperature equations, which constitutes the katabatic source. Hence, it is not possible to adequatly model a katabatic flow (z_j O(1m)) using the MOST, especially with a vertical resolution of the order of magnitude of z_j. The aim of the this PhD work is thus to improve the current understanding and modelling capacity of the katabatic winds.Since data sets of turbulent-katabatic-flow measurements are still scarce, a new field campain was carried out on a steep slope (20-40 deg): the west face of the Grand-Colon mountain (Belledonne ridge, French Alps). The experimental setup was mainly composed of a 6m mast with four sonic-anemometer levels (1, 2, 4 and 6m) to measure the turbulence on both sides of the katabatic jet. The spectral analysis shows the hight sensitivity of the local flow to external perturbations, even when these are weak. The hight-frequency subrange shows a classical behaviour (energy-injection frequency, inertial subrange), but the spectra of the intermediate and low-frequency subranges are less typical: turbulent perturbations with an energy of the order of magnitude of the local injection are present. A specific cospectra behaviour of the katabatic flows is shown: negative and positive cross-correlations overlap gradually, increasing z. The MOST fails in representing the observed flow and a surface-flux alternative estimation is succesfully used to describe the friction vellocity.The 1D surface model of ISBA (Météo-France) is modified to model katabatic flows. The model is firstly validated with a standard calibrated Prandtl model (with variable eddy difusivity). Secondly, the field data are modelled both with a prescribed effective diffusivity (from data) and using the 1.5-order turbulence scheme. The mean velocity and temperature fields are well reproduced, but it appears that the model is over-diffusive (which generates excessive fluxes), even when an adapted mixing-length is used.Realistic 3D LES simulations (Meso-NH, Météo-France) are computed with high resolution to model the field data. Spatial flow variabilities over sloping terrain are finely represented, but are biased, mainly due to the using of MOST for the surface boundary counditions. The using of MOST shifts the start of the katabatic source detection by the atmospheric model to a height of 2 m, while the katabatic source reaches its maximum at the surface. Analytical katabatic models (of the Prandtl type, which could be easily used to feed surface boundary counditions) need an "apriori" definition of the eddy and heat diffusivities. Currently, the general definition of these diffusivities is only possible by the use of turbulent models that include closures. The coupling of the previously-presented 1D surface model (validated off-line) is suggested to overcome the lack of physics description in the classic surface boundary counditions. Preliminary work on this coupling is developed and perspective solutions are proposed

Observation et modélisation de couche limite atmosphérique stable en relief complexe : le processus turbulent d'écoulement catabatique

The stable atmospheric boundary layer, particularly in complex terrain, is not yet fully understood and it is thus still inadequately modelled. A surface cooling of a sloping terrain generates katabatic wind due to local density increase. This flow behaves as a wall-bounded turbulent jet, often simply modelled by a local balance between the buoyancy force and the turbulent friction. In mountainous regions, the wind maximum is typically observed at a height (z_j) of 1-10 m above the ground. The wall-bounded jet is responsible for a momentum-flux sign change and a heat-flux variability close to the ground. Those turbulent-flux variabilities are fully conflicting with the aplicability of the Monin-Obukhov similarity theory (MOST), which is nevertheless universally used in the atmospheric models to provide the surface boundary condition. If the MOST is already questionable for the very stable cases, it is obviously not valide over sloping surfaces because it neglectes the coupling of the wind and temperature equations, which constitutes the katabatic source. Hence, it is not possible to adequatly model a katabatic flow (z_j O(1m)) using the MOST, especially with a vertical resolution of the order of magnitude of z_j. The aim of the this PhD work is thus to improve the current understanding and modelling capacity of the katabatic winds.Since data sets of turbulent-katabatic-flow measurements are still scarce, a new field campain was carried out on a steep slope (20-40 deg): the west face of the Grand-Colon mountain (Belledonne ridge, French Alps). The experimental setup was mainly composed of a 6m mast with four sonic-anemometer levels (1, 2, 4 and 6m) to measure the turbulence on both sides of the katabatic jet. The spectral analysis shows the hight sensitivity of the local flow to external perturbations, even when these are weak. The hight-frequency subrange shows a classical behaviour (energy-injection frequency, inertial subrange), but the spectra of the intermediate and low-frequency subranges are less typical: turbulent perturbations with an energy of the order of magnitude of the local injection are present. A specific cospectra behaviour of the katabatic flows is shown: negative and positive cross-correlations overlap gradually, increasing z. The MOST fails in representing the observed flow and a surface-flux alternative estimation is succesfully used to describe the friction vellocity.The 1D surface model of ISBA (Météo-France) is modified to model katabatic flows. The model is firstly validated with a standard calibrated Prandtl model (with variable eddy difusivity). Secondly, the field data are modelled both with a prescribed effective diffusivity (from data) and using the 1.5-order turbulence scheme. The mean velocity and temperature fields are well reproduced, but it appears that the model is over-diffusive (which generates excessive fluxes), even when an adapted mixing-length is used.Realistic 3D LES simulations (Meso-NH, Météo-France) are computed with high resolution to model the field data. Spatial flow variabilities over sloping terrain are finely represented, but are biased, mainly due to the using of MOST for the surface boundary counditions. The using of MOST shifts the start of the katabatic source detection by the atmospheric model to a height of 2 m, while the katabatic source reaches its maximum at the surface. Analytical katabatic models (of the Prandtl type, which could be easily used to feed surface boundary counditions) need an "apriori" definition of the eddy and heat diffusivities. Currently, the general definition of these diffusivities is only possible by the use of turbulent models that include closures. The coupling of the previously-presented 1D surface model (validated off-line) is suggested to overcome the lack of physics description in the classic surface boundary counditions. Preliminary work on this coupling is developed and perspective solutions are proposed.La couche-limite atmosphérique turbulente stable, particulièrement en zone de relief, n'est pas totalement comprise. Elle est, donc, mal représentée par les modèles atmosphériques. En présence de pente et d'un refroidissement du sol, l'augmentation locale de masse volumique génère un écoulement catabatique. En région de montagne, le maximum de vent est généralement enregistré à une hauteur (z_j) de 1-10 m. Le jet de paroi engendre un changement de signe du flux de qdm ainsi qu'une variabilité du flux de chaleur sensible proche du sol. Ces variabilités de flux turbulents contredisent l'applicabilité de la théorie des similitudes de Monin-Obukhov (TSMO), pourtant utilisée de manière universelle dans les modèles atmosphériques. Si la TSMO est discutable pour les cas très stables, c'est en présence de pente qu'elle n'est naturellement plus valide puisqu'elle néglige le couplage entre les équations de vitesse et de température. Il est donc impossible de représenter correctement un écoulement catabatique (z_j O(1m)) par un modèle utilisant la TSMO et avec une résolution verticale de l'ordre de la hauteur du maximum de vent. L'objectif du travail de thèse est d'apporter une contribution dans la compréhension et la modélisation de ces écoulements.Afin de compléter les observations peu nombreuses, une campagne de mesure a été réalisée sur une pente raide (20-40 deg) : la pente ouest du Grand Colon (chaîne de Belledonne, Alpes). Les analyses spectrales témoignent de la sensibilité de l'écoulement local aux perturbations externes, même faibles. Les caractéristiques turbulentes classiques sont observées à haute-fréquence alors que des comportements moins standards sont observés aux fréquences intermédiaires ou basses et expliqués par la présence de perturbations turbulentes d'énergie du même ordre de grandeur que l'injection locale. Les cospectres montrent un comportement propre aux écoulements catabatiques: recouvrement progressif selon z des corrélations croisées 0. La TSMO est mise en défaut pour l'écoulement observé et une solution alternative est utilisée pour estimer les flux en surface, permettant une bonne description de la vitesse de frottement.Le modèle 1D de surface de ISBA (Météo-France) est modifié pour répondre à la modélisation des écoulements catabatiques. Dans un premier temps, le modèle est validé sur un cas standard: en comparaison avec un modèle de Prandtl adapté. Dans un second temps, les données in-situ sont modélisées, d'abord en fournissant des profils de diffusivités effectives puis en utilisant un modèle modifié de turbulence d'ordre 1.5. Les modélisations 1D représentent correctement les champs moyens de vitesse et température mais montrent cependant des comportements trop diffusifs. Le modèle de longueur de mélange est principalement remis en cause, y compris en utilisant des paramétrisations adaptées.Des simulations LES 3D réalistes (Meso-NH, Météo-France) sont effectuées à haute résolution pour représenter le cas d'étude. Ces modélisations représentent finement les variabilités spatiales de l'écoulement catabatique. Cependant, des biais sont engendrés principalement par l'utilisation de la TSMO en condition aux limites de surface. Malgré la forte résolution spatiale, l'utilisation de la TSMO repousse à seulement z=2 m la perception des termes sources de l'écoulement catabatique par le modèle, alors que la source de l'écoulement atteint son maximum précisément en surface. Les modèles analytiques d'écoulement catabatique (de type Prandtl, qui pourraient aisément être intégrés en conditions aux limites) nécessitent de connaître "a priori" les profils de diffusivité. Ceci implique l'utilisation d'un modèle de turbulence. Le couplage du modèle 1D de surface (précédemment modifié et validé "off-line") est donc proposé pour répondre au manque de description de la physique par les CaL classiques de surface. Le travail préliminaire du couplage est présenté et des solutions sont proposées en perspective

Large Eddy Simulation at high resolution for a stably stratified atmospheric boundary layer on an alpine slope overlooking Grenoble.

International audienc

Parametrizing the mesoscale enhancement of oceanic surface turbulent fluxes: A physical–statistical approach

International audienceThe mesoscale enhancement of surface turbulent fluxes at the air–sea interface is driven by the mesoscale surface wind-speed variability, especially the gustiness velocity and the mesoscale wind-speed standard variation. This study proposes a parametrization of these two variables. A large dataset based on the operational 2.5-km AROME convection-permitting model is used in a coarse-graining framework, to quantify various quantities that are subgrid at the scale of a 100-km resolution global circulation model grid cell. This provides a learning dataset to help build the parametrization. The analysis of two case studies of intense wind-speed mesoscale variability, combined with a literature review, provides a physically based set of 12 potential predictors, accounting for the convection activity and the large-scale dynamics. The least absolute shrinkage and selection operator then frames a penalized multivariate linear regression approach to identify the most relevant predictors objectively. Five predictors are selected for predicting the gustiness velocity: the updraft mass flux at the lifting condensation level, the density-current spreading velocity, the large-scale horizontal shear and divergence, and the large-scale wind speed. The parametrization of the mesoscale wind-speed standard deviation requires an additional predictor, namely the cold-pool object aggregation index. The proposed parametrization performs significantly better than the previously published parametrizations and is able to capture 80, 99, and 93% of the mesoscale enhancement of the momentum, sensible heat, and latent heat fluxes, respectively. From the perspective of a global circulation model implementation, in which some predictors may be unavailable, simpler versions of the parametrization, that is, involving fewer predictors, are also discussed

Large-Eddy Simulation of a Katabatic Jet along a Convexly Curved Slope. Part I: Statistical Results

International audienc

Turbulence spectral-characteristics of pure katabatic flow over a steep Alpine slope

International audienc

Role of the Görtler instability on the turbulent mixing in a katabatic flow along a curved slope for stably stratified atmospheric boundary layer.

International audienc

Turbulence spectral-characteristics of pure katabatic flow over a steep Alpine slope

International audienc

Turbulence anisotropy in katabatic flows. From Antarctica plateau to alpine slopes.

International audienc