An analytic similarity theory for the planetary boundary layer stabilized by surface buoyancy

An analytic solution for a steady, horizontally homogeneous boundary layer with rotation, % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaadAgaaaa!38AA! \ f \] , and surface friction velocity, û_*, subjected to surface buoyancy characterized by Obukhov length L, is proposed as follows. Nondimensional variables are % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeA7a6jabg2 % da9iaadAgacaWG6bGaai4laiabeE7aOnaaBaaaleaacqGHxiIkaeqa % aOGaamyDamaaBaaaleaacqGHxiIkaeqaaOGaaiilaiqadwhagaqcai % abg2da9iabeE7aOnaaBaaaleaacqGHxiIkaeqaaOGabmyvayaajaGa % ai4laiqadwhagaqcamaaBaaaleaacqGHxiIkaeqaaOGaaiilaiqads % fagaqcaiabg2da9iqbes8a0zaajaGaai4laiaadwhadaWgaaWcbaGa % ey4fIOcabeaakiqadwhagaqcamaaBaaaleaacqGHxiIkcaGGSaaabe % aaaaa!5587! \ \zeta = fz/\eta _ * u_ * ,\hat u = \eta _ * \hat U/\hat u_ * ,\hat T = \hat \tau /u_ * \hat u_{ * ,} \] , where carets denote complex (vector) quantities; Û is the mean velocity; % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiqbes8a0zaaja% aaaa!3994!\\hat \tau \]is the kinematic turbulent stress; and % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeE7aOnaaBa % aaleaacqGHxiIkaeqaaOGaeyypa0JaaiikaiaaigdacqGHRaWkcqaH % +oaEdaWgaaWcbaGaamOtaaqabaGccaWG1bWaaSbaaSqaaiabgEHiQa % qabaGccaGGVaGaamOuamaaBaaaleaacaWGJbaabeaakiaadAgacaWG % mbGaaiykamaaCaaaleqabaGaeyOeI0IaaGymaiaac+cacaaIYaaaaa % aa!4B1F! \ \eta _ * = (1 + \xi _N u_ * /R_c fL)^{ - 1/2} \]is a stability parameter. The constant % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabe67a4naaBa % aaleaacaWGobaabeaaaaa!3A81! \\xi _N \] is the ratio of the maximum mixing length(% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaamaaBaaaleaaca% WGTbaabeaaaaa!38DD!\_m \]) to the PBL depth, % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaadwhadaWgaa % WcbaGaey4fIOcabeaakiaac+cacaWGMbaaaa!3B7C! \ u_ * /f \] , for neutrally stable conditions; and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaadkfadaWgaa% WcbaGaam4yaaqabaaaaa!39AA!\R_c\](the critical flux Richardson number) is the ratio % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaadYgadaWgaa % WcbaGaamyBaaqabaGccaGGVaGaamitaaaa!3B5C! \ l_m /L \] under highly stable conditions. Profiles of stress and velocity in the ocean (<0) are given by % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaamaaxacabaGabm % yDayaajaGaeyypa0ZaaiqaaqaabeqaaiabgkHiTiaadMgacqaH0oaz % caWGLbWaaWbaaSqabeaacqaH0oazcqaH2oGEaaGccaqGGaGaaeiiai % aabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGa % aeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccaca % qGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaa % bccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaae % iiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqG % GaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabc % cacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeii % aiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGa % GaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabcca % caqGGaGaaeiiaiaabccacaqGGaGaeqOTdONaeyizImQaeyOeI0Iaeq % OVdG3aaSbaaSqaaiaad6eaaeqaaaGcbaGaeyOeI0IaamyAaiabes7a % KjaadwgadaahaaWcbeqaaiabes7aKjabe67a4naaBaaameaacaWGob % aabeaaaaGccqGHsisldaWcaaqaaiabeE7aOnaaBaaaleaacaGGQaaa % beaaaOqaaiaadUgaaaWaamWaaeaaciGGSbGaaiOBamaalaaabaWaaq % WaaeaacqaH2oGEaiaawEa7caGLiWoaaeaacqaH+oaEdaWgaaWcbaGa % amOtaaqabaaaaOGaey4kaSIaaiikaiabes7aKjabgkHiTiaadggaca % GGPaGaaiikaiabeA7a6jabgUcaRiabe67a4naaBaaaleaacaWGobaa % beaakiaacMcacqGHsisldaWcaaqaaiaadggaaeaacaaIYaaaaiabes % 7aKjaacIcacqaH2oGEdaahaaWcbeqaaiaaikdaaaGccqGHsislcqaH % +oaEdaqhaaWcbaGaamOtaaqaaiaaikdaaaGccaGGPaaacaGLBbGaay % zxaaGaaeiiaiaabccacaqGGaGaaeiiaiabeA7a6naaBaaaleaacaaI % WaaabeaakiabgwMiZkabeA7a6jabg6da+iabgkHiTiabe67a4naaBa % aaleaacaWGobaabeaaaaGccaGL7baaaSqabKazbaiabaGabmivayaa % jaGaeyypa0JaamyzamaaCaaajqMaacqabeaacaWGPbGaeqiTdqMaeq % OTdOhaaaaaaaa!C5AA! \ \mathop {\hat u = \left\{ \begin{array}{l} - i\delta e^{\delta \zeta } {\rm{ }}\zeta \le - \xi _N \\ - i\delta e^{\delta \xi _N } - \frac{{\eta _* }}{k}\left {\ln \frac{{\left| \zeta \right|}}{{\xi _N }} + (\delta - a)(\zeta + \xi _N ) - \frac{a}{2}\delta \end{array} \right.}\limits^{\hat T = e^{i\delta \zeta } } \] where % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l % b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr % 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabes7aKjabg2 % da9maabmaabaGaamyAaiaac+cacaWGRbGaeqOVdG3aaSbaaSqaaiaa % d6eaaeqaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacaaIXaGaai4lai % aaikdaaaGccaGG7aGaamyyaiabg2da9iabeE7aOnaaBaaaleaacqGH % xiIkaeqaaOGaaiikaiaaigdacaGGVaGaeqOVdG3aaSbaaSqaaiaad6 % eaaeqaaOGaey4kaSIaamyDamaaBaaaleaacqGHxiIkaeqaaOGaai4l % aiaadAgacaWGmbGaamOuamaaBaaaleaacaWGJbaabeaakiaacMcaca % GGOaGaaGymaiabgkHiTiabeE7aOnaaBaaaleaacqGHxiIkaeqaaOGa % aiykaiaacUdaaaa!5CB6! \ \delta = \left( {i/k\xi _N } \right)^{1/2} ;a = \eta _ * (1/\xi _N + u_ * /fLR_c )(1 - \eta _ * ); \] and ₀ is the nondimensional surface roughness. The constants are% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaadkfadaWgaa% WcbaGaam4yaaqabaaaaa!39AA!\R_c \]= 0.2 and% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabe67a4naaBa% aaleaacaWGobaabeaaaaa!3A81!\\xi _N \]= 0.052. The solutions for the atmosphere are similar except û is the nondimensional velocity The model produces satisfactory predictions of geostrophic drag and near-surface current (wind) profiles under stable stratification.