On using mixed-layer transport parameterizations with radiometric surface temperature for computing regional scale sensible heat flux

Recent mixed-layer formulations for computing large-scale surface energy fluxes under daytime convective conditions do not require the estimation of surface-layer parameters, such as the roughness lengths for momentum and heat. This greatly simplifies approaches using operational satellite measurements of surface temperature for computing the surface energy balance at regional scales because the surface roughness parameters are not well known for many landscapes. The utility of such mixed-layer formulations is tested using data from several recent multidisciplinary field experiments (HAPEX-MOBILHY, FIFE and Monsoon 90). The results indicate that specific mixed-layer formulations adequately simulate surface sensible heat fluxes in the grassland and shrubland sites. However, use of the original values of proposed empirical coefficients for the forested site yield poor results. This is probably due to the fact that the forested site has significantly different surface geometry and associated distribution of temperature among the surface components (especially the relative importance of soil background temperatures) compared to the other sites. Therefore, the relationship between aerodynamic and radiometric surface temperature may have greatly differed between the forested site and the other locations. However, differences in aerodynamic roughness between the experimental sites were not correlated with changes required in the values of the coefficients. Instead, a two-source model which makes the distinction between aerodynamic and radiative temperature is proposed, as a means to determine which surface properties significantly affect the magnitude of the mixed-layer coefficients.     

A Simple Parametrization for the Concentration Variance Dissipation in a Lagrangian Single-Particle Model

A new quasi-analytical mixed-layer model is formulated describing the evolution of the convective atmospheric boundary layer (ABL) during cold-air outbreaks (CAO) over polar oceans downstream of the marginal sea-ice zones. The new model is superior to previous ones since it predicts not only temperature and mixed-layer height but also the height-averaged horizontal wind components. Results of the mixed-layer model are compared with dropsonde and aircraft observations carried out during several CAOs over the Fram Strait and also with results of a 3D non-hydrostatic (NH3D) model. It is shown that the mixed-layer model reproduces well the observed ABL height, temperature, low-level baroclinicity and its influence on the ABL wind speed. The mixed-layer model underestimates the observed ABL temperature only by about 10 %, most likely due to the neglect of condensation and subsidence. The comparison of the mixed-layer and NH3D model results shows good agreement with respect to wind speed including the formation of wind-speed maxima close to the ice edge. It is concluded that baroclinicity within the ABL governs the structure of the wind field while the baroclinicity above the ABL is important in reproducing the wind speed. It is shown that the baroclinicity in the ABL is strongest close to the ice edge and slowly decays further downwind. Analytical solutions demonstrate that the \(\mathrm{e}\)-folding distance of this decay is the same as for the decay of the difference between the surface temperature of open water and of the mixed-layer temperature. This distance characterizing cold-air mass transformation ranges from 450 to 850 km for high-latitude CAOs.     

Application of a mixed-layer model to bare soil surfaces

A model that couples the surface energy balance equation, a surface hydraulic resistance equation, and the force-restore soil temperature model to a mixed-layer model of the planetary boundary layer is described. The mixed layer is separated from the soil by a relatively thin surface layer and is overlain by a stable free atmosphere with prescribed profiles of potential temperature and water vapour density. The model is in reasonably good agreement with daytime micrometeorological measurements made at a wet bare site at Agassiz, British Columbia, and a desert site at Pampa de La Joya, Peru. The sensitivity of the mixed-layer model to conditions in the free atmosphere, to the parameters describing the growth of the mixed layer, and to surface roughness lengths, surface hydraulic resistance, and windspeed is examined.     

Climate variability in the southern Indian Ocean as revealed by self-organizing maps

Using a non-linear statistical analysis called “self-organizing maps”, the interannual sea surface temperature (SST) variations in the southern Indian Ocean are investigated. The SST anomalies during austral summer from 1951 to 2006 are classified into nine types with differences in the position of positive and negative SST anomaly poles. To investigate the evolution of these SST anomaly poles, heat budget analysis of mixed-layer using outputs from an ocean general circulation model is conducted. The warming of the mixed-layer by the climatological shortwave radiation is enhanced (suppressed) as a result of negative (positive) mixed-layer thickness anomaly over the positive (negative) SST anomaly pole. This contribution from shortwave radiation is most dominant in the growth of SST anomalies. In contrast to the results reported so far, the contribution from latent heat flux anomaly is not so important. The discrepancy in the analysis is explained by the modulation in the contribution from the climatological heat flux by the interannual mixed-layer depth anomaly that was neglected in the past studies.     

Summer boundary-layer height at the plateau site of Dome’C,antarctica

Measurements of the mean and turbulent structure of the planetary boundary layer using a sodar and a sonic anemometer, and radiative measurements using a radiometer, were carried out in the summer of 1999–2000 at the Antarctic plateau station of Dome C during a two-month period. At Dome C strong ground-based inversions dominate for most of the year. However, in spite of the low surface temperatures (between −50 and −20 °C), and the surface always covered by snow and ice, a regular daytime boundary-layer evolution, similar to that observed at mid-latitudes, was observed during summertime. The mixed-layer height generally reaches 200–300 m at 1300–1400 LST in high summer (late December, early January); late in the summer (end of January to February), as the solar elevation decreases, it reduces to 100–200 m. A comparison between the mixed-layer height estimated from sodar measurements and that calculated using a mixed-layer growth model shows a rather satisfactory agreement if we assign a value of 0.01–0.02 m s⁻¹ to the subsidence velocity at the top of the mixed layer, and a value of 0.003–0.004 K m⁻¹ to the potential temperature gradient above the mixed layer.     

An advective mixed-layer model for heat and moisture incorporating an analytic expression for moisture entrainment

A slab mixed-layer model with zero-order entrainment for both temperature and humidity is developed in order to examine the relative magnitude of advective and turbulence flux convergence effects. The model formulation provides an analytic function for the ratio of surface-layer to entrainment-layer humidity flux. Model results are compared with measured mixed-layer properties over one day at a coastal location. It is concluded that the model is highly successful at simulating the mixed-layer depth, and mean mixed-layer humidity. It is suggested that a first-order model may be more appropriate for the latter half of the day when the mixed-layer depth is decreasing due to the dominance of advection over vertical turbulence flux convergence.     

A ROUTINE FOR THE CALCULATION OF THE TIME-DEPENDENT HEIGHT OF THE ATMOSPHERIC BOUNDARY LAYER FROM SURFACE-LAYER PARAMETERS

A simple routine has been implemented to deduce the 24-hourevolution of the height of the atmospheric boundary layer. This uses a reduced data set of surface-layer parameters, as obtained for examplefrom surface automatic stations.The routine is based on the combination and fitting of the three alreadyexistent models for the evolution of the convectiveboundary layer, the stable boundary layer, and the surface inversionlayer.Hourly values of temperature, friction velocity and potentialtemperature scale (or sensible heat flux) in the surface layer need onlyto be supplied as input data. The lapse rate at the top of the daytime mixed layer is derived fromthe calculated surface inversion profile at sunrise, so that only a roughevaluation of the lapse rate in the free atmosphere remains to be given.The sensitivity of the mixed-layer height is expected to be not verystrong with respect to this last parameter (final part of the growth). The routine has shown satisfactory performances when compared withsodar measurements, working with only a rough average estimate of thefree atmosphere lapse rate.     

A model of internal boundary-layer development

A slab model of the boundary layer was used to study the dynamics of the internal boundary layer associated with changes in surface temperature. The usual numerical procedure involving finite differences was avoided by solving the governing equations in a Lagrangian framework. The results of the modelling study showed that mixed-layer growth was enhanced by: (a) an increase in surface roughness; (b) an increase in the surface temperature change; and (c) a decrease in the horizontal velocity. It was found that the vertical velocity induced by variations in the horizontal velocity could play an important role in controlling the expansion of the mixed layer.The second part of the study involved the formulation of a model by simplifying the governing equations. The analytical solution obtained from the model compared favourably with the results of the numerical model. Furthermore, the analytical expression for the mixed-layer height was virtually identical to that presented by Raynor et al. (1974) to fit their observational data.     

A mixed-layer model for regional evaporation

A simple mixed-layer model is developed to describe evaporation into a convective planetary boundary layer (PBL). The model comprises volume budget equations for temperature and humidity, equations to describe transport through the surface layer which is treated as part of the lower boundary, and equations to describe entrainment at the top of the PBL. The ground surface is modelled as a canopy resistance. The model was integrated with canopy resistance, surface-layer resistance and available energy, (R _n – G), input as given functions of time, and the simulated PBL was allowed to grow into an atmosphere with known temperature and humidity profiles.Two variants of the mixed-layer model were tested using data from the KNMI tower site at Cabauw in the Netherlands. These variants differed only in the formulation of entrainment: one used a formulation developed by Driedonks (1982) while the other was a simpler formulation. Simulated evaporation agreed very well with observations irrespective of which entrainment formulation was used, despite discrepancies between simulated and observed PBL height growth which were sometimes quite large for the simpler formulation. Sensitivity analysis of the model confirms that good PBL height-growth predictions are not always a prerequisite for good evaporation predictions.     

Climate Comparisons and Change Projections for the Northwest Atlantic from Six CMIP5 Models

Abstract

Key physical variables for the Northwest Atlantic (NWA) are examined in the “historical” and two future Representative Concentration Pathway (RCP) simulations of six Earth System Models (ESMs) available through Phase 5 of the Climate Model Intercomparison Project (CMIP5). The variables are air temperature, sea-ice concentration, surface and subsurface ocean temperature and salinity, and ocean mixed-layer depth. Comparison of the historical simulations with observations indicates that the models provide a good qualitative and approximate quantitative representation of many of the large-scale climatological features in the NWA (e.g., annual cycles and spatial patterns). However, the models represent the detailed structure of some important NWA ocean and ice features poorly, such that caution is needed in the use of their projected future changes. Monthly “climate change” fields between the bidecades 1986–2005 and 2046–2065 are described, using ensemble statistics of the changes across the six ESMs. The results point to warmer air temperatures everywhere, warmer surface ocean temperatures in most areas, reduced sea-ice extent and, in most areas, reduced surface salinities and mixed-layer depths. However, the magnitudes of the inter-model differences in the projected changes are comparable to those of the ensemble-mean changes in many cases, such that robust quantitative projections are generally not possible for the NWA.     

Feedbacks in the Land-Surface and Mixed-Layer Energy Budgets

A mixed-layer model of the surface energy budget and the planetary boundary layer (PBL) is developed, based on the prognostic equations for soil temperature, mixed layer potential temperature and specific humidity and the growth and abrupt collapse of the PBL. Detailed parameterizations of the longwave radiative fluxes are included. The feedbacks in the uncoupled (i.e. surface energy budget with non-responding PBL) and coupled land surface and atmospheric mixed-layer energy budgets are examined. A simplified, time continuous, version of the model, in which the specific humidity budget is the balance of evapotranspiration and dry-air entrainment, and the PBL height is given by the lifted condensation level, is shown to be in good agreement with the complete model. By forcing the simplified model with daily mean rather than periodic solar radiation, an equilibrium model state is achieved where the fluxes are in close agreement with the daily mean fluxes corresponding to the periodic forcing. The model also agrees favorably with measurements from the FIFE field experiment. Feedbacks are examined using the equilibrium model state. The uncoupled and coupled model sensitivities with respect to the minimal stomatal resistance and the atmospheric specific humidity not only differ in magnitude, but in sign as well. This results puts into question the extent to which uncoupled land-surface models that are forced with atmospheric variables may be used in sensitivity studies.     

Ocean temperature and salinity components of the Madden–Julian oscillation observed by Argo floats

New diagnostics of the Madden–Julian oscillation (MJO) cycle in ocean temperature and, for the first time, salinity are presented. The MJO composites are based on 4 years of gridded Argo float data from 2003 to 2006, and extend from the surface to 1,400 m depth in the tropical Indian and Pacific Oceans. The MJO surface salinity anomalies are consistent with precipitation minus evaporation fluxes in the Indian Ocean, and with anomalous zonal advection in the Pacific. The Argo sea surface temperature and thermocline depth anomalies are consistent with previous studies using other data sets. The near-surface density changes due to salinity are comparable to, and partially offset, those due to temperature, emphasising the importance of including salinity as well as temperature changes in mixed-layer modelling of tropical intraseasonal processes. The MJO-forced equatorial Kelvin wave that propagates along the thermocline in the Pacific extends down into the deep ocean, to at least 1,400 m. Coherent, statistically significant, MJO temperature and salinity anomalies are also present in the deep Indian Ocean.     

Variability of the Mixed-Layer Height Over Mexico City

The diurnal and seasonal variability of the mixed-layer height in urban areas has implications for ground-level air pollution and the meteorological conditions. Measurements of the backscatter of light pulses with a commercial lidar system were performed for a continuous period of almost six years between 2011 and 2016 in the southern part of Mexico City. The profiles were temporally and vertically smoothed, clouds were filtered out, and the mixed-layer height was determined with an ad hoc treatment of both the filtered and unfiltered profiles. The results are in agreement when compared with values of mixed-layer height reconstructed from, (i) radiosonde data, and (ii) surface and vertical column densities of a trace gas. The daily maxima of the mean mixed-layer height reach values \(> 3\hbox { km}\) above ground level in the months of March–April, and are clearly lower (\(< 2.7\hbox { km}\)) during the colder months from September–December. Mean daily minima are typically observed at 0700 local time (UTC ? 6h), and are lowest during the winter months with values on average below 500 m. The data presented here show an anti-correlation between high-pollution episodes and the height of the mixed layer. The growth rate of the convective mixed-layer height has a seasonal behaviour, which is characterized together with the mixed-layer-height anomalies. A clear residual layer is evident from the backscattered signals recorded in days with specific atmospheric conditions, but also from the cloud-filtered mean diurnal profiles. The occasional presence of a residual layer results in an overestimation of the reported mixed-layer height during the night and early morning hours.     

Impact of the Horizontal Heat Flux in the Mixed Layer on an Extreme Heat Event in North China:A Case Study

Extreme heat over the North China Plain is typically induced by anomalous descending flows associated with anticyclonic circulation anomalies. However, an extreme heat event that happened in the North China Plain region on 12–13 July 2015,with maximum temperature higher than 40℃ at some stations, was characterized by only a weak simultaneous appearance of an anomalous anticyclone and descending flow, suggesting that some other factor(s) may have induced this heat event. In this study, we used the forecast data produced by the Beijing Rapid Updated Cycling operational forecast system, which predicted the heat event well, to investigate the formation mechanism of this extreme heat event. We calculated the cumulative heat in the mixed-layer air column of North China to represent the change in surface air temperature. The cumulative heat was composed of sensible heat flux from the ground surface and the horizontal heat flux convergence. The results indicated that the horizontal heat flux in the mixed layer played a crucial role in the temporal and spatial distribution of high temperatures.The horizontal heat flux was found to be induced by distinct distributions of air temperatures and horizontal winds at low levels during the two days, implying a complexity of the low-level atmosphere in causing the extreme heat.     

Understanding Madden-Julian-Induced sea surface temperature variations in the North Western Australian Basin

The strongest large-scale intraseasonal (30–110 day) sea surface temperature (SST) variations in austral summer in the tropics are found in the eastern Indian Ocean between Australia and Indonesia (North-Western Australian Basin, or NWAB). TMI and Argo observations indicate that the temperature signal (std. ~0.4 °C) is most prominent within the top 20 m. This temperature signal appears as a standing oscillation with a 40–50 day timescale within the NWAB, associated with ~40 Wm^?2 net heat fluxes (primarily shortwave and latent) and ~0.02 Nm^?2 wind stress perturbations. This signal is largely related to the Madden-Julian Oscillation. A slab ocean model with climatological observed mixed-layer depth and an ocean general circulation model both accurately reproduce the observed intraseasonal SST oscillations in the NWAB. Both indicate that most of the intraseasonal SST variations in the NWAB in austral winter are related to surface heat flux forcing, and that intraseasonal SST variations are largest in austral summer because the mixed-layer is shallow (~20 m) and thus more responsive during that season. The general circulation model indicates that entrainment cooling plays little role in intraseasonal SST variations. The larger intraseasonal SST variations in the NWAB as compared to the widely-studied thermocline-ridge of the Indian Ocean region is explained by the larger convective and air-sea heat flux perturbations in the NWAB.     

CO2 climate sensitivity and snow-sea-ice albedo parameterization in an atmospheric GCM coupled to a mixed-layer ocean model

The snow-sea-ice albedo parameterization in an atmospheric general circulation model (GCM), coupled to a simple mixed-layer ocean and run with an annual cycle of solar forcing, is altered from a version of the same model described by Washington and Meehl (1984). The model with the revised formulation is run to equilibrium for 1 × CO₂ and 2 × CO₂ experiments. The 1 ×CO₂ (control) simulation produces a global mean climate about 1° warmer than the original version, and sea-ice extent is reduced. The model with the altered parameterization displays heightened sensitivity in the global means, but the geographical patterns of climate change due to increased carbon dioxide (CO₂) are qualitatively similar. The magnitude of the climate change is affected, not only in areas directly influenced by snow and ice changes but also in other regions of the globe, including the tropics where sea-surface temperature, evaporation, and precipitation over the oceans are greater. With the less-sensitive formulation, the global mean surface air temperature increase is 3.5 °C, and the increase of global mean precipitation is 7.12%. The revised formulation produces a globally averaged surface air temperature increase of 4.04 °C and a precipitation increase of 7.25%, as well as greater warming of the upper tropical troposphere. Sensitivity of surface hydrology is qualitatively similar between the two cases with the larger-magnitude changes in the revised snow and ice-albedo scheme experiment. Variability of surface air temperature in the model is comparable to observations in most areas except at high latitudes during winter. In those regions, temporal variation of the sea-ice margin and fluctuations of snow cover dependent on the snow-ice-albedo formulation contribute to larger-than-observed temperature variability. This study highlights an uncertainty associated with results from current climate GCMs that use highly parameterized snow-sea-ice albedo schemes with simple mixed-layer ocean models.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

人为热源和城市绿化对冬季山谷城市边界层结构影响的数值模拟

为了探讨人为热源和城市绿化对城市边界层的影响,本文在RAMS模式中引入周期性日变化的人为热源和人工改变下垫面状况,初步模拟分析了人为热源和城市绿化对城市边界层结构的影响.结果表明:人为热源具有明显的增温效应,增强了城市的湍流交换,白天增加了大气不稳定度,促进了混合层的发展,夜间降低了大气稳定度,减弱了城市夜间逆温.城市绿化减小了地表反照率,增加了到达地面的净辐射,模拟期间土壤冻结,增加的净辐射其中一部分以感热的形式来加热大气;绿化后地气之间的湍流交换增强,增加了大气不稳定度,减弱了白天高空逆温;本文还讨论了不同绿化布局对白天高空逆温的影响以及人为热源和城市绿化之间的非线性相互作用.     

An applied model for the height of the daytime mixed layer and the entrainment zone

A model is presented for the height of the mixed layer and the depth of the entrainment zone under near-neutral and unstable atmospheric conditions. It is based on the zero-order mixed-layer height model of Batchvarova and Gryning (1991) and the parameterization of the entrainment zone depth proposed by Gryning and Batchvarova (1994). However, most zero-order slab type models of mixed-layer height may be applied. The use of the model requires only information on those meteorological parameters that are needed in operational applications of ordinary zero-order slab type models of mixed-layer height: friction velocity, kinematic heat flux near the ground and potential temperature gradient in the free atmosphere above the entrainment zone. When information is available on the horizontal divergence of the large-scale flow field, the model also takes into account the effect of subsidence, although this is usually neglected in operational models of mixed-layer height owing to lack of data. Model performance is tested using data from the CIRCE experiment.     

Can local linear stochastic theory explain sea surface temperature and salinity variability?

 Sea surface temperature (SST) and salinity (SSS) time series from four ocean weather stations and data from an integration of the GFDL coupled ocean-atmosphere model are analyzed to test the applicability of local linear stochastic theory to the mixed-layer ocean. According to this theory, mixed-layer variability away from coasts and fronts can be explained as a ‘red noise’ response to the ‘white noise’ forcing by atmospheric disturbances. At one weather station, Papa (northeast Pacific), this stochastic theory can be applied to both salinity and temperature, explaining the relative redness of the SSS spectrum. Similar results hold for a model grid point adjacent to Papa, where the relationships between atmospheric energy and water fluxes and actual changes in SST and SSS are what is expected from local linear stochastic theory. At the other weather stations, this theory cannot adequately explain mixed-layer variability. Two oceanic processes must be taken into account: at Panulirus (near Bermuda), mososcale eddies enhance the observed variability at high frequencies. At Mike and India (North Atlantic), variations in SST and SSS advection, indicated by the coherence and equal persistence of SST and SSS anomalies, contribute to much of the low frequency variability in the model and observations. To achieve a global perspective, TOPEX altimeter data and model results are used to identify regions of the ocean where these mechanisms of variability are important. Where mesoscale eddies are as energetic as at Panulirus, indicated by the TOPEX global distribution of sea level variability, one would expect enhanced variability on short time scales. In regions exhibiting signatures of variability similar to Mike and India, variations in SST and SSS advection should dominate at low frequencies. According to the model, this mode of variability is found in the circumpolar ocean and the northern North Atlantic, where it is associated with the irregular oscillations of the model’s thermohaline circulation. Received: 11 March 1996 / Accepted: 6 September 1996     

Low-frequency variability and CO2 transient climate change. Part 3. Intermonthly and interannual variability

Components of interannual, intermonthly, and total monthly variability of lower troposphere temperature are calculated from a global coupled ocean-atmosphere general circulation model (GCM) (referred to as the coupled model), from the same atmospheric model coupled to a nondynamic mixedlayer ocean (referred to as the mixed-layer model), and from microwave sounding unit (MSU) satellite data. The coupled model produces most features of intermonthly and interannual variability compared to the MSU data, but with somewhat reduced amplitude in the extratropics and increased variability in the tropical western Pacific and tropical Atlantic. The relatively short 14-year period of record of the MSU data precludes definitive conclusions about variability in the observed system at longer time scales (e.g., decadal or longer). Different 14-year periods from the coupled model show variability on those longer time scales that were noted in Part 1 of this series. The relative contributions of intermonthly and interannual variability that make up the total monthly variability are similar between the coupled model and the MSU data, suggesting that similar mechanisms are at work in both the model and observed system. These include El Niño-Southern Oscillation (ENSO)-type interannual variability in the tropics, Madden-Julian Oscillation (MJO) type intermonthly variability in the tropics, and blocking-type intermonthly variability in the extratropics. Manifestations of all of these features have been noted in various versions of the model. Significant changes of variability noted in the coupled model with doubled carbon dioxide differ from those in our mixed-layer model and earlier studies with mixed-layer models. In particular, in our mixed-layer model intermonthly and interannual variability changes are similar with a mixture of regional increases and decreases, but with mainly decreases in the zonal mean from about 20°S to 60°N and near 60°S. In the coupled model, intermonthly and interannual changes of variability with doubled CO₂ show mostly increases of tropical interannual variability and decreases of intermonthly variability near 60°N. These changes in the tropics are related to changes in ENSO, the south Asian monsoon, and other regional hydrological regimes, while the alterations near 60°N are likely associated with changes in blocking activity. These results point to the important contribution from ENSO seen in the coupled model and the MSU data that are not present in the mixed-layer model.