Predicting thermal vulnerability of stream and river ecosystems to climate change

We use a predictive model of mean summer stream temperature to assess the vulnerability of USA streams to thermal alteration associated with climate change. The model uses air temperature and watershed features (e.g., watershed area and slope) from 569 US Geological Survey sites in the conterminous USA to predict stream temperatures. We assess the model for predicting climate-related variation in stream temperature by comparing observed and predicted historical stream temperature changes. Analysis of covariance confirms that observed and predicted changes in stream temperatures respond similarly to historical changes in air temperature. When applied to spatially-downscaled future air temperature projections (A2 emission scenario), the model predicts mean warming of 2.2 °C for the conterminous USA by 2100. Stream temperatures are most responsive to climate changes in the Cascade and Appalachian Mountains and least responsive in the southeastern USA. We then use random forests to conduct an empirical sensitivity analysis to identify those stream features most strongly associated with both observed historical and predicted future changes in summer stream temperatures. Larger changes in stream temperature are associated with warmer future air temperatures, greater air temperature changes, and larger watershed areas. Smaller changes in stream temperature are predicted for streams with high initial rates of heat loss associated with longwave radiation and evaporation, and greater base-flow index values. These models provide important insight into the potential extent of stream temperature warming at a near-continental scale and why some streams will likely be more vulnerable to climate change than others.     

Projecting impacts of climate change on surface water temperatures of a large subalpine lake: Lake Tahoe, USA

Predicted increases in atmospheric CO₂ concentration are expected to cause increases in air temperatures in many regions around the world, and this will likely lead to increases in the surface water temperatures of aquatic ecosystems in these regions. Using daily air and littoral water temperature data collected from Lake Tahoe, a large sub-alpine lake located in the Sierra Nevada mountains (USA), we developed and tested an empirical approach for constructing models designed to estimate site-specific daily surface water temperatures from daily air temperature projections generated from a regional climate model. We used cluster analysis to identify thermally distinct groups among sampled sites within the lake and then developed and independently validated a set of linked regression models designed to estimate daily water temperatures for each spatially distinct thermal group using daily air temperature data. When daily air temperatures projections, generated for 2080–2099 by a regional climate model, were used as input to these group models, projected increases in summer surface water temperatures of as much as 3 °C were projected. This study demonstrates an empirical approach for generating models capable of using daily air temperature projections from established climate models to project site specific impacts on littoral surface waters within large limnetic ecosystems.     

Observations of the convective plume of a lake under cold-air advective conditions

Moderating effects of Lake Apopka, Florida on downwind surface temperatures were evaluated under cold-air advective conditions. Point temperature measurements north and south of the lake and data obtained from a thermal scanner flown at 1.6 km indicate that surface temperatures directly downwind may be higher than surrounding surface temperatures by as much as 5 °C under conditions of moderate winds (~4 m s^–1). No substantial temperature effects were observed with surface wind speed less than 1 m s^–1. Fluxes of sensible and latent heat from Lake Apopka were calculated from measurements of lake temperature, net radiation, relative humidity and air temperature above the lake. Bulk transfer coefficients and the Bowen ratio were calculated and found to be in agreement with reported data for non-advective conditions.IFAS Journal Series No. 1006.     

Relative effects of climate change and wildfires on stream temperatures: a simulation modeling approach in a Rocky Mountain watershed

Freshwater ecosystems are warming globally from the direct effects of climate change on air temperature and hydrology and the indirect effects on near-stream vegetation. In fire-prone landscapes, vegetative change may be especially rapid and cause significant local stream temperature increases but the importance of these increases relative to broader changes associated with air temperature and hydrology are not well understood. We linked a spatially explicit landscape fire and vegetation model (FireBGCv2) to an empirical regression equation that predicted daily stream temperatures to explore how climate change and its impacts on fire might affect stream thermal conditions across a partially forested, mountainous landscape in the western U.S. We used the model to understand the roles that wildfire and management actions such as fuel reduction and fire suppression could play in mitigating stream thermal responses to climate change. Results indicate that air temperature increases associated with future climates could account for a much larger proportion of stream temperature increases (as much as 90 % at a basin scale) than wildfire. Similarly, land management scenarios that limited wildfire prevalence had negligible effects on future stream temperature increases. These patterns emerged at broader spatial scales because wildfires typically affected only a subset of a stream’s network. However, at finer spatial and temporal scales stream temperatures were sensitive to wildfire. Although wildfires will continue to cause local, short-term effects on stream temperatures, managers of aquatic systems may need to find other solutions to cope with the larger impact from climate change on future stream warming that involves adapting to the increases while developing broad strategies for riparian vegetation restoration.     

Stream temperature sensitivity to climate warming in California’s Sierra Nevada: impacts to coldwater habitat

Water temperature influences the distribution, abundance, and health of aquatic organisms in stream ecosystems, so understanding the impacts of climate warming on stream temperature will help guide management and restoration. This study assesses climate warming impacts on stream temperatures in California’s west-slope Sierra Nevada watersheds, and explores stream temperature modeling at the mesoscale. We used natural flow hydrology to isolate climate induced changes from those of water operations and land use changes. A 21 year time series of weekly streamflow estimates from WEAP21, a spatially explicit rainfall-runoff model were passed to RTEMP, an equilibrium temperature model, to estimate stream temperatures. Air temperature was uniformly increased by 2°C, 4°C, and 6°C as a sensitivity analysis to bracket the range of likely outcomes for stream temperatures. Other meteorological conditions, including precipitation, were unchanged from historical values. Raising air temperature affects precipitation partitioning into snowpack, runoff, and snowmelt in WEAP21, which change runoff volume and timing as well as stream temperatures. Overall, stream temperatures increased by an average of 1.6°C for each 2°C rise in air temperature, and increased most during spring and at middle elevations. Viable coldwater habitat shifted to higher elevations and will likely be reduced in California. Thermal heterogeneity existed within and between basins, with the high elevations of the southern Sierra Nevada and the Feather River watershed most resilient to climate warming. The regional equilibrium temperature modeling approach used here is well suited for climate change analysis because it incorporates mechanistic heat exchange, is not overly data or computationally intensive, and can highlight which watersheds are less vulnerable to climate warming. Understanding potential changes to stream temperatures from climate warming will affect how fish and wildlife are managed, and should be incorporated into modeling studies, restoration assessments, and licensing operations of hydropower facilities to best estimate future conditions and achieve desired outcomes.     

Tropical and extratropical cyclone damages under climate change

This paper provides the first quantitative synthesis of the rapidly growing literature on future tropical and extratropical cyclone damages under climate change. We estimate a probability distribution for the predicted impact of changes in global surface air temperatures on future storm damages, using an ensemble of 478 estimates of the temperature-damage relationship from nineteen studies. Our analysis produces three main empirical results. First, we find strong but not conclusive support for the hypothesis that climate change will cause damages from tropical cyclones and wind storms to increase, with most models predicting higher future storm damages due to climate change. Second, there is substantial variation in projected changes in losses across regions. Potential changes in damages are greatest in the North Atlantic basin, where the multi-model average predicts that a 2.5 °C increase in global surface air temperature would cause hurricane damages to increase by 63 %. The ensemble predictions for Western North Pacific tropical cyclones and European wind storms (extratropical cyclones) are +28 % and +23 %, respectively. Finally, our analysis shows that existing models of storm damages under climate change generate a wide range of predictions, ranging from moderate decreases to very large increases in losses.     

Role of Climate in the Modern Condition of Lake Victoria

Summary Comparison of historical and recent climatic data for Lake Victoria, a great lake in the high elevation tropics of East Africa, demonstrates changes in surface temperature, air moisture, atmospheric transparency, and wind shear from conditions 35 years ago. The changes appear to be part of a global change in climate conditions of the high elevation tropics. The physical and ecological conditions of Lake Victoria appear responsive to lake-atmosphere interactions through mechanisms of wind stress and surface heat fluxes. Lake temperatures, mixing regimes, oxygen levels, and primary production have changed in parallel to air temperature, humidity, atmospheric transparency, and wind speeds between about 1960 and the present. Data indicate strong coupling between meteorological components of lake heat balance and both biological and chemical conditions. Many features of the apparent modern eutrophication of Lake Victoria may have been accelerated or exaggerated by the climate effects. The analysis suggests a mechanism for feedback of climate on lake condition which, if general, might provide a modern analog for periodic changes reported in the fossil diatom community of the lake over the past 10 millenia. Received April 20, 1997 Revised January 1, 1998     

Projections and downscaling of 21st century temperatures, precipitation, radiative fluxes and winds for the Southwestern US, with focus on Lake Tahoe

Recent projections of global climate changes in response to increasing greenhouse-gas concentrations in the atmosphere include warming in the Southwestern US and, especially, in the vicinity of Lake Tahoe of from about +3°C to +6°C by end of century and changes in precipitation on the order of 5–10 % increases or (more commonly) decreases, depending on the climate model considered. Along with these basic changes, other climate variables like solar insolation, downwelling (longwave) radiant heat, and winds may change. Together these climate changes may result in changes in the hydrology of the Tahoe basin and potential changes in lake overturning and ecological regimes. Current climate projections, however, are generally spatially too coarse (with grid cells separated by 1 to 2° latitude and longitude) for direct use in assessments of the vulnerabilities of the much smaller Tahoe basin. Thus, daily temperatures, precipitation, winds, and downward radiation fluxes from selected global projections have been downscaled by a statistical method called the constructed-analogues method onto 10 to 12 km grids over the Southwest and especially over Lake Tahoe. Precipitation, solar insolation and winds over the Tahoe basin change only moderately (and with indeterminate signs) in the downscaled projections, whereas temperatures and downward longwave fluxes increase along with imposed increases in global greenhouse-gas concentrations.     

Sensitivity of climate change projections to uncertainties in the estimates of observed changes in deep-ocean heat content

The MIT 2D climate model is used to make probabilistic projections for changes in global mean surface temperature and for thermosteric sea level rise under a variety of forcing scenarios. The uncertainties in climate sensitivity and rate of heat uptake by the deep ocean are quantified by using the probability distributions derived from observed twentieth century temperature changes. The impact on climate change projections of using the smallest and largest estimates of twentieth century deep ocean warming is explored. The impact is large in the case of global mean thermosteric sea level rise. In the MIT reference (“business as usual”) scenario the median rise by 2100 is 27 and 43 cm in the respective cases. The impact on increases in global mean surface air temperature is more modest, 4.9 and 3.9 C in the two respective cases, because of the correlation between climate sensitivity and ocean heat uptake required by twentieth century surface and upper air temperature changes. The results are also compared with the projections made by the IPCC AR4’s multi-model ensemble for several of the SRES scenarios. The multi-model projections are more consistent with the MIT projections based on the largest estimate of ocean warming. However, the range for the rate of heat uptake by the ocean suggested by the lowest estimate of ocean warming is more consistent with the range suggested by the twentieth century changes in surface and upper air temperatures, combined with the expert prior for climate sensitivity.     

Simulated effects of land immersion on regional arid climate: a case study of the pre-Saharan playa of Chott el-Jerid (south of Tunisia)

The potential effects on the regional climate induced by partially immersing the arid pre-Saharan playa basin of Chott el-Jerid (south of Tunisia) are investigated by comparing two multi-year (1991–2011) sets of numerical simulations each consisting of ten-member ensemble and performed using the WRF regional climate model. The first WRF ensemble is performed under current land use and land cover, while the second is carried out after introducing a virtual large and shallow surface water reservoir (a lake) in Chott el-Jerid. The most pronounced effects generated by the artificial lake are circumscribed over its surface and slightly spread downwind to the other parts of the Chott. The lake has a clear moderating effect on near-surface air temperatures by increasing (decreasing) the wintertime (summertime) air temperatures. Sensible heat fluxes are remarkably increased in winter and decreased in summer over the lake following the temperature gradient between the lake surface and the overlying atmosphere. Latent heat fluxes, moisture convergence, and water vapor mixing ratio are increased over the lake throughout the year, especially in winter. The lake also induces domain-wide decreased (increased) surface pressures and land (lake) breeze circulation in winter (summer). Simulated rainfall amounts are most increased over the lake in winter likely because of an enhanced atmospheric instability, while they slightly decrease in summer.     

The Warming of the Tibetan Plateau in Response to Transient and Stabilized 2.0°C/1.5°C Global Warming Targets

As "the third pole", the Tibetan Plateau (TP) is sensitive to climate forcing and has experienced rapid warming in recent decades. This study analyzes annual and seasonal near-surface air temperature changes on the TP in response to transient and stabilized 2.0°C/1.5°C global warming targets based on simulations of the Community Earth System Model (CESM). Elevation-dependent warming (EDW) with faster warming at higher elevations is predicted. A surface energy budget analysis is adopted to uncover the mechanisms responsible for the temperature changes. Our results indicate a clear amplified warming on the TP with positive EDW in 2.0°C/1.5°C warmer futures, especially in the cold season. Mean TP warming relative to the reference period (1961–90) is dominated by an enhanced downward longwave radiation flux, while the variations in surface albedo shape the detailed pattern of EDW. For the same global warming level, the temperature changes under transient scenarios are ~0.2°C higher than those under stabilized scenarios, and the characteristics of EDW are broadly similar for both scenarios. These differences can be primarily attributed to the combined effects of differential downward longwave radiation, cloud radiative forcing, and surface sensible and latent heat fluxes. These findings contribute to a more detailed understanding of regional climate on the TP in response to the long-term climate goals of the Paris Agreement and highlight the differences between transient and stabilized warming scenarios.     

A role for icebergs in the 8.2 ka climate event

We investigate the potential role of icebergs in the 8.2 ka climate event, using a coupled climate model equipped with an iceberg component. First, we evaluate the effect of a large iceberg discharge originating from the decaying Laurentide ice sheet on ocean circulation, compared to a release of an identical volume of freshwater alone. Our results show that, on top of the freshwater effect, a large iceberg discharge facilitates sea-ice growth as a result of lower sea-surface temperatures induced by latent heat of melting. This causes an 8% increased sea-ice cover, 5% stronger reduction in North Atlantic Deep Water production and 1°C lower temperature in Greenland. Second, we use the model to investigate the effect of a hypothetical two-stage lake drainage, which is suggested by several investigators to have triggered the 8.2 ka climate event. To account for the final collapse of the ice-dam holding the Laurentide Lakes we accompany the secondary freshwater pulse in one scenario with a fast 5-year iceberg discharge and in a second scenario with a slow 100-year iceberg discharge. Our experiments show that a two-stage lake drainage accompanied by the collapsing ice-dam could explain the anomalies observed around the 8.2 ka climate event in various climate records. In addition, they advocate a potential role for icebergs in the 8.2 ka climate event and illustrate the importance of latent heat of melting in the simulation of climate events that involve icebergs. Our two-stage lake drainage experiments provide a framework in the discussion of two-stage lake drainage and ice sheet collapse.     

The Warming of Lake Tahoe

Summary We investigated the effects of climate variability on the thermal structure of Lake Tahoe, California-Nevada, 1970–2002, and with principal components analysis and step-wise multiple regression, related the volume-weighed average lake temperature to trends in climate. We then used a 1-dimensional hydrodynamic model to show that the observed trends in the climatic forcing variables can reasonably explain the observed changes in the lake. Between 1970 and 2002, the volume-weighted mean temperature of the lake increased at an average rate of 0.015 ^∘C yr⁻¹. Trends in the climatic drivers include 1) upward trends in maximum and minimum daily air temperature at Tahoe City; and 2) a slight upward trend in downward long-wave radiation. Changes in the thermal structure of the lake include 1) a long-term warming trend, with the highest rates near the surface and at 400 m; 2) an increase in the resistance of the lake to mixing and stratification, as measured by the Schmidt Stability and Birge Work; 3) a trend toward decreasing depth of the October thermocline. The long-term changes in the thermal structure of Lake Tahoe may interact with and exacerbate the well-documented trends in the lake's clarity and primary productivity.     

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.     

Development and validation of a regional coupled atmosphere lake model for the Caspian Sea Basin

We present a validation analysis of a regional climate model coupled to a distributed one dimensional (1D) lake model for the Caspian Sea Basin. Two model grid spacings are tested, 50 and 20 km, the simulation period is 1989–2008 and the lateral boundary conditions are from the ERA-Interim reanalysis of observations. The model is validated against atmospheric as well as lake variables. The model performance in reproducing precipitation and temperature mean seasonal climatology, seasonal cycles and interannual variability is generally good, with the model results being mostly within the observational uncertainty range. The model appears to overestimate cloudiness and underestimate surface radiation, although a large observational uncertainty is found in these variables. The 1D distributed lake model (run at each grid point of the lake area) reproduces the observed lake-average sea surface temperature (SST), although differences compared to observations are found in the spatial structure of the SST, most likely as a result of the absence of 3 dimensional lake water circulations. The evolution of lake ice cover and near surface wind over the lake area is also reproduced by the model reasonably well. Improvements resulting from the increase of resolution from 50 to 20 km are most significant in the lake model. Overall the performance of the coupled regional climate—1D lake model system appears to be of sufficient quality for application to climate change scenario simulations over the Caspian Sea Basin.     

Impact of Secular Climate Change on the Thermal Structure of a Large Temperate Central European Lake

Strong climate-related secular trends are apparent in a 52-yr long (1947–1998) uninterrupted series of monthly temperature profiles fromLake Zurich, a large, deep (136 m), temperate lake on the Swiss Plateau. Decadal mean water temperatures have undergone a secular increase at all depths, reflecting the high degree of regional warming that occurred in the European Alpine area during the 20th century. From the 1950s to the 1990s, high warming rates ( 0.24 K per decade) in the uppermost 20 m of the lake (i.e., the epi/metalimnion) combined with lower warming rates ( 0.13 K per decade) below 20 m (i.e., in the hypolimnion), have resulted in a20% increase in thermal stability and a consequent extension of 2–3 weeksin the stratification period. In common with many other parts of the world, 20th-century climate change on the Swiss Plateau has involved a steep secular increase in daily minimum (nighttime) air temperatures, but not in daily maximum (daytime) air temperatures. With respect to both secular change and decadal-scale variability, the temporal structure of the temperature of the surface mixed layer of Lake Zurich faithfully reflects that of the regional daily minimum air temperature, but not that of the daily maximum. The processes responsible for longer-term changes in the temperature structure of the lake therefore act during the night, presumably by suppressing nighttime convective cooling of the surface mixed layer. Application of a one-box heat exchange model suggests that the observed secular changes in thermal structure are due to shifts in the nighttime rate of emission of infrared radiation from the atmosphere and in the nighttime rates of latent and sensible heat exchange at the air-water interface. The increase in hypolimnetic temperatures is mainly a result of the increased prevalence of warm winters in Europe.     

An estimate of equilibrium sensitivity of global terrestrial carbon cycle using NCAR CCSM4

Increasing concentrations of atmospheric CO₂ influence climate, terrestrial biosphere productivity and ecosystem carbon storage through its radiative, physiological and fertilization effects. In this paper, we quantify these effects for a doubling of CO₂ using a low resolution configuration of the coupled model NCAR CCSM4. In contrast to previous coupled climate-carbon modeling studies, we focus on the near-equilibrium response of the terrestrial carbon cycle. For a doubling of CO₂, the radiative effect on the physical climate system causes global mean surface air temperature to increase by 2.14 K, whereas the physiological and fertilization on the land biosphere effects cause a warming of 0.22 K, suggesting that these later effects increase global warming by about 10 % as found in many recent studies. The CO₂-fertilization leads to total ecosystem carbon gain of 371 Gt-C (28 %) while the radiative effect causes a loss of 131 Gt-C (~10 %) indicating that climate warming damps the fertilization-induced carbon uptake over land. Our model-based estimate for the maximum potential terrestrial carbon uptake resulting from a doubling of atmospheric CO₂ concentration (285–570 ppm) is only 242 Gt-C. This highlights the limited storage capacity of the terrestrial carbon reservoir. We also find that the terrestrial carbon storage sensitivity to changes in CO₂ and temperature have been estimated to be lower in previous transient simulations because of lags in the climate-carbon system. Our model simulations indicate that the time scale of terrestrial carbon cycle response is greater than 500 years for CO₂-fertilization and about 200 years for temperature perturbations. We also find that dynamic changes in vegetation amplify the terrestrial carbon storage sensitivity relative to a static vegetation case: because of changes in tree cover, changes in total ecosystem carbon for CO₂-direct and climate effects are amplified by 88 and 72 %, respectively, in simulations with dynamic vegetation when compared to static vegetation simulations.     

The vertical distribution of climate forcings and feedbacks from the surface to top of atmosphere

The radiative forcings and feedbacks that determine Earth’s climate sensitivity are typically defined at the top-of-atmosphere (TOA) or tropopause, yet climate sensitivity itself refers to a change in temperature at the surface. In this paper, we describe how TOA radiative perturbations translate into surface temperature changes. It is shown using first principles that radiation changes at the TOA can be equated with the change in energy stored by the oceans and land surface. This ocean and land heat uptake in turn involves an adjustment of the surface radiative and non-radiative energy fluxes, with the latter being comprised of the turbulent exchange of latent and sensible heat between the surface and atmosphere. We employ the radiative kernel technique to decompose TOA radiative feedbacks in the IPCC Fourth Assessment Report climate models into components associated with changes in radiative heating of the atmosphere and of the surface. (We consider the equilibrium response of atmosphere-mixed layer ocean models subjected to an instantaneous doubling of atmospheric CO₂). It is shown that most feedbacks, i.e., the temperature, water vapor and cloud feedbacks, (as well as CO₂ forcing) affect primarily the turbulent energy exchange at the surface rather than the radiative energy exchange. Specifically, the temperature feedback increases the surface turbulent (radiative) energy loss by 2.87 W m^?2 K^?1 (0.60 W m^?2 K^?1) in the multimodel mean; the water vapor feedback decreases the surface turbulent energy loss by 1.07 W m^?2 K^?1 and increases the surface radiative heating by 0.89 W m^?2 K^?1; and the cloud feedback decreases both the turbulent energy loss and the radiative heating at the surface by 0.43 and 0.24 W m^?2 K^?1, respectively. Since changes to the surface turbulent energy exchange are dominated in the global mean sense by changes in surface evaporation, these results serve to highlight the fundamental importance of the global water cycle to Earth’s climate sensitivity.     

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.     

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.