Representations of the atmospheric boundary layer in global climate models

Last updated

Representations of the atmospheric boundary layer in global climate models play a role in simulations of past, present, and future climates. Representing the atmospheric boundary layer (ABL) within global climate models (GCMs) are difficult due to differences in surface type, scale mismatch between physical processes affecting the ABL and scales at which GCMs are run, and difficulties in measuring different physical processes within the ABL. Various parameterization techniques described below attempt to address the difficulty in ABL representations within GCMs.

Contents

What is the ABL?

Idealized ABL structure Atmospheric boundary layer.svg
Idealized ABL structure

The ABL is the lowest part of the Earth's troposphere, loosely about the altitude zone 0 km to 1.5 km. The ABL is the only part of the troposphere directly affected by daily cycled contact with the Earth's surface, so the ABL is directly affected by forcings originating at the surface. [1] Such forcings include: heat flux, moisture flux, convection, friction, pollutant emission, and topographically modified flow. Response times to these forcings typically are an hour or less. [1]

Why is modeling of the ABL so difficult?

Surface forcings must be accounted for in GCMs in order to have accurate simulations of the Earth's climate. Unfortunately, difficulty in simulating these forcings arise for a number of reasons. First, the Earth's surface is not uniform. It consists (broadly) of land, water, and ice, and each surface interacts differently with the atmosphere. Secondly, the time and spatial scales for forcings between the Earth's surface and atmosphere occur on scales much less than the scales at which GCMs are run. Third, many of these processes are difficult to measure directly. So, these forcings need to be parameterized as a result. [2]

In general, the Earth is divided into a 3-D grid. Prognostic equations for each relevant physical process are run for each grid point. From there, grid point values for each variable are interpolated into each grid cell, which can then be analyzed. [2]

Schematic describing GCMs AtmosphericModelSchematic.png
Schematic describing GCMs

Grid resolution for GCMs varies considerably, on the order of 1 to 5 degrees (approximately 110 to 550 km for latitude, up to 110 km for longitude) in the horizontal, and 10 levels in the vertical. As grid resolution is made finer, the amount of computational time needed to run the model increases exponentially, because there are many more grid points comprising the model. Furthermore, even with the highest resolution of 1 degree, the spatial scales of many of the physical processes incorporated in the model are still far smaller than the model resolution. [2]

GCM cloud parameterizations

Clouds and convection often arise from small-scale processes that occur within the ABL. Additionally, clouds and convection help tie together the ABL with the free atmosphere, as convection helps grow the ABL. Furthermore, when the environment is sufficiently unstable, convection may help wash away the temperature inversion that caps the ABL. Also, "the convective motions associated with clouds produce important fluxes of mass, momentum, heat, and moisture". [2] The scales at which these fluxes are set come about are usually much smaller than GCM grids. However, these fluxes are often greater than those of the synoptic flow. Parameterizations of clouds and convection aim to address the scaling differences between GCM grids and cloud/convective scales. [2]

GCM cloud parameterizations account for at least two cloud types: convective clouds and large-scale supersaturation clouds. "Large-scale supersaturation clouds occur when the relative humidity in a grid box at some model level exceeds a critical value". [2] One way of accounting for large-scale supersaturation clouds is by setting the critical relative humidity value to 80%, assigning grid boxes with relative humidity values at or above 80% as cloud covered. Another way to account for large-scale supersaturation clouds is to benchmark clouds via temperature variability, in that wherever the temperature "causes the relative humidity to reach 100% is cloud covered". [2]

Moist adiabatic adjustment scheme

Convective clouds follow one of three general parameterization schemes. The first scheme is moist adiabatic adjustment. The main advantage of this method is that it is simple, in that if the moist adiabatic lapse rate is exceeded, moisture and heat within the vertical layer are adjusted so that air within the layer is saturated. [2] [3] In turn the lapse rate is then adjusted to the moist adiabatic lapse rate, conserving energy, raining out excess moisture while not transporting momentum. The disadvantage of this scheme is that it over-estimates convection, as constrains an entire grid box to act in a convective manner, unrealistic in the real atmosphere. [2]

Kuo scheme

The second scheme is the Kuo parameterization, and is more complex than the moist adiabatic adjustment scheme. [2] [3] This scheme is concerned with large-scale moisture convergence as the moisture source for convection. A weakness of this scheme is convective heating, in that this parameterization scheme assumes that convection heats by mixing cloud and environmental air, instead of heating produced by subsidence between clouds. [2]

Arakawa-Schubert scheme

The Arakawa-Schubert scheme is the third scheme, and is the most comprehensive. [2] [3] This scheme models cumulus cloud/environment interactions, entrainment and detrainment of air/moisture, cloud downdrafts, and subsidence outside of clouds. Additionally, this scheme assumes quasi-equilibrium, with cloud dissipation at a "rate sufficient to keep the atmosphere near equilibrium in the face of large-scale destabilization." [2]

Turbulence closure

Most atmospheric turbulence occurs within the ABL, whereas the free atmosphere is largely non-turbulent. So, in order to have a proper accounting for motions within the ABL, GCMs must properly account for turbulence within the ABL. To do this, GCMs must have some way of achieving turbulence closure.

Turbulent processes can be roughly combined into three categories: heat flux, moisture flux, and momentum flux. [2] When calculating these fluxes, one ends up with more unknowns than equations, which means those equations cannot be solved directly. In order to calculate turbulent fluxes and close the equations, one must make parameterizations for the higher order terms. The sections below describe methods for parameterizing turbulent fluxes and turbulence closure.

Parameterizations

Local and non-local closure

There are two methods of parameterizing turbulent fluxes. The first is local closure. Local closure ties the unknown turbulent quantity at a specific point in space to values and gradients of known quantities at the same point. Additionally, local closure likens turbulent transport to molecular diffusion, and is usually first or second order. [1]

The second method of parameterizing turbulent fluxes is nonlocal closure. Turbulence does not just depend on local values and gradients due to the superposition of many individual eddies. Unlike local closure, nonlocal closure links unknown turbulent quantities to known quantities at many points in space. [1]

K-theory

K-theory (eddy diffusivity/viscosity theory) is a form of local closure, and is the main first order closure scheme within the surface layer. K-theory follows a similar concept as molecular viscosity, in that the turbulent flux of a quantity is proportional to its spatial gradient, with K as the eddy viscosity/diffusivity. K-theory is powerful because the flux of quantities can be directly tied into the gradient of the mean quantities with height, multiplied by a value K. The idea behind K-theory is that turbulence eliminates gradients, creating a net flow of quantities down their gradients. [1] So K is positive.

Turbulent fluxGradient

The eddy diffusivity/viscosity equations are shown in the table above. Km represents eddy viscosity, while KH and KW represent eddy diffusivity of heat and moisture, respectively. τx and τy represent Reynolds stress (momentum flux) in the x and y directions, Hv represents the turbulent heat flux, and E represents the turbulent moisture flux. ρ is the density of air, u', v', and w' are velocity perturbations, and θ' and q' are potential temperature and moisture perturbations. Other key points regarding K: [1]

Diagram of eddy diffusivity Diagram of eddy diffusivity.jpg
Diagram of eddy diffusivity

K-theory is best applied in statically neutral environments, with predominantly mechanically generated turbulence, though it has been applied to statically stable environments. [1] The figure to the right demonstrates how temperature flux works in a stable environment. If a parcel is displaced upward in a stable environment, w' is positive 0, and θ' is greater than the mean θ. So w'θ' is greater than zero. Conversely, when a parcel is displaced downward, w' is negative, and θ' is negative.

Depending on the vertical resolution of the model, K-theory may be used to account for the various fluxes within the ABL. [2] [3] Additionally, K-theory can be applied in unstable conditions within the mixed layer if a counter-gradient term is applied, which accounts for buoyancy-driven turbulence.

Bulk aerodynamic formulae

Bulk aerodynamic formulae Bulk aerodynamic formulas.jpg
Bulk aerodynamic formulae

As an alternative to K-theory, ABL fluxes can be accounted for through the use of bulk aerodynamic formulae. [2]

The bulk aerodynamic formulae use Cd, Ch, and Cq, where Cd is the drag coefficient, Ch is the heat exchange (transfer) coefficient, and Cq is the moisture exchange (transfer) coefficient. [3] Each of these coefficients are functions of the known variables, such as the average wind speed at 10m, and average potential temperature and moisture at the surface and 10m. Additionally, each of these coefficients varies based on stability. From the known variables the exchange coefficients can be easily calculated, and thus, the fluxes can be easily calculated as well. The image below shows the equations used for calculating the exchange coefficients and fluxes, with u* as the frictional velocity. [1]

Related Research Articles

<span class="mw-page-title-main">Fluid dynamics</span> Aspects of fluid mechanics involving flow

In physics, physical chemistry and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluids—liquids and gases. It has several subdisciplines, including aerodynamics and hydrodynamics. Fluid dynamics has a wide range of applications, including calculating forces and moments on aircraft, determining the mass flow rate of petroleum through pipelines, predicting weather patterns, understanding nebulae in interstellar space and modelling fission weapon detonation.

In fluid dynamics, turbulence or turbulent flow is fluid motion characterized by chaotic changes in pressure and flow velocity. It is in contrast to a laminar flow, which occurs when a fluid flows in parallel layers, with no disruption between those layers.

<span class="mw-page-title-main">General circulation model</span> Type of climate model

A general circulation model (GCM) is a type of climate model. It employs a mathematical model of the general circulation of a planetary atmosphere or ocean. It uses the Navier–Stokes equations on a rotating sphere with thermodynamic terms for various energy sources. These equations are the basis for computer programs used to simulate the Earth's atmosphere or oceans. Atmospheric and oceanic GCMs are key components along with sea ice and land-surface components.

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

<span class="mw-page-title-main">Surface layer</span> Layer of a turbulent fluid affected by interaction with a surface

The surface layer is the layer of a turbulent fluid most affected by interaction with a solid surface or the surface separating a gas and a liquid where the characteristics of the turbulence depend on distance from the interface. Surface layers are characterized by large normal gradients of tangential velocity and large concentration gradients of any substances transported to or from the interface.

<span class="mw-page-title-main">Cloud physics</span> Study of the physical processes in atmospheric clouds

Cloud physics is the study of the physical processes that lead to the formation, growth and precipitation of atmospheric clouds. These aerosols are found in the troposphere, stratosphere, and mesosphere, which collectively make up the greatest part of the homosphere. Clouds consist of microscopic droplets of liquid water, tiny crystals of ice, or both, along with microscopic particles of dust, smoke, or other matter, known as condensation nuclei. Cloud droplets initially form by the condensation of water vapor onto condensation nuclei when the supersaturation of air exceeds a critical value according to Köhler theory. Cloud condensation nuclei are necessary for cloud droplets formation because of the Kelvin effect, which describes the change in saturation vapor pressure due to a curved surface. At small radii, the amount of supersaturation needed for condensation to occur is so large, that it does not happen naturally. Raoult's law describes how the vapor pressure is dependent on the amount of solute in a solution. At high concentrations, when the cloud droplets are small, the supersaturation required is smaller than without the presence of a nucleus.

<span class="mw-page-title-main">Planetary boundary layer</span> Lowest part of the atmosphere directly influenced by contact with the planetary surface

In meteorology, the planetary boundary layer (PBL), also known as the atmospheric boundary layer (ABL) or peplosphere, is the lowest part of the atmosphere and its behaviour is directly influenced by its contact with a planetary surface. On Earth it usually responds to changes in surface radiative forcing in an hour or less. In this layer physical quantities such as flow velocity, temperature, and moisture display rapid fluctuations (turbulence) and vertical mixing is strong. Above the PBL is the "free atmosphere", where the wind is approximately geostrophic, while within the PBL the wind is affected by surface drag and turns across the isobars.

<span class="mw-page-title-main">Large eddy simulation</span>

Large eddy simulation (LES) is a mathematical model for turbulence used in computational fluid dynamics. It was initially proposed in 1963 by Joseph Smagorinsky to simulate atmospheric air currents, and first explored by Deardorff (1970). LES is currently applied in a wide variety of engineering applications, including combustion, acoustics, and simulations of the atmospheric boundary layer.

Parameterization in a weather or climate model is a method of replacing processes that are too small-scale or complex to be physically represented in the model by a simplified process. This can be contrasted with other processes—e.g., large-scale flow of the atmosphere—that are explicitly resolved within the models. Associated with these parameterizations are various parameters used in the simplified processes. Examples include the descent rate of raindrops, convective clouds, simplifications of the atmospheric radiative transfer on the basis of atmospheric radiative transfer codes, and cloud microphysics. Radiative parameterizations are important to both atmospheric and oceanic modeling alike. Atmospheric emissions from different sources within individual grid boxes also need to be parameterized to determine their impact on air quality.

<span class="mw-page-title-main">Turbulence modeling</span> Use of mathematical models to simulate turbulent flow

In fluid dynamics, turbulence modeling is the construction and use of a mathematical model to predict the effects of turbulence. Turbulent flows are commonplace in most real-life scenarios, including the flow of blood through the cardiovascular system, the airflow over an aircraft wing, the re-entry of space vehicles, besides others. In spite of decades of research, there is no analytical theory to predict the evolution of these turbulent flows. The equations governing turbulent flows can only be solved directly for simple cases of flow. For most real-life turbulent flows, CFD simulations use turbulent models to predict the evolution of turbulence. These turbulence models are simplified constitutive equations that predict the statistical evolution of turbulent flows.

In fluid dynamics, turbulence kinetic energy (TKE) is the mean kinetic energy per unit mass associated with eddies in turbulent flow. Physically, the turbulence kinetic energy is characterised by measured root-mean-square (RMS) velocity fluctuations. In the Reynolds-averaged Navier Stokes equations, the turbulence kinetic energy can be calculated based on the closure method, i.e. a turbulence model.

<span class="mw-page-title-main">Entrainment (meteorology)</span>

Entrainment is a phenomenon of the atmosphere which occurs when a turbulent flow captures a non-turbulent flow. It is typically used to refer to the capture of a wind flow of high moisture content, or in the case of tropical cyclones, the capture of drier air.

<span class="mw-page-title-main">Eddy diffusion</span> Mixing of fluids due to eddy currents

In fluid dynamics, eddy diffusion, eddy dispersion, or turbulent diffusion is a process by which fluid substances mix together due to eddy motion. These eddies can vary widely in size, from subtropical ocean gyres down to the small Kolmogorov microscales, and occur as a result of turbulence. The theory of eddy diffusion was first developed by Sir Geoffrey Ingram Taylor.

The MEMO model is a Eulerian non-hydrostatic prognostic mesoscale model for wind-flow simulation. It was developed by the Aristotle University of Thessaloniki in collaboration with the Universität Karlsruhe. The MEMO Model together with the photochemical dispersion model MARS are the two core models of the European zooming model (EZM). This model belongs to the family of models designed for describing atmospheric transport phenomena in the local-to-regional scale, frequently referred to as mesoscale air pollution models.

<span class="mw-page-title-main">Horizontal convective rolls</span>

Horizontal convective rolls, also known as horizontal roll vortices or cloud streets, are long rolls of counter-rotating air that are oriented approximately parallel to the ground in the planetary boundary layer. Although horizontal convective rolls, also known as cloud streets, have been clearly seen in satellite photographs for the last 30 years, their development is poorly understood, due to a lack of observational data. From the ground, they appear as rows of cumulus or cumulus-type clouds aligned parallel to the low-level wind. Research has shown these eddies to be significant to the vertical transport of momentum, heat, moisture, and air pollutants within the boundary layer. Cloud streets are usually more or less straight; rarely, cloud streets assume paisley patterns when the wind driving the clouds encounters an obstacle. Those cloud formations are known as von Kármán vortex streets.

Tropical convective clouds play an important part in the Earth's climate system. Convection and release of latent heat transports energy from the surface into the upper atmosphere. Clouds have a higher albedo than the underlying ocean, which causes more incoming solar radiation to be reflected back to space. Since the tops of tropical systems are much cooler than the surface of the Earth, the presence of high convective clouds cools the climate system.

In engineering, physics, and chemistry, the study of transport phenomena concerns the exchange of mass, energy, charge, momentum and angular momentum between observed and studied systems. While it draws from fields as diverse as continuum mechanics and thermodynamics, it places a heavy emphasis on the commonalities between the topics covered. Mass, momentum, and heat transport all share a very similar mathematical framework, and the parallels between them are exploited in the study of transport phenomena to draw deep mathematical connections that often provide very useful tools in the analysis of one field that are directly derived from the others.

Ocean general circulation models (OGCMs) are a particular kind of general circulation model to describe physical and thermodynamical processes in oceans. The oceanic general circulation is defined as the horizontal space scale and time scale larger than mesoscale. They depict oceans using a three-dimensional grid that include active thermodynamics and hence are most directly applicable to climate studies. They are the most advanced tools currently available for simulating the response of the global ocean system to increasing greenhouse gas concentrations. A hierarchy of OGCMs have been developed that include varying degrees of spatial coverage, resolution, geographical realism, process detail, etc.

Monin–Obukhov (M–O) similarity theory describes the non-dimensionalized mean flow and mean temperature in the surface layer under non-neutral conditions as a function of the dimensionless height parameter, named after Russian scientists A. S. Monin and A. M. Obukhov. Similarity theory is an empirical method that describes universal relationships between non-dimensionalized variables of fluids based on the Buckingham π theorem. Similarity theory is extensively used in boundary layer meteorology since relations in turbulent processes are not always resolvable from first principles.

The convective planetary boundary layer (CPBL), also known as the daytime planetary boundary layer, is the part of the lower troposphere most directly affected by solar heating of the earth's surface.

References

  1. 1 2 3 4 5 6 7 8 Roland B. Stull, ed. (1988). An introduction to boundary layer meteorology. Vol. 13. Springer. pp. 2, 200, 204, 208, 261–269.
  2. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Hartmann, Dennis L. (1994). Global physical climatology. Vol. 56. Academic press. pp. 257–258, 260–263.
  3. 1 2 3 4 5 Kalnay, Eugenia (2003). Atmospheric modeling, data assimilation and predictability. Cambridge University Press. pp. 130–135.