Characterising and visualizing spatio-temporal patterns in hourly precipitation records

We develop new techniques to summarise and visualise spatial patterns of coincidence in weather events such as more or less heavy precipitation at a network of meteorological stations. The cosine similarity measure, which has a simple probabilistic interpretation for vectors of binary data, is generalised to characterise spatial dependencies of events that may reach different stations with a variable time lag. More specifically, we reduce such patterns into three parameters (dominant time lag, maximum cross-similarity, and window-maximum similarity) that can easily be computed for each pair of stations in a network. Furthermore, we visualise such three-parameter summaries by using colour-coded maps of dependencies to a given reference station and distance-decay plots for the entire network. Applications to hourly precipitation data from a network of 93 stations in Sweden illustrate how this method can be used to explore spatial patterns in the temporal synchrony of precipitation events.     

Sensitivity of North Atlantic subpolar gyre and overturning to stratification-dependent mixing: response to global warming

We use a reduced complexity climate model with a three-dimensional ocean component and realistic topography to investigate the effect of stratification-dependent mixing on the sensitivity of the North Atlantic subpolar gyre (SPG), and the Atlantic meridional overturning circulation (AMOC), to idealized CO₂ increase and peaking scenarios. The vertical diffusivity of the ocean interior is parameterized as κ ∼ N ^−α, where N is the local buoyancy frequency. For all parameter values 0 ≤ α ≤ 3, we find the SPG, and subsequently the AMOC, to weaken in response to increasing CO₂ concentrations. The weakening is significantly stronger for α ≥ α_cr ≈ 1.5. Depending on the value of α, two separate model states develop. These states remain different after the CO₂ concentration is stabilized, and in some cases even after the CO₂ concentration has been decreased again to the pre-industrial level. This behaviour is explained by a positive feedback between stratification and mixing anomalies in the Nordic Seas, causing a persistent weakening of the SPG.     

Available energy and energy balance closure at four coniferous forest sites across Europe

The available energy (AE), driving the turbulent fluxes of sensible heat and latent heat at the earth surface, was estimated at four partly complex coniferous forest sites across Europe (Tharandt, Germany; Ritten/Renon, Italy; Wetzstein, Germany; Norunda, Sweden). Existing data of net radiation were used as well as storage change rates calculated from temperature and humidity measurements to finally calculate the AE of all forest sites with uncertainty bounds. Data of the advection experiments MORE II (Tharandt) and ADVEX (Renon, Wetzstein, Norunda) served as the main basis. On-site data for referencing and cross-checking of the available energy were limited. Applied cross checks for net radiation (modelling, referencing to nearby stations and ratio of net radiation to global radiation) did not reveal relevant uncertainties. Heat storage of sensible heat J _H, latent heat J _E, heat storage of biomass J _veg and heat storage due to photosynthesis J _C were of minor importance during day but of some importance during night, where J _veg turned out to be the most important one. Comparisons of calculated storage terms (J _E, J _H) at different towers of one site showed good agreement indicating that storage change calculated at a single point is representative for the whole canopy at sites with moderate heterogeneity. The uncertainty in AE was assessed on the basis of literature values and the results of the applied cross checks for net radiation. The absolute mean uncertainty of AE was estimated to be between 41 and 52 W m^?2 (10–11 W m^?2 for the sum of the storage terms J and soil heat flux G) during mid-day (approximately 12% of AE). At night, the absolute mean uncertainty of AE varied from 20 to about 30 W m^?2 (approximately 6 W m^?2 for J plus G) resulting in large relative uncertainties as AE itself is small. An inspection of the energy balance showed an improvement of closure when storage terms were included and that the imbalance cannot be attributed to the uncertainties in AE alone.     

An evaluation of neutral and convective planetary boundary-layer parameterizations relative to large eddy simulations

This paper compares a number of one-dimensional closure models for the planetary boundary layer (PBL) that are currently in use in large-scale atmospheric models. Using the results of a large-eddy simulation (LES) model as the standard of comparison, the PBL models are evaluated over a range of stratifications from free convective to neutral and a range of surface shear stresses. Capping inversion strengths for the convective cases range from weakly to strongly capped. Six prototypical PBL models are evaluated in this study, which focuses on the accuracy of the boundary-layer fluxes of momentum, heat, and two passive scalars. One scalar mimics humidity and the other is a top-down scalar entrained into the boundary layer from above. A set of measures based on the layer-averaged differences of these fluxes from the LES solutions is developed. In addition to the methodological framework and suite of LES solutions, the main result of the evaluation is the recognition that all of the examined PBL parameterizations have difficulty reproducing the entrainment at the top of the PBL, as given by the LES, in most parameter regimes. Some of the PBL models are relatively accurate in their entrainment flux in a subset of parameter regimes. The sensitivity of the PBL models to vertical resolution is explored, and substantive differences are observed in the performance of the PBL models, relative to LES, at low resolution typical of large scale atmospheric models.     

Pine forest microclimate simulation using different diffusivities

Proper understanding of, e.g., evaporation from a forest requires an understanding of its microclimate. A well established, steady-state model was used to simulate microclimate and evaporation of a sparse pine forest in central Sweden. Model input included profiles of turbulent diffusivity, boundary-layer resistance, stomatal resistance, wind speed, net and global radiation and needle area density. Momentum balance, energy balance and exponentially decreasing diffusivities were used to study the sensitivity of the evaporation rates and of the temperature and humidity profiles. Model output proved to be unreliable when measured temperature and humidity at the bottom of the stand were used instead of a measured ground heat flux as the lower boundary condition. Energy balance diffusivity was usually larger than momentum balance diffusivity at the canopy top but decreased rapidly to a minimum at approximately the height where the momentum balance diffusivity had its maximum. Energy balance diffusivity commonly showed a secondary maximum below the height of the maximum needle area density. Profiles of Richardson number showed that thermal effects became important just below the canopy top. Bluff-body effects distinguished the energy balance from the momentum balance diffusivity and both were subject to shelter effects. Total evaporation was not very sensitive to the choice of diffusivity when soil heat flux was given as the lower boundary condition.     

Morphodynamic effects of riparian vegetation growth after stream restoration

The prediction of the morphological evolution of renaturalized streams is important for the success of restoration projects. Riparian vegetation is a key component of the riverine landscape and is therefore essential for the natural rehabilitation of rivers. This complicates the design of morphological interventions, since riparian vegetation is influenced by and influences the river dynamics. Morphodynamic models, useful tools for project planning, should therefore include the interaction between vegetation, water flow and sediment processes. Most restoration projects are carried out in USA and Europe, where rivers are highly intervened and where the climate is temperate and vegetation shows a clear seasonal cycle. Taking into account seasonal variations might therefore be relevant for the prediction of the river morphological adaptation. This study investigates the morphodynamic effects of riparian vegetation on a re‐meandered lowland stream in the Netherlands, the Lunterse Beek. The work includes the analysis of field data covering 5 years and numerical modelling. The results allow assessment of the performance of a modelling tool in predicting the morphological evolution of the stream and the relevance of including the seasonal variations of vegetation in the computations. After the establishment of herbaceous plants on its banks, the Lunterse Beek did not show any further changes in channel alignment. This is here attributed to the stabilizing effects of plant roots together with the small size of the stream. It is expected that the morphological restoration of similarly small streams may result in important initial morphological adaptation followed by negligible changes after full vegetation establishment. Copyright © 2018 John Wiley & Sons, Ltd.     

Observations on the association between mercury and organic matter dissolved in natural waters

Dissolved mercury in estuarine waters from the Mississippi Delta and Florida Everglades is associated with dissolved organic matter which has the properties of fulvic matter found in soils. Ultrafiltration of water samples demonstrated that mercury and dissolved organic carbon are selectively enriched in the < 500 molecular size cut-off fraction. A decrease in high molecular weight dissolved organic matter with increasing salinity in the Everglades exerts a partial control on the mercury content of these estuarine waters.     

Epithermal Gold Mineralization in the Trenggalek District,East Java,Indonesia

Gold‐mineralized quartz veins at the Trenggalek district of the Southern Mountains Range in East Java, Indonesia, are hosted by Oligo‐Miocene volcaniclastic and volcanic rocks, and are distributed close to andesitic plugs in the northern prospects (Dalangturu, Suruh, Jati, Gregah, Jombok, Salak, and Kojan) and the southern prospects (Sentul and Buluroto). The plugs are subalkaline tholeiitic basaltic‐andesite to calc‐alkaline andesite in composition. ⁴⁰Ar–³⁹Ar dating of a quartz‐adularia vein at the Dalangturu prospect yielded an age of 16.29 ± 0.56 Ma (2σ), and a crystal tuff of a limestone‐pyroclastic rock sequence at the southwest of the Dalangturu prospect was determined as 15.6 ± 0.5 Ma (2σ). Statistic overlap of ages suggests that the gold mineralization in the northern prospects took place in a shallow marine to subaerial transitional environment. Hydrothermal alteration of the host rocks is characterized by the replacement of quartz, illite and adularia. Quartz veins in surface outcrops are up to 50 cm wide in the northern prospects and up to 3 m wide in the southern prospects, showing a banded or brecciated texture, and are composed of quartz, adularia, carbonates with pyrite, electrum, sphalerite, galena, and polybasite. Gold contents of quartz veins are positively correlated with Ag, Zn, Pb, and Cu contents in both the northern and southern prospects. The quartz veins at the Jati, Gregah, and Sentul prospects have relatively lower gold‐silver ratios (Ag/Au = 23.2) compared to those at the Kojan, Dalangturu, Salak, and Suruh prospects (Ag/Au = 66.8). The quartz veins at the Dalangturu prospect are relatively rich in base metal sulfides. Ag/(Au+Ag) ratios of electrum in the Dalangturu prospect range from 45.2 to 65.0 at%, and FeS contents of sphalerite range from 1.2 to 6.4 mol%. Fluid inclusion microthermometry indicates ore‐forming temperatures of 190–200°C and 220–230°C at the Sentul and Kojan prospects, respectively. Widely variable vapor/liquid ratio of fluid inclusions indicates that fluid boiling took place within the hydrothermal system at the Sentul prospect. Salinities of ore‐fluids range from 0 to 0.7 wt% (av. 0.4 wt% NaCl equiv.) and from 0.5 to 1.4 wt% (av. 0.9 wt%) for the Sentul and Kojan prospects, respectively. The boiling of hydrothermal fluid was one of the gold deposition mechanisms in the Sentul prospect.     

Benthic Phosphorus Dynamics in the Gulf of Finland, Baltic Sea

Benthic fluxes of soluble reactive phosphorus (SRP) and dissolved inorganic carbon (DIC) were measured in situ using autonomous landers in the Gulf of Finland in the Baltic Sea, on four expeditions between 2002 and 2005. These measurements together with model estimates of bottom water oxygen conditions were used to compute the magnitude of the yearly integrated benthic SRP flux (also called internal phosphorus load). The yearly integrated benthic SRP flux was found to be almost 10 times larger than the external (river and land sources) phosphorus load. The average SRP flux was 1.25?±?0.56?mmol?m^?2?d^?1 on anoxic bottoms, and ?0.01?±?0.08?mmol?m^?2?d^?1 on oxic bottoms. The bottom water oxygen conditions determined whether the SRP flux was in a high or low regime, and degradation of organic matter (as estimated from benthic DIC fluxes) correlated positively with SRP fluxes on anoxic bottoms. From this correlation, we estimated a potential increase in phosphorus flux of 0.69?±?0.26?mmol?m^?2?d^?1 from presently oxic bottoms, if they would turn anoxic. An almost full annual data set of in situ bottom water oxygen measurements showed high variability of oxygen concentration. Because of this, an estimate of the time which the sediments were exposed to oxygenated overlying bottom water was computed using a coupled thermohydrodynamic ocean?Csea and ecosystem model. Total phosphorus burial rates were calculated from vertical profiles of total phosphorus in sediment and sediment accumulation rates. Recycling and burial efficiencies for phosphorus of 97 and 3%, respectively, were estimated for anoxic accumulation bottoms from a benthic mass balance, which was based on the measured effluxes and burial rates.