Refine
Year of publication
Document Type
- Article (39)
- Book (1)
- Conference Proceeding (1)
- Preprint (1)
- Review (1)
Has Fulltext
- yes (43)
Is part of the Bibliography
- no (43)
Keywords
- COSMO-CLM (5)
- precipitation (4)
- Clausius–Clapeyron scaling (3)
- convective storms (3)
- tracking (3)
- Europe (2)
- Odra (2)
- convection-permitting simulation (2)
- homogeneity (2)
- transfer entropy (2)
Influence of sea surface roughness length parameterization on Mistral and Tramontane simulations
(2016)
The Mistral and Tramontane are mesoscale winds in southern France and above the Western Mediterranean Sea. They are phenomena well suited for studying channeling effects as well as atmosphere–land/ocean processes. This sensitivity study deals with the influence of the sea surface roughness length parameterizations on simulated Mistral and Tramontane wind speed and wind direction. Several simulations with the regional climate model COSMO-CLM were performed for the year 2005 with varying values for the Charnock parameter α. Above the western Mediterranean area, the simulated wind speed and wind direction pattern on Mistral days changes depending on the parameterization used. Higher values of α lead to lower simulated wind speeds. In areas, where the simulated wind speed does not change much, a counterclockwise rotation of the simulated wind direction is observed.
Evaluation of radiation components in a global freshwater model with station-based observations
(2016)
In many hydrological models, the amount of evapotranspired water is calculated using the potential evapotranspiration (PET) approach. The main driver of several PET approaches is net radiation, whose downward components are usually obtained from meteorological input data, whereas the upward components are calculated by the model itself. Thus, uncertainties can be large due to both the input data and model assumptions. In this study, we compare the radiation components of the WaterGAP Global Hydrology Model, driven by two meteorological input datasets and two radiation setups from ERA-Interim reanalysis. We assess the performance with respect to monthly observations provided by the Baseline Surface Radiation Network (BSRN) and the Global Energy Balance Archive (GEBA). The assessment is done for the global land area and specifically for energy/water limited regions. The results indicate that there is no optimal radiation input throughout the model variants, but standard meteorological input datasets perform better than those directly obtained by ERA-Interim reanalysis for the key variable net radiation. The low number of observations for some radiation components, as well as the scale mismatch between station observations and 0.5° × 0.5° grid cell size, limits the assessment.
Projections from coarse-grid global circulation models are not suitable for regional estimates of water balance or trends of extreme precipitation and temperature, especially not in complex terrain. Thus, downscaling of global to regionally resolved projections is necessary to provide input to integrated water resources management approaches for river basins like the Upper Danube River Basin (UDRB) and the Upper Brahmaputra River Basin (UBRB).
This paper discusses the application of the regional climate model COSMO-CLM as a dynamical downscaling tool. To provide accurate data the COSMO-CLM model output was post-processed by statistical means. This downscaling chain performs well in the baseline period 1971 to 2000. However, COSMO-CLM performs better in the UDRB than in the UBRB because of a longer application experience and a less complex climate in Europe.
Different climate change scenarios were downscaled for the time period 1960–2100. The projections show an increase of temperature in both basins and for all seasons. The values are generally higher in the UBRB with the highest values occurring in the region of the Tibetan Plateau. Annual precipitation shows no substantial change. However, seasonal amounts show clear trends, for instance an increasing amount of spring precipitation in the UDRB. Again, the largest trends for different precipitation statistics are projected in the region of the Tibetan Plateau. Here, the projections show up to 50% longer dry periods in the months June to September with a simultaneous increase of about 10% for the maximum amount of precipitation on five consecutive days. For the Assam region in India, the projections also show an increase of 25% in the number of consecutive dry days during the monsoon season leading to prolonged monsoon breaks.
Recently, new soil data maps were developed, which include vertical soil properties like soil type. Implementing those into a multilayer Soil-Vegetation-Atmosphere-Transfer (SVAT) scheme, discontinuities in the water content occur at the interface between dissimilar soils. Therefore, care must be taken in solving the Richards equation for calculating vertical soil water fluxes. We solve a modified form of the mixed (soil water and soil matric potential based) Richards equation by subtracting the equilibrium state of soil matrix potential ψE from the hydraulic potential ψh. The sensitivity of the modified equation is tested under idealized conditions. The paper will show that the modified equation can handle with discontinuities in soil water content at the interface of layered soils.
Towards the goal to understand the role of land-surface processes over the Indian sub-continent, a series of soil-moisture sensitivity simulations have been performed using a non-hydrostatic regional climate model COSMO-CLM. The experiments were driven by the lateral boundary conditions provided by the ERA-Interim (ECMWF) reanalysis. The simulation results show that the pre-monsoonal soil moisture has a significant influence on the monsoonal precipitation. Both, positive and negative soil-moisture precipitation (S-P) feedback processes are of importance. The negative S-P feedback process is especially influential in the western and the northern parts of India.
Mistral and tramontane wind speed and wind direction patterns in regional climate simulations
(2016)
The Mistral and Tramontane are important wind phenomena that occur over southern France and the northwestern Mediterranean Sea. Both winds travel through constricting valleys before flowing out towards the Mediterranean Sea. The Mistral and Tramontane are thus interesting phenomena, and represent an opportunity to study channeling effects, as well as the interactions between the atmosphere and land/ocean surfaces. This study investigates Mistral and Tramontane simulations using five regional climate models with grid spacing of about 50 km and smaller. All simulations are driven by ERA-Interim reanalysis data. Spatial patterns of surface wind, as well as wind development and error propagation along the wind tracks from inland France to offshore during Mistral and Tramontane events, are presented and discussed. To disentangle the results from large-scale error sources in Mistral and Tramontane simulations, only days with well simulated large-scale sea level pressure field patterns are evaluated. Comparisons with the observations show that the large-scale pressure patterns are well simulated by the considered models, but the orographic modifications to the wind systems are not well simulated by the coarse-grid simulations (with a grid spacing of about 50 km), and are reproduced slightly better by the higher resolution simulations. On days with Mistral and/or Tramontane events, most simulations underestimate (by 13 % on average) the wind speed over the Mediterranean Sea. This effect is strongest at the lateral borders of the main flow—the flow width is underestimated. All simulations of this study show a clockwise wind direction bias over the sea during Mistral and Tramontane events. Simulations with smaller grid spacing show smaller biases than their coarse-grid counterparts.
Spatial interpolation of rain gauge data is important in forcing of hydrological simulations or evaluation of weather predictions, for example. This paper investigates the application of statistical distance, like one minus common variance of observation time series, between data sites instead of geographical distance in interpolation. Here, as a typical representative of interpolation methods the inverse distance weighting interpolation is applied and the test data is daily precipitation observed in Austria. Choosing statistical distance instead of geographical distance in interpolation of available coarse network observations to sites of a denser network, which is not reporting for the interpolation date, yields more robust interpolation results. The most distinct performance enhancement is in or close to mountainous terrain. Therefore, application of statistical distance in the inverse distance weighting interpolation or in similar methods can parsimoniously densify the currently available observation network. Additionally, the success further motivates search for conceptual rain-orography interaction models as components of spatial rain interpolation algorithms in mountainous terrain.
This study presents a method for adjusting long-term climate data records (CDRs) for the integrated use with near-real-time data using the example of surface incoming solar irradiance (SIS). Recently, a 23-year long (1983–2005) continuous SIS CDR has been generated based on the visible channel (0.45–1 μm) of the MVIRI radiometers onboard the geostationary Meteosat First Generation Platform. The CDR is available from the EUMETSAT Satellite Application Facility on Climate Monitoring (CM SAF). Here, it is assessed whether a homogeneous extension of the SIS CDR to the present is possible with operationally generated surface radiation data provided by CM SAF using the SEVIRI and GERB instruments onboard the Meteosat Second Generation satellites. Three extended CM SAF SIS CDR versions consisting of MVIRI-derived SIS (1983–2005) and three different SIS products derived from the SEVIRI and GERB instruments onboard the MSG satellites (2006 onwards) were tested. A procedure to detect shift inhomogeneities in the extended data record (1983–present) was applied that combines the Standard Normal Homogeneity Test (SNHT) and a penalized maximal T-test with visual inspection. Shift detection was done by comparing the SIS time series with the ground stations mean, in accordance with statistical significance. Several stations of the Baseline Surface Radiation Network (BSRN) and about 50 stations of the Global Energy Balance Archive (GEBA) over Europe were used as the ground-based reference. The analysis indicates several breaks in the data record between 1987 and 1994 probably due to artefacts in the raw data and instrument failures. After 2005 the MVIRI radiometer was replaced by the narrow-band SEVIRI and the broadband GERB radiometers and a new retrieval algorithm was applied. This induces significant challenges for the homogenisation across the satellite generations. Homogenisation is performed by applying a mean-shift correction depending on the shift size of any segment between two break points to the last segment (2006–present). Corrections are applied to the most significant breaks that can be related to satellite changes. This study focuses on the European region, but the methods can be generalized to other regions. To account for seasonal dependence of the mean-shifts the correction was performed independently for each calendar month. In comparison to the ground-based reference the homogenised data record shows an improvement over the original data record in terms of anomaly correlation and bias. In general the method can also be applied for the adjustment of satellite datasets addressing other variables to bridge the gap between CDRs and near-real-time data.
So-called medicanes (Mediterranean hurricanes) are meso-scale, marine, and warm-core Mediterranean cyclones that exhibit some similarities to tropical cyclones. The strong cyclonic winds associated with medicanes threaten the highly populated coastal areas around the Mediterranean basin. To reduce the risk of casualties and overall negative impacts, it is important to improve the understanding of medicanes with the use of numerical models. In this study, we employ an atmospheric limited-area model (COSMO-CLM) coupled with a one-dimensional ocean model (1-D NEMO-MED12) to simulate medicanes. The aim of this study is to assess the robustness of the coupled model in simulating these extreme events. For this purpose, 11 historical medicane events are simulated using the atmosphere-only model, COSMO-CLM, and coupled model, with different setups (horizontal atmospheric grid-spacings of 0.44°, 0.22°, and 0.08°; with/without spectral nudging, and an ocean grid-spacing of 1/12°). The results show that at high-resolution, the coupled model is able to not only simulate most of medicane events but also improve the track length, core temperature, and wind speed of simulated medicanes compared to the atmosphere-only simulations. The results suggest that the coupled model is more proficient for systemic and detailed studies of historical medicane events, and that this model can be an effective tool for future projections.
A satellite-based climate record of monthly mean surface solar irradiance (SIS) is investigated with regard to possible inhomogeneities in time. The data record is provided by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Climate Monitoring (CM SAF) for the period of 1983 to 2005, covering a disk area between ±70° in latitude and longitude. The Standard Normal Homogeneity Test (SNHT) and two other homogeneity tests are applied with and without the use of reference SIS data (from the Baseline Surface Radiation Network (BSRN) and from the ECMWF (European Centre for Medium-Range Weather Forecasts) ERA -Interim reanalysis. The focus is on the detection of break-like inhomogeneities, which may occur due to satellite or SIS retrieval algorithm changes. In comparison with the few suitable BSRN SIS observation series with limited extension in time (no data before 1992), the CM SAF SIS time series do not show significant inhomogeneities, even though slight discrepancies in the surface measurements appear. The investigation of the full CM SAF SIS domain reveal inhomogeneities related to most of the documented satellite and retrieval changes, but only for relatively small domain fractions (especially in mountainous desert-like areas in Africa). In these regions the retrieval algorithm is not capable of adjusting for the changes of the satellite instruments. For other areas, e.g., Europe, no such breaks in the time series are found. We conclude that the CM SAF SIS data record has to be further assessed and regionally homogenized before climate trend investigations can be conducted.
In this study, two different methods were applied to derive daily and monthly sunshine duration based on high-resolution satellite products provided by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Climate Monitoring using data from Meteosat Second Generation (MSG) SEVIRI (Spinning Enhanced Visible and Infrared Imager). The satellite products were either hourly cloud type or hourly surface incoming direct radiation. The satellite sunshine duration estimates were not found to be significantly different using the native 15-minute temporal resolution of SEVIRI. The satellite-based sunshine duration products give additional spatial information over the European continent compared with equivalent in situ-based products. An evaluation of the satellite sunshine duration by product intercomparison and against station measurements was carried out to determine their accuracy. The satellite data were found to be within ±1 h/day compared to high-quality Baseline Surface Radiation Network or surface synoptic observations (SYNOP) station measurements. The satellite-based products differ more over the oceans than over land, mainly because of the treatment of fractional clouds in the cloud type-based sunshine duration product. This paper presents the methods used to derive the satellite sunshine duration products and the performance of the different retrievals. The main benefits and disadvantages compared to station-based products are also discussed.
The frequency of extreme events has changed, having a direct impact on human lives. Regional climate models help us to predict these regional climate changes. This work presents an atmosphere–ocean coupled regional climate system model (RCSM; with the atmospheric component COSMO-CLM and the ocean component NEMO) over the European domain, including three marginal seas: the Mediterranean, North, and Baltic Sea. To test the model, we evaluate a simulation of more than 100 years (1900–2009) with a spatial grid resolution of about 25 km. The simulation was nested into a coupled global simulation with the model MPI-ESM in a low-resolution configuration, whose ocean temperature and salinity were nudged to the ocean–ice component of the MPI-ESM forced with the NOAA 20th Century Reanalysis (20CR). The evaluation shows the robustness of the RCSM and discusses the added value by the coupled marginal seas over an atmosphere-only simulation. The coupled system is stable for the complete 20th century and provides a better representation of extreme temperatures compared to the atmosphere-only model. The produced long-term dataset will help us to better understand the processes leading to meteorological and climate extremes.
A twentieth century-long coupled atmosphere-ocean regional climate simulation with COSMO-CLM (Consortium for Small-Scale Modeling, Climate Limited-area Model) and NEMO (Nucleus for European Modelling of the Ocean) is studied here to evaluate the added value of coupled marginal seas over continental regions. The interactive coupling of the marginal seas, namely the Mediterranean, the North and the Baltic Seas, to the atmosphere in the European region gives a comprehensive modelling system. It is expected to be able to describe the climatological features of this geographically complex area even more precisely than an atmosphere-only climate model. The investigated variables are precipitation and 2 m temperature. Sensitivity studies are used to assess the impact of SST (sea surface temperature) changes over land areas. The different SST values affect the continental precipitation more than the 2 m temperature. The simulated variables are compared to the CRU (Climatic Research Unit) observational data, and also to the HOAPS/GPCC (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data, Global Precipitation Climatology Centre) data. In the coupled simulation, added skill is found primarily during winter over the eastern part of Europe. Our analysis shows that, over this region, the coupled system is dryer than the uncoupled system, both in terms of precipitation and soil moisture, which means a decrease in the bias of the system. Thus, the coupling improves the simulation of precipitation over the eastern part of Europe, due to cooler SST values and in consequence, drier soil.
Often in climate system studies, linear and symmetric statistical measures are applied to quantify interactions among subsystems or variables. However, they do not allow identification of the driving and responding subsystems. Therefore, in this study, we aimed to apply asymmetric measures from information theory: the axiomatically proposed transfer entropy and the first principle-based information flow to detect and quantify climate interactions. As their estimations are challenging, we initially tested nonparametric estimators like transfer entropy (TE)-binning, TE-kernel, and TE k-nearest neighbor and parametric estimators like TE-linear and information flow (IF)-linear with idealized two-dimensional test cases along with their sensitivity on sample size. Thereafter, we experimentally applied these methods to the Lorenz-96 model and to two real climate phenomena, i.e., (1) the Indo-Pacific Ocean coupling and (2) North Atlantic Oscillation (NAO)–European air temperature coupling. As expected, the linear estimators work for linear systems but fail for strongly nonlinear systems. The TE-kernel and TE k-nearest neighbor estimators are reliable for linear and nonlinear systems. Nevertheless, the nonparametric methods are sensitive to parameter selection and sample size. Thus, this work proposes a composite use of the TE-kernel and TE k-nearest neighbor estimators along with parameter testing for consistent results. The revealed information exchange in Lorenz-96 is dominated by the slow subsystem component. For real climate phenomena, expected bidirectional information exchange between the Indian and Pacific SSTs was detected. Furthermore, expected information exchange from NAO to European air temperature was detected, but also unexpected reversal information exchange. The latter might hint to a hidden process driving both the NAO and European temperatures. Hence, the limitations, availability of time series length and the system at hand must be taken into account before drawing any conclusions from TE and IF-linear estimations.
Moisture sources of heavy precipitation in Central Europe in synoptic situations with Vb-cyclones
(2022)
During the past century, several extreme summer floods in Central Europe were associated with so-called Vb-cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. The processes intensifying the precipitation in synoptic situations with Vb-cyclones in the Danube, Elbe, and Odra catchments are only partially understood. Our study aims to investigate these processes with Lagrangian moisture-source diagnostics for 16 selected Vb-events. Moreover, we analyse the characteristics of typical moisture source regions during 1107 Vb-events from 1901 to 2010 based on ERA-20C reanalysis dynamically downscaled with COSMO-CLM+NEMO. We observe moisture contributions by various source regions highlighting the complex dynamical interplay of different air masses leading to moisture convergence in synoptic situations with Vb-cyclones. Overall, up to 80% of the precipitation originates from the European continent, indicating the importance of continental moisture recycling, especially within the respective river catchment. Other major moisture uptake regions are the North Sea, the Baltic Sea, the North Atlantic, and for a few events the Black Sea. Remarkably, anomalies in these oceanic source regions show no connection to precipitation amounts in synoptic situations with Vb-cyclones. In contrast, the Vb-cyclones with the highest precipitation are associated with anomalously high evaporation in the Mediterranean Sea, even though the Mediterranean Sea is only a minor moisture source region on average. Interestingly, the evaporation anomalies are not connected with sea-surface temperature but with wind-speed anomalies (Spearman’s rank correlation coefficient R≈0.7, significant with p<0.01) indicating mainly dynamically driven evaporation. The particular role of the Mediterranean Sea hints towards possible importance of Mediterranean moisture for the early-stage intensification of Vb-cyclones and the pre-moistening of the continental uptake regions upstream of the target catchments.
The goal of limited area models (LAMs) is to downscale coarse-gridded general circulation model output to represent small-scale features of weather and climate. The LAM needs information from the driving coarse-gridded model passing through its lateral boundaries. The treatment of this information transfer causes inconsistencies between driving and nested models and, subsequently, issues in regional weather and climate simulations. This work examines errors arising from choices taken by the modeler (temporal update frequency of boundary data, spatial resolution jump, and numerical lateral boundary formulation) systematically in an idealized simulation environment. So-called Big-Brother Experiments were performed with the LAM COSMO-CLM (0.11° grid spacing). A baroclinic wave in a zonal channel was simulated over flat terrain with and without a Gaussian hill. The results reveal that the quality of the driving data, here represented by simulations only differing from the LAM simulations by reduced spatial resolution, dominates the performance of the nested model. Consequently, at the simulated mesoscale, the performance of the nested small-scale model simulations is weakly sensitive to the numerical lateral boundary formulation (Davies relaxation or the newly implemented, computationally less demanding Mesinger Eta-model formulation). The performance sensitivity to boundary update frequency and resolution jump is small when at least 6-hourly updates and a resolution jump factor of maximally six is used. Gaussian hill LAM simulations illustrated the strength of downscaling; they can represent small-scale features missing in the coarse-scale driving simulations. In the idealized simulation experiments, spectral nudging is not advisable as it imprints the driving models deficits on the nested simulation.
Background:Aedes aegypti is a potential vector for several arboviruses including dengue and Zika viruses. The species seems to be restricted to subtropical/tropical habitats and has difficulties in establishing permanent populations in southern Europe, probably due to constraints during the winter season. The aim of this study was to systematically analyze the cold tolerance (CT) of Ae. aegypti in its most cold-resistant life stage, the eggs.
Methods: The CT of Ae. aegypti eggs was compared with that of Ae. albopictus which is well established in large parts of Europe. By systematically studying the literature (meta-analysis), we recognized that CT has been rarely tested in Ae. aegypti eggs, but eggs can survive at zero and sub-zero temperatures for certain exposure periods. To overcome potential bias from experimental differences between studies, we then conducted species comparisons using a harmonized high-resolution CT measuring method. From subtropical populations of the same origin, the survival (hatching in %) and emergence of adults of both species were measured after zero and sub-zero temperature exposures for up to 9 days (3 °C, 0 °C and − 2 °C: ≤ 9 days; − 6 °C: ≤ 2 days).
Results: Our data show that Ae. aegypti eggs can survive low and sub-zero temperatures for a short time period similar to or even better than those of Ae. albopictus. Moreover, after short sub-zero exposures of eggs of both species, individuals still developed into viable adults (Ae. aegypti: 3 adults emerged after 6 days at − 2 °C, Ae. albopictus: 1 adult emerged after 1 day at − 6 °C).
Conclusions: Thus, both the literature and the present experimental data indicate that a cold winter may not be the preventing factor for the re-establishment of the dengue vector Ae. aegypti in southern Europe.
Container-breeding Aedes spp. (Diptera: Culicidae) mosquitoes can be surveilled at low cost using ovitraps. Hence, this method is a preferred monitoring approach of dengue vectors in low-resource settings. The ovitraps consist of a cup filled with water and an oviposition substrate for female mosquitoes. The attractiveness of the substrates for female mosquitoes can greatly differ due to differences in texture, color, and smell of the materials used. We compare four oviposition substrates, which are all low priced, easy to transport, and easy to purchase, to maximize the success of Aedes egg sampling. Sampled egg material is often reared to adulthood for further taxonomic identification and transported to (international) laboratories for specialized vector research. Here we introduce a transport technique for sampled eggs. In addition, we explored the impact of international transport by means of a bilateral hatching experiment in Nepal, the country of origin, and in Germany, in a laboratory specialized in ecophysiological research. The best low-cost oviposition substrate for the dengue vectors Aedes albopictus (Skuse) and Aedes aegypti (L.) was found to be a white cotton sheet. The introduced transport technique of sampled eggs is easy to build from laboratory and household materials and ensures good transport conditions (i.e., temperature and relative humidity). Even under good temperature (17.4–31.0 °C) and humidity conditions (58.9–94.2%), hatching success of eggs was found to be reduced after international transport to Germany when compared to the hatching success of eggs in Nepal. We postulate that air pressure during international transport may have reduced the hatching success and strongly recommend pressure-regulated transport boxes for egg transport via airplane. As the proposed operation procedure is useful in assisting the monitoring of Ae. albopictus and Ae. aegypti in low-resource settings, Aedes researchers are encouraged to follow it for the sampling and transport of Aedes eggs.
Several past summer floods in Central Europe were associated with so-called Vb‑cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. This study illustrates the usefulness of the parametric transfer entropy measure TE‑linear in investigating heavy Vb‑cyclone precipitation events in the Odra catchment (Poland). With the application of the TE‑linear approach, we confirm the impact of the Mediterranean Sea on precipitation intensification. Moreover, we also detect significant information exchange to Vb‑cyclone precipitation from evaporation over the European continent along the typical Vb‑cyclone pathway. Thus, the Mediterranean Sea could enhance the Vb‑cyclone precipitation by pre-moistening continental moisture source regions that contribute to precipitation downstream in the investigated catchments. Overall, the transfer entropy approach with the measure TE‑linear proved to be computationally effective and complementary to traditional methods such as Lagrangian and Eulerian diagnostics.
In the last decade, the Climate Limited-area Modeling Community (CLM-Community) has contributed to the Coordinated Regional Climate Downscaling Experiment (CORDEX) with an extensive set of regional climate simulations. Using several versions of the COSMO-CLM-Community model, ERA-Interim reanalysis and eight global climate models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) were dynamically downscaled with horizontal grid spacings of 0.44∘ (∼ 50 km), 0.22∘ (∼ 25 km), and 0.11∘ (∼ 12 km) over the CORDEX domains Europe, South Asia, East Asia, Australasia, and Africa. This major effort resulted in 80 regional climate simulations publicly available through the Earth System Grid Federation (ESGF) web portals for use in impact studies and climate scenario assessments. Here we review the production of these simulations and assess their results in terms of mean near-surface temperature and precipitation to aid the future design of the COSMO-CLM model simulations. It is found that a domain-specific parameter tuning is beneficial, while increasing horizontal model resolution (from 50 to 25 or 12 km grid spacing) alone does not always improve the performance of the simulation. Moreover, the COSMO-CLM performance depends on the driving data. This is generally more important than the dependence on horizontal resolution, model version, and configuration. Our results emphasize the importance of performing regional climate projections in a coordinated way, where guidance from both the global (GCM) and regional (RCM) climate modeling communities is needed to increase the reliability of the GCM–RCM modeling chain.
In the last decade, the Climate Limited-area Modeling (CLM) Community has contributed to the Coordinated Regional Climate Downscaling Experiment (CORDEX) with an extensive set of regional climate simulations. Using several versions of the COSMO-CLM community model, ERA-Interim reanalysis and eight Global Climate Models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) were dynamically downscaled with horizontal grid spacings of 0.44◦(∼50 km), 0.22◦ (∼25 km) and 0.11◦ (∼12 km) over the CORDEX domains Europe, South Asia, East Asia, Australasia and Africa. This major effort resulted in 80 regional climate simulations publicly available through the Earth System Grid Federation (ESGF) web portals for use in impact studies and climate scenario assessments. Here we review the production of these simulations and assess their results in terms of mean near-surface temperature and precipitation to aid the future design of the COSMO-CLM model simulations. It is found that a domain-specific parameter tuning is beneficial, while increasing horizontal model resolution (from 50 to 25 or 12 km grid spacing) alone does not always improve the performance of the simulation. Moreover, the COSMO-CLM performance depends on the driving data. This is generally more important than the dependence on horizontal resolution, model version and configuration. Our results emphasize the importance of performing regional climate projections in a coordinated way, where guidance from both the global (GCM) and regional (RCM) climate modelling communities is needed to increase the reliability of the GCM-RCM modelling chain.
Im Rahmen einer Zusammenarbeit zwischen der Thüringer Landesanstalt für Umwelt und Geologie und der Goethe-Universität Frankfurt fand in Kooperation mit dem Deutschen Wetterdienst (DWD) eine umfassende Studie zum konvektiven Unwetterpotential über Thüringen statt. Unwetterereignisse, die durch konvektive Prozesse in der Atmosphäre verursacht werden, besitzen ein nicht unerhebliches Schadenspotential, obwohl sie oftmals eine räumlich eng begrenzte Ausdehnung aufweisen. Aufgrund ihrer Charakteristik ist sowohl die Vorhersage solcher Ereignisse, als auch eine vollständige, systematische Erfassung für eine detaillierte Auswertung längerer Zeitreihen noch immer eine Herausforderung. Zusätzliches Interesse besteht in der Abschätzung der durch den Klimawandel abhängigen Entwicklung des zukünftigen Gefährdungspotentials konvektiver Unwetter. Für eine gezielte Untersuchung des Themenkomplexes ist eine Vielzahl unterschiedlicher Daten und Methoden verwendet worden. Mit Hilfe von Fernerkundungsdatensätzen wird ein räumlich differenziertes Gefährdungspotential über Thüringen nachgewiesen. Bedingt durch das Relief ist das Auftreten von Konvektion am häufigsten und intensivsten über dem südlichen Thüringer Wald und dessen Ostrand zu beobachten, während Nordthüringen eine deutlich geringere Aktivität solcher Unwetterereignisse aufweist. Eine Abschätzung mittels globaler Klimamodelle und daraus abgeleiteten Wetterlagen zeigt unter Berücksichtigung des RCP8.5 Klimaszenarios für die nahe Zukunft (2016-2045) eine Zunahme des Gefährdungspotentials durch konvektive Unwetter. Aufgrund des Anstiegs feuchter Wetterlagen (49 % auf 82 %) erhöht sich die Zunahme der Gefährdung für den Zeitraum 2071-2100 noch deutlicher. Im Vergleich zu diesem statistischen Ansatz nimmt die projizierte Gefährdung durch extreme Ereignisse erheblich zu (Faktor 6), wenn die Ergebnisse expliziter Simulationen konvektiver Ereignisse mit einem regionalen Klimamodell (mit horizontaler Gitterdistanz von 1 km) und eine Zunahme der Tage mit konvektiven Extremereignissen berücksichtigt werden. Ein Anstieg der Gefährdung durch konvektive Unwetter in der Zukunft ist wahrscheinlich. Eine Quantifizierung bleibt jedoch unsicher.
Convective shower characteristics simulated with the convection-permitting climate model COSMO-CLM
(2019)
This paper evaluates convective precipitation as simulated by the convection-permitting climate model (CPM) Consortium for Small-Scale Modeling in climate mode (COSMO-CLM) (with 2.8 km grid-spacing) over Germany in the period 2001–2015. Characteristics of simulated convective precipitation objects like lifetime, area, mean intensity, and total precipitation are compared to characteristics observed by weather radar. For this purpose, a tracking algorithm was applied to simulated and observed precipitation with 5-min temporal resolution. The total amount of convective precipitation is well simulated, with a small overestimation of 2%. However, the simulation underestimates convective activity, represented by the number of convective objects, by 33%. This underestimation is especially pronounced in the lowlands of Northern Germany, whereas the simulation matches observations well in the mountainous areas of Southern Germany. The underestimation of activity is compensated by an overestimation of the simulated lifetime of convective objects. The observed mean intensity, maximum intensity, and area of precipitation objects increase with their lifetime showing the spectrum of convective storms ranging from short-living single-cell storms to long-living organized convection like supercells or squall lines. The CPM is capable of reproducing the lifetime dependence of these characteristics but shows a weaker increase in mean intensity with lifetime resulting in an especially pronounced underestimation (up to 25%) of mean precipitation intensity of long-living, extreme events. This limitation of the CPM is not identifiable by classical evaluation techniques using rain gauges. The simulation can reproduce the general increase of the highest percentiles of cell area, total precipitation, and mean intensity with temperature but fails to reproduce the increase of lifetime. The scaling rates of mean intensity and total precipitation resemble observed rates only in parts of the temperature range. The results suggest that the evaluation of coarse-grained (e.g., hourly) precipitation fields is insufficient for revealing challenges in convection-permitting simulations.
Convection-permitting models (CPMs) have proven their usefulness in representing precipitation on a sub-daily scale. However, investigations on sub-hourly scales are still lacking, even though these are the scales for which showers exhibit the most variability. A Lagrangian approach is implemented here to evaluate the representation of showers in a CPM, using the limited-area climate model COSMO-CLM. This approach consists of tracking 5‑min precipitation fields to retrieve different features of showers (e.g., temporal pattern, horizontal speed, lifetime). In total, 312 cases are simulated at a resolution of 0.01 ° over Central Germany, and among these cases, 78 are evaluated against a radar dataset. The model is able to represent most observed features for different types of convective cells. In addition, the CPM reproduced well the observed relationship between the precipitation characteristics and temperature indicating that the COSMO-CLM model is sophisticated enough to represent the climatological features of showers.
Observed weather and projected climate change suggest an increase in the transmission of vector-borne diseases (VBDs) in the Hindu Kush Himalayan (HKH) region. In this study, we systematically explore the literature for empiric associations between the climate variables and specific VBDs and their vectors in the HKH region. We conducted a systematic synthesis of the published literature on climate variables, VBDs and vectors in the HKH region until the 8th of December 2020. The majority of studies show significant positive associations of VBDs with climatic factors, such as temperature, precipitation, relative humidity, etc. This systematic review allowed us to identify the most significant variables to be considered for evidence-based trend estimates of the effects of climate change on VBDs and their vectors in the HKH region. This evidence-based trend was set into the context of climate change as well as the observed expansion of VBDs and disease vectors in the HKH region. The geographic range of VBDs expanded into previously considered non-endemic areas of highlands (mountains) in the HKH region. Based on scarce, but clear evidence of a positive relationship of most climate variables and VBDs and the observed climatic changes, we strongly recommend an expansion of vector control and surveillance programmes in areas of the HKH region that were previously considered to be non-endemic.
Lightning climate change projections show large uncertainties caused by limited empirical knowledge and strong assumptions inherent to coarse-grid climate modeling. This study addresses the latter issue by implementing and applying the lightning potential index parameterization (LPI) into a fine-grid convection-permitting regional climate model (CPM). This setup takes advantage of the explicit representation of deep convection in CPMs and allows for process-oriented LPI inputs such as vertical velocity within convective cells and coexistence of microphysical hydrometeor types, which are known to contribute to charge separation mechanisms. The LPI output is compared to output from a simpler flash rate parameterization, namely the CAPE × PREC parameterization, applied in a non-CPM on a coarser grid. The LPI’s implementation into the regional climate model COSMO-CLM successfully reproduces the observed lightning climatology, including its latitudinal gradient, its daily and hourly probability distributions, and its diurnal and annual cycles. Besides, the simulated temperature dependence of lightning reflects the observed dependency. The LPI outperforms the CAPE × PREC parameterization in all applied diagnostics. Based on this satisfactory evaluation, we used the LPI to a climate change projection under the RCP8.5 scenario. For the domain under investigation centered over Germany, the LPI projects a decrease of 4.8% in flash rate by the end of the century, in opposition to a projected increase of 17.4% as projected using the CAPE × PREC parameterization. The future decrease of LPI occurs mostly during the summer afternoons and is related to (i) a change in convection occurrence and (ii) changes in the microphysical mixing. The two parameterizations differ because of different convection occurrences in the CPM and non-CPM and because of changes in the microphysical mixing, which is only represented in the LPI lightning parameterization.
Spatial interpolation of precipitation data is uncertain. How important is this uncertainty and how can it be considered in evaluation of high-resolution probabilistic precipitation forecasts? These questions are discussed by experimental evaluation of the COSMO consortium's limited-area ensemble prediction system COSMO-LEPS. The applied performance measure is the often used Brier skill score (BSS). The observational references in the evaluation are (a) analyzed rain gauge data by ordinary Kriging and (b) ensembles of interpolated rain gauge data by stochastic simulation. This permits the consideration of either a deterministic reference (the event is observed or not with 100% certainty) or a probabilistic reference that makes allowance for uncertainties in spatial averaging. The evaluation experiments show that the evaluation uncertainties are substantial even for the large area (41 300 km2) of Switzerland with a mean rain gauge distance as good as 7 km: the one- to three-day precipitation forecasts have skill decreasing with forecast lead time but the one- and two-day forecast performances differ not significantly.
Driven by globalization, urbanization and climate change, the distribution range of invasive vector species has expanded to previously colder ecoregions. To reduce health-threatening impacts on humans, insect vectors are extensively studied. Population genomics can reveal the genomic basis of adaptation and help to identify emerging trends of vector expansion. By applying whole genome analyses and genotype-environment associations to populations of the main dengue vector Aedes aegypti, sampled along an altitudinal gradient in Nepal (200–1300 m), we identify putatively adaptive traits and describe the species' genomic footprint of climate adaptation to colder ecoregions. We found two differentiated clusters with significantly different allele frequencies in genes associated to climate adaptation between the highland population (1300 m) and all other lowland populations (≤800 m). We revealed nonsynonymous mutations in 13 of the candidate genes associated to either altitude, precipitation or cold tolerance and identified an isolation-by-environment differentiation pattern. Other than the expected gradual differentiation along the altitudinal gradient, our results reveal a distinct genomic differentiation of the highland population. Local high-altitude adaptation could be one explanation of the population's phenotypic cold tolerance. Carrying alleles relevant for survival under colder climate increases the likelihood of this highland population to a worldwide expansion into other colder ecoregions.
Extreme convective precipitation is expected to increase with global warming. However, the rate of increase and the understanding of contributing processes remain highly uncertain. We investigated characteristics of convective rain cells like area, intensity, and lifetime as simulated by a convection-permitting climate model in the area of Germany under historical (1976–2005) and future (end-of-century, RCP8.5 scenario) conditions. To this end, a tracking algorithm was applied to 5-min precipitation output. While the number of convective cells is virtually similar under historical and future conditions, there are more intense and larger cells in the future. This yields an increase in hourly precipitation extremes, although mean precipitation decreases. The relative change in the frequency distributions of area, intensity, and precipitation sum per cell is highest for the most extreme percentiles, suggesting that extreme events intensify the most. Furthermore, we investigated the temperature and moisture scaling of cell characteristics. The temperature scaling drops off at high temperatures, with a shift in drop-off towards higher temperatures in the future, allowing for higher peak values. In contrast, dew point temperature scaling shows consistent rates across the whole dew point range. Cell characteristics scale at varying rates, either below (mean intensity), at about (maximum intensity and area), or above (precipitation sum) the Clausius–Clapeyron rate. Thus, the widely investigated extreme precipitation scaling at fixed locations is a complex product of the scaling of different cell characteristics. The dew point scaling rates and absolute values of the scaling curves in historical and future conditions are closest for the highest percentiles. Therefore, near-surface humidity provides a good predictor for the upper limit of for example, maximum intensity and total precipitation of individual convective cells. However, the frequency distribution of the number of cells depending on dew point temperature changes in the future, preventing statistical inference of extreme precipitation from near-surface humidity.
Convective rain cell properties and the resulting precipitation scaling in a warm-temperate climate
(2022)
Convective precipitation events have been shown to intensify at rates exceeding the Clausius–Clapeyron rate (CC rate) of ca. 7% K−1 under current climate conditions. In this study, we relate atmospheric variables (low-level dew point temperature, convective available potential energy, and vertical wind shear), which are regarded as ingredients for severe deep convection, to properties of convective rain cells (cell area, maximum precipitation intensity, lifetime, precipitation sum, and cell speed). The rain cell properties are obtained from a rain gauge-adjusted radar dataset in a mid-latitude region, which is characterized by a temperate climate with warm summers (Germany). Different Lagrangian cell properties scale with dew point temperature at varying rates. While the maximum precipitation intensity of cells scales consistently at the CC rate, the area and precipitation sum per cell scale at varying rates above the CC rate. We show that this super-CC scaling is caused by a covarying increase of convective available potential energy with dew point temperature. Wind shear increases the precipitation sum per cell mainly by increasing the spatial cell extent. From a Eulerian point of view, this increase is partly compensated by a higher cell velocity, which leads to Eulerian precipitation scaling rates close to and slightly above the CC rate. Thus, Eulerian scaling rates of convective precipitation are modulated by convective available potential energy and vertical wind shear, making it unlikely that present scaling rates can be applied to future climate conditions. Furthermore, we show that cells that cause heavy precipitation at fixed locations occur at low vertical wind shear and, thus, move relatively slowly compared to typical cells.
Background: Driven by globalization, urbanization and climate change, the distribution range of invasive vector species has expanded to previously colder ecoregions. To reduce health-threatening impacts on humans, insect vectors are extensively studied. Population genomics can reveal the genomic basis of adaptation and help to identify emerging trends of vector expansion.
Results: By applying whole genome analyses and genotype-environment associations to populations of the main dengue vector Ae. aegypti, sampled along an altitudinal temperature gradient in Nepal (200- 1300m), we identify adaptive traits and describe the species’ genomic footprint of climate adaptation to colder ecoregions. We found two clusters of differentiation with significantly different allele frequencies in genes associated to climate adaptation between the highland population (1300m) and all other lowland populations (≤ 800 m). We revealed non-synonymous mutations in 13 of the candidate genes associated to either altitude, precipitation or cold tolerance and identified an isolation-by-environment differentiation pattern.
Conclusion: Other than the expected gradual differentiation along the altitudinal gradient, our results reveal a distinct genomic differentiation of the highland population. This finding either indicates a differential invasion history to Nepal or local high-altitude adaptation explaining the population’s phenotypic cold tolerance. In any case, this highland population can be assumed to carry pre-adapted alleles relevant for the species’ invasion into colder ecoregions worldwide that way expanding their climate niche.
Seasonal forecasting systems still have difficulties predicting temperature over continental regions, while their performance is better over some maritime regions. On the other hand, the land surface is a substantial source of (sub-)seasonal predictability. A crucial land surface component in focus here is the snow cover, which stores water and modulates the surface radiation balance. This paper’s goal is to attribute snow cover seasonal forecasting biases and lack of skill to either initialization or parameterization errors. For this purpose, we compare the snow representation in five seasonal forecasting systems (from DWD, ECMWF, Météo-France, CMCC, and ECCC) and their performances in predicting snow and 2-m temperature over a Siberian region against ERA5 reanalysis and station data. Although all systems use similar atmospheric and land initialization approaches and data, their snow and temperature biases differ in sign and amplitude. Too-large initial snow biases persist over the forecast period, delaying and prolonging the melting phase. The simplest snow scheme (used in DWD’s system) shows too-early and fast melting in spring. However, systems including multi-layer snow schemes (Météo-France and CMCC) do not necessarily perform better. Both initialization and parameterization are causes of snow biases, but, depending on the system, one can be more dominant.