Boreal forests are an important source of trace gases and atmospheric aerosols, as well as a crucial carbon sink. As such, they form a strongly interconnected coupled system with the atmosphere. The SMEAR II station is located in a boreal Scots pine forest in Hyytiala, Finland, and has over 25 years of continuous measurements of atmospheric and ecosystem variables. In this study, we analyse the seasonal variations of trace gases, atmospheric aerosols, greenhouse gases, and meteorological variables, measured at the SMEAR II sta-tion during the past two and a half decades. Several ecosystem and atmospheric variables show seasonal correlations with each other, which suggests seasonal interactions within the climate system that links together ecosystem processes, greenhouse gases, trace gases and atmospheric aerosols. For instance, increased global radiation in summer increases air temperature and consequently affects the plant phenology, which promotes the ecosystem carbon exchange and biogenic volatile organic compound (biogenic VOC) release. This further affects the ambient concentrations of highly oxygenated organic molecules (HOMs) as well as the formation and growth of atmospheric organic aerosols. Organic aerosols subsequently influence aerosol optical properties and, through increased scattering, have the potential to cool the climate. We also discuss the impacts of the warm and dry summers of 2010 and 2018 on the studied variables. For these years, we find a higher-than-average ecosystem primary production especially in June, leading to an increased VOC flux from the forest. The increased VOC flux in turn leads to higher HOM and secondary aerosol concentration in the atmosphere. The latter increases light scattering by atmospheric aero-sol particles and thus leads to climate cooling. The results obtained in this study improve our understanding of how boreal forests respond to climate change.

In this article, we present temporally dynamic population distribution data from the Helsinki Metropolitan Area, Finland, at the level of 250 m by 250 m statistical grid cells. An hourly population distribution dataset is provided for regular workdays (Mon – Thu), Saturdays and Sundays. The data are based on aggregated mobile phone data collected by the biggest mobile network operator in Finland. Mobile phone data are assigned to statistical grid cells using an advanced dasymetric interpolation method based on ancillary data about land cover, buildings and a time use survey. The dataset is validated by comparing population register data from Statistics Finland for night hours and a daytime workplace registry. The resulting 24-hour population data can be used to reveal the temporal dynamics of the city, and examine population variations relevant to spatial accessibility analyses, crisis management, planning and beyond.

'Precision seismology' encompasses a set of methods which use differential measurements of time-delays to estimate the relative locations of earthquakes and explosions. Delay-times estimated from signal correlations often allow far more accurate estimates of one event location relative to another than is possible using classical hypocentre determination techniques. Many different algorithms and software implementations have been developed and different assumptions and procedures can often result in significant variability between different relative event location estimates. We present a Ground Truth (GT) dataset of 55 military surface explosions in northern Finland in 2007 that all took place within 300 m of each other. The explosions were recorded with a high signal-to-noise ratio to distances of about 2 degrees, and the exceptional waveform similarity between the signals from the different explosions allows for accurate correlation-based time-delay measurements. With exact coordinates for the explosions, we are able to assess the fidelity of relative location estimates made using any location algorithm or implementation. Applying double-difference calculations using two different 1-D velocity models for the region results in hypocentre-to-hypocentre distances which are too short and it is clear that the wavefield leaving the source region is more complicated than predicted by the models. Using the GT event coordinates, we are able to measure the slowness vectors associated with each outgoing ray from the source region. We demonstrate that, had such corrections been available, a significant improvement in the relative location estimates would have resulted. In practice we would of course need to solve for event hypocentres and slowness corrections simultaneously, and significant work will be needed to upgrade relative location algorithms to accommodate uncertainty in the form of the outgoing wavefield. We present this data set, together with GT coordinates, raw waveforms for all events on six regional stations, and tables of time-delay measurements, as a reference benchmark by which relative location algorithms and software can be evaluated.

Phase equilibria modeling is a powerful petrological tool to address both forward and inverse geological problems over a broad range of crustal and upper mantle conditions of pressure (P), temperature (T), composition (X), and redox (f(o2)). The development of thermodynamic databases, relatively realistic activity-composition (a-X) relations for solids, melts and fluids, pressure-volume-temperature (PVT) equations of state (EOS), and efficient numerical algorithms represent an inflection point in our ability to understand the nexus between tectonics and petrogenesis. While developed-and typically applied in isolation-by either metamorphic or igneous petrologists, some of the published thermodynamic models have overlapping P-T-X calibration ranges, which enables comparisons of model outcomes for similar conditions within the range of applicability. In this paper, we systematically compare the results of two such models that are routinely used for calculating phase equilibria in melt-bearing systems: rhyolite-MELTS (Gualda et al. 2012; Ghiorso and Gualda 2015) and the metabasite set of Green et al. (2016) using the thermodynamic database ds62 (Holland and Powell 2011) (hereafter denoted as "HPx-mb16"). We selected a N-MORB composition and modeled closed system equilibrium phase relations as a function of temperature at 0.25 and 1 GPa for N-MORB with 0.5 and 4 wt% H2O. Our results show that phase relations exhibit some key differences that, in some instances, impact geological inferences. For example, clinopyroxene and plagioclase stabilities are expanded to higher temperatures in HPx-mb16 compared to predictions from rhyolite-MELTS. Orthopyroxene and olivine are stable in greater proportions and at wider temperature ranges in rhyolite-MELTS compared to HPx-mb16. Importantly, HPx-mb16 predicts amphibole in all runs presented here, whereas amphibole is only predicted at high-P-high-H2O (1 GPa and 4 wt% H2O) in rhyolite-MELTS, and in lesser amounts. Garnet stability is systematically expanded at higher temperatures, and the proportion is greater in rhyolite-MELTS. In addition to phase assemblage differences, phase compositions may differ. For example, plagioclase anorthite content is systematically higher in HPx-mb16 (for the same set of conditions), whereas garnet Mg# is higher in rhyolite-MELTS. Calculated amphibole compositions are substantially different between the two models as well. Liquid compositions also show important differences. High-T liquids are generally similar in SiO2 contents but diverge at lower temperatures; in these cases, HIPx-mb16 liquids are SiO2-depleted compared to those produced by rhyolite-MELTS. Liquids are also systematically and substantially more mafic in HPx-mb16, and alumina and the alkali concentrations are relatively different and show different trends as a function of temperature at constant pressure. Overall, liquid compositions show the greatest differences near the solidus. Differences in modal abundances of phases and liquid compositions influence liquid trace-element signatures, and these differences can affect geological interpretations. Finally, a comparison between melting experiments of basaltic bulk composition and both thermodynamic models shows that rhyolite-MELTS better reproduces the higher temperature experiments, whereas HPx-mb16 better reproduces the lower temperature experiments. We discuss these and other similarities and differences to highlight the strengths and limitations of each model and to recognize that modeling results have important implications for interpretations of geologic processes. We recognize that our results are informed by a small subset of calculations over a limited range of conditions-but encourage further comparisons over a wider range of conditions and compositions.

Empirical evidence demonstrates that lakes and reservoirs are warming across the globe. Consequently, there is an increased need to project future changes in lake thermal structure and resulting changes in lake biogeochemistry in order to plan for the likely impacts. Previous studies of the impacts of climate change on lakes have often relied on a single model forced with limited scenario-driven projections of future climate for a relatively small number of lakes. As a result, our understanding of the effects of climate change on lakes is fragmentary, based on scattered studies using different data sources and modelling protocols, and mainly focused on individual lakes or lake regions. This has precluded identification of the main impacts of climate change on lakes at global and regional scales and has likely contributed to the lack of lake water quality considerations in policy-relevant documents, such as the Assessment Reports of the Intergovernmental Panel on Climate Change (IPCC). Here, we describe a simulation protocol developed by the Lake Sector of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) for simulating climate change impacts on lakes using an ensemble of lake models and climate change scenarios for ISIMIP phases 2 and 3. The protocol prescribes lake simulations driven by climate forcing from gridded observations and different Earth system models under various representative greenhouse gas concentration pathways (RCPs), all consistently bias-corrected on a 0.5 degrees x 0.5 degrees global grid. In ISIMIP phase 2, 11 lake models were forced with these data to project the thermal structure of 62 well-studied lakes where data were available for calibration under historical conditions, and using uncalibrated models for 17 500 lakes defined for all global grid cells containing lakes. In ISIMIP phase 3, this approach was expanded to consider more lakes, more models, and more processes. The ISIMIP Lake Sector is the largest international effort to project future water temperature, thermal structure, and ice phenology of lakes at local and global scales and paves the way for future simulations of the impacts of climate change on water quality and biogeochemistry in lakes.

The large-scale shifts in the salinity of the Baltic Sea over the Holocene are well understood and have been comprehensively documented using sedimentary proxy records. More recent work has focused on understanding how past salinity fluctuations have affected other ecological parameters (e.g. primary productivity, nutrient content) of the Baltic basin, and salinity changes over key events and over short time scales are still not well understood. The International Ocean Drilling Program Expedition 347 cored the Baltic basin in order to collect basin-wide environmental records through a glacial-interglacial cycle. Site M0059 is located in the Little Belt between the Baltic Sea and the Atlantic Ocean. A composite splice section from Site M0059 was analysed at a decadal resolution to study changes in salinity, nutrient conditions and other surface water column parameters based on changes in diatom assemblages and on quantitative diatom-based salinity inferences. A mesotrophic slightly brackish assemblage is seen in the lowermost analysed depths, corresponding to 7800-7500 cal. a BP. An increase in salinity and nutrient content of the water column leads into a meso-eutrophic brackish phase. The observed salinity increase is rapid, lasting from 7500 to 7150 cal. a BP. Subsequently, the Little Belt becomes oligotrophic and is dominated by tychopelagic diatoms from c. 7100 to c. 3900 cal. a BP. This interval contains some of the highest salinities observed followed by diatom assemblages similar to those of the Northern Atlantic Ocean, composed primarily of cosmopolitan open oceanmarine diatoms. A return to tychopelagic productivity is seen from 3850 to 980 cal. a BP. Anthropogenic eutrophication is detected in the last 300 years of the record, which intensifies in the uppermost sediments. These results represent the first decadally resolved record in the region and provide new insight into the transition to a brackish basin and subsequent ecological development.