Colloidal meteorological processes in the formation of precipitation

This paper is the edited translation of the paper by Walter Findeisen “Die kolloidmeteorologischen Vorgänge bei der Niederschlagsbildung” (Colloidal meteorological processes in the formation of precipitation) that was published 1938 in the Meteorologische Zeitschrift 55, 121–133

Establishing the skill of climate field reconstruction techniques for precipitation with pseudoproxy experiments

This study aims at assessing the skill of several climate field reconstruction techniques (CFR) to reconstruct past precipitation over continental Europe and the Mediterranean at seasonal time scales over the last two millennia from proxy records. A number of pseudoproxy experiments are performed within the virtual reality ofa regional paleoclimate simulation at 45 km resolution to analyse different aspects of reconstruction skill. Canonical Correlation Analysis (CCA), two versions of an Analog Method (AM) and Bayesian hierarchical modeling (BHM) are applied to reconstruct precipitation from a synthetic network of pseudoproxies that are contaminated with various types of noise. The skill of the derived reconstructions is assessed through comparison with precipitation simulated by the regional climate model. Unlike BHM, CCA systematically underestimates the variance. The AM can be adjusted to overcome this shortcoming, presenting an intermediate behaviour between the two aforementioned techniques. However, a trade-off between reconstruction-target correlations and reconstructed variance is the drawback of all CFR techniques. CCA (BHM) presents the largest (lowest) skill in preserving the temporal evolution, whereas the AM can be tuned to reproduce better correlation at the expense of losing variance. While BHM has been shown to perform well for temperatures, it relies heavily on prescribed spatial correlation lengths. While this assumption is valid for temperature, it is hardly warranted for precipitation. In general, none of the methods outperforms the other. All experiments agree that a dense and regularly distributed proxy network is required to reconstruct precipitation accurately, reflecting its high spatial and temporal variability. This is especially true in summer, when a specifically short de-correlation distance from the proxy location is caused by localised summertime convective precipitation events.

A Parallel Solver for the Time-Periodic Navier–Stokes Equations

We investigate parallel algorithms for the solution of the Navier–Stokes equations in space-time. For periodic solutions, the discretized problem can be written as a large non-linear system of equations. This system of equations is solved by a Newton iteration. The Newton correction is computed using a preconditioned GMRES solver. The parallel performance of the algorithm is illustrated.

Synoptic drivers of 400 years of summer temperature and precipitation variability on Mt. Olympus, Greece

The Mediterranean region has been identified as a global warming hotspot, where future climate impacts are expected to have significant consequences on societal and ecosystem well-being. To put ongoing trends of summer climate into the context of past natural variability, we reconstructed climate from maximum latewood density (MXD) measurements of Pinus heldreichii (1521–2010) and latewood width (LWW) of Pinus nigra (1617–2010) on Mt. Olympus, Greece. Previous research in the northeastern Mediterranean has primarily focused on inter-annual variability, omitting any low-frequency trends. The present study utilizes methods capable of retaining climatically driven long-term behavior of tree growth. The LWW chronology corresponds closely to early summer moisture variability (May–July, r = 0.65, p < 0.001, 1950–2010), whereas the MXD-chronology relates mainly to late summer warmth (July–September, r = 0.64, p < 0.001; 1899–2010). The chronologies show opposing patterns of decadal variability over the twentieth century (r = −0.68, p < 0.001) and confirm the importance of the summer North Atlantic Oscillation (sNAO) for summer climate in the northeastern Mediterranean, with positive sNAO phases inducing cold anomalies and enhanced cloudiness and precipitation. The combined reconstructions document the late twentieth—early twenty-first century warming and drying trend, but indicate generally drier early summer and cooler late summer conditions in the period ~1700–1900 CE. Our findings suggest a potential decoupling between twentieth century atmospheric circulation patterns and pre-industrial climate variability. Furthermore, the range of natural climate variability stretches beyond summer moisture availabilityobserved in recent decades and thus lends credibility to the significant drying trends projected for this region in current Earth System Model simulations.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

Climatology of Vb cyclones, physical mechanisms and their impact on extreme precipitation over Central Europe

Cyclones, which develop over the western Mediterranean and move northeastward are a major source of extreme weather and known to be responsible for heavy precipitation over the northern side of the Alpine range and Central Europe. As the relevant processes triggering these so-called Vb events and their impact on extreme precipitation are not yet fully understood, this study focuses on gaining insight into the dynamics of past events. For this, a cyclone detection and tracking tool is applied to the ERA-Interim reanalysis (1979–2013) to identify prominent Vb situations. Precipitation in the ERA-Interim and the E-OBS data sets is used to evaluate case-to-case precipitation amounts and to assess consistency between the two data sets. Both data sets exhibit high variability in precipitation amounts among different Vb events. While only 23 % of all Vb events are associated with extreme precipitation, around 15 % of all extreme precipitation days (99 percentile) over the northern Alpine region and Central Europe are induced by Vb events, although Vb cyclones are rare events (2.3 per year). To obtain a better understanding of the variability within Vb events, the analysis of the 10 heaviest and lowest precipitation Vb events reveals noticeable differences in the state of the atmosphere. These differences are most pronounced in the geopotential height and potential vorticity field, indicating a much stronger cyclone for heavy precipitation events. The related differences in wind direction are responsible for the moisture transport around the Alps and the orographical lifting along the northern slopes of the Alps. These effects are the main reasons for a disastrous outcome of Vb events, and consequently are absent in the Vb events associated with low precipitation. Hence, our results point out that heavy precipitation related to Vb events is mainly related to large-scale dynamics rather than to thermodynamic processes.

Temperature and precipitation effects on δ13C depth profiles in SOM under temperate beech forests

Enrichment of 13C in SOM with soil depth is related to interacting processes influenced by temperature and precipitation. Our objectives were to derive climate effects on patterns of vertical δ13C values of soil organic matter (SOM) while minimizing the effect of confounding variables. We investigated vertical changes in δ13C values of SOM in 1-cm depth intervals in silvicultural mature beech (Fagus sylvatica L.) forest ecosystems in northern Rhineland-Palatinate across gradients of MAT (7.9 to 9.7 °C mean annual temperature) and MAP (607 to 1085 mm mean annual precipitation) in winter 2011. Forest stands (n = 10) were chosen based on data sets provided by the Rhineland-Palatinate Forest Administration so that variations in these gradients occurred while other environmental factors like physico-chemical soil properties, tree species, stand age, exposition and precipitation (for the temperature gradient) or temperature (for the precipitation gradient) did not differ among study sites. From litter down to the mineral soil at 10 cm depth, soil organic carbon (SOC) content decreased (47.5 ± SE 0.1% to 2.5 ± 0.1%) while the δ13C values increased (− 29.4 ± 0.1‰ to − 26.1 ± 0.1‰). Litter of sites under higher MAP/lower MAT had lower δ13C values which was in line with literature data on climate driven plant physiological process. To compare the dimension of the vertical 13C enrichment, δ13C values were regressed linearly against log-transformed carbon contents yielding absolute values of these slopes (beta). Beta values ranged between 0.6 and 4.5 (range of r from − 0.7 to − 1.0; p < 0.01). Due to an assumed decay continuum and similar variations of δ13C values in litter and in 10 cm depth, we conclude that effects on isotope composition in the Oi layer continue vertically and therefore, δ13C values in litter do not solely control beta values. Beta values decreased with increasing MAT (r = − 0.83; p < 0.05). Reduced soil moisture and therefore both, reduced microbial activity and reduced downward transport of microbial cycled DOM (=13C enriched) might be responsible for less pronounced δ13C depth profiles in case of high temperatures. Greater C:N ratios (lower degradability) of the litter under higher temperatures likely contributed to these depth trends. Beta values increased with increasing MAP (r = 0.73; p < 0.05). We found decreasing C:N ratios in the mineral soil that possibly indicates higher decomposition under higher precipitation. Exclusion of the organic layers from linear regressions indicated a stronger impact of MAP on the development of δ13C depth profiles. Our results confirm temperature and precipitation effects on δ13C depth profiles and indicate stronger 13C enrichment under lower MAT/higher MAP. Therefore, time series of vertical δ13C depth profiles might provide insights into climate change effects.

Flood triggering in Switzerland: the role of daily to monthly preceding precipitation

Determining the role of different precipitation periods for peak discharge generation is crucial for both projecting future changes in flood probability and for short- and medium-range flood forecasting. In this study, catchment-averaged daily precipitation time series are analyzed prior to annual peak discharge events (floods) in Switzerland. The high number of floods considered – more than 4000 events from 101 catchments have been analyzed – allows to derive significant information about the role of antecedent precipitation for peak discharge generation. Based on the analysis of precipitation times series, a new separation of flood-related precipitation periods is proposed: (i) the period 0 to 1 day before flood days, when the maximum flood-triggering precipitation rates are generally observed, (ii) the period 2 to 3 days before flood days, when longer-lasting synoptic situations generate "significantly higher than normal" precipitation amounts, and (iii) the period from 4 days to 1 month before flood days when previous wet episodes may have already preconditioned the catchment. The novelty of this study lies in the separation of antecedent precipitation into the precursor antecedent precipitation (4 days before floods or earlier, called PRE-AP) and the short range precipitation (0 to 3 days before floods, a period when precipitation is often driven by one persistent weather situation like e.g., a stationary low-pressure system). A precise separation of "antecedent" and "peak-triggering" precipitation is not attempted. Instead, the strict definition of antecedent precipitation periods permits a direct comparison of all catchments. The precipitation accumulating 0 to 3 days before an event is the most relevant for floods in Switzerland. PRE-AP precipitation has only a weak and region-specific influence on flood probability. Floods were significantly more frequent after wet PRE-AP periods only in the Jura Mountains, in the western and eastern Swiss plateau, and at the outlet of large lakes. As a general rule, wet PRE-AP periods enhance the flood probability in catchments with gentle topography, high infiltration rates, and large storage capacity (karstic cavities, deep soils, large reservoirs). In contrast, floods were significantly less frequent after wet PRE-AP periods in glacial catchments because of reduced melt. For the majority of catchments however, no significant correlation between precipitation amounts and flood occurrences is found when the last 3 days before floods are omitted in the precipitation amounts. Moreover, the PRE-AP was not higher for extreme floods than for annual floods with a high frequency and was very close to climatology for all floods. The fact that floods are not significantly more frequent nor more intense after wet PRE-AP is a clear indicator of a short discharge memory of Pre-Alpine, Alpine and South Alpine Swiss catchments. Our study poses the question whether the impact of long-term precursory precipitation for floods in such catchments is not overestimated in the general perception. The results suggest that the consideration of a 3–4 days precipitation period should be sufficient to represent (understand, reconstruct, model, project) Swiss Alpine floods.

Triple isotope (δD, δ17O, δ18O) study on precipitation, drip water and speleothem fluid inclusions for a Western Central European cave (NW Switzerland)

Deuterium (δD) and oxygen (δ18O) isotopes are powerful tracers of the hydrological cycle and have been extensively used for paleoclimate reconstructions as they can provide information on past precipitation, temperature and atmospheric circulation. More recently, the use of δ17O excess derived from precise measurement of δ17O and δ18O gives new and additional insights in tracing the hydrological cycle whereas uncertainties surround this proxy. However, 17O excess could provide additional information on the atmospheric conditions at the moisture source as well as about fractionations associated with transport and site processes. In this paper we trace water stable isotopes (δD,δ17O and δ18O) along their path from precipitation to cave drip water and finally to speleothem fluid inclusions for Milandre cave in northwestern Switzerland. A two year-long daily resolved precipitation isotope record close to the cave site is compared to collected cave drip water (3 months average resolution) and fluid inclusions of modern and Holocene stalagmites. Amount weighted mean δD,δ18O and δ17O are -71.0‰, -9.9‰, -5.2‰ for precipitation, -60.3‰, -8.7‰, -4.6‰ for cave drip water and -61.3‰, -8.3‰, -4.7‰ for recent fluid inclusions respectively. Second order parameters have also been derived in precipitation and drip water and present similar values with 18 per meg for 17O excess whereas d-excess is 1.5‰ more negative in drip water. Furthermore, the atmospheric signal is shifted towards enriched values in the drip water and fluid inclusions (Δ of ~ + 10‰ for δD). The isotopic composition of cave drip water exhibits a weak seasonal signal which is shifted by around 8 - 10 months (groundwater residence time) when compared to the precipitation. Moreover, we carried out the first δ17O measurement in speleothem fluid inclusions, as well as the first comparison of the δ17 O behaviour from the meteoric water to the fluid inclusions entrapment in speleothems. This study on precipitation, drip water and fluid inclusions will be used as a speleothem proxy calibration for Milandre cave in order to reconstruct paleotemperatures and moisture source variations for Western Central Europe.

Reconciling Precipitation with Runoff: Observed Hydrological Change in the Midlatitudes

Century-long observed gridded land precipitation datasets are a cornerstone of hydrometeorological research. But recent work has suggested that observed Northern Hemisphere midlatitude (NHML) land mean precipitation does not show evidence of an expected negative response to mid-twentieth-century aerosol forcing. Utilizing observed river discharges, the observed runoff is calculated and compared with observed land precipitation. The results show a near-zero twentieth-century trendinobserved NHML landmean runoff,in contrast to the significant positive trend in observed NHML land mean precipitation. However, precipitation and runoff share common interannual and decadal variability. An obvious split, or breakpoint, is found in the NHML land mean runoff–precipitation relationship in the 1930s. Using runoff simulated by six land surface models (LSMs), which are driven by the observed precipitation dataset, such breakpoints are absent. These findings support previous hypotheses that inhomogeneities exist in the early-twentieth-century NHML land mean precipitation record. Adjusting the observed precipitation record according to the observed runoff record largely accounts for the departure of the observed precipitation response from that predicted given the real-world aerosol forcing estimate, more than halving the discrepancy from about 6 to around 2 W m 22. Consideration of complementary observed runoff adds support to the suggestion that NHML-wide early-twentieth-century precipitation observations are unsuitable for climate change studies. The agreement between precipitation and runoff over Europe, however, is excellent, supporting the use of whole-twentieth-century observed precipitation datasets here.