Moving Beyond Monin–Obukhov Similarity Theory in Modelling Wind-Speed Profiles in the Lower Atmospheric Boundary Layer under Stable Stratification

Monin–Obukhov similarity theory (MOST) is commonly used to model the wind-speed profile at altitudes relevant to wind-power production (e.g. 10–200 m). Though reasonably accurate for unstable to weakly stable stratification, this approach becomes less accurate under increasingly stable stratification, largely due to the constant-flux surface layer assumed by MOST becoming shallower than the altitude range of interest. Furthermore, above the surface layer, the Coriolis force has a considerable influence on the wind-speed profile (in particular in the formation of low-level jets) that cannot be modelled using similarity theory. Our goal is to compare the accuracy of alternative extrapolation models that are more physically appropriate above the surface layer. Using data from the 213-m Cabauw meteorological tower in the Netherlands between July 2007 and June 2008, it is shown that MOST is accurate only at low altitudes and low stability, and breaks down at high altitudes and high stability. Local similarity is generally more accurate than MOST across all altitudes and stabilities, though the model requires turbulent flux data at multiple altitudes that is generally impractical. In contrast, a two-layer MOST–Ekman model is found to be comparable to the other models at low stability ranges and considerably more accurate in the high stability range, while requiring only a measure of surface stability and the geostrophic wind.     

Oceanic origin of southeast tropical Atlantic biases

Most coupled general circulation models suffer from a prominent warm sea surface temperature bias in the southeast tropical Atlantic Ocean off the coast of Africa. The origin of the bias is not understood and remains highly controversial. Previous studies suggest that the origin of the bias stems from systematic errors of atmospheric models in simulating surface heat flux and coastal wind, or poorly simulated coastal upwelling. In this study, we show, using different reanalysis and observational data sets combined with a set of eddy-resolving regional ocean model simulations, that systematic errors in ocean models also make a significant contribution to the bias problem. In particular (1) the strong warm bias at the Angola-Benguela front that is maintained by the local wind and the convergence of Angola and Benguela Currents is caused by an overshooting of the Angola Current in ocean models and (2) the alongshore warm bias to the south of the front is caused by ocean model deficiencies in simulating the sharp thermocline along the Angola coast, which is linked to biases in the equatorial thermocline, and the complex circulation system within the Benguela upwelling zone.     

A comparison of the vegetation response to rainfall in the Sahel and East Africa,using normalized difference vegetation index from NOAA AVHRR

This article presents the results of a study of the relationship between rainfall and Normalized Difference Vegetation Index (NDVI) in East Africa and the Sahel. Monthly data for the years 1982 to 1985 have been analyzed. We have evaluated NDVI-rainfall relationships by vegetation type, using the major formations described by White (1983). In the article, a comparison of the differential response of vegetation growth to rainfall in the two study regions is emphasized. The most important conclusions of our research are as follows:

1. The spatial patterns of annually-integrated NDVI closely reflect mean annual rainfall.
2. There is a good relationship between rainfall variations and NDVI on seasonal and interannual time scales for areas where mean annual rainfall ranges from approximately 200 to 1200 mm.
3. In most cases, NDVI is best correlated with the rainfall total for the concurrent plus two antecedent months; the correlation is better in the Sahel than in East Africa.
4. The ratios of NDVI to rainfall are considerably higher in East Africa than in the Sahel.
5. Mean annually-integrated NDVI is linearly related to mean annual rainfall in the Sahel. In East Africa the relationship is approximately log-linear; above some threshold value of rainfall, NDVI values level off and vary minimally with rainfall.

Two possible explanations of this last conclusion are suggested: above this threshold, rainfall is no longer the limiting factor in vegetation growth and/or NDVI is not a good indicator of vegetation growth. The latter is a likely possibility since NDVI directly reflects photosynthetic activity and becomes a poor indicator of biomass (i.e., growth) as high canopy densities are reached. The NDVI-rainfall relationship for East Africa is markedly similar to the relationship between NDVI and Leaf Area Index demonstrated by Sellers (1985) and Asrar et al. (1984).     

Atmospheric teleconnection mechanisms of extratropical North Atlantic SST influence on Sahel rainfall

Extratropical North Atlantic cooling has been tied to droughts over the Sahel in both paleoclimate observations and modeling studies. This study, which uses an atmospheric general circulation model (GCM) coupled to a slab ocean model that simulates this connection, explores the hypothesis that the extratropical North Atlantic cooling causes the Sahel droughts via an atmospheric teleconnection mediated by tropospheric cooling. The drying is also produced in a regional climate model simulation of the Sahel when reductions in air temperature (and associated geopotential height and humidity changes) from the GCM simulation are imposed as the lateral boundary conditions. This latter simulation explicitly demonstrates the central role of tropospheric cooling in mediating the atmospheric teleconnection from extratropical North Atlantic cooling. Diagnostic analyses are applied to the GCM simulation to infer teleconnection mechanisms. An analysis of top of atmosphere radiative flux changes diagnosed with a radiative kernel technique shows that extratropical North Atlantic cooling is augmented by a positive low cloud feedback and advected downstream, cooling Europe and North Africa. The cooling over North Africa is further amplified by a reduced greenhouse effect from decreased atmospheric specific humidity. A moisture budget analysis shows that the direct moisture effect and monsoon weakening, both tied to the ambient cooling and resulting circulation changes, and feedbacks by vertical circulation and evaporation augment the rainfall reduction. Cooling over the Tropical North Atlantic in response to the prescribed extratropical cooling also augments the Sahel drying. Taken together, they suggest a thermodynamic pathway for the teleconnection. The teleconnection may also be applicable to understanding the North Atlantic influence on Sahel rainfall over the twentieth century.     

The eastward shift of the Walker Circulation in response to global warming and its relationship to ENSO variability

This study investigates the global warming response of the Walker Circulation and the other zonal circulation cells (represented by the zonal stream function), in CMIP3 and CMIP5 climate models. The changes in the mean state are presented as well as the changes in the modes of variability. The mean zonal circulation weakens in the multi model ensembles nearly everywhere along the equator under both the RCP4.5 and SRES A1B scenarios. Over the Pacific the Walker Circulation also shows a significant eastward shift. These changes in the mean circulation are very similar to the leading mode of interannual variability in the tropical zonal circulation cells, which is dominated by El Niño Southern Oscillation variability. During an El Niño event the circulation weakens and the rising branch over the Maritime Continent shifts to the east in comparison to neutral conditions (vice versa for a La Niña event). Two-thirds of the global warming forced trend of the Walker Circulation can be explained by a long-term trend in this interannual variability pattern, i.e. a shift towards more El Niño-like conditions in the multi-model mean under global warming. Further, interannual variability in the zonal circulation exhibits an asymmetry between El Niño and La Niña events. El Niño anomalies are located more to the east compared with La Niña anomalies. Consistent with this asymmetry we find a shift to the east of the dominant mode of variability of zonal stream function under global warming. All these results vary among the individual models, but the multi model ensembles of CMIP3 and CMIP5 show in nearly all aspects very similar results, which underline the robustness of these results. The observed data (ERA Interim reanalysis) from 1979 to 2012 shows a westward shift and strengthening of the Walker Circulation. This is opposite to what the results in the CMIP models reveal. However, 75 % of the trend of the Walker Circulation can again be explained by a shift of the dominant mode of variability, but here towards more La Niña-like conditions. Thus in both climate change projections and observations the long-term trends of the Walker Circulation seem to follow to a large part the pre-existing dominant mode of internal variability.     

Validation of boundary layer parameters of climate model REMO: estimation of leaf area index from NOAA-AVHRR data for the Baltimos region

The preparation of time- and space-dependent input surface parameters for the climate model REMO was one task of the Baltimos project “Validation of Boundary Layer Parameters and Extension of Boundary conditions of Climate Model REMO”. The leaf area index (LAI) is one of these parameters. It is used in REMO as defined value per month for each land-use class with a defined seasonal trend during the year. Since 1982 at the Institute of Meteorology of the Free University Berlin, a high-resolved AVHRR data set of the NOAA satellite has been available (1/100 degree, approximately 1?×?1 km at nadir in a geographic coordination system) (Koslowsky 1996). The vegetation periods of the years 1997 until 2001 were selected from the dataset to estimate the LAI within the Baltimos region on the basis of an algorithm by Sellers et al. (J Climate 9:706–737, 1996) and a modified United States Geological Survey (USGS) land-use classification. The calculated high-resolved NOAA LAI values were converted to the 1/6 degree grid of the REMO climate model. Then, they were compared to the fixed LAI values, which are used in the model.     

The atmospheric boundary layer structure over the open and ice-covered Baltic Sea: in situ measurements compared to simulations with the regional model REMO

The regional model REMO, which is the atmospheric component of the coupled atmosphere–ice–ocean–land climate model system BALTIMOS, is tested with respect to its ability to simulate the atmospheric boundary layer over the open and ice-covered Baltic Sea. REMO simulations are compared to ship, radiosonde, and aircraft observations taken during eight field experiments. The main results of the comparisons are: (1) The sharpness and strength of the temperature inversion are underestimated by REMO. Over open water, this is connected with an overestimation of cloud coverage and moisture content above the inversion. (2) The vertical temperature stratification in the lowest 200 m over sea ice is too stable. (3) The horizontal inhomogeneity of sea ice concentration as observed by aircraft could not be properly represented by the prescribed ice concentration in REMO; large differences in the surface heat fluxes arise especially under cold-air advection conditions. The results of the comparisons suggest a reconsideration of the parameterization of subgrid-scale vertical exchange both under unstable und stable conditions.     

Assessment of uncertainties in the response of the African monsoon precipitation to land use change simulated by a regional model

Land use and land cover (LULC) over Africa have changed substantially over the last 60 years and this change has been proposed to affect monsoon circulation and precipitation. This study examines the uncertainties of model simulated response in the African monsoon system and Sahel precipitation due to LULC change using a set of regional model simulations with different combinations of land surface and cumulus parameterization schemes. Although the magnitude of the response covers a broad range of values, most of the simulations show a decline in Sahel precipitation due to the expansion of pasture and croplands at the expense of trees and shrubs and an increase in surface air temperature. The relationship between the model responses to LULC change and the climatologists of the control simulations is also examined. Simulations that are climatologically too dry or too wet compared to observations and reanalyses have weak response to land use change because they are in moisture or energy limited regimes respectively. The ones that lie in between have stronger response to the LULC changes, showing a more significant role in land–atmosphere interactions. Much of the change in precipitation is related to changes in circulation, particularly to the response of the intensity and latitudinal position of the African Easterly Jet, which varies with the changes in meridional surface temperature gradients. The study highlights the need for measurements of the surface fluxes across the meridional cross-section of the Sahel to evaluate models and thereby allowing human impacts such as land use change on the monsoon to be projected more realistically.     

The Effect of Surface Heterogeneity on the Structure Parameters of Temperature and Specific Humidity: A Large-Eddy Simulation Case Study for the LITFASS-2003 Experiment

We conduct a high-resolution large-eddy simulation (LES) case study in order to investigate the effects of surface heterogeneity on the (local) structure parameters of potential temperature \(C_T^2\) and specific humidity \(C_q^2\) in the convective boundary layer (CBL). The kilometre-scale heterogeneous land-use distribution as observed during the LITFASS-2003 experiment was prescribed at the surface of the LES model in order to simulate a realistic CBL development from the early morning until early afternoon. The surface patches are irregularly distributed and represent different land-use types that exhibit different roughness conditions as well as near-surface fluxes of sensible and latent heat. In the analysis, particular attention is given to the Monin–Obukhov similarity theory (MOST) relationships and local free convection (LFC) scaling for structure parameters in the surface layer, relating \(C_T^2\) and \(C_q^2\) to the surface fluxes of sensible and latent heat, respectively. Moreover we study possible effects of surface heterogeneity on scintillometer measurements that are usually performed in the surface layer. The LES data show that the local structure parameters reflect the surface heterogeneity pattern up to heights of 100–200 m. The assumption of a blending height, i.e. the height above the surface where the surface heterogeneity pattern is no longer visible in the structure parameters, is studied by means of a two-dimensional correlation analysis. We show that no such blending height is found at typical heights of scintillometer measurements for the studied case. Moreover, \(C_q^2\) does not follow MOST, which is ascribed to the entrainment of dry air at the top of the boundary layer. The application of MOST and LFC scaling to elevated \(C_T^2\) data still gives reliable estimates of the surface sensible heat flux. We show, however, that this flux, derived from scintillometer data, is only representative of the footprint area of the scintillometer, whose size depends strongly on the synoptic conditions.     

Non-stationary extreme value analysis in a changing climate

This paper introduces a framework for estimating stationary and non-stationary return levels, return periods, and risks of climatic extremes using Bayesian inference. This framework is implemented in the Non-stationary Extreme Value Analysis (NEVA) software package, explicitly designed to facilitate analysis of extremes in the geosciences. In a Bayesian approach, NEVA estimates the extreme value parameters with a Differential Evolution Markov Chain (DE-MC) approach for global optimization over the parameter space. NEVA includes posterior probability intervals (uncertainty bounds) of estimated return levels through Bayesian inference, with its inherent advantages in uncertainty quantification. The software presents the results of non-stationary extreme value analysis using various exceedance probability methods. We evaluate both stationary and non-stationary components of the package for a case study consisting of annual temperature maxima for a gridded global temperature dataset. The results show that NEVA can reliably describe extremes and their return levels.     

Projecting date palm distribution in Iran under climate change using topography,physicochemical soil properties,soil taxonomy,land use,and climate data

This study set out to model potential date palm distribution under current and future climate scenarios using an emission scenario, in conjunction with two different global climate models (GCMs): CSIRO-Mk3.0 (CS), and MIROC-H (MR), and to refine results based on suitability under four nonclimatic parameters. Areas containing suitable physicochemical soil properties and suitable soil taxonomy, together with land slopes of less than 10° and suitable land uses for date palm (Phoenix dactylifera) were selected as appropriate refining tools to ensure the CLIMEX results were accurate and robust. Results showed that large regions of Iran are projected as likely to become climatically suitable for date palm cultivation based on the projected scenarios for the years 2030, 2050, 2070, and 2100. The study also showed CLIMEX outputs merit refinement by nonclimatic parameters and that the incremental introduction of each additional parameter decreased the disagreement between GCMs. Furthermore, the study indicated that the least amount of disagreement in terms of areas conducive to date palm cultivation resulted from CS and MR GCMs when the locations of suitable physicochemical soil properties and soil taxonomy were used as refinement tools.     

TEST OF A TROPICAL LIMITED AREA NUMERICAL PREDICTION MODEL INCLUDING EFFECT OF REAL TOPOGRAPHY

In this paper, a tropical limited numerical prediction p.e. model in sigma coordinate is developed. The predicted variables are deviations from a rest reference atmosphere. This transformation is of benefit in reducing truncation error and guarantees computational stability after the very steep real topography is introduced into the model. The numerical tests with this model show that the large-scale obstacle of the Tibetan Plateau is not negligible in forecasting tropical systems like the South China Sea typhoon.     

Exploring spatiotemporal patterns and physical controls of soil moisture at various spatial scales

Soil moisture variability of various spatial scales is analyzed based on empirical orthogonal function (EOF) method using soil moisture datasets with various spatial resolutions: 1 km eco-hydrological model simulation, 0.25° passive microwave (Advanced Microwave Scanning Radiometer for the Earth Observing System, AMSR-E) dataset, and 0.5° land surface model simulation from Climate Predictor Center (CPC). All three datasets generate EOFs that explain similar variances with those generated from in situ observations from agro-meteorological network. Using AMSR-E product and eco-hydrological model simulation, it is found that the primary spatial pattern of soil moisture obtained from watershed scale has a strong connection to topographic attributes, followed by soil texture and rainfall variability, as suggested by the correlation between the primary EOF mode (EOF1) of soil moisture and landscape attributes. However, the EOF analysis of both AMSR-E and CPC datasets at regional scale reaches the conclusion that soil texture indices, such as sand and clay content, is of higher importance to soil moisture EOF1 spatial pattern (explaining 61 % variance) than topography is. Furthermore, correlation between soil moisture EOF1 and soil property is higher in spring than in autumn, which indicates that soil water-holding and drainage capabilities are more important under dry conditions. At national scale, the combined effects of topographic feature and soil property are clearly exhibited in EOF1. The study results reveal that different emphases should be placed on accurate acquisition of landscape attributes for soil moisture estimation according to various spatial scales.     

Global and regional trends in greenhouse gas emissions from livestock

Following IPCC guidelines (IPCC 2006), we estimate greenhouse gas emissions related to livestock in 237 countries and 11 livestock categories during the period 1961–2010. We find that in 2010 emissions of methane and nitrous oxide related to livestock worldwide represented approximately 9 % of total greenhouse gas (GHG) emissions. Global GHG emissions from livestock increased by 51 % during the analyzed period, mostly due to strong growth of emissions in developing (Non-Annex I) countries (+117 %). In contrast, developed country (Annex I) emissions decreased (?23 %). Beef and dairy cattle are the largest source of livestock emissions (74 % of global livestock emissions). Since developed countries tend to have lower CO₂-equivalent GHG emissions per unit GDP and per quantity of product generated in the livestock sector, the amount of wealth generated per unit GHG emitted from the livestock sector can be increased by improving both livestock farming practices in developing countries and the overall state of economic development. Our results reveal important details of how livestock production and associated GHG emissions have occurred in time and space. Discrepancies with higher tiers, demonstrate the value of more detailed analyses, and discourage over interpretation of smaller-scale trends in the Tier 1 results, but do not undermine the value of global Tier 1 analysis.     

Cause of severe droughts in Southwest China during 1951–2010

The cause of severe droughts over the Southwest China (SWC) during the local dry season is investigated based on the station rainfall data and the National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis data during 1951–2010. The droughts are in general consistent with local anomalous descent in the middle troposphere. The diagnosis of the vertical motion (omega) equation indicates that the local descent is primarily maintained by the anomalous cold temperature advection processes. Both the advection of anomalous temperature by mean wind and the advection of mean temperature by anomalous wind contribute to maintaining the anomalous descent over the SWC region. A composite analysis shows that the circulation anomaly over SWC is induced by remote forcing from the tropical Pacific and North Atlantic Oceans. During La Niña years, enhanced heating over the Maritime Continent induces anomalous downward motion over SWC through the connection of local Hadley circulation. Adiabatic warming associated with the downward motion helps to set up and maintain the local anomalous anticyclone. Another possible route is through the North Atlantic-Asia teleconnection, in which downstream Rossby wave energy propagation plays a crucial role. A negative-phase North Atlantic Oscillation may trigger a large-scale wave train pattern that induces an anomalous anticyclone over the subtropical Asia and promotes the dry condition over SWC.     

The gas-phase reaction of ClONO2 with hydrogen chloride

The gas-phase reaction of ClONO₂ with HCl was investigated using two large-volume environmental chambers with analysis by in situ long pathlength Fourier transform infrared absorption spectroscopy. In these chambers the reaction was observed to proceed, at least in part, by heterogenous routes, and an upper limit to the rate constant for the homogeneous gas-phase reaction of geneous routes, and an upper limit to the rate constant for the homogeneous gas-phase reaction of $$k\left( {{\text{ClONO}}_{\text{2}} + {\text{HCl}}} \right) < 1.5 \times 10^{ - 19} {\text{ cm}}^{\text{3}} {\text{ molecule}}^{{\text{ - 1}}} {\text{ s}}^{{\text{ - 1}}}$$ Was derived at 298±2K. Assuming that this room-temperature upper limit to the rate constant is applicable to stratospheric temperatures, this homogeneous gas-phase reaction can be estimated to be of negligible importance as a ClONO₂ loss process in the stratosphere.     

Impact of projected climate change on hydrologic regime of the Upper Paraguay River basin

We present an assessment of climate change impacts on the hydrologic regime of the 600,000 km² Upper Paraguay River basin, located in central South America based on predictions of 20 Atmospheric/Ocean General Circulation Models (AOGCMs). We considered two climate change scenarios from the Intergovernmental Panel on Climate Change (IPCC) and two 30-years time intervals centered at 2030 and 2070. Projected temperature and precipitation anomalies estimated by the AOGCMs for the study site are spatially downscaled. Time series of projected temperature and precipitation were estimated using the delta change approach. These time series were used as input to a detailed coupled hydrologic-hydraulic model aiming to estimate projected streamflow in climate change scenarios at several control points in the basin. Results show that impacts on streamflow are highly dependent on the AOGCM used to obtain the climate predictions. Patterns of temperature increase persist over the entire year for almost all AOGCMs resulting in an increase in the evapotranspiration rate of the hydrological model. The precipitation anomalies show large dispersion, being projected as either an increase or decrease in precipitation rates. Based on these inputs, results from the coupled hydrologic-hydraulic model show nearly one half of projections as increasing river discharge, and other half as decreasing river discharge. If the mean or median of the predictions is considered, no discernible change in river discharge should be expected, despite the dispersion among results of the AOGCMs that reached +/?10 % in the short horizon and +/? 20 % in the long horizon, at several control points.     

Simulation of the variations of the temperature of droplets in freezing precipitation

Carried out is the modeling of the process of temperature variations in a droplet and possible ice coating under synoptic conditions observed in the Perm krai on December 14, 2010 that are characterized with the freezing rain fall at the temperature inversion in the atmospheric boundary layer. The phase transition on the moving water-ice interface, the curvature of the phase transition boundary, and variations of heat exchange on the boundary between the ice coating and atmosphere are taken into account. Computed is the thickness of the ice crust formed on fine droplets under conditions that provide the initiation of freezing precipitation. The evolution of the rate of the front of water crystallization on the droplet surface is defined and explained.