Understanding the spatial differences in terrestrial water storage variations in the Tibetan Plateau from 2002 to 2016

Climate change has been driving terrestrial water storage variations in the high mountains of Asia in the recent decades. This study is based on Gravity Recovery and Climate Experiment (GRACE) data to analyse spatial and temporal variations in terrestrial water storage (TWS) across the Tibetan Plateau (TP) from April 2002 to December 2016. Regional averaged TWS anomaly has increased by 0.20 mm/month (p?<?0.01) during the 2002–2012 period, but decreased by ??0.68 mm/month (p?<?0.01) since 2012. The seasonal variations in TWS anomalies also showed a decreasing trend from May 2012 to December 2016. TWS variations in the TP also showed significant spatial differences, which were decreasing in southern TP but increasing in the Inner TP. And a declining trend was clearly evident in the seasonal variability of TWS anomalies in the south TP (about ??30 to ??55 mm/a), but increasing in the inner TP (about 10–35 mm/a). Meanwhile, this study links temperature/precipitation changes, glacial retreat and lake area expansion to explain the spatial differences in TWS. Results indicated that precipitation increases and lake area expansion drove increasing TWS in the Inner TP during the 2002–2016 period, but temperature increases and glacial retreat drove decreasing TWS in southern TP.     

Evaluating observed and projected future climate changes for the Arctic using the K?ppen-Trewartha climate classification

The ecosystems in the Arctic region are known to be very sensitive to climate changes. The accelerated warming for the past several decades has profoundly influenced the lives of the native populations and ecosystems in the Arctic. Given that the K?ppen-Trewartha (K-T) climate classification is based on reliable variations of land-surface types (especially vegetation), this study used the K-T scheme to evaluate climate changes and their impact on vegetation for the Arctic (north of 50°N) by analyzing observations as well as model simulations for the period 1900–2099. The models include 16 fully coupled global climate models from the Intergovernmental Panel on Climate Change Fourth Assessment. By the end of this century, the annual-mean surface temperature averaged over Arctic land regions is projected to increase by 3.1, 4.6 and 5.3°C under the Special Report on Emissions Scenario (SRES) B1, A1b, and A2 emission scenarios, respectively. Increasing temperature favors a northward expansion of temperate climate (i.e., Dc and Do in the K-T classification) and boreal oceanic climate (i.e., Eo) types into areas previously covered by boreal continental climate (i.e., Ec) and tundra; and tundra into areas occupied by permanent ice. The tundra region is projected to shrink by ?1.86?×?10⁶?km² (?33.0%) in B1, ?2.4?×?10⁶?km² (?42.6%) in A1b, and ?2.5?×?10⁶?km² (?44.2%) in A2 scenarios by the end of this century. The Ec climate type retreats at least 5° poleward of its present location, resulting in ?18.9, ?30.2, and ?37.1% declines in areal coverage under the B1, A1b and A2 scenarios, respectively. The temperate climate types (Dc and Do) advance and take over the area previously covered by Ec. The area covered by Dc climate expands by 4.61?×?10⁶?km² (84.6%) in B1, 6.88?×?10⁶?km² (126.4%) in A1b, and 8.16?×?10⁶?km² (149.6%) in A2 scenarios. The projected redistributions of K-T climate types also differ regionally. In northern Europe and Alaska, the warming may cause more rapid expansion of temperate climate types. Overall, the climate types in 25, 39.1, and 45% of the entire Arctic region are projected to change by the end of this century under the B1, A1b, and A2 scenarios, respectively. Because the K-T climate classification was constructed on the basis of vegetation types, and each K-T climate type is closely associated with certain prevalent vegetation species, the projected large shift in climate types suggests extensive broad-scale redistribution of prevalent ecoregions in the Arctic.     

Potential spread of recently naturalised plants in New Zealand under climate change

Climate change and biological invasions are major causes of biodiversity loss and may also have synergistic effects, such as range shifts of invaders due to changing climate. Bioclimatic models provide an important tool to assess how the threat of invasive species may change with altered temperature and precipitation regimes. In this study, potential distributions of three recently naturalised plant species in New Zealand are modelled (Archontophoenix cunninghamiana, Psidium guajava and Schefflera actinophylla), using four different general circulation models (CCCMA-CGCM3, CSIRO-Mk3.0, GFDL-CM2.0 and UKMO-HADCM3) with two emission scenarios (A2 and B1) each. Based on a maximum entropy approach, models were trained on global data using a small set of uncorrelated predictors. The models were projected to the country of interest, using climate models that had been statistically downscaled to New Zealand, in order to obtain high resolution predictions. This study provides evidence of the potential range expansion of these species, with potentially suitable habitat increasing by as much as 169 % (A. cunninghamiana; with up to 115,805 km² of suitable habitat), 133 % (P. guajava; 164,450 km²) and 208 % (S. actinophylla; 31,257 km²) by the end of the century compared to the currently suitable habitat. The results show that while predictions vary depending on the chosen climate scenario, there is remarkable consistency amongst most climate models within the same emission scenario, with overlaps in areas of predicted presence ranging between 81 % and 99.5 % (excluding CSIRO-Mk3.0). By having a better understanding of how climate change will affect distribution of invasive plants, appropriate management measures can be taken.     

Evidence of climate change impact upon glaciers’ recession within the Italian Alps

We present evidence of climate change impact upon recent changes of glaciers within Lombardy region, in Northern Italy. We illustrate the recent area evolution of a set of 249 glaciers in the area using three surface area records for 1991, 1999 and 2003. The 1999 and 2003 surface area data are processed by combining glacier limits manually digitized upon registered color orthophotos and differential GPS (DGPS) glaciers’ surveys. Glaciers’ area was 117.4?km² in 1991, 104.7?km² in 1999, and to 92.4?km² 2003, with a 21% reduction. Glaciers smaller than 1?km² accounted for 53% of the total loss in area (13.1?km² during 1991–2003). The area change rate was higher lately, with ca. 11.7 % reduction during 1999–2003. We split Alps and fore Alps of Lombardy into six mountain groups, and we separately investigate relative area variations. We use climate series from local stations within each group to assess climate change during a 30-year window (1976–2005). We focus upon temperature and snow cover depth at thaw, known to impact glaciers’ changes. We compare local year-round temperature anomalies against global ones to evidence enhanced warming within this area, and we investigate the correlation of our target climate variables against NAO. Eventually, we highlight the link between the rate of change of our climate variables to the observed scaling of area loss against glaciers’ size, showing that in rapidly warming areas glaciers’ size affects less relative melting.     

The response of lake levels and areas to climatic change

The levels and areas of lakes, particularly closed lakes (those without outlet) are known to be sensitive indicators of changes in climate, and so in principle they can play a role in monitoring current and future climatic changes. In this paper we derive solutions to the water balance equation giving the response of the level and area of closed lakes to steps, spikes and sinusoidal variations in aridity, and show that such a lake acts approximately as a simple low pass filter having a characteristic equilibrium response timeτ _e. We also review the similar response of open lakes to changes in runoff. It is shown howτ _e depends on lake geomorphology and climate, and for a number of lakes we calculateτ _e values, ranging from 1.5 to 350 yr for closed lakes, and from 2.9 days to 2 yr for open lakes. We propose methods involving Fourier analysis for inverting closed lake level or area records to obtain plots of a simple aridity index C. It is shown that for all ≈ 200 of the world's large (? 100 km²) closed lakes, satellite remote sensing of lake levels and areas is currently sensitive enough to monitor variations in C (e.g. in average basin precipitation) of order 1% to 10% on the time-scale of years to decades.     

Modelling climate change effects on the spatial distribution of mountain permafrost at three sites in northwest Canada

Spatial models of present-day mountain permafrost probability were perturbed to examine potential climate change impacts. Mean annual air temperature (MAAT) changes were simulated by adjusting elevation in the models, and cloud cover changes were examined by altering the partitioning of direct beam and diffuse radiation within the calculation for potential incoming solar radiation (PISR). The effects of changes in MAAT on equilibrium permafrost distribution proved to be more important than those due to cloud cover. Under a ?2 K scenario (approximating Little Ice Age conditions), permafrost expanded into an additional 22?C43% of the study areas as zonal boundaries descended by 155?C290 m K^???1. Under warming scenarios, permafrost probabilities progressively declined and zonal boundaries rose in elevation. A MAAT change of +5 K, caused two of the areas to become essentially permafrost-free. The absolute values of these predictions were affected up to ±10% when lapse rates were altered by ±1.5 K km^???1 but patterns and trends were maintained. A higher proportion of diffuse radiation (greater cloud cover) produced increases in permafrost extent of only 2?C4% while decreases in the diffuse radiation fraction had an equal but opposite effect. Notwithstanding the small change in overall extent, permafrost probabilities on steep south-facing slopes were significantly impacted by the altered partitioning. Combined temperature and PISR partitioning scenarios produced essentially additive results, but the impact of changes in the latter declined as MAAT increased. The modelling illustrated that mountain permafrost in the discontinuous zone is sensitive spatially to long-term climate change and identified those areas where changes may already be underway following recent atmospheric warming.     

Physical and chemical impacts of a major storm on a temperate lake: a taste of things to come?

Extreme weather can have a substantial influence on lakes and is expected to become more frequent with climate change. We explored the influence of one particular extreme event, Storm Ophelia, on the physical and chemical environment of England’s largest lake, Windermere. We found that the substantial influence of Ophelia on meteorological conditions at Windermere, in particular wind speed, resulted in a 25-fold increase (relative to the study-period average) in the wind energy flux at the lake-air interface. Following Ophelia, there was a short-lived mixing event in which the Schmidt stability decreased by over 100 Jm^?2 and the thermocline deepened by over 10 m during a 12-h period. As a result of changes to the strength of stratification, Ophelia also changed the internal seiche regime of Windermere with the dominant seiche period increasing from ~?17 h pre-storm to ~?21 h post-storm. Following Ophelia, there was an upwelling of cold and low-oxygenated waters at the southern-end of the lake. This had a substantial influence on the main outflow of Windermere, the River Leven, where dissolved oxygen concentrations decreased by ~?48%, from 9.3 to 4.8 mg L^?1, while at the mid-lake monitoring station in Windermere, it decreased by only ~?3%. This study illustrates that the response of a lake to extreme weather can cause important effects downstream, the influence of which may not be evident at the lake surface. To understand the impact of future extreme events fully, the whole lake and downstream-river system need to be studied together.     

Glacial mass balance changes in the Karakoram and Himalaya based on CMIP5 multi-model climate projections

The impact of future climate change on the glaciers in the Karakoram and Himalaya (KH) is investigated using CMIP5 multi-model temperature and precipitation projections, and a relationship between glacial accumulation-area ratio and mass balance developed for the region based on the last 30 to 40 years of observational data. We estimate that the current glacial mass balance (year 2000) for the entire KH region is -6.6?±?1 Gta^?1, which decreases about sixfold to -35?±?2 Gta^?1 by the 2080s under the high emission scenario of RCP8.5. However, under the low emission scenario of RCP2.6 the glacial mass loss only doubles to -12?±?2 Gta^?1 by the 2080s. We also find that 10.6 and 27 % of the glaciers could face ‘eventual disappearance’ by the end of the century under RCP2.6 and RCP8.5 respectively, underscoring the threat to water resources under high emission scenarios.     

Estimation of Nocturnal Area Fluxes of Carbon Cycle Gases in Saint Petersburg Suburbs

Carbon dioxide, methane, and carbon monoxide are the carbon cycle gases, the data on their emissions are needed when monitoring air pollution and developing methods for reducing anthropogenic emissions to the atmosphere and for climate forecasting. The estimates of nocturnal area fluxes for CO₂, CH₄, and CO presented for a suburb of Saint Petersburg (Peterhof) are obtained using the box model and continuous observations of concentration of these gases. The mean values of CH₄, CO₂, and CO fluxes estimated for Peterhof for 2014–2015 are 44 ± 27, 6100 ± 4000, and 90 ± 100 t/(km² year), respectively. The intensity of the CO area flux has pronounced seasonal variations characterized by the maximum of ~(160 ± 120) t/(km² year) in November—February and by the minimum of ~(30 ± 20) t/(km² year) in June-July. The analysis of the ratio of CO/CO₂ fluxes identified the main types of anthropogenic sources typical of Peterhof: motor transport, natural gas combustion, and the use of wood stoves for the heating of private low-rise buildings (in the cold season).     

Quantitative reconstruction of precipitation and runoff during MIS 5a,MIS 3a,and Holocene,arid China

Marine oxygen isotope stage 5a (MIS 5a), MIS 3a, and Holocene were highlighted periods in paleoclimate studies. Many scientists have published a great number of studies in this regard, but they paid more attention to qualitative research, and there was often a lack of quantitative data. In this paper, based on chronological evidence from a paleolake in arid China, MIS 5a, MIS 3a, and Holocene lake area, the precipitation of the drainage area and the runoff of the inflowing rivers of the lake were reconstructed with ArcGIS spatial analysis software and the improved water and energy balance model which was calibrated by modern meteorological and hydrological data in the Shiyang River drainage basin. The results showed that the paleolake areas were 1824, 1124, and 628 km² for MIS 5a, MIS 3a, and Holocene; meanwhile, the paleoprecipitation and runoff were 293.992–297.433, 271.105–274.294, and 249.431–252.373 mm and 29.103 × 10⁸–29.496 × 10⁸, 18.810 × 10⁸–18.959 × 10⁸, and 10.637 × 10⁸–10.777 × 10⁸ mm, respectively. The quantitative data can help us not only strengthen the understanding of paleoclimatic characteristics but also recognize the complexity and diversity of the climate system.     

The impacts of climate change on river flood risk at the global scale

This paper presents an assessment of the implications of climate change for global river flood risk. It is based on the estimation of flood frequency relationships at a grid resolution of 0.5?×?0.5°, using a global hydrological model with climate scenarios derived from 21 climate models, together with projections of future population. Four indicators of the flood hazard are calculated; change in the magnitude and return period of flood peaks, flood-prone population and cropland exposed to substantial change in flood frequency, and a generalised measure of regional flood risk based on combining frequency curves with generic flood damage functions. Under one climate model, emissions and socioeconomic scenario (HadCM3 and SRES A1b), in 2050 the current 100-year flood would occur at least twice as frequently across 40 % of the globe, approximately 450 million flood-prone people and 430 thousand km² of flood-prone cropland would be exposed to a doubling of flood frequency, and global flood risk would increase by approximately 187 % over the risk in 2050 in the absence of climate change. There is strong regional variability (most adverse impacts would be in Asia), and considerable variability between climate models. In 2050, the range in increased exposure across 21 climate models under SRES A1b is 31–450 million people and 59 to 430 thousand km² of cropland, and the change in risk varies between ?9 and +376 %. The paper presents impacts by region, and also presents relationships between change in global mean surface temperature and impacts on the global flood hazard. There are a number of caveats with the analysis; it is based on one global hydrological model only, the climate scenarios are constructed using pattern-scaling, and the precise impacts are sensitive to some of the assumptions in the definition and application.     

Changing inland lakes responding to climate warming in Northeastern Tibetan Plateau

The main portion of Tibetan Plateau has experienced statistically significant warming over the past 50 years, especially in cold seasons. This paper aims to identify and characterize the dynamics of inland lakes that located in the hinterland of Tibetan Plateau responding to climate change. We compared satellite imageries in late 1970s and early 1990s with recent to inventory and track changes in lakes after three decades of rising temperatures in the region. It showed warm and dry trend in climate with significant accelerated increasing annual mean temperature over the last 30 years, however, decreasing periodically annual precipitation and no obvious trend in potential evapotranspiration during the same period. Our analysis indicated widespread declines in inland lake??s abundance and area in the whole origin of the Yellow River and southeastern origin of the Yangtze River. In contrast, the western and northern origin of the Yangtze River revealed completely reverse change. The regional lake surface area decreased by 11,499 ha or 1.72% from the late 1970s to the early 1990s, and increased by 6,866 ha or 1.04% from the early 1990s to 2004. Shrinking inland lakes may become a common feature in the discontinuous permafrost regions as a consequence of warming climate and thawing permafrost. Furthermore, obvious expanding were found in continuous permafrost regions due to climate warming and glacier retreating. The results may provide information for the scientific recognition of the responding events to the climate change recorded by the inland lakes.     

Vegetation patterns in South America associated with rising CO2: uncertainties related to sea surface temperatures

Spatially precise forecasts of the impacts of climate change on the distribution of major vegetation types are essential for the implementation of effective conservation and land use policy. However, existing studies frequently omit major sources of climate variability that can significantly increase the uncertainty of projections. In this study we demonstrate how different predictions for sea surface temperature (SST) for the first half of the twenty-first century increase the uncertainty associated with forecasts of the future distribution of major ecosystems in South America. This is demonstrated through a numerical experiment using a coupled climate–vegetation model (CCM3-IBIS) for IPCC emission scenario A2 that incorporates the SST data from ten different models. The study reveals an increasing uncertainty in the ability to forecast future vegetation patterns, such that by 2050 the simulation is unable to robustly forecast the vegetation cover in an area equivalent to 28 % in South America (5?×?10⁶ km²). The future of the central and northeastern regions of Brazil is especially uncertain, with outcomes, ranging from savanna, and open shrubland to grassland. Recognizing and managing such uncertainty should be a priority for decision makers.     

Simulations of the Impact of Lakes on Local and Regional Climate Over the Tibetan Plateau

ABSTRACT

Because of the high elevation and complex topography of the Tibetan Plateau (TP), the role of lakes in the climate system over the Tibetan Plateau is not well understood. For this study, we investigated the impact of lake processes on local and regional climate using the Weather Research and Forecasting (WRF) model, which includes a one-dimensional physically based lake model. The first simulation with the WRF model was performed for the TP over the 2000–2010 period, and the second was carried out during the same period but with the lakes filled with nearby land-use types. Results with the lake simulation show that the model captures the spatial and temporal patterns of annual mean precipitation and temperature well over the TP. Through comparison of the two simulations, we found that the TP lakes mainly cool the near-surface air, inducing a decreasing sensible heat flux for the entire year. Meanwhile, stronger evaporation produced by the lakes is found in the fall. During the summer, the cooling effect of the lakes decreases precipitation in the surrounding area and generates anomalous circulation patterns. In conclusion, the TP lakes cool the near-surface atmosphere most of the time, weaken the sensible heat flux, and strengthen the latent heat flux, resulting in changes in mesoscale precipitation and regional-scale circulation.     

The radiative forcing benefits of “cool roof” construction in California: quantifying the climate impacts of building albedo modification

Exploiting surface albedo change has been proposed as a form of geoengineering to reduce the heating effect of anthropogenic increases in greenhouse gases (GHGs). Recent modeling experiments have projected significant negative radiative forcing from large-scale implementation of albedo reduction technologies (“cool” roofs and pavements). This paper complements such model studies with measurement-based calculations of the direct radiation balance impacts of replacement of conventional roofing with “cool” roof materials in California. This analysis uses, as a case study, the required changes to commercial buildings embodied in California’s building energy efficiency regulations, representing a total of 4300 ha of roof area distributed over 16 climate zones. The estimated statewide mean radiative forcing per 0.01 increase in albedo (here labeled RF01) is ?1.38 W/m². The resulting unit-roof-area mean annual radiative forcing impact of this regulation is ?44.2 W/m². This forcing is computed to counteract the positive radiative forcing of ambient atmospheric CO₂ at a rate of about 41 kg for each square meter of roof. Aggregated over the 4300 ha of cool roof estimated built in the first decade after adoption of the State regulation, this is comparable to removing about 1.76 million metric tons (MMT) of CO₂ from the atmosphere. The point radiation data used in this study also provide perspective on the spatial variability of cool roof radiative forcing in California, with individual climate zone effectiveness ranging from ?37 to ?59 W/m² of roof. These “bottom-up” calculations validate the estimates reported for published “top down” modeling, highlight the large spatial diversity of the effects of albedo change within even a limited geographical area, and offer a potential methodology for regulatory agencies to account for the climate effects of “cool” roofing in addition to its well-known energy efficiency benefits.     

Nitrous oxide emissions from gasohol, ethanol and CNG light duty vehicles

A sample of 21 light duty vehicles powered by Otto cycle engines were tested on a chassis dynamometer to measure the exhaust emissions of nitrous oxide (N₂O). The tests were performed at the Vehicle Emission Laboratory of CETESB (Environmental Company of the State of Sao Paulo) using the US-FTP-75 (Federal Test Procedure) driving cycle. The sample tested included passenger cars running on three types of fuels used in Brazil: gasohol, ethanol and CNG. The measurement of N₂O was made using two methods: Non Dispersive InfraRed (NDIR) analyzer and Fourier Transform InfraRed spectroscopy (FTIR). Measurements of regulated pollutants were also made in order to establish correlations between N₂O and NO_x. The average N₂O emission factors obtained by the NDIR method was 78?±?41?mg.km^?1 for vehicles running with gasohol, 73?±?45?mg.km^?1 for ethanol vehicles and 171?±?69?mg.km^?1 for CNG vehicles. Seventeen results using the FTIR method were also obtained. For gasohol vehicles the results showed a good agreement between the two methods, with an average emission factor of 68?±?41?mg.km^?1. The FTIR measurement results of N₂O for ethanol and CNG vehicles were much lower than those obtained by the NDIR method. The emission factors were 17?±?10?mg.km^?1 and 33?±?17?mg.km^?1, respectively, possibly because of the interference of water vapor (present at a higher concentration in the exhaust gases of these vehicles) on measurements by the NDIR method.     

Oscillatory dynamics do not mask linear trends in the timing of ice breakup for Northern Hemisphere lakes from 1855 to 2004

Our analyses partition the relative influence of progressive climate change and large-scale climate drivers that can be associated with the Quasi-Biennial Oscillation (QBO), El Niño Southern Oscillation (ENSO), North Atlantic Oscillation (NAO), solar sunspot cycle, and multi-decadal oscillations on lake ice breakup dates for thirteen Northern Hemisphere lakes. Oscillatory dynamics explain 26 % of the total variance in the time series compared with 15 % for linear trends, leaving 60 % unexplained and likely attributable, in part, to local weather. Significant oscillatory dynamics include frequencies in 2–3 year periods (9.4 % of the total variance), 3–6 year periods (8.2 %), 10–12 year periods (1.6 %) and various multidecadal periods (0.4–1.3 %). All 13 study lakes, although widely scattered in the Northern Hemisphere, had similar oscillatory dynamics and linear trends, emphasizing that global processes influence lake ice breakup locally. We illustrate that while quasi-periodic dynamics associated with large-scale climate drivers are important, they do not mask the clear evidence for progressive climate change.     

Range expansion and redefinition of a crop-raiding rodent associated with global warming and temperature increase

Since the 1970s, the crested porcupine Hystrix cristata has shown a marked range expansion in Italy. A web page has been created to collect occurrences of this species to monitor its distribution redefinition. Thus, aims of this work were (i) to identify the main predictors promoting the distribution of this large rodent in Italy and (ii) to predict its potential expansion under future climate change scenarios. A total of 1674 locations were used for this analysis, i.e., all those collected through the web page, with the exception of recently introduced populations (Sardinia, Western Liguria and Province of Varese). The current distribution of the crested porcupine covers a total of 135,177 km², as estimated through ensemble predictions. Future climate change scenarios for 2050 and 2070 show that a further range expansion by this species would occur up to 225,576 km², mainly towards areas where the species was historically absent. The increase of isothermality (i.e., the ratio between the mean diurnal and the annual temperature range) and the mean temperature of the driest months would help crested porcupines to reach high altitudes, e.g., in the Alps. In mountain habitats, the ongoing global warming is shifting the distribution of European forests to high elevations, thus potentially providing porcupines with suitable habitats. A reduction in snow cover and the snow period at ground level would remove an important barrier to the range expansion of the crested porcupine in Italy, and thus facilitate digging and food search by this large rodent. Despite being protected at national and international levels, the crested porcupine is reported to be an introduced species to Italy and, therefore, monitoring its range expansion is required. Furthermore, there are complaints about crop damage in agriculture ecosystems, and the species is still widely poached, thus additional management practices are required. Thus, given the conservation interest of this large rodent, an integrated and constantly updated monitoring system that sustains an addressed set of decision-making tools is recommended.     

Impact of 3-D topography on surface radiation budget over the Tibetan Plateau

The 3-D complex topography effect on the surface solar radiative budget over the Tibetan Plateau is investigated by means of a parameterization approach on the basis of “exact” 3-D Monte Carlo photon tracing simulations, which use 90 m topography data as building blocks. Using a demonstrative grid size of 10?×?10 km², we show that differences in downward surface solar fluxes for a clear sky without aerosols between the 3-D model and the conventional plane-parallel radiative transfer scheme are substantial, on the order of 200 W/m² at shaded or sunward slopes. Deviations in the reflected fluxes of the direct solar beam amount to about +100 W/m² over snow-covered areas, which would lead to an enhanced snowmelt if the 3-D topography effects had been accounted for in current climate models. We further demonstrate that the entire Tibetan Plateau would receive more solar flux by about 14 W/m², if its 3-D mountain structure was included in the calculations, which would result in larger sensible and latent heat transfer from the surface to the atmosphere.     

Simulation of surface energy fluxes using high-resolution non-hydrostatic simulations and comparisons with measurements for the LITFASS-2003 experiment

Land-surface heterogeneity effects on the subgrid scale of regional climate and numerical weather prediction models are of vital interest for the energy and mass exchange between the surface and the atmospheric boundary layer. High-resolution numerical model simulations can be used to quantify these effects, and are a tool used to obtain area-averaged surface fluxes over heterogeneous land surfaces. We present high-resolution model simulations for the LITFASS area near Berlin during the LITFASS-2003 experiment, which were carried out using the non-hydrostatic model FOOT3DK of the University of Köln with horizontal resolutions of 1 km and 250 m. The LITFASS-2003 experimental dataset is used for comparison. The screen level quantities show good quality for the simulated pressure, temperature, humidity and wind speed and direction. Averaged over the four week experimental period, simulated surface energy fluxes at land stations show a small bias for the turbulent heat fluxes and an underestimation of the net radiation caused by excessive cloudiness in the simulations. For eight selected days with low cloud amounts, the net radiation bias is close to zero, but the sensible heat flux shows a strong positive bias. Large differences are found for latent heat fluxes over a lake, which are partly due to local effects on the measurements, but an additional problem seems to be the overestimation of the turbulent exchange under stable conditions in the daytime internal boundary layer over the lake. In the area average over the LITFASS area of 20 ×  20 km², again a strong positive bias of 70 W m^?2 for the sensible heat is present. For the low soil moisture conditions during June 2003, the simulation of the turbulent heat fluxes is sensitive to variations in the soil type and its hydrological properties. Under these conditions, the supply of ground water to the lowest soil layer should be accounted for. Different area-averaging methods are tested. The experimental set-up of the LITFASS-2003 experiment is found to be well suited for the computation of area-averaged turbulent heat fluxes.