ALBERT

Description: The concept of cloud radiative forcing (CRF) is commonly applied to quantify the impact of clouds on the surface radiative energy budget (REB). In the Arctic, specific radiative interactions between microphysical and macrophysical properties of clouds and the surface strongly modify the warming or cooling effect of clouds, complicating the estimate of CRF obtained from observations or models. Clouds tend to increase the broadband surface albedo over snow or sea ice surfaces compared to cloud-free conditions. However, this effect is not adequately considered in the derivation of CRF in the Arctic so far. Therefore, we have quantified the effects caused by surface-albedo–cloud interactions over highly reflective snow or sea ice surfaces on the CRF using radiative transfer simulations and below-cloud airborne observations above the heterogeneous springtime marginal sea ice zone (MIZ) during the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) campaign. The impact of a modified surface albedo in the presence of clouds, as compared to cloud-free conditions, and its dependence on cloud optical thickness is found to be relevant for the estimation of the shortwave CRF. A method is proposed to consider this surface albedo effect on CRF estimates by continuously retrieving the cloud-free surface albedo from observations under cloudy conditions, using an available snow and ice albedo parameterization. Using ACLOUD data reveals that the estimated average shortwave cooling by clouds almost doubles over snow- and ice-covered surfaces (−62 W m−2 instead of −32 W m−2), if surface-albedo–cloud interactions are considered. As a result, the observed total (shortwave plus longwave) CRF shifted from a warming effect to an almost neutral one. Concerning the seasonal cycle of the surface albedo, it is demonstrated that this effect enhances shortwave cooling in periods when snow dominates the surface and potentially weakens the cooling by optically thin clouds during the summertime melting season. These findings suggest that the surface-albedo–cloud interaction should be considered in global climate models and in long-term studies to obtain a realistic estimate of the shortwave CRF to quantify the role of clouds in Arctic amplification.

Description: The evaluation of potential landslides in mountain areas is a very complex process. Currently, event understanding is scarce due to information limitations. Identifying the whole chain of events is not a straightforward task, and the impacts of mass-wasting processes depend on the conditions downstream of the origin. In this paper, we present an example that illustrates the complexities in the evaluation of the chain of events that may lead to a natural disaster. On 16 December 2017, a landslide occurred in the Yelcho mountain range (southern Chile). In that event, 7 million m3 of rocks and soil fell on the Yelcho glacier, depositing 2 million m3 on the glacier terminal, and the rest continued downstream, triggering a mudflow that hit Villa Santa Lucía in Chilean Patagonia and killing 22 people. The complex event was anticipated in the region by the National Geological and Mining Survey (Sernageomin in Spanish). However, the effects of the terrain characteristics along the run-out area were more significant than anticipated. In this work, we evaluate the conditions that enabled the mudflow that hit Villa Santa Lucía. We used the information generated by Sernageomin's professionals after the mudflow. We carried out geotechnical tests to characterize the soil. We simulated the mudflow using two hydrodynamic programs (r.avaflow and Flo-2D) that can handle the rheology of the water–soil mixture. Our results indicate that the soil is classified as volcanic pumices. This type of soil can be susceptible to the collapse of the structure when subjected to shearing (molding), flowing as a viscous liquid. From the numerical modeling, we concluded that r.avaflow performs better than Flo-2D. The mudflow was satisfactorily simulated using a water content in the mixture ranging from 30 % to 40 %, which would have required a source of about 3 million m3 of water. Coupling the simulations and the soil tests that we performed, we estimated that in the area scoured by the mudflow, there were probably around 2 800 000 m3 of water within the soil. Therefore, the conditions of the valley were crucial to enhancing the impacts of the landslide. This result is relevant because it highlights the importance of evaluating the complete chain of events to map hazards. We suggest that in future hazard mapping, geotechnical studies in combination with hydrodynamic simulation should be included, in particular when human lives are at risk.

Description: Production and reduction of nitrous oxide (N2O) by soil denitrifiers influence atmospheric concentrations of this potent greenhouse gas. Accurate projections of the net N2O flux have three key uncertainties: (1) short- vs. long-term responses to warming, (2) interactions among soil horizons, and (3) temperature responses of different steps in the denitrification pathway. We addressed these uncertainties by sampling soil from a boreal forest climate transect encompassing a 5.2 ∘C difference in the mean annual temperature and incubating the soil horizons in isolation and together at three ecologically relevant temperatures in conditions that promote denitrification. Both short-term exposure to warmer temperatures and long-term exposure to a warmer climate increased N2O emissions from organic and mineral soils; an isotopic tracer suggested that an increase in N2O production was more important than a decline in N2O reduction. Short-term warming promoted the reduction of organic horizon-derived N2O by mineral soil when these horizons were incubated together. The abundance of nirS (a precursor gene for N2O production) was not sensitive to temperature, whereas that of nosZ clade I (a gene for N2O reduction) decreased with short-term warming in both horizons and was higher from a warmer climate. These results suggest a decoupling of gene abundance and process rates in these soils that differs across horizons and timescales. In spite of these variations, our results suggest a consistent, positive response of denitrifier-mediated net N2O efflux rates to temperature across timescales in these boreal forests. Our work also highlights the importance of understanding cross-horizon N2O fluxes for developing a predictive understanding of net N2O efflux from soils.

Description: Wetlands are the largest and most uncertain natural sources of atmospheric methane (CH4). Several process-based models have been developed to quantify the magnitude and estimate spatial and temporal variations in CH4 emissions from global wetlands. Reliable models are required to estimate global wetland CH4 emissions. This study aimed to test two process-based models, CH4MODwetland and Terrestrial Ecosystem Model (TEM), against the CH4 flux measurements of marsh, swamp, peatland and coastal wetland sites across the world; specifically, model accuracy and generality were evaluated for different wetland types and in different continents, and then the global CH4 emissions from 2000 to 2010 were estimated. Both models showed similar high correlations with the observed seasonal/annual total CH4 emissions, and the regression of the observed versus computed total seasonal/annual CH4 emissions resulted in R2 values of 0.81 and 0.68 for CH4MODwetland and TEM, respectively. The CH4MODwetland produced accurate predictions for marshes, peatlands, swamps and coastal wetlands, with model efficiency (EF) values of 0.22, 0.52, 0.13 and 0.72, respectively. TEM produced good predictions for peatlands and swamps, with EF values of 0.69 and 0.74, respectively, but it could not accurately simulate marshes and coastal wetlands (EF 

Description: As knowledge about the cirrus clouds in the lower stratosphere is limited, reliable long-term measurements are needed to assess their characteristics, radiative impact and important role in upper troposphere and lower stratosphere (UTLS) chemistry. We used 6 years (2006–2012) of Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) measurements to investigate the global and seasonal distribution of stratospheric cirrus clouds and compared the MIPAS results with results derived from the latest version (V4.x) of the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) data. For the identification of stratospheric cirrus clouds, precise information on both the cloud top height (CTH) and the tropopause height is crucial. Here, we used lapse rate tropopause heights estimated from the ERA-Interim global reanalysis. Considering the uncertainties of the tropopause heights and the vertical sampling grid, we define CTHs more than 0.5 km above the tropopause as stratospheric for CALIPSO data. For MIPAS data, we took into account the coarser vertical sampling grid and the broad field of view so that we considered cirrus CTHs detected more than 0.75 km above the tropopause as stratospheric. Further sensitivity tests were conducted to rule out sampling artefacts in MIPAS data. The global distribution of stratospheric cirrus clouds was derived from night-time measurements because of the higher detection sensitivity of CALIPSO. In both data sets, MIPAS and CALIPSO, the stratospheric cirrus cloud occurrence frequencies are significantly higher in the tropics than in the extra-tropics. Tropical hotspots of stratospheric cirrus clouds associated with deep convection are located over equatorial Africa, South and Southeast Asia, the western Pacific, and South America. Stratospheric cirrus clouds were more often detected in December–February (15 %) than June–August (8 %) in the tropics (±20∘). At northern and southern middle latitudes (40–60∘), MIPAS observed about twice as many stratospheric cirrus clouds (occurrence frequencies of 4 %–5 % for MIPAS rather than about 2 % for CALIPSO). We attribute more frequent observations of stratospheric cirrus clouds with MIPAS to the higher detection sensitivity of the instrument to optically thin clouds. In contrast to the difference between daytime and night-time occurrence frequencies of stratospheric cirrus clouds by a factor of about 2 in zonal means in the tropics (4 % and 10 %, respectively) and at middle latitudes for CALIPSO data, there is little diurnal cycle in MIPAS data, in which the difference of occurrence frequencies in the tropics is about 1 percentage point in zonal mean and about 0.5 percentage point at middle latitudes. The difference between CALIPSO day and night measurements can also be attributed to their differences in detection sensitivity. Future work should focus on better understanding the origin of the stratospheric cirrus clouds and their impact on radiative forcing and climate.

Description: Air quality on our planet has been changing in particular since the industrial revolution (1750s) because of anthropogenic emissions. It is becoming increasingly important to realize air cleanliness, since clean air is as valuable a resource as clean water. A global standard for quantifying the level of air cleanliness is swiftly required, and we defined a novel concept, namely the Clean aIr Index (CII). The CII is a simple index defined by the normalization of the amount of a set of individual air pollutants. A CII value of 1 indicates completely clean air (no air pollutants), and 0 indicates the presence of air pollutants that meet the numerical environmental criteria for the normalization. In this time, the air pollutants used in the CII were taken from the Air Quality Guidelines (AQG) set by the World Health Organization (WHO), namely O3, particulate matters, NO2, and SO2. We chose Japan as a study area to evaluate CII because of the following reasons: (i) accurate validation data, as the in situ observation sites of the Atmospheric Environmental Regional Observation System (AEROS) provide highly accurate values of air pollutant amounts, and (ii) fixed numerical criteria from the Japanese Environmental Quality Standards (JEQS) as directed by the Ministry of the Environment (MOE) of Japan. We quantified air cleanliness in terms of the CII for the all 1896 municipalities in Japan and used data from Seoul and Beijing to evaluate Japanese air cleanliness. The amount of each air pollutant was calculated using a model that combined the Weather Research and Forecasting (WRF) and Community Multiscale Air Quality (CMAQ) models for 1 April 2014 to 31 March 2017. The CII values calculated by the WRF–CMAQ model and the AEROS measurements showed good agreement. The mean of the correlation coefficient for the CII values of 498 municipalities where the AEROS measurements operated was 0.66±0.05, which was higher than that of the Air Quality Index (AQI) of 0.57±0.06. The CII values averaged for the study period were 0.67, 0.52, and 0.24 in Tokyo, Seoul, and Beijing, respectively; thus, the air in Tokyo was 1.5 and 2.3 times cleaner, i.e., lower amounts of air pollutants, than the air in Seoul and Beijing, respectively. The average CII value for the all Japanese municipalities was 0.72 over the study period. The extremely clean air, CII ≈0.90, occurred in the southern remote islands of Tokyo and to the west of the Pacific coast, i.e., Kochi, Mie, and Wakayama prefectures during summer, with the transport of clean air from the ocean. We presented the top 100 clean air cities in Japan as one example of an application of CII in society. We confirmed that the CII enabled the quantitative evaluation of air cleanliness. The CII can be useful and valuable in various scenarios such as encouraging sightseeing and migration, investment and insurance business, and city planning. The CII is a simple and fair index that can be applied to all nations.

Description: Calculating a multi-model mean, a commonly used method for ensemble averaging, assumes model independence and equal model skill. Sharing of model components amongst families of models and research centres, conflated by growing ensemble size, means model independence cannot be assumed and is hard to quantify. We present a methodology to produce a weighted-model ensemble projection, accounting for model performance and model independence. Model weights are calculated by comparing model hindcasts to a selection of metrics chosen for their physical relevance to the process or phenomena of interest. This weighting methodology is applied to the Chemistry–Climate Model Initiative (CCMI) ensemble to investigate Antarctic ozone depletion and subsequent recovery. The weighted mean projects an ozone recovery to 1980 levels, by 2056 with a 95 % confidence interval (2052–2060), 4 years earlier than the most recent study. Perfect-model testing and out-of-sample testing validate the results and show a greater projective skill than a standard multi-model mean. Interestingly, the construction of a weighted mean also provides insight into model performance and dependence between the models. This weighting methodology is robust to both model and metric choices and therefore has potential applications throughout the climate and chemistry–climate modelling communities.

Description: Connections between vegetation and soil thermal dynamics are critical for estimating the vulnerability of permafrost to thaw with continued climate warming and vegetation changes. The interplay of complex biophysical processes results in a highly heterogeneous soil temperature distribution on small spatial scales. Moreover, the link between topsoil temperature and active layer thickness remains poorly constrained. Sixty-eight temperature loggers were installed at 1–3 cm depth to record the distribution of topsoil temperatures at the Trail Valley Creek study site in the northwestern Canadian Arctic. The measurements were distributed across six different vegetation types characteristic for this landscape. Two years of topsoil temperature data were analysed statistically to identify temporal and spatial characteristics and their relationship to vegetation, snow cover, and active layer thickness. The mean annual topsoil temperature varied between −3.7 and 0.1 ∘C within 0.5 km2. The observed variation can, to a large degree, be explained by variation in snow cover. Differences in snow depth are strongly related with vegetation type and show complex associations with late-summer thaw depth. While cold winter soil temperature is associated with deep active layers in the following summer for lichen and dwarf shrub tundra, we observed the opposite beneath tall shrubs and tussocks. In contrast to winter observations, summer topsoil temperature is similar below all vegetation types with an average summer topsoil temperature difference of less than 1 ∘C. Moreover, there is no significant relationship between summer soil temperature or cumulative positive degree days and active layer thickness. Altogether, our results demonstrate the high spatial variability of topsoil temperature and active layer thickness even within specific vegetation types. Given that vegetation type defines the direction of the relationship between topsoil temperature and active layer thickness in winter and summer, estimates of permafrost vulnerability based on remote sensing or model results will need to incorporate complex local feedback mechanisms of vegetation change and permafrost thaw.

Description: Evapotranspiration (ET) from tropical forests serves as a critical moisture source for regional and global climate cycles. However, the magnitude, seasonality, and interannual variability of ET in the Congo Basin remain poorly constrained due to a scarcity of direct observations, despite the Congo being the second-largest river basin in the world and containing a vast region of tropical forest. In this study, we applied a water balance model to an array of remotely sensed and in situ datasets to produce monthly, basin-wide ET estimates spanning April 2002 to November 2016. Data sources include water storage changes estimated from the Gravity Recovery and Climate Experiment (GRACE) satellites, in situ measurements of river discharge, and precipitation from several remotely sensed and gauge-based sources. An optimal precipitation dataset was determined as a weighted average of interpolated data by Nicholson et al. (2018), Climate Hazards InfraRed Precipitation with Station data version 2 (CHIRPS2) , and the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks–Climate Data Record product (PERSIANN-CDR), with the relative weights based on the error magnitudes of each dataset as determined by triple collocation. The resulting water-balance-derived ET (ETwb) features a long-term average that is consistent with previous studies (117.2±3.5 cm yr−1) but displays greater seasonal and interannual variability than seven global ET products. The seasonal cycle of ETwb generally tracks that of precipitation over the basin, with the exception that ETwb is greater in March–April–May (MAM) than in the relatively wetter September–October–November (SON) periods. This pattern appears to be driven by seasonal variations in the diffuse photosynthetically active radiation (PAR) fraction, net radiation (Rn), and soil water availability. From 2002 to 2016, Rn, PAR, and vapor-pressure deficit (VPD) all increased significantly within the Congo Basin; however, no corresponding trend occurred in ETwb. We hypothesize that the stability of ETwb over the study period despite sunnier and less humid conditions may be due to increasing atmospheric CO2 concentrations that offset the impacts of rising VPD and irradiance on stomatal water use efficiency (WUE).

Description: A range of future climate scenarios are projected for high atmospheric CO2 concentrations, given uncertainties over future human actions as well as potential environmental and climatic feedbacks. The geological record offers an opportunity to understand climate system response to a range of forcings and feedbacks which operate over multiple temporal and spatial scales. Here, we examine a single interglacial during the late Pliocene (KM5c, ca. 3.205±0.01 Ma) when atmospheric CO2 exceeded pre-industrial concentrations, but were similar to today and to the lowest emission scenarios for this century. As orbital forcing and continental configurations were almost identical to today, we are able to focus on equilibrium climate system response to modern and near-future CO2. Using proxy data from 32 sites, we demonstrate that global mean sea-surface temperatures were warmer than pre-industrial values, by ∼2.3 ∘C for the combined proxy data (foraminifera Mg∕Ca and alkenones), or by ∼3.2–3.4 ∘C (alkenones only). Compared to the pre-industrial period, reduced meridional gradients and enhanced warming in the North Atlantic are consistently reconstructed. There is broad agreement between data and models at the global scale, with regional differences reflecting ocean circulation and/or proxy signals. An uneven distribution of proxy data in time and space does, however, add uncertainty to our anomaly calculations. The reconstructed global mean sea-surface temperature anomaly for KM5c is warmer than all but three of the PlioMIP2 model outputs, and the reconstructed North Atlantic data tend to align with the warmest KM5c model values. Our results demonstrate that even under low-CO2 emission scenarios, surface ocean warming may be expected to exceed model projections and will be accentuated in the higher latitudes.

Description: The meridional heat transport (MHT) in the North Atlantic is critically important to climate variability and the global overturning circulation. A wide range of ocean processes contribute to North Atlantic MHT, ranging from basin-scale overturning and gyre motions to mesoscale instabilities (such as eddies). However, previous analyses of “eddy” MHT in the region have mostly focused on the contributions of time-variable velocity and temperature, rather than considering the association of MHT with distinct spatial scales within the basin. In this study, a zonal spatial-scale decomposition separates large-scale from mesoscale velocity and temperature contributions to MHT, in order to characterize the physical processes driving MHT. Using this approach, we found that the mesoscale contributions to the time-mean and interannual/decadal (ID) variability of MHT in the latitude range 39–45∘ N are larger than large-scale horizontal contributions, though smaller than the overturning contributions. Considering the 40∘ N transect as a case study, large-scale ID variability is mostly generated close to the western boundary. In contrast, most ID MHT variability associated with mesoscales originates in two distinct regions: a western boundary region (70–60∘ W) associated with 1- to 4-year interannual variations and an interior region (50–35∘ W) associated with decadal variations. Surface eddy kinetic energy is not a reliable indicator of high MHT episodes, but the large-scale meridional temperature gradient is an important factor, by influencing the local temperature variance as well as the local correlation of velocity and temperature. Most of the mesoscale contribution to MHT at 40∘ N is associated with transient and propagating processes, but stationary mesoscale structures explain most of the mesoscale MHT south of the Gulf Stream separation, highlighting the differences between the temporal and spatial decomposition of meridional temperature fluxes.

Description: The aerosol size distribution and cloud condensation nuclei (CCN) number concentration were measured using a wide-range particle spectrometer (WPS) and a cloud condensation nuclei counter (CCNC) on Mt. Tian from 31 July to 9 September, 2019. Combined with meteorological data, distribution characteristics of aerosol size and CCN and their influencing factors were analyzed. The results indicated that the mean aerosol number concentration was 5475.6 ± 5636.5 cm−3. The mean CCN concentrations were 183.7 ± 114.5 cm−3, 729.8 ± 376.1 cm−3, 1630.5 ± 980.5 cm−3, 2162.5 ± 1345.3 cm−3, and 2575.7 ± 1632.9 cm−3 at supersaturation levels of 0.1%, 0.2%, 0.4%, 0.6%, and 0.8%, respectively. The aerosol number size distribution is unimodal, and the dominant particle size is 30–60 nm. Affected by the height of the boundary layer and the valley wind, the diurnal variation in aerosol number concentration shows a unimodal distribution with a peak at 17:00, and the CCN number concentration showed a bimodal distribution with peaks at 18:00 and 21:00. The particle size distribution and supersaturation have a major impact on the activation of the aerosol into CCN. At 0.1% supersaturation (S), the 300–500 nm particles are most likely to activate to CCN. Particles of 100–300 nm are most easily activated at 0.2% (S), while particles of 60–80 nm are most likely activated at high supersaturation (≥0.4%). The concentrations of aerosol and CCN are higher in the northerly wind. Ambient relative humidity (RH) has little relationship with the aerosol activation under high supersaturation. According to N = CSk fitting the CCN spectrum, C = 3297 and k = 0.90 on Mt. Tian, characteristic of the clean continental type.

Description: The atmospheric boundary layer height is important for constraining air pollution and meteorological models. This study attempted to validate the MODIS-estimated atmospheric boundary layer height (ABLH), and variation in the ABLH in Uganda was evaluated. The ABLH was estimated from MODIS data using the mixing ratio profile gradient method and compared to the ABLH estimated from radiosonde data using three different methods. Unlike in studies in other regions of the world, correlations between ABLH estimated using MODIS and radiosonde data were weak, implying limited usefulness of MODIS data for determining ABLH. However, the diurnal variation in MODIS-derived ABLH and particulate matter (PM10) was consistent with the expected inverse relationship between PM10 mass concentration and ABLH, and the mean MODIS-derived ABLH values were significantly lower during wet seasons than dry seasons, as expected.

Description: Accurate reference spectroscopic information for the water molecule from the microwave to the near-ultraviolet is of paramount importance in atmospheric research. A semi-empirical potential energy surface for the ground electronic state of H216O has been created by refining almost 4000 experimentally determined energy levels. These states extend into regions with large values of rotational and vibrational excitation. For all states considered in our refinement procedure, which extend to 37 000 cm−1 and J=20 (total angular momentum), the average root-mean-square deviation is approximately 0.05 cm−1. This potential energy surface offers significant improvements when compared to recent models by accurately predicting states possessing high values of J. This feature will offer significant improvements in calculated line positions for high-temperature spectra where transitions between high J states become more prominent. Combining this potential with the latest dipole moment surface for water vapour, a line list has been calculated which extends reliably to 37 000 cm−1. Obtaining reliable results in the ultraviolet is of special importance as it is a challenging spectral region for the water molecule both experimentally and theoretically. Comparisons are made against several experimental sources of cross sections in the near-ultraviolet and discrepancies are observed. In the near-ultraviolet our calculations are in agreement with recent atmospheric retrievals and the upper limit obtained using broadband spectroscopy by Wilson et al. (2016, p. 194), but they do not support recent suggestions of very strong absorption in this region.

Description: In our previous studies, we investigated the performance of a novel pneumatic planar magnetic separator (PMS) for the dry beneficiation of a selected magnetite ore. In the present study, we have extended the studies on the PMS with the focus on investigating how various PMS processing flowsheet configurations influence its performance. The outcomes were subsequently compared with those of a Davis tube recovery (DTR) tester. The study demonstrated that the use of PMS in the dry beneficiation of magnetite ores is feasible, and operating the PMS in different flowsheet configurations positively influences the magnetite concentrate grade and purity. Finally, the study showed that the PMS performance compares well with that of DTR and can potentially replace DTR in operations that are carried out in arid regions.

Description: Geospatial data is urgently needed in decision-making processes to achieve Sustainable Development Goals (SDGs) at global, national, regional and local scales. While the advancement of geo-technologies to obtain or produce geospatial data has become faster and more affordable, many countries in the global south still experience a geospatial data scarcity at the rural level due to complex geographical terrains, weak coordination among institutions and a lack of knowledge and technologies to produce visualised geospatial data like maps. We proposed a collaborative spatial learning framework that integrates the spatial knowledge of stakeholders to obtain geospatial data. By conducting participatory mapping workshops in three villages in the Deli Serdang district in Indonesia, we tested the framework in terms of facilitating communication and collaboration of the village stakeholders while also supporting knowledge co-production and social learning among them. Satellite images were used in digital and non-digital mapping workshops to support village stakeholders to produce proper village maps while fulfilling the SDGs’ emphasis to make geospatial data available through a participatory approach.

Description: In this study, the flotation behavior of malachite was investigated using hydrophobic talc nanoparticles (TNs) as collectors. To improve the floatability of TN-deposited malachite, various experimental parameters were systematically investigated. We found that the floatability sharply increased as the size of the TNs decreased. The floatability of malachite was enhanced in the presence of smaller TNs, since higher amounts of smaller TNs were deposited on the surface of the malachite, thus rendering the surface more hydrophobic. Moreover, the floatability of the TN-deposited malachite increased as the pH decreased, likely due to the more favorable interaction between TNs and malachite by means of electrostatic attraction. Furthermore, the floatability became more enhanced as the TN concentration increased, likely associated with increases in the amount of TNs deposited on the surface of the malachite, thus enhancing the floatability by altering the hydrophobicity of the surface. Our findings suggest that the application of natural hydrophobic TNs as collectors in malachite flotation should be introduced as a new concept.

Description: Various studies investigated the fate of evaporation and the origin of precipitation. The more recent studies among them were often carried out with the help of numerical moisture tracking. Many research questions could be answered within this context, such as dependencies of atmospheric moisture transfers between different regions, impacts of land cover changes on the hydrological cycle, sustainability-related questions, and questions regarding the seasonal and interannual variability of precipitation. In order to facilitate future applications, global datasets on the fate of evaporation and the sources of precipitation are needed. Since most studies are on a regional level and focus more on the sources of precipitation, the goal of this study is to provide a readily available global dataset on the fate of evaporation for a fine-meshed grid of source and receptor cells. The dataset was created through a global run of the numerical moisture tracking model Water Accounting Model-2layers (WAM-2layers) and focused on the fate of land evaporation. The tracking was conducted on a 1.5∘×1.5∘ grid and was based on reanalysis data from the ERA-Interim database. Climatic input data were incorporated in 3- to 6-hourly time steps and represent the time period from 2001 to 2018. Atmospheric moisture was tracked forward in time and the geographical borders of the model were located at ±79.5∘ latitude. As a result of the model run, the annual, the monthly and the interannual average fate of evaporation were determined for 8684 land grid cells (all land cells except those located within Greenland and Antarctica) and provided via source–receptor matrices. The gained dataset was complemented via an aggregation to country and basin scales in order to highlight possible usages for areas of interest larger than grid cells. This resulted in data for 265 countries and 8223 basins. Finally, five types of source–receptor matrices for average moisture transfers were chosen to build the core of the dataset: land grid cell to grid cell, country to grid cell, basin to grid cell, country to country, basin to basin. The dataset is, to our knowledge, the first ready-to-download dataset providing the overall fate of evaporation for land cells of a global fine-meshed grid in monthly resolution. At the same time, information on the sources of precipitation can be extracted from it. It could be used for investigations into average annual, seasonal, and interannual sink and source regions of atmospheric moisture from land masses for most of the regions in the world and shows various application possibilities for studying interactions between people and water, such as land cover changes or human water consumption patterns. The dataset is accessible under https://doi.org/10.1594/PANGAEA.908705 (Link et al., 2019a) and comes along with example scripts for reading and plotting the data.

Description: Protected areas are the backbone of biodiversity conservation but are fixed in space and vulnerable to anthropogenic climate change. Myanmar is exceptionally rich in biodiversity but has a small protected area system. This study aimed to assess the potential vulnerability of this system to climate change. In the absence of good biodiversity data, we used a spatial modeling approach based on a statistically derived bioclimatic stratification (the Global Environmental Stratification, GEnS) to understand the spatial implications of projected climate change for Myanmar’s protected area system by 2050 and 2070. Nine bioclimatic zones and 41 strata were recognized in Myanmar, but their representation in the protected area system varied greatly, with the driest zones especially underrepresented. Under climate change, most zones will shift upslope, with some protected areas projected to change entirely to a new bioclimate. Potential impacts on biodiversity include mountaintop extinctions of species endemic to isolated peaks, loss of climate specialists from small protected areas and those with little elevational range, and woody encroachment into savannas and open forests as a result of both climate change and rising atmospheric CO2. Myanmar needs larger, better connected, and more representative protected areas, but political, social, and economic problems make this difficult.

Description: Minerals from mantle xenoliths in the Zapolyarnaya pipe in the Upper Muna field, Russia and from mineral separates from other large diamondiferous kimberlite pipes in this field (Deimos, Novinka and Komsomolskaya-Magnitnaya) were studied with EPMA and LA-ICP-MS. All pipes contain very high proportions of sub-calcic garnets. Zapolyarnaya contains mainly dunitic xenoliths with veinlets of garnets, phlogopites and Fe-rich pyroxenes similar in composition to those from sheared peridotites. PT estimates for the clinopyroxenes trace the convective inflection of the geotherm (40–45 mW·m−2) to 8 GPa, inflected at 6 GPa and overlapping with PT estimates for ilmenites derived from protokimberlites. The Upper Muna mantle lithosphere includes dunite channels from 8 to 2 GPa, which were favorable for melt movement. The primary layering deduced from the fluctuations of CaO in garnets was smoothed by the refertilization events, which formed additional pyroxenes. Clinopyroxenes from the Novinka and Komsomolskaya-Magnitnaya pipes show a more linear geotherm and three branches in the P-Fe# plot from the lithosphere base to the Moho, suggesting several episodes of pervasive melt percolation. Clinopyroxenes from Zapolyarnaya are divided into four groups according to thermobarometry and trace element patterns, which show a stepwise increase of REE and incompatible elements. Lower pressure groups including dunitic garnets have elevated REE with peaks in Rb, Th, Nb, Sr, Zr, and U, suggesting mixing of the parental protokimberlitic melts with partially melted metasomatic veins of ancient subduction origin. At least two stages of melt percolation formed the inclined PT paths: (1) an ancient garnet semi-advective geotherm (35–45 mW·m−2) formed by volatile-rich melts during the major late Archean event of lithosphere growth; and (2) a hotter megacrystic PT path (Cpx-Ilm) formed by feeding systems for kimberlite eruptions (40–45 mW·m−2). Ilmenite PT estimates trace three separate PT trajectories, suggesting a multistage process associated with metasomatism and formation of the Cpx-Phl veinlets in dunites. Heating associated with intrusions of protokimberlite caused reactivation of the mantle metasomatites rich in H2O and alkali metals and possibly favored the growth of large megacrystalline diamonds.

Description: People spend most of their time in indoors and, as a result, indoor air quality has become an issue of increasing concern. Due to the use of coal and heavy transportation in Beijing, China, concentrations of polycyclic aromatic hydrocarbons (PAHs) bound to PM2.5 have risen and caused concerns about health risk, both outdoors and indoors. This study carried out quantitative investigation of PM2.5-bound PAHs in middle school classrooms and estimated the health risk to adolescents. According to the results, indoor PM2.5 concentrations ranged from 20.9 μg/m3 to 257.6 μg/m3, indoor PAH concentrations ranged from 8.0 ng/m3 to 83.0 ng/m3, and both were statistically correlated with outdoor concentrations. Results of diagnostic ratios (DR) and the PMF (positive matrix factorization) model indicated that coal combustion was the main source of PAHs in the classroom environment. The average value of incremental lifetime cancer risk (ILCR) was estimated to be 1.49 × 10−6, which indicated a potential health risk to students according to USEPA standards. Predictions showed that by 2021–2022, the risk will be reduced to an acceptable level. Results of this study could provide useful information for air pollution control in Beijing and proposing targeted solution against indoor air pollution.

Description: An experimental product study of the reactions of furfural with the main tropospheric oxidants (Cl, OH and NO3) has been carried out using a Fourier Transform Infrared spectrophotometer (FTIR) and a gas chromatograph–mass spectrometer with a time of flight detector (GC–TOFMS). The main gas-phase products detected were 5-chloro-2(5H)-furanone, maleic anhydride, 2-nitrofuran and CO. Molar yields were quantified for the detected products in these reactions, thus suggesting the existence of nongaseous products that could not be observed with the analytical techniques employed. The formation of Secondary Organic Aerosol (SOA) from the oxidation of furfural with Cl atoms, OH, NO3 and ozone was investigated in a smog chamber in the absence of inorganic seed aerosols. The experimental results show the formation of ultrafine particles (less than 1 µm in diameter) for all of the studied reactions except for the nitrate radical. Given their small size, these ultrafine particles (

Description: Precise exposure assessment of air pollutants is crucial in epidemiologic studies to ensure valid estimates of health effects. We conducted a longitudinal study to evaluate the role of air quality monitoring (AQM) measurements and high-resolution modeling outcomes focusing on nitrogen dioxide (NO2) exposure and atopic dermatitis (AD). A total of 128 young children with AD in Seoul Metropolitan Area, Korea, were recruited as a panel. We estimated the participants’ exposure to NO2 for four months, from 1 April through 31 July 2014 based on (1) monitored levels from 60 AQM stations located at varying distances from residential areas (AQM station-based NO2, AQM-NO2) and (2) estimates from a community multi-scale air quality (CMAQ) modeling system with a high-resolution (1 × 1 km) (CMAQ-NO2). We then compared the effect of AQM-NO2 on AD symptoms with that of CMAQ-NO2. The average distance between the participants’ residences and the nearest AQM station was 2.03 ± 1.06 km, ranging from 0.28 km to 5.73 km. Based on AQM-NO2, the AD symptoms increased by 10.28% (95% confidence interval (CI): 3.24, 17.79) with an increase of 10 ppb of NO2. The effect estimates of CMAQ-NO2 were similar to those of AQM-NO2 when assessed in patients living within 3 km from the nearest AQM station. Even within 1 km, the CI estimate obtained from the CMAQ was much narrower than from AQM (44.18–49.54 vs. 7.02–64.75). However, the association of AQM-NO2 with AD symptoms of patients living beyond 3 km was not positive, whereas that of CMAQ-NO2 maintained positive. In conclusion, exposure to ambient NO2 is significantly associated with aggravation of AD symptoms in young children. In addition, our study suggests that exposure assessment of NO2 using measurement data obtained from monitoring stations far from residential locations can lead to misclassification bias.

Description: Reducing plastic pollution in rivers, lakes, and oceans is beneficial to aquatic animals and human livelihood. To achieve this, reliable observations of the abundance, spatiotemporal variation, and composition of plastics in aquatic ecosystems are crucial. Current efforts mainly focus on collecting data on the open ocean, on beaches and coastlines, and in river systems. Urban areas are the main source of plastic leakage into the natural environment, yet data on plastic pollution in urban water systems are scarce. In this paper, we present a simple method for plastic hotspot mapping in urban water systems. Through visual observations, macroplastic abundance and polymer categories are determined. Due to its simplicity, this method is suitable for citizen science data collection. A first application in the Dutch cities of Leiden and Wageningen showed similar mean plastic densities (111–133 items/km canal) and composition (75–80% soft plastics), but different spatial distributions. These observations emphasize the importance of long-term data collection to further understand and quantify spatiotemporal variations of plastics in urban water systems. In turn, this will support improved estimates of the contribution of urban areas to the plastic pollution of rivers and oceans.

Description: Black carbon (BC), organic carbon (OC), and total carbon (TC) in snow are important for their climatic and cryospheric effects. They are also part of the global carbon cycle. Atmospheric black and organic carbon (including brown carbon) may deposit and darken snow surfaces. Currently, there are no standardized methods for sampling, filtering, and analysis protocols to detect carbon in snow. Here, we describe our current methods and protocols to detect carbon in seasonal snow using the OCEC thermal optical method, a European standard for atmospheric elemental carbon (EC). We analyzed snow collected within and around the urban background SMEARIII (Station for Measuring Ecosystem-Atmosphere Relations) at Kumpula (60° N) and the Arctic GAW (Global Atmospheric Watch) station at Sodankylä (67° N). The median BC, OC, and TC in snow samples (ntot = 30) in Kumpula were 1118, 5279, and 6396 ppb, and in Sodankylä, they were 19, 1751, and 629 ppb. Laboratory experiments showed that error due to carbon attached to a sampling bag (n = 11) was

Description: Optimum management of irrigated crops in regions with shallow saline groundwater requires a careful balance between application of irrigation water and upward movement of salinity from the groundwater. Few field-validated surrogate models are available to aid in the management of irrigation water under shallow groundwater conditions. The objective of this research is to develop a model that can aid in the management using a minimum of input data that are field validated. In this paper a 2-year field experiment was carried out in the Hetao irrigation district in Inner Mongolia, China, and a physically based integrated surrogate model for arid irrigated areas with shallow groundwater was developed and validated with the collected field data. The integrated model that links crop growth with available water and salinity in the vadose zone is called Evaluation of the Performance of Irrigated Crops and Soils (EPICS). EPICS recognizes that field capacity is reached when the matric potential is equal to the height above the groundwater table and thus not by a limiting hydraulic conductivity. In the field experiment, soil moisture contents and soil salt conductivity at five depths in the top 100 cm, groundwater depth, crop height, and leaf area index were measured in 2017 and 2018. The field results were used for calibration and validation of EPICS. Simulated and observed data fitted generally well during both calibration and validation. The EPICS model that can predict crop growth, soil water, groundwater depth, and soil salinity can aid in optimizing water management in irrigation districts with shallow aquifers.

Description: Recently, it has been established that interplanetary coronal mass ejections (ICMEs) can dramatically affect both trapped electron fluxes in the outer radiation belt and precipitating electron fluxes lost from the belt into the atmosphere. Precipitating electron fluxes and energies can vary over a range of timescales during these events. These variations depend on the initial energy and location of the electron population and the ICME characteristics and structures. One important factor controlling electron dynamics is the magnetic field orientation within the ejecta that is an integral part of the ICME. In this study, we examine Van Allen Probes (RBSPs) and Polar Orbiting Environmental Satellites (POESs) data to explore trapped and precipitating electron fluxes during two ICMEs. The ejecta in the selected ICMEs have magnetic cloud characteristics that exhibit the opposite sense of the rotation of the north–south magnetic field component (BZ). RBSP data are used to study trapped electron fluxes in situ, while POES data are used for electron fluxes precipitating into the upper atmosphere. The trapped and precipitating electron fluxes are qualitatively analysed to understand their variation in relation to each other and to the magnetic cloud rotation during these events. Inner magnetospheric wave activity was also estimated using RBSP and Geostationary Operational Environmental Satellite (GOES) data. In each event, the largest changes in the location and magnitude of both the trapped and precipitating electron fluxes occurred during the southward portion of the magnetic cloud. Significant changes also occurred during the end of the sheath and at the sheath–ejecta boundary for the cloud with south to north magnetic field rotation, while the ICME with north to south rotation had significant changes at the end boundary of the cloud. The sense of rotation of BZ and its profile also clearly affects the coherence of the trapped and/or precipitating flux changes, timing of variations with respect to the ICME structures, and flux magnitude of different electron populations. The differing electron responses could therefore imply partly different dominant acceleration mechanisms acting on the outer radiation belt electron populations as a result of opposite magnetic cloud rotation.

Description: In this paper, we quantify the CO2 and N2O emissions from denitrification over the Amazonian wetlands. The study concerns the entire Amazonian wetland ecosystem with a specific focus on three floodplain (FP) locations: the Branco FP, the Madeira FP and the FP alongside the Amazon River. We adapted a simple denitrification model to the case of tropical wetlands and forced it by open water surface extent products from the Soil Moisture and Ocean Salinity (SMOS) satellite. A priori model parameters were provided by in situ observations and gauging stations from the HYBAM Observatory. Our results show that the denitrification and the trace gas emissions present a strong cyclic pattern linked to the inundation processes that can be divided into three distinct phases: activation, stabilization and deactivation. We quantify the average yearly denitrification and associated emissions of CO2 and N2O over the entire watershed at 17.8 kgN ha−1 yr−1, 0.37 gC-CO2 m−2 yr−1 and 0.18 gN-N2O m−2 yr−1 respectively for the period 2011–2015. When compared to local observations, it was found that the CO2 emissions accounted for 0.01 % of the integrated ecosystem, which emphasizes the fact that minor changes to the land cover may induce strong impacts on the Amazonian carbon budget. Our results are consistent with the state of the art of global nitrogen models with a positive bias of 28 %. When compared to other wetlands in different pedoclimatic environments we found that the Amazonian wetlands have similar emissions of N2O with the Congo tropical wetlands and lower emissions than the temperate and tropical anthropogenic wetlands of the Garonne (France), the Rhine (Europe) and south-eastern Asia rice paddies. In summary our paper shows that a data-model-based approach can be successfully applied to quantify N2O and CO2 fluxes associated with denitrification over the Amazon basin. In the future, the use of higher-resolution remote sensing products from sensor fusion or new sensors like the Surface Water and Ocean Topography (SWOT) mission will permit the transposition of the approach to other large-scale watersheds in tropical environments.

Description: In the framework of the EU Copernicus programme, the European Centre for Medium-Range Weather Forecasts (ECMWF) on behalf of the Joint Research Centre (JRC) is forecasting daily fire weather indices using its medium-range ensemble prediction system. The use of weather forecasts in place of local observations can extend early warnings by up to 1–2 weeks, allowing for greater proactive coordination of resource-sharing and mobilization within and across countries. Using 1 year of pre-operational service in 2017 and the Fire Weather Index (FWI), here we assess the capability of the system globally and analyse in detail three major events in Chile, Portugal and California. The analysis shows that the skill provided by the ensemble forecast system extends to more than 10 d when compared to the use of mean climate, making a case for extending the forecast range to the sub-seasonal to seasonal timescale. However, accurate FWI prediction does not translate into accuracy in the forecast of fire activity globally. Indeed, when all fires detected in 2017 are considered, including agricultural- and human-induced burning, high FWI values only occur in 50 % of the cases and are limited to the Boreal regions. Nevertheless for very large events which were driven by weather conditions, FWI forecasts provide advance warning that could be instrumental in setting up management and containment strategies.

Description: North African dust reaches the southeastern United States every summer. Size-resolved dust mass measurements taken in Miami, Florida, indicate that more than one-half of the surface dust mass concentrations reside in particles with geometric diameters less than 2.1 µm, while vertical profiles of micropulse lidar depolarization ratios show dust reaching above 4 km during pronounced events. These observations are compared to the representation of dust in the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) aerosol reanalysis and closely related Goddard Earth Observing System model version 5 (GEOS-5) Forward Processing (FP) aerosol product, both of which assimilate satellite-derived aerosol optical depths using a similar protocol and inputs. These capture the day-to-day variability in aerosol optical depth well, in a comparison to an independent sun-photometer-derived aerosol optical depth dataset. Most of the modeled dust mass resides in diameters between 2 and 6 µm, in contrast to the measurements. Model-specified mass extinction efficiencies equate light extinction with approximately 3 times as much aerosol mass, in this size range, compared to the measured dust sizes. GEOS-5 FP surface-layer sea salt mass concentrations greatly exceed observed values, despite realistic winds and relative humidities. In combination, these observations help explain why, despite realistic total aerosol optical depths, (1) free-tropospheric model volume extinction coefficients are lower than those retrieved from the micro-pulse lidar, suggesting too-low model dust loadings in the free troposphere, and (2) model dust mass concentrations near the surface can be higher than those measured. The modeled vertical distribution of dust, when captured, is reasonable. Large, aspherical particles exceeding the modeled dust sizes are also occasionally present, but dust particles with diameters exceeding 10 µm contribute little to the measured total dust mass concentrations after such long-range transport. Remaining uncertainties warrant a further integrated assessment to confirm this study's interpretations.

Description: Anaerobic nitrate-dependent Fe(II) oxidation (NDFeO) is widespread in various aquatic environments and plays a major role in iron and nitrogen redox dynamics. However, evidence for truly enzymatic, autotrophic NDFeO remains limited, with alternative explanations involving the coupling of heterotrophic denitrification with the abiotic oxidation of structurally bound or aqueous Fe(II) by reactive intermediate nitrogen (N) species (chemodenitrification). The extent to which chemodenitrification is caused (or enhanced) by ex vivo surface catalytic effects has not been directly tested to date. To determine whether the presence of either an Fe(II)-bearing mineral or dead biomass (DB) catalyses chemodenitrification, two different sets of anoxic batch experiments were conducted: 2 mM Fe(II) was added to a low-phosphate medium, resulting in the precipitation of vivianite (Fe3(PO4)2), to which 2 mM nitrite (NO2-) was later added, with or without an autoclaved cell suspension (∼1.96×108 cells mL−1) of Shewanella oneidensis MR-1. Concentrations of nitrite (NO2-), nitrous oxide (N2O), and iron (Fe2+, Fetot) were monitored over time in both set-ups to assess the impact of Fe(II) minerals and/or DB as catalysts of chemodenitrification. In addition, the natural-abundance isotope ratios of NO2- and N2O (δ15N and δ18O) were analysed to constrain the associated isotope effects. Up to 90 % of the Fe(II) was oxidized in the presence of DB, whereas only ∼65 % of the Fe(II) was oxidized under mineral-only conditions, suggesting an overall lower reactivity of the mineral-only set-up. Similarly, the average NO2- reduction rate in the mineral-only experiments (0.004±0.003 mmol L−1 d−1) was much lower than in the experiments with both mineral and DB (0.053±0.013 mmol L−1 d−1), as was N2O production (204.02±60.29 nmol L−1 d−1). The N2O yield per mole NO2- reduced was higher in the mineral-only set-ups (4 %) than in the experiments with DB (1 %), suggesting the catalysis-dependent differential formation of NO. N-NO2- isotope ratio measurements indicated a clear difference between both experimental conditions: in contrast to the marked 15N isotope enrichment during active NO2- reduction (15εNO2=+10.3 ‰) observed in the presence of DB, NO2- loss in the mineral-only experiments exhibited only a small N isotope effect (

Description: The Mediterranean (MED) Basin is a climate change hotspot that has seen drying and a pronounced increase in heatwaves over the last century. At the same time, it is experiencing increased heavy precipitation during wintertime cold spells. Understanding and quantifying the risks from compound events over the MED is paramount for present and future disaster risk reduction measures. Here, we apply a novel method to study compound events based on dynamical systems theory and analyse compound temperature and precipitation events over the MED from 1979 to 2018. The dynamical systems analysis quantifies the strength of the coupling between different atmospheric variables over the MED. Further, we consider compound warm–dry anomalies in summer and cold–wet anomalies in winter. Our results show that these warm–dry and cold–wet compound days are associated with large values of the temperature–precipitation coupling parameter of the dynamical systems analysis. This indicates that there is a strong interaction between temperature and precipitation during compound events. In winter, we find no significant trend in the coupling between temperature and precipitation. However in summer, we find a significant upward trend which is likely driven by a stronger coupling during warm and dry days. Thermodynamic processes associated with long-term MED warming can best explain the trend, which intensifies compound warm–dry events.

Description: Recent advances in particle sensor technologies have led to an increased development and utilization of low-cost, compact, particulate matter (PM) monitors. These devices can be deployed in dense monitoring networks, enabling an improved characterization of the spatiotemporal variability in ambient levels and exposure. However, the reliability of their measurements is an important prerequisite, necessitating rigorous performance evaluation and calibration in comparison to reference-grade instrumentation. In this study, field evaluation of Purple Air PA-II devices (low-cost PM sensors) is performed in two urban environments and across three seasons in Greece, in comparison to different types of reference instruments. Measurements were conducted in Athens (the largest city in Greece with nearly four-million inhabitants) for five months spanning over the summer of 2019 and winter/spring of 2020 and in Ioannina, a medium-sized city in northwestern Greece (100,000 inhabitants) during winter/spring 2019–2020. The PM2.5 sensor output correlates strongly with reference measurements (R2 = 0.87 against a beta attenuation monitor and R2 = 0.98 against an optical reference-grade monitor). Deviations in the sensor-reference agreement are identified as mainly related to elevated coarse particle concentrations and high ambient relative humidity. Simple and multiple regression models are tested to compensate for these biases, drastically improving the sensor’s response. Large decreases in sensor error are observed after implementation of models, leading to mean absolute percentage errors of 0.18 and 0.12 for the Athens and Ioannina datasets, respectively. Overall, a quality-controlled and robustly evaluated low-cost network can be an integral component for air quality monitoring in a smart city. Case studies are presented along this line, where a network of PA-II devices is used to monitor the air quality deterioration during a peri-urban forest fire event affecting the area of Athens and during extreme wintertime smog events in Ioannina, related to wood burning for residential heating.

Description: The article is devoted to a reconstruction of the sedimentation processes in Kamyshovoye Lake (the Kaliningrad Region, Russia) during the Late Glacial and Holocene. The results of the geochemical analysis of Kamyshovoye Lake’s bottom sediments, accompanied by statistical processing and detailed radiocarbon dating, are presented. It was established that a high proportion of mineral matter dominated in the intervals between 15,000 and 11,400 and between 1400 and 600 cal y BP; enrichment with carbonates was noted between 11,400 and 5200 cal y BP and during the past 600 years; and a high percentage of organic matter was recorded between 7800 and 600 cal y BP. We conclude that the increase in mineral matter was influenced by such factors as reduced vegetation cover due to natural and anthropogenic processes, aeolian transfer, and dead-ice melting during the Late Glacial. The increase in carbonate matter was mainly associated with humidity and the reduction conditions of the lake ecosystem. Organogenic matter content was affected by the autochthonic (biological) productivity of the lake, which directly depends on more favorable climatic conditions.

Description: Fertilized agricultural soils serve as a primary source of anthropogenic N2O emissions. In South Africa, there is a paucity of data on N2O emissions from fertilized, irrigated dairy-pasture systems and emission factors (EF) associated with the amount of N applied. A first study aiming to quantify direct N2O emissions and associated EFs of intensive pasture-based dairy systems in sub-Sahara Africa was conducted in South Africa. Field trials were conducted to evaluate fertilizer rates (0, 220, 440, 660, and 880 kg N ha−1 year−1) on N2O emissions from irrigated kikuyu–perennial ryegrass (Pennisetum clandestinum–Lolium perenne) pastures. The static chamber method was used to collect weekly N2O samples for one year. The highest daily N2O fluxes occurred in spring (0.99 kg ha−1 day−1) and summer (1.52 kg ha−1 day−1). Accumulated N2O emissions ranged between 2.45 and 15.5 kg N2O-N ha−1 year−1 and EFs for mineral fertilizers applied had an average of 0.9%. Nitrogen in yielded herbage varied between 582 and 900 kg N ha−1. There was no positive effect on growth of pasture herbage from adding N at high rates. The relationship between N balance and annual N2O emissions was exponential, which indicated that excessive fertilization of N will add directly to N2O emissions from the pastures. Results from this study could update South Africa’s greenhouse gas inventory more accurately to facilitate Tier 3 estimates.

Description: Achieving sustainable zero-waste and carbon neutral solutions that contribute to a circular economy is critically important for the long-term prosperity and continuity of traditional carbon-based energy industries. The Estonian oil shale (OS) sector is an example where such solutions are more than welcome. The combustion of OS generates a continuous flow of ashes destined to landfills. In this study, the technical feasibility of producing monolith building materials incorporating different OS ashes from Estonia was evaluated. Three binder systems were studied: self-cementation of the ashes, ceramic sintering in clay brick production and accelerated carbonation of OS ash (OSA) compacts. Results showed that most of the OSAs studied have low self-cementitious properties and these properties were affected by ash fineness and mineralogical composition. In case of clay bricks, OSA addition resulted in a higher porosity and improved insulation properties. The carbonated OSA compacts showed promising compressive strength. Accelerated carbonation of compacted samples was found to be the most promising way for the future utilization of OSAs as sustainable zero-waste and carbon neutral solution.

Description: This paper presents a review and analysis of large-scale air convection tests and the establishment of intrinsic permeability in coarse open-graded materials. Natural air convection can make a significant contribution to heat transfer during cooling periods. In seasonally freezing environments this can result in excessive frost penetration and subsequent frost-related problems. Intrinsic permeability largely defines the onset of convective heat transfer in granular materials. Conventional methods for measuring intrinsic permeability cannot be applied to very coarse materials. Large-scale laboratory experiments on natural air convection can serve as an alternative method for determining this crucial parameter. This paper gives an overview of four different experimental test setups for measuring natural air convection, all differing in physical shape, boundary conditions and heat flux/temperature measurement devices. Comparison between these is difficult because the air convection pattern can differ and in some cases the shape and number of convection cells cannot be validated. Most of the studies available in the literature use theoretical equations to approximate intrinsic permeability. A method based on the analytical Nu-Ra number relationship is employed to establish the values of intrinsic permeability. Tests that provide enough data to enable the use of the Nu-Ra relationship are very limited. The overall results show a reasonable correlation between experiment-based intrinsic permeability and theoretical approximation. However, several issues must be addressed: first, differences may exist between the intrinsic permeability of natural and of crushed materials due to the shape effect. Second, the method used is in theory valid only for two-dimensional air convection within a square enclosure heated from below. Yet the results show that this method could be extended to other conditions with a certain degree of confidence. Third, a good estimate of intrinsic permeability is possible only with accurate experimental measurement.

Description: Anhand von sedimentologischen und geländemorphologischen Untersuchungen wird die Abschmelzgeschichte des südöstlichen Chiemsee-Gletschers beschrieben. Mit dem Trockenfallen der Bad Adelholzen-Erlstätter Rinne im Verlaufe des Spätwürm entwickelt sich aus dem Abschmelzen des Eislappens in der Grabenstätter Bucht eine sich ständig tiefer legende konzentrische Abfolge von zunächst peripheren Entwässerungsrinnen, wobei die ältesten Rinnen dieser Phase bei Chieming, die jüngeren dann entsprechend weiter im Süden, in die zentripetale Richtung umschwenken. Die Entstehung des Tüttensee-Komplexes ist im Kontext dieser Entwicklung zu sehen. Er ist das Ergebnis der glazifluvialen und glazilakustrinen Sedimentation im Einflussbereich des sukzessiven Eisabbaus in der Grabenstätter Bucht in Kombination mit einer Toteisbildung im Bereich des heutigen Tüttensees. Dafür sprechen die stufenartige Abfolge der beschriebenen peripheren Abflussrinnen mit ihren immer tiefer liegenden Abflussniveaus, die Höhengleichheit von drei dieser Rinnen mit den Tüttensee-Terrassen sowie die für die jeweilige Terrassenentstehung typische glazifluviale bzw. delta-artige Sedimentstruktur und -reife. Dieses Ergebnis stellt ein Korrektiv zur Hypothese des Chiemgau-Impakts dar, wonach der Tüttensee ein Impaktkrater sein soll. Da diese nun falsifizierte Annahme vor allem im deutschsprachigen Raum von zahlreichen Medien propagiert wird, ist der folgende Artikel auf Deutsch verfasst, um einer breiten Leserschaft zugänglich zu sein.

Description: Temperature appears to be the best predictor of species composition of planktonic foraminifera communities, making it possible to use their fossil assemblages to reconstruct sea surface temperature (SST) variation in the past. However, the role of other environmental factors potentially modulating the spatial and vertical distribution of planktonic foraminifera species is poorly understood. This is especially relevant for environmental factors affecting the subsurface habitat. If such factors play a role, changes in the abundance of subsurface-dwelling species may not solely reflect SST variation. In order to constrain the effect of subsurface parameters on species composition, we here characterize the vertical distribution of living planktonic foraminifera community across an east–west transect through the subtropical South Atlantic Ocean, where SST variability was small, but the subsurface water mass structure changed dramatically. Four planktonic foraminifera communities could be identified across the top 700 m of the transect. Gyre and Agulhas Leakage surface faunas were predominantly composed of Globigerinoides ruber, Globigerinoides tenellus, Trilobatus sacculifer, Globoturborotalita rubescens, Globigerinella calida, Tenuitella iota, and Globigerinita glutinata, and these only differed in terms of relative abundances (community composition). Upwelling fauna was dominated by Neogloboquadrina pachyderma, Neogloboquadrina incompta, Globorotalia crassaformis, and Globorotalia inflata. Thermocline fauna was dominated by Tenuitella fleisheri, Globorotalia truncatulinoides, and Globorotalia scitula in the west and by G. scitula only in the east. The largest part of the standing stock was consistently found in the surface layer, but SST was not the main predictor of species composition either for the depth-integrated fauna across the stations or at individual depth layers. Instead, we identified a pattern of vertical stacking of different parameters controlling species composition, reflecting different aspects of the pelagic habitat. Whereas productivity appears to dominate in the mixed layer (0–60 m), physical properties (temperature, salinity) become important at intermediate depths and in the subsurface, a complex combination of factors including oxygen concentration is required to explain the assemblage composition. These results indicate that the seemingly straightforward relationship between assemblage composition and SST in sedimentary assemblages reflects vertically and seasonally integrated processes that are only indirectly linked to SST. It also implies that fossil assemblages of planktonic foraminifera should also contain a signature of subsurface processes, which could be used for paleoceanographic reconstructions.

Description: A remarkable exposure of tubular authigenic carbonates was found on the seafloor in the Dongsha area of the South China Sea (SCS). The tubular carbonates, around 2–3 cm in diameter and usually less than 10 cm in length, represent broken fragments of once-larger pipes that now protrude from muddy sediments. The morphology, carbon and oxygen stable isotope compositions, and trace and rare earth element contents of the carbonates were analyzed to decipher the mode of carbonate formation. The tubular carbonates exhibit a dark brown coating of iron and manganese hydrous oxides, indicating prolonged exposure to oxic bottom waters. The carbonate content of the micritic pipes falls between 12.5 and 67.3 wt.% with an average of 42.0 wt.%, suggesting formation within the sediment. This inference is supported by trace and rare earth element patterns including a moderate enrichment of middle rare earth elements. Low δ13C values (as low as −50.3‰, Vienna Pee Dee Belemnite (VPDB)) suggest that carbonate precipitation was induced by the anaerobic oxidation of methane. The unusually positive δ18O values of the carbonates (as high as +5.3‰, VPDB) are believed to reflect the destabilization of locally abundant gas hydrate. Taken together, it is suggested that pipe formation was initiated by sediment-dwelling organisms, such as crustaceans or bivalves. The burrows subsequently acted as conduits for upward fluid migration. The lithification of the sediment directly surrounding the conduits and the partial filling of the conduits with carbonate cement resulted in the formation of tubular carbonates. Turbidity currents, sediment slumps, or the vigorous emission of fluids probably induced the fragmentation of tubular carbonates within the sediment. The carbonate fragments had been further subjected to winnowing by bottom currents. This study provides insight into the interaction of megafauna burrowing with fluid migration and carbonate formation at hydrocarbon seeps, highlighting the role of bottom currents and mass wasting on the formation of fragmented tubular carbonates.

Description: In November 2019, the Supreme Court of India issued a notification to all the states in the National Capital Region of Delhi to install smog towers for clean air and allocated INR 36 crores (~USD 5.2 million) for a pilot. Can we vacuum our air pollution problem using smog towers? The short answer is “no”. Atmospheric science defines the air pollution problem as (a) a dynamic situation where the air is moving at various speeds with no boundaries and (b) a complex mixture of chemical compounds constantly forming and transforming into other compounds. With no boundaries, it is unscientific to assume that one can trap air, clean it, and release into the same atmosphere simultaneously. In this paper, we outline the basics of atmospheric science to describe why the idea of vacuuming outdoor air pollution is unrealistic, and the long view on air quality management in Indian cities.

Description: The problem of wave identification is formulated as applied to the results of measurements of the temperature and the density of the neutral atmosphere in the range height 90–120 km by the artificial periodic irregularities (APIs) technique. The technique is based on the resonant scattering of radio waves by artificial periodic irregularities of the ionospheric plasma emerging in the field of a standing wave arising from the interference of the incident and reflected waves from the ionosphere. APIs were created using SURA heating facility (named as SURA experiment). The acoustic wave theory is reformulated on the base of data which can be observed in the given experimental setup. The basic system of equations is reduced so that it accounts only upward and downward directed waves, ignoring entropy mode. The algorithm of wave identification based on usage of dynamic projection operators for such a reduced case is proposed and explicit form of projection operators is derived. Its application to finite number dataset via Discrete Fourier Transform (DFT) is described and results of its application to the DFT-transformed set of experimental observation of the temperature and density perturbations are presented. The result yields hybrid amplitudes, that allow us to calculate energy of the directed waves that enter the observed superposition. The problem of entropy mode detection is discussed, the corresponding projecting operators for the full evolution system are built and a way to apply the method to quantification of it is proposed.

Description: Light-use efficiency defines the ability of primary producers to convert sunlight energy to primary production and is computed as the ratio between the gross primary production and the intercepted photosynthetic active radiation. While this measure has been applied broadly within terrestrial ecology to investigate habitat resource-use efficiency, it remains underused within the aquatic realm. This report provides a conceptual framework to compute hourly and daily light-use efficiency using underwater O2 eddy covariance, a recent technological development that produces habitat-scale rates of primary production under unaltered in situ conditions. The analysis, tested on two benthic flux datasets, documents that hourly light-use efficiency may approach the theoretical limit of 0.125 O2 per photon under low-light conditions, but it decreases rapidly towards the middle of the day and is typically 10-fold lower on a 24 h basis. Overall, light-use efficiency provides a useful measure of habitat functioning and facilitates site comparison in time and space.

Description: This study is focused on a barian titanian phlogopite found in an alkaline ultramafic dyke transecting Mesozoic limestones of the Gargano Promontory (Apulia, Italy). The rock containing the barian titanian phlogopite, an olivine-clinopyroxene-rich lamprophyre with nepheline and free of feldspars, has been classified as monchiquite. The present study combines chemical analyses, single crystal X-ray diffraction and Raman spectroscopy. Chemical variations suggest that the entry of Ba into the phlogopite structure can be explained by the exchange Ba + Al = K + Si. The crystal structure refinement indicates that the Ti uptake is consistent with the Ti–oxy exchange mechanism. The structural parameters associated with the oxy substitution mechanism are extremely enhanced and rarely reported in natural phlogopite: (a) displacement of M2 cation toward the O4 site (~0.7); (b) M2 octahedron bond-length distortion (~2.5); (c) very short c cell parameter (~10.14 Å). Raman analysis showed most prominent features in the 800–200 cm−1 region with the strongest peaks occurring at 773 and 735 cm−1. Only a weak, broad band was observed to occur in the OH-stretching region. As concerns the origin of the barian titanian phlogopite, the rock textural features clearly indicate that it crystallized from pockets of the interstitial melt. Here, Ba and Ti enrichment took place after major crystallization of olivine under fast-cooling conditions, close to the dyke margin.

Description: The lithobiont community encrusting an early Pleistocene palaeocliff cropping out north of Augusta (SE Sicily, Italy) was investigated based on field observations and laboratory inspection of two rocky samples. Bryozoans, serpulids, brachiopods and bivalves encrusted part of the exposed surfaces that were bored mostly by clionaid sponges. Bryozoans, with at least 25 species detected on the rocky samples, are the most diversified skeletonized lithobionts also accounting for the highest number of colonies/specimens and highest coverage. Brachiopods, with the only species Novocrania anomala and a few but large cemented valves, cover wide surfaces. Serpulids, with two species identified on the sampled rocks and further two on the outcrop, were intermediate. A multiphase colonization is present, including a final epilithobiont community locally formed on eroded surfaces exposing a network of pervasive borings. The co-occurrence of very sciaphilic species having circalittoral to bathyal distributions suggests that the studied community thrived on a rocky substratum located near or at the shelf break, probably belonging to the shelf break (or RL) biocoenosis, also in agreement with observations on the fossil content of neighboring marly sediments. The observed relationships among colonizers largely represent mere superimpositions, and real interactions are not enough to state species competitiveness.

Description: Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96∘ S, 14.35∘ W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d−1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.

Description: Information on global cropland distribution and agricultural production is critical for the world's agricultural monitoring and food security. We present datasets of cropland extent and agricultural production in a two-paper series of a cultivated planet in 2010. In the first part, we propose a new Self-adapting Statistics Allocation Model (SASAM) to develop the global map of cropland distribution. SASAM is based on the fusion of multiple existing cropland maps and multilevel statistics of the cropland area, which is independent of training samples. First, cropland area statistics are used to rank the input cropland maps, and then a scoring table is built to indicate the agreement among the input datasets. Secondly, statistics are allocated adaptively to the pixels with higher agreement scores until the cumulative cropland area is close to the statistics. The multilevel allocation results are then integrated to obtain the extent of cropland. We applied SASAM to produce a global cropland synergy map with a 500 m spatial resolution for circa 2010. The accuracy assessments show that the synergy map has higher accuracy than the input datasets and better consistency with the cropland statistics. The synergy cropland map is available via an open-data repository (https://doi.org/10.7910/DVN/ZWSFAA; Lu et al., 2020). This new cropland map has been used as an essential input to the Spatial Production Allocation Model (SPAM) for producing the global dataset of agricultural production for circa 2010, which is described in the second part of the two-paper series.

Description: The western European kinematic evolution results from the opening of the western Neotethys and the Atlantic oceans since the late Paleozoic and the Mesozoic. Geological evidence shows that the Iberian domain recorded the propagation of these two oceanic systems well and is therefore a key to significantly advancing our understanding of the regional plate reconstructions. The late-Permian–Triassic Iberian rift basins have accommodated extension, but this tectonic stage is often neglected in most plate kinematic models, leading to the overestimation of the movements between Iberia and Europe during the subsequent Mesozoic (Early Cretaceous) rift phase. By compiling existing seismic profiles and geological constraints along the North Atlantic margins, including well data over Iberia, as well as recently published kinematic and paleogeographic reconstructions, we propose a coherent kinematic model of Iberia that accounts for both the Neotethyan and Atlantic evolutions. Our model shows that the Europe–Iberia plate boundary was a domain of distributed and oblique extension made of two rift systems in the Pyrenees and in the Iberian intra-continental basins. It differs from standard models that consider left-lateral strike-slip movement localized only in the northern Pyrenees in introducing a significant strike-slip movement south of the Ebro block. At a larger scale it emphasizes the role played by the late-Permian–Triassic rift and magmatism, as well as strike-slip faulting in the evolution of the western Neotethys Ocean and their control on the development of the Atlantic rift.

Description: This paper addresses the question of how much uncertainties in CO2 fluxes over Australia can be reduced by assimilation of total-column carbon dioxide retrievals from the Orbiting Carbon Observatory-2 (OCO-2) satellite instrument. We apply a four-dimensional variational data assimilation system, based around the Community Multiscale Air Quality (CMAQ) transport-dispersion model. We ran a series of observing system simulation experiments to estimate posterior error statistics of optimized monthly-mean CO2 fluxes in Australia. Our assimilations were run with a horizontal grid resolution of 81 km using OCO-2 data for 2015. Based on four representative months, we find that the integrated flux uncertainty for Australia is reduced from 0.52 to 0.13 Pg C yr−1. Uncertainty reductions of up to 90 % were found at grid-point resolution over productive ecosystems. Our sensitivity experiments show that the choice of the correlation structure in the prior error covariance plays a large role in distributing information from the observations. We also found that biases in the observations would significantly impact the inverted fluxes and could contaminate the final results of the inversion. Biases in prior fluxes are generally removed by the inversion system. Biases in the boundary conditions have a significant impact on retrieved fluxes, but this can be mitigated by including boundary conditions in our retrieved parameters. In general, results from our idealized experiments suggest that flux inversions at this unusually fine scale will yield useful information on the carbon cycle at continental and finer scales.

Description: In this paper, we present the first multiyear time series of atmospheric ammonia (NH3) ground-based measurements in the Paris region (Créteil, 48.79∘ N, 2.44∘ E, France) retrieved with the midresolution “Observations of the Atmosphere by Solar absorption Infrared Spectroscopy” (OASIS) ground-based Fourier transform infrared solar observatory. Located in an urban region, OASIS has previously been used for monitoring air quality (tropospheric ozone and carbon monoxide) thanks to its specific column sensitivity across the whole troposphere down to the atmospheric boundary layer. A total of 4920 measurements of atmospheric total columns of ammonia have been obtained from 2009 to 2017, with uncertainties ranging from 20 % to 35 %, and have been compared with NH3 concentrations derived from the Infrared Atmospheric Sounding Interferometer (IASI). OASIS ground-based measurements show significant interannual and seasonal variabilities of atmospheric ammonia. NH3 total columns over the Paris megacity (12 million people) vary seasonally by 2 orders of magnitude from approximately 0.1×1016 molec. cm−2 in winter to 10×1016 molec. cm−2 for spring peaks, probably due to springtime spreading of fertilizers on surrounding croplands.

Description: In order to track progress towards the global climate targets, the parties that signed the Paris Climate Agreement will regularly report their anthropogenic carbon dioxide (CO2) emissions based on energy statistics and CO2 emission factors. Independent evaluation of this self-reporting system is a fast-growing research topic. Here, we study the value of satellite observations of the column CO2 concentrations to estimate CO2 anthropogenic emissions with 5 years of the Orbiting Carbon Observatory-2 (OCO-2) retrievals over and around China. With the detailed information of emission source locations and the local wind, we successfully observe CO2 plumes from 46 cities and industrial regions over China and quantify their CO2 emissions from the OCO-2 observations, which add up to a total of 1.3 Gt CO2 yr−1 that accounts for approximately 13 % of mainland China's annual emissions. The number of cities whose emissions are constrained by OCO-2 here is 3 to 10 times larger than in previous studies that only focused on large cities and power plants in different locations around the world. Our satellite-based emission estimates are broadly consistent with the independent values from China's detailed emission inventory MEIC but are more different from those of two widely used global gridded emission datasets (i.e., EDGAR and ODIAC), especially for the emission estimates for the individual cities. These results demonstrate some skill in the satellite-based emission quantification for isolated source clusters with the OCO-2, despite the sparse sampling of this instrument not designed for this purpose. This skill can be improved by future satellite missions that will have a denser spatial sampling of surface emitting areas, which will come soon in the early 2020s.

Description: Violent tornadoes are rare in Europe but they can have devastating effects. Damage associated with individual tornadoes can reach several billion euros and they have caused hundreds of fatalities. The tornado risk varies considerably over Europe, but so far only a few national maps of tornado risk and one Europe-wide map exist. We show several different ways to create quantitative maps of tornado occurrence rates as follows: Kernel smoothing of observations, climatologies of convective parameters from reanalysis, output of a logistic regression model to link convective parameters with observed tornadoes, orography-dependent climatologies and finally the population-bias corrected tornado occurrence rates from the Risk Management Solutions (RMS) Europe Severe Convective Storm Model. We discuss advantages and disadvantages of each approach and compare the results. While the climatologies created from the individual methods show a lot of qualitative similarities, we advocate to combine the methods to achieve the most reliable quantitative climatology.

Description: We develop a new large-scale hydrological and water resources model, the Community Water Model (CWatM), which can simulate hydrology both globally and regionally at different resolutions from 30 arcmin to 30 arcsec at daily time steps. CWatM is open source in the Python programming environment and has a modular structure. It uses global, freely available data in the netCDF4 file format for reading, storage, and production of data in a compact way. CWatM includes general surface and groundwater hydrological processes but also takes into account human activities, such as water use and reservoir regulation, by calculating water demands, water use, and return flows. Reservoirs and lakes are included in the model scheme. CWatM is used in the framework of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP), which compares global model outputs. The flexible model structure allows for dynamic interaction with hydro-economic and water quality models for the assessment and evaluation of water management options. Furthermore, the novelty of CWatM is its combination of state-of-the-art hydrological modeling, modular programming, an online user manual and automatic source code documentation, global and regional assessments at different spatial resolutions, and a potential community to add to, change, and expand the open-source project. CWatM also strives to build a community learning environment which is able to freely use an open-source hydrological model and flexible coupling possibilities to other sectoral models, such as energy and agriculture.

Description: More people appreciate the importance of global sustainability [...]

Description: Retrievals of ice cloud properties require accurate estimates of ice particle mass. Empirical mass–dimensional (m–D) relationships in the form m = a D b are widely used and usually universally applied across the complete range of particle sizes. For the first time, the dependence of a and b coefficients in m–D relationships on median mass diameter (Dmm) is studied. Using combined cloud microphysical data collected during the Olympic Mountains Experiment and coincident observations from Airborne Precipitation Radar Third Generation, Dmm-dependent (a, b) coefficients are derived and represented as surfaces of equally plausible solutions determined by some tolerance in the chi-squared difference χ 2 that minimizes the difference between observed and retrieved radar reflectivity. Robust dependences of a and b on Dmm are shown with both parameters significantly decreasing with Dmm, leading to smaller effective densities for larger Dmm ranges. A universally applied constant m–D relationship overestimates the mass of large aggregates when Dmm is between 3–6 mm and temperatures are between −15–0 °C. Multiple m–D relations should be applied for different Dmm ranges in retrievals and simulations to account for the variability of particle sizes that are responsible for the mass and thus for the variability of particle shapes and densities.

Description: Roadside vegetation can affect the performance of installed road lighting. We demonstrate a workflow in which a car-mounted measurement system is used to assess the light-obstructing effect of roadside vegetation. The mobile mapping system (MMS) includes a panoramic camera system, laser scanner, inertial measurement unit, and satellite positioning system. The workflow and the measurement system were applied to a road section of Munkkiniemenranta, Helsinki, Finland, in 2015 and 2019. The relative luminance distribution on a road surface and the obstructing vegetation were measured before and after roadside vegetation pruning applying a luminance-calibrated mobile mapping system. The difference between the two measurements is presented, and the opportunities provided by the mobile 3D luminance measurement system are discussed.

Description: Density-dependent population growth regulates long-term urban expansion and shapes distinctive socioeconomic trends. Despite a marked heterogeneity in the spatial distribution of the resident population, Mediterranean European countries are considered more homogeneous than countries in other European regions as far as settlement structure and processes of metropolitan growth are concerned. However, rising socioeconomic inequalities among Southern European regions reflect latent demographic and territorial transformations that require further investigation. An integrated assessment of the spatio-temporal distribution of resident populations in more than 1000 municipalities (1961–2011) was carried out in this study to characterize density-dependent processes of metropolitan growth in Greece. Using geographically weighted regressions, the results of our study identified distinctive local relationships between population density and growth rates over time. Our results demonstrate that demographic growth rates were non-linearly correlated with other variables, such as population density, with positive and negative impacts during the first (1961–1971) and the last (2001–2011) observation decade, respectively. These findings outline a progressive shift over time from density-dependent processes of population growth, reflecting a rapid development of large metropolitan regions (Athens, Thessaloniki) in the 1960s, to density-dependent processes more evident in medium-sized cities and accessible rural regions in the 2000s. Density-independent processes of population growth have been detected in the intermediate study period (1971–2001). This work finally discusses how a long-term analysis of demographic growth, testing for density-dependent mechanisms, may clarify the intrinsic role of population concentration and dispersion in different phases of the metropolitan cycle in Mediterranean Europe.

Description: The complexity and high sensitivity of proteins to environmental factors give rise to a multitude of variables, which affect the stabilization mechanisms in protein foams. Interfacial and foaming properties of proteins have been widely studied, but the reported unique effect of pH, which can be of great interest to applications, has been investigated to a lesser extent. In this paper, we focus on the impact of pH on the stability of black foam films and corresponding foams obtained from solutions of a model globular protein—the whey β-lactoglobulin (BLG). Foam stability was analyzed utilizing three characteristic parameters (deviation time, transition time and half-lifetime) for monitoring the foam decay, while foam film stability was measured in terms of the critical disjoining pressure of film rupture. We attempt to explain correlations between the macroscopic properties of a foam system and those of its major building blocks (foam films and interfaces), and thus, to identify structure-property relationships in foam. Good correlations were found between the stabilities of black foam films and foams, while relations to the properties of adsorption layers appeared to be intricate. That is because pH-dependent interfacial properties of proteins usually exhibit an extremum around the isoelectric point (pI), but the stability of BLG foam films increases with increasing pH (3–7), which is well reflected in the foam stability. We discuss the possible reasons behind these intriguingly different behaviors on the basis of pH-induced changes in the molecular properties of BLG, which seem to be determining the mechanism of film rupture at the critical disjoining pressure.

Description: The oxidation state and local atomic environment of admixtures of In, Cu, and Ag in synthetic sphalerite crystals were determined by X-ray absorption spectroscopy (XAS). The sphalerite crystals doped with In, Cu, Ag, In–Cu, and In–Ag were synthesized utilizing gas transport, salt flux, and dry synthesis techniques. Oxidation states of dopants were determined using X-ray absorption near edge structure (XANES) technique. The local atomic structure was studied by X-ray absorption fine structure spectroscopy (EXAFS). The spectra were recorded at Zn, In, Ag, and Cu K-edges. In all studied samples, In was in the 3+ oxidation state and replaced Zn in the structure of sphalerite, which occurs with the expansion of the nearest coordination shells due to the large In ionic radius. In the presence of In, the oxidation state of Cu and Ag is 1+, and both metals can form an isomorphous solid solution where they substitute for Zn according to the coupled substitution scheme 2Zn2+ ↔ Me+ + In3+. Moreover, Ag K-edges EXAFS spectra fitting, combined with the results obtained for In- and Au-bearing sphalerite shows that the Me-S distances in the first coordination shell in the solid solution state are correlated with the ionic radii and increase in the order of Cu 〈 Ag 〈 Au. The distortion of the atomic structure increases in the same order. The distant (second and third) coordination shells of Cu and Ag in sphalerite are split into two subshells, and the splitting is more pronounced for Ag. Analysis of the EXAFS spectra, coupled with the results of DFT (Density Function Theory) simulations, showed that the In–In and Me+–In3+ clustering is absent when the metals are present in the sphalerite solid solution. Therefore, all studied admixtures (In, Cu, Ag), as well as Au, are randomly distributed in the matrix of sphalerite, where the concentration of the elements in the “invisible” form can reach a few tens wt.%.

Description: In urban planning and transportation management, the centrality characteristics of urban streets are vital measures to consider. Centrality can help in understanding the structural properties of dense traffic networks that affect both human life and activity in cities. Many cities classify urban streets to provide stakeholders with a group of street guidelines for possible new rehabilitation such as sidewalks, curbs, and setbacks. Transportation research always considers street networks as a connection between different urban areas. The street functionality classification defines the role of each element of the urban street network (USN). Some potential factors such as land use mix, accessible service, design goal, and administrators’ policies can affect the movement pattern of urban travelers. In this study, nine centrality measures are used to classify the urban roads in four cities evaluating the structural importance of street segments. In our work, a Stacked Denoising Autoencoder (SDAE) predicts a street’s functionality, then logistic regression is used as a classifier. Our proposed classifier can differentiate between four different classes adopted from the U.S. Department of Transportation (USDT): principal arterial road, minor arterial road, collector road, and local road. The SDAE-based model showed that regular grid configurations with repeated patterns are more influential in forming the functionality of road networks compared to those with less regularity in their spatial structure.

Description: In the past decades, severe heat waves have frequently occurred in many parts of the world. These conspicuous heat waves exerted terrible influences on human health, society, the economy, agriculture, the ecosystem and so on. Based on observed daily temperatures in China, an integrated index of heat waves and extreme-temperature days was established involving the frequency, duration, intensity and scale of these events across large cities in China. Heat waves and extreme-temperature days showed an increasing trend in most regions except northwest China from 1955 to 2014. After the late 1980s, the increasing trend was more obvious than the decades before. The cities in the middle and lower reaches of the Yangtze River were threatened by the most serious heat events in the past 60 years, especially Chongqing and Changsha. Due to the subtropical monsoon climate and special terrain, Chongqing experienced the most heat events in a long period of time. In particular, there was obvious fluctuation of hot years in 31 cities, which did not continuously rise with global warming; 21 cities mainly located in the eastern and southern regions of China had an obvious rising trend; eight cities had a clear declining trend which was mainly distributed in the western and northern regions of China; and there were no extreme-temperature days in Kunming and Lhasa in the past 60 years. The study revealed an obvious differentiation of heat events for 31 cities under climate change; heat threat in most cities is increasing but declining or remaining unchanged in the other cities. The trend is likely to intensify with global warming.

Description: Understanding the behavior of halogens (Cl, Br, and I) in subduction zones is critical to constrain the geochemical cycle of these volatiles and associated trace metals, as well as to quantify the halogen fluxes to the atmosphere via volcanic degassing. Here, the partitioning of bromine between coexisting aqueous fluids and hydrous granitic melts and its speciation in slab-derived fluids have been investigated in situ up to 840 ∘C and 2.2 GPa by synchrotron x-ray fluorescence (SXRF) and x-ray absorption spectroscopy (XAS) in diamond anvil cells. The partition coefficients DBrf/m range from ∼2 to ∼15, with an average value of 6.7±3.6 (1σ) over the whole pressure–temperature (P–T) range, indicating a moderate Br enrichment in aqueous fluids, in agreement with previous work. Extended x-ray-absorption fine-structure (EXAFS) analysis further evidences a gradual evolution of Br speciation from hydrated Br ions [Br(H2O)6]− in slab dehydration fluids to more complex structures involving both Na ions and water molecules, [BrNax(H2O)y], in hydrous silicate melts and supercritical fluids released at greater depth (〉 200 km). In denser fluids (ρ 〉 1.5 g cm−3) containing 60 wt % dissolved alkali–silicates and in hydrous Na2Si2O5 melts (10 wt % H2O), Br is found to be in a “salt-like” structure involving the six nearest Na ions and several next-nearest O neighbors that are either from water molecules and/or the silicate network. Bromine (and likely chlorine and iodine) complexing with alkalis is thus an efficient mechanism for the mobilization and transport of halogens by hydrous silicate melts and silica-rich supercritical fluids. Our results suggest that both shallow dehydration fluids and deeper silicate-bearing fluids efficiently remove halogens from the slab in the sub-arc region, thus favoring an efficient transfer of halogens across subduction zones.

Description: The total Antarctic sea ice extent (SIE) experiences a distinct annual cycle, peaking in September and reaching its minimum in February. In this paper we propose a mathematical and statistical decomposition of this temporal variation in SIE. Each component is interpretable and, when combined, gives a complete picture of the variation in the sea ice. We consider timescales varying from the instantaneous and not previously defined to the multi-decadal curvilinear trend, the longest. Because our representation is daily, these timescales of variability give precise information about the timing and rates of advance and retreat of the ice and may be used to diagnose physical contributors to variability in the sea ice. We define a number of annual cycles each capturing different components of variation, especially the yearly amplitude and phase that are major contributors to SIE variation. Using daily sea ice concentration data, we show that our proposed invariant annual cycle explains 29 % more of the variation in daily SIE than the traditional method. The proposed annual cycle that incorporates amplitude and phase variation explains 77 % more variation than the traditional method. The variation in phase explains more of the variability in SIE than the amplitude. Using our methodology, we show that the anomalous decay of sea ice in 2016 was associated largely with a change of phase rather than amplitude. We show that the long term trend in Antarctic sea ice extent is strongly curvilinear and the reported positive linear trend is small and dependent strongly on a positive trend that began around 2011 and continued until 2016.

Description: Aerosol–cloud interactions are the largest source of uncertainty in quantifying anthropogenic radiative forcing. The large uncertainty is, in part, due to the difficulty of predicting cloud microphysical parameters, such as the cloud droplet number concentration (Nd). Even though rigorous first-principle approaches exist to calculate Nd, the cloud and aerosol research community also relies on empirical approaches such as relating Nd to aerosol mass concentration. Here we analyze relationships between Nd and cloud water chemical composition, in addition to the effect of environmental factors on the degree of the relationships. Warm, marine, stratocumulus clouds off the California coast were sampled throughout four summer campaigns between 2011 and 2016. A total of 385 cloud water samples were collected and analyzed for 80 chemical species. Single- and multispecies log–log linear regressions were performed to predict Nd using chemical composition. Single-species regressions reveal that the species that best predicts Nd is total sulfate (Radj2=0.40). Multispecies regressions reveal that adding more species does not necessarily produce a better model, as six or more species yield regressions that are statistically insignificant. A commonality among the multispecies regressions that produce the highest correlation with Nd was that most included sulfate (either total or non-sea-salt), an ocean emissions tracer (such as sodium), and an organic tracer (such as oxalate). Binning the data according to turbulence, smoke influence, and in-cloud height allowed for examination of the effect of these environmental factors on the composition–Nd correlation. Accounting for turbulence, quantified as the standard deviation of vertical wind speed, showed that the correlation between Nd with both total sulfate and sodium increased at higher turbulence conditions, consistent with turbulence promoting the mixing between ocean surface and cloud base. Considering the influence of smoke significantly improved the correlation with Nd for two biomass burning tracer species in the study region, specifically oxalate and iron. When binning by in-cloud height, non-sea-salt sulfate and sodium correlated best with Nd at cloud top, whereas iron and oxalate correlated best with Nd at cloud base.

Description: To study the influence of complex terrain with different scales on the structure of near-surface turbulence, the turbulence observational data from Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) were analyzed. SACOL is located in typical Loess Plateau topography. The terrain around the site varies greatly with the direction. Representative data from the Northeast, Southeast, Southwest, and Northwest were selected to study the structure characteristics of the near surface turbulence. The complex topography within the flux footprint enhances the vertical scale of turbulence σw and thereby increases the vertical–longitudinal ratio of turbulence Ar. While the turbulent horizontal scale σu is also significantly affected by the distant terrain beyond the flux footprint. The upwind terrain undulation increases σu and reduces Ar. Affected by the complex terrain, the ratio of the spectrum of the vertical velocity to that of the longitudinal velocity, Sw(n)/Su(n), is far less than 4/3 in the southwest direction, and the turbulence is significantly anisotropic.

Description: Providing recommendations in cold start situations is one of the most challenging problems for collaborative filtering based recommender systems (RSs). Although user social context information has largely contributed to the cold start problem, most of the RSs still suffer from the lack of initial social links for newcomers. For this study, we are going to address this issue using a proposed user similarity detection engine (USDE). Utilizing users’ personal smart devices enables the proposed USDE to automatically extract real-world social interactions between users. Moreover, the proposed USDE uses user clustering algorithm that includes contextual information for identifying similar users based on their profiles. The dynamically updated contextual information for the user profiles helps with user similarity clustering and provides more personalized recommendations. The proposed RS is evaluated using movie recommendations as a case study. The results show that the proposed RS can improve the accuracy and personalization level of recommendations as compared to two other widely applied collaborative filtering RSs. In addition, the performance of the USDE is evaluated in different scenarios. The conducted experimental results on USDE show that the proposed USDE outperforms widely applied similarity measures in cold start and data sparsity situations.

Description: Temperature, H2O, and O3 profiles, as well as CO2, N2O, CH4, chlorofluorocarbon-12 (CFC-12), and sea surface temperature (SST) scalar anomalies are computed using a clear subset of AIRS observations over ocean for the first 16 years of NASA's Earth-Observing Satellite (EOS) Aqua Atmospheric Infrared Sounder (AIRS) operation. The AIRS Level-1c radiances are averaged over 16 d and 40 equal-area zonal bins and then converted to brightness temperature anomalies. Geophysical anomalies are retrieved from the brightness temperature anomalies using a relatively standard optimal estimation approach. The CO2, N2O, CH4, and CFC-12 anomalies are derived by applying a vertically uniform multiplicative shift to each gas in order to obtain an estimate for the gas mixing ratio. The minor-gas anomalies are compared to the National Oceanic and Atmospheric Administration (NOAA) Earth System Research Laboratory (ESRL) in situ values and used to estimate the radiometric stability of the AIRS radiances. Similarly, the retrieved SST anomalies are compared to the SST values used in the ERA-Interim reanalysis and to NOAA's Optimum Interpolation SST (OISST) product. These intercomparisons strongly suggest that many AIRS channels are stable to better than 0.02 to 0.03 K per decade, well below climate trend levels, indicating that the AIRS blackbody is not drifting. However, detailed examination of the anomaly retrieval residuals (observed – computed) shows various small unphysical shifts that correspond to AIRS hardware events (shutdowns, etc.). Some examples are given highlighting how the AIRS radiance stability could be improved, especially for channels sensitive to N2O and CH4. The AIRS shortwave channels exhibit larger drifts that make them unsuitable for climate trending, and they are avoided in this work. The AIRS Level 2 surface temperature retrievals only use shortwave channels. We summarize how these shortwave drifts impacts recently published comparisons of AIRS surface temperature trends to other surface climatologies.

Description: The CloudRoots field experiment was designed to obtain a comprehensive observational dataset that includes soil, plant, and atmospheric variables to investigate the interaction between a heterogeneous land surface and its overlying atmospheric boundary layer at the sub-hourly and sub-kilometre scale. Our findings demonstrate the need to include measurements at leaf level to better understand the relations between stomatal aperture and evapotranspiration (ET) during the growing season at the diurnal scale. Based on these observations, we obtain accurate parameters for the mechanistic representation of photosynthesis and stomatal aperture. Once the new parameters are implemented, the model reproduces the stomatal leaf conductance and the leaf-level photosynthesis satisfactorily. At the canopy scale, we find a consistent diurnal pattern on the contributions of plant transpiration and soil evaporation using different measurement techniques. From highly resolved vertical profile measurements of carbon dioxide (CO2) and other state variables, we infer a profile of the CO2 assimilation in the canopy with non-linear variations with height. Observations taken with a laser scintillometer allow us to quantify the non-steadiness of the surface turbulent fluxes during the rapid changes driven by perturbation of photosynthetically active radiation by cloud flecks. More specifically, we find 2 min delays between the cloud radiation perturbation and ET. To study the relevance of advection and surface heterogeneity for the land–atmosphere interaction, we employ a coupled surface–atmospheric conceptual model that integrates the surface and upper-air observations made at different scales from leaf to the landscape. At the landscape scale, we calculate a composite sensible heat flux by weighting measured fluxes with two different land use categories, which is consistent with the diurnal evolution of the boundary layer depth. Using sun-induced fluorescence measurements, we also quantify the spatial variability of ET and find large variations at the sub-kilometre scale around the CloudRoots site. Our study shows that throughout the entire growing season, the wide variations in stomatal opening and photosynthesis lead to large diurnal variations of plant transpiration at the leaf, plant, canopy, and landscape scales. Integrating different advanced instrumental techniques with modelling also enables us to determine variations of ET that depend on the scale where the measurement were taken and on the plant growing stage.

Description: Most available verification metrics for ensemble forecasts focus on univariate quantities. That is, they assess whether the ensemble provides an adequate representation of the forecast uncertainty about the quantity of interest at a particular location and time. For spatially indexed ensemble forecasts, however, it is also important that forecast fields reproduce the spatial structure of the observed field and represent the uncertainty about spatial properties such as the size of the area for which heavy precipitation, high winds, critical fire weather conditions, etc., are expected. In this article we study the properties of the fraction of threshold exceedance (FTE) histogram, a new diagnostic tool designed for spatially indexed ensemble forecast fields. Defined as the fraction of grid points where a prescribed threshold is exceeded, the FTE is calculated for the verification field and separately for each ensemble member. It yields a projection of a – possibly high-dimensional – multivariate quantity onto a univariate quantity that can be studied with standard tools like verification rank histograms. This projection is appealing since it reflects a spatial property that is intuitive and directly relevant in applications, though it is not obvious whether the FTE is sufficiently sensitive to misrepresentation of spatial structure in the ensemble. In a comprehensive simulation study we find that departures from uniformity of the FTE histograms can indeed be related to forecast ensembles with biased spatial variability and that these histograms detect shortcomings in the spatial structure of ensemble forecast fields that are not obvious by eye. For demonstration, FTE histograms are applied in the context of spatially downscaled ensemble precipitation forecast fields from NOAA's Global Ensemble Forecast System.

Description: On 22 March 2020, Zagreb was struck by an M5.5 earthquake that had been expected for more than 100 years and revealed all the failures in the construction of residential buildings in the Croatian capital, especially those built in the first half of the 20th century. Because of that, extensive seismological, geological, geodetic and structural engineering surveys were conducted immediately after the main shock. This study provides descriptions of damage, specifying the building performances and their correlation with the local soil characteristics, i.e., seismic motion amplification. Co-seismic vertical ground displacement was estimated, and the most affected area is identified according to Sentinel-1 interferometric wide-swath data. Finally, preliminary 3D structural modeling of the earthquake sequence was performed, and two major faults were modeled using inverse distance weight (IDW) interpolation of the grouped hypocenters. The first-order assessment of seismic amplification (due to site conditions) in the Zagreb area for the M5.5 earthquake shows that ground motions of approximately 0.16–0.19 g were amplified at least twice. The observed co-seismic deformation (based on Sentinel-1A IW SLC images) implies an approximately 3 cm uplift of the epicentral area that covers approximately 20 km2. Based on the preliminary spatial and temporal analyses of the Zagreb 2020 earthquake sequence, the main shock and the first aftershocks evidently occurred in the subsurface of the Medvednica Mountains along a deep-seated southeast-dipping thrust fault, recognized as the primary (master) fault. The co-seismic rupture propagated along the thrust towards northwest during the first half-hour of the earthquake sequence, which can be clearly seen from the time-lapse visualization. The preliminary results strongly support one of the debated models of the active tectonic setting of the Medvednica Mountains and will contribute to a better assessment of the seismic hazard for the wider Zagreb area.

Description: Micro, small, and medium-sized enterprises (MSMEs) are an essential part of economies at the national, regional, and local levels. Understanding the determinants of the development of this sector is interesting not only for researchers but also for local governments to support the development of this sector. This paper analyses micro, small, and medium enterprises at the gmina (local) level in one region, the Kujawsko-Pomorskie voivodship (NUTS2) in Poland. The authors use multivariate linear regression, spatial econometrics, and classification trees to model the influence of different factors on the number of enterprises relative to population size. The authors found that the most crucial factor in all cases, independently of the method used, is the local government’s revenue from personal income tax per capita. This finding, together with the lack of significance of variables related to the distance to technological parks or economic zones, indicates that the enterprises in the region produce mainly for local consumption and lack innovativeness. The authors also examined the influence of spatial context on the number of enterprises. The most important factor seems to be the percentage of built-up areas, but there are also others, depending on the model type; again, this confirms the local character of the activity of micro, small, and medium enterprises in the region. Variables representing the spatial context can explain the relative number of enterprises with coefficient of determination (R2) between 0.30 and 0.45, which shows that this context played a relatively significant role in the development of the MSME sector in the region. On the other hand, the econometric models (that include the neighborhood) are only significant (improving R2) for medium enterprises, which means that medium enterprises expand their activity beyond the local range.

Description: In this paper, the performance of three machine-learning methods for predicting short-term evolution and for reproducing the long-term statistics of a multiscale spatiotemporal Lorenz 96 system is examined. The methods are an echo state network (ESN, which is a type of reservoir computing; hereafter RC–ESN), a deep feed-forward artificial neural network (ANN), and a recurrent neural network (RNN) with long short-term memory (LSTM; hereafter RNN–LSTM). This Lorenz 96 system has three tiers of nonlinearly interacting variables representing slow/large-scale (X), intermediate (Y), and fast/small-scale (Z) processes. For training or testing, only X is available; Y and Z are never known or used. We show that RC–ESN substantially outperforms ANN and RNN–LSTM for short-term predictions, e.g., accurately forecasting the chaotic trajectories for hundreds of numerical solver's time steps equivalent to several Lyapunov timescales. The RNN–LSTM outperforms ANN, and both methods show some prediction skills too. Furthermore, even after losing the trajectory, data predicted by RC–ESN and RNN–LSTM have probability density functions (pdf's) that closely match the true pdf – even at the tails. The pdf of the data predicted using ANN, however, deviates from the true pdf. Implications, caveats, and applications to data-driven and data-assisted surrogate modeling of complex nonlinear dynamical systems, such as weather and climate, are discussed.

Description: In the western Mediterranean Sea, Levantine intermediate waters (LIW), which circulate below the surface productive zone, progressively accumulate nutrients along their pathway from the Tyrrhenian Sea to the Algerian Basin. This study addresses the role played by diffusion in the nutrient enrichment of the LIW, a process particularly relevant inside step-layer structures extending down to deep waters – structures known as thermohaline staircases. Profiling float observations confirmed that staircases develop over epicentral regions confined in large-scale circulation features and maintained by saltier LIW inflows on the periphery. Thanks to a high profiling frequency over the 4-year period 2013–2017, float observations reveal the temporal continuity of the layering patterns encountered during the cruise PEACETIME and document the evolution of layer properties by about +0.06 ∘C in temperature and +0.02 in salinity. In the Algerian Basin, the analysis of in situ lateral density ratios untangled double-diffusive convection as a driver of thermohaline changes inside epicentral regions and isopycnal diffusion as a driver of heat and salt exchanges with the surrounding sources. In the Tyrrhenian Sea, the nitrate flux across thermohaline staircases, as opposed to the downward salt flux, contributes up to 25 % of the total nitrate pool supplied to the LIW by vertical transfer. Overall, however, the nutrient enrichment of the LIW is driven mostly by other sources, coastal or atmospheric, as well as by inputs advected from the Algerian Basin.

Description: Geologic dating methods for the most part do not directly measure ages. Instead, interpreting a geochemical observation as a geologically useful parameter – an age or a rate – requires an interpretive middle layer of calculations and supporting data sets. These are the subject of active research and evolve rapidly, so any synoptic analysis requires repeated recalculation of large numbers of ages from a growing data set of raw observations, using a constantly improving calculation method. Many important applications of geochronology involve regional or global analyses of large and growing data sets, so this characteristic is an obstacle to progress in these applications. This paper describes the ICE-D (Informal Cosmogenic-Nuclide Exposure-age Database) database project, a prototype computational infrastructure for dealing with this obstacle in one geochronological application – cosmogenic-nuclide exposure dating – that aims to enable visualization or analysis of diverse data sets by making middle-layer calculations dynamic and transparent to the user. An important aspect of this concept is that it is designed as a forward-looking research tool rather than a backward-looking archive: only observational data (which do not become obsolete) are stored, and derived data (which become obsolete as soon as the middle-layer calculations are improved) are not stored but instead calculated dynamically at the time data are needed by an analysis application. This minimizes “lock-in” effects associated with archiving derived results subject to rapid obsolescence and allows assimilation of both new observational data and improvements to middle-layer calculations without creating additional overhead at the level of the analysis application.

Description: Moss biomonitoring technique was used for a heavy-metal pollution study in Macedonia in the framework of the International Cooperative Program on Effects of Air Pollution on Natural Vegetation and Crops (UNECE IPC Vegetation). Moss samples (n = 72) were collected during the summers of 2002, 2005, and 2010. The contents of 41 elements were determined by neutron activation analysis, atomic absorption spectrometry, and inductively coupled plasma atomic emission spectrometry. Using factor and cluster analyses, three geogenic factors were determined (Factor 1, including Al, As, Co, Cs, Fe, Hf, Na, Rb, Sc, Ta, Th, Ti, U, V, Zr, and rare-earth elements–RE; Factor 4 with Ba, K, and Sr; and Factor 5 with Br and I), one anthropogenic factor (Factor 2, including Cd, Pb, Sb, and Zn), and one geogenic-anthropogenic factor (Factor 3, including Cr and Ni). The highest anthropogenic impact of heavy metal to the air pollution in the country was from the ferronickel smelter near Kavadraci (Ni and Cr), the lead and zinc mines in the vicinity of Makedonska Kamenica, Probištip, and Kriva Palanka in the eastern part of the country (Cd, Pb, and Zn), and the former lead and zinc smelter plant in Veles. Beside the anthropogenic influences, the lithology and the composition of the soil also play an important role in the distribution of the elements.

Description: The invasion of aquifer microbial communities by aboveground microorganisms, a phenomenon known as community coalescence, is likely to be exacerbated in groundwaters fed by stormwater infiltration systems (SISs). Here, the incidence of this increased connectivity with upslope soils and impermeabilized surfaces was assessed through a meta-analysis of 16S rRNA gene libraries. Specifically, DNA sequences encoding 16S rRNA V5-V6 regions from free-living and attached aquifer bacteria (i.e., water and biofilm samples) were analysed upstream and downstream of a SIS and compared with those from bacterial communities from watershed runoffs and surface sediments from the SIS detention and infiltration basins. Significant bacterial transfers were inferred by the SourceTracker Bayesian approach, with 23 % to 57 % of the aquifer bacterial biofilms being composed of taxa from aboveground sediments and urban runoffs. Sediments from the detention basin were found more significant contributors of taxa involved in the buildup of these biofilms than soils from the infiltration basin. Inferred taxa among the coalesced biofilm community were predicted to be high in hydrocarbon degraders such as Sphingobium and Nocardia. The 16S rRNA-based bacterial community structure of the downstream-SIS aquifer waters showed lower coalescence with aboveground taxa (8 % to 38 %) than those of biofilms and higher numbers of taxa predicted to be involved in the N and S cycles. A DNA marker named tpm enabled the tracking of bacterial species from 24 genera including Pseudomonas, Aeromonas and Xanthomonas, among these communities. Several tpm sequence types were found to be shared between the aboveground and aquifer samples. Reads related to Pseudomonas were allocated to 50 species, of which 16 were found in the aquifer samples. Several of these aquifer species were found to be involved in denitrification but also hydrocarbon degradation (P. aeruginosa, P. putida and P. fluorescens). Some tpm sequence types allocated to P. umsongensis and P. chengduensis were found to be enriched among the tpm-harbouring bacteria, respectively, of the aquifer biofilms and waters. Reads related to Aeromonas were allocated to 11 species, but only those from A. caviae were recovered aboveground and in the aquifer samples. Some tpm sequence types of the X. axonopodis phytopathogen  were recorded in higher proportions among the tpm-harbouring bacteria of the aquifer waters than in the aboveground samples. A significant coalescence of microbial communities from an urban watershed with those of an aquifer was thus observed, and recent aquifer biofilms were found to be significantly colonized by runoff-opportunistic taxa able to use urban C sources from aboveground compartments.

Description: The Metropolitan Region of Sao Paulo (MRSP) is one of the main regions of Brazil that in recent years has shown an increase in the number of days with heat waves, mainly affecting the health of the most sensitive populations, such as the elderly. In this study, we identified the heat waves in the MRSP using three different definitions regarding the maximum daily temperature threshold. To analyze the impact of heat waves on elderly mortality, we used distributed lag nonlinear models (dlnm) and we quantified the heat wave-related excess mortality of elderly people from 1985 to 2005 and made projections for the near future (2030 to 2050) and the distant future (2079–2099) under the climate change scenarios RCP4.5 and RCP8.5 (RCP: Representative Concentration Paths). An important aspect of this research is that for the projections we take into account two assumptions: non-adaptation and adaptation to the future climate. Our projections show that the heat wave-related excess of elderly mortality will increase in the future, being highest when we consider no adaptation, mainly from cardiovascular diseases in women (up to 587 deaths per 100,000 inhabitants per year). This study can be used for public policies to implement preventive and adaptive measures in the MRSP.

Description: Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics.

Description: Landfast sea ice (fast ice) attached to Antarctic (near-)coastal elements is a critical component of the local physical and ecological systems. Through its direct coupling with the atmosphere and ocean, fast-ice properties are also a potential indicator of processes related to a changing climate. However, in situ fast-ice observations in Antarctica are extremely sparse because of logistical challenges and harsh environmental conditions. Since 2010, a monitoring program observing the seasonal evolution of fast ice in Atka Bay has been conducted as part of the Antarctic Fast Ice Network (AFIN). The bay is located on the northeastern edge of Ekström Ice Shelf in the eastern Weddell Sea, close to the German wintering station Neumayer III. A number of sampling sites have been regularly revisited each year between annual ice formation and breakup to obtain a continuous record of sea-ice and sub-ice platelet-layer thickness, as well as snow depth and freeboard across the bay. Here, we present the time series of these measurements over the last 9 years. Combining them with observations from the nearby Neumayer III meteorological observatory as well as auxiliary satellite images enables us to relate the seasonal and interannual fast-ice cycle to the factors that influence their evolution. On average, the annual consolidated fast-ice thickness at the end of the growth season is about 2 m, with a loose platelet layer of 4 m thickness beneath and 0.70 m thick snow on top. Results highlight the predominately seasonal character of the fast-ice regime in Atka Bay without a significant interannual trend in any of the observed variables over the 9-year observation period. Also, no changes are evident when comparing with sporadic measurements in the 1980s and 1990s. It is shown that strong easterly winds in the area govern the year-round snow distribution and also trigger the breakup of fast ice in the bay during summer months. Due to the substantial snow accumulation on the fast ice, a characteristic feature is frequent negative freeboard, associated flooding of the snow–ice interface, and a likely subsequent snow ice formation. The buoyant platelet layer beneath negates the snow weight to some extent, but snow thermodynamics is identified as the main driver of the energy and mass budgets for the fast-ice cover in Atka Bay. The new knowledge of the seasonal and interannual variability of fast-ice properties from the present study helps to improve our understanding of interactions between atmosphere, fast ice, ocean, and ice shelves in one of the key regions of Antarctica and calls for intensified multidisciplinary studies in this region.

Description: Chemistry plays an indispensable role in investigations of the atmosphere; however, many climate models either ignore or greatly simplify atmospheric chemistry, limiting both their accuracy and their scope. We present the development and evaluation of the online global atmospheric chemical model BCC-GEOS-Chem v1.0, coupling the GEOS-Chem chemical transport model (CTM) as an atmospheric chemistry component in the Beijing Climate Center atmospheric general circulation model (BCC-AGCM). The GEOS-Chem atmospheric chemistry component includes detailed tropospheric HOx–NOx–volatile organic compounds–ozone–bromine–aerosol chemistry and online dry and wet deposition schemes. We then demonstrate the new capabilities of BCC-GEOS-Chem v1.0 relative to the base BCC-AGCM model through a 3-year (2012–2014) simulation with anthropogenic emissions from the Community Emissions Data System (CEDS) used in the Coupled Model Intercomparison Project Phase 6 (CMIP6). The model captures well the spatial distributions and seasonal variations in tropospheric ozone, with seasonal mean biases of 0.4–2.2 ppbv at 700–400 hPa compared to satellite observations and within 10 ppbv at the surface to 500 hPa compared to global ozonesonde observations. The model has larger high-ozone biases over the tropics which we attribute to an overestimate of ozone chemical production. It underestimates ozone in the upper troposphere which is likely due either to the use of a simplified stratospheric ozone scheme or to biases in estimated stratosphere–troposphere exchange dynamics. The model diagnoses the global tropospheric ozone burden, OH concentration, and methane chemical lifetime to be 336 Tg, 1.16×106 molecule cm−3, and 8.3 years, respectively, which is consistent with recent multimodel assessments. The spatiotemporal distributions of NO2, CO, SO2, CH2O, and aerosol optical depth are generally in agreement with satellite observations. The development of BCC-GEOS-Chem v1.0 represents an important step for the development of fully coupled earth system models (ESMs) in China.

Description: We report the results of a multiproxy study that combines structural analysis of a fracture–stylolite network and isotopic characterization of calcite vein cements and/or fault coating. Together with new paleopiezometric and radiometric constraints on burial evolution and deformation timing, these results provide a first-order picture of the regional fluid systems and pathways that were present during the main stages of contraction in the Tuscan Nappe and Umbria–Marche Apennine Ridge (northern Apennines). We reconstruct four steps of deformation at the scale of the belt: burial-related stylolitization, Apenninic-related layer-parallel shortening with a contraction trending NE–SW, local extension related to folding, and late-stage fold tightening under a contraction still striking NE–SW. We combine the paleopiezometric inversion of the roughness of sedimentary stylolites – that constrains the range of burial depth of strata prior to layer-parallel shortening – with burial models and U–Pb absolute dating of fault coatings in order to determine the timing of development of mesostructures. In the western part of the ridge, layer-parallel shortening started in Langhian time (∼15 Ma), and then folding started at Tortonian time (∼8 Ma); late-stage fold tightening started by the early Pliocene (∼5 Ma) and likely lasted until recent/modern extension occurred (∼3 Ma onward). The textural and geochemical (δ18O, δ13C, Δ47CO2 and 87Sr∕86Sr) study of calcite vein cements and fault coatings reveals that most of the fluids involved in the belt during deformation either are local or flowed laterally from the same reservoir. However, the western edge of the ridge recorded pulses of eastward migration of hydrothermal fluids (〉140 ∘C), driven by the tectonic contraction and by the difference in structural style of the subsurface between the eastern Tuscan Nappe and the Umbria–Marche Apennine Ridge.

Description: We developed the WRF-GC model, an online coupling of the Weather Research and Forecasting (WRF) mesoscale meteorological model and the GEOS-Chem atmospheric chemistry model, for regional atmospheric chemistry and air quality modeling. WRF and GEOS-Chem are both open-source community models. WRF-GC offers regional modellers access to the latest GEOS-Chem chemical module, which is state of the science, well documented, traceable, benchmarked, actively developed by a large international user base, and centrally managed by a dedicated support team. At the same time, WRF-GC enables GEOS-Chem users to perform high-resolution forecasts and hindcasts for any region and time of interest. WRF-GC uses unmodified copies of WRF and GEOS-Chem from their respective sources; the coupling structure allows future versions of either one of the two parent models to be integrated into WRF-GC with relative ease. Within WRF-GC, the physical and chemical state variables are managed in distributed memory and translated between WRF and GEOS-Chem by the WRF-GC coupler at runtime. We used the WRF-GC model to simulate surface PM2.5 concentrations over China during 22 to 27 January 2015 and compared the results to surface observations and the outcomes from a GEOS-Chem Classic nested-China simulation. Both models were able to reproduce the observed spatiotemporal variations of regional PM2.5, but the WRF-GC model (r=0.68, bias =29 %) reproduced the observed daily PM2.5 concentrations over eastern China better than the GEOS-Chem Classic model did (r=0.72, bias =55 %). This was because the WRF-GC simulation, nudged with surface and upper-level meteorological observations, was able to better represent the pollution meteorology during the study period. The WRF-GC model is parallelized across computational cores and scales well on massively parallel architectures. In our tests where the two models were similarly configured, the WRF-GC simulation was 3 times more efficient than the GEOS-Chem Classic nested-grid simulation due to the efficient transport algorithm and the Message Passing Interface (MPI)-based parallelization provided by the WRF software framework. WRF-GC v1.0 supports one-way coupling only, using WRF-simulated meteorological fields to drive GEOS-Chem with no chemical feedbacks. The development of two-way coupling capabilities, i.e., the ability to simulate radiative and microphysical feedbacks of chemistry to meteorology, is under way. The WRF-GC model is open source and freely available from http://wrf.geos-chem.org (last access: 10 July 2020).

Description: Mixed-method bicultural research in Aotearoa New Zealand, including the weaving of Indigenous and other knowledge, is emerging within many academic disciplines. However, mātauranga Māori (the knowledge, culture, values, and world view of the Indigenous peoples of Aotearoa New Zealand) and Te Ao Māori (the Māori world) is poorly represented within geomorphological investigations. Here, we review international efforts to include Indigenous knowledge in geologic and geomorphic studies and provide an overview of the current state of mātauranga Māori within research endeavours in Aotearoa New Zealand. We review three theoretical frameworks (i.e. methodologies) for including mātauranga Māori in research projects and three models (i.e. methods) for including Māori values within research. We identify direct benefits to geomorphology and discuss how these frameworks and models can be adapted for use with Indigenous knowledge systems outside of Aotearoa New Zealand. The aim of this review is to encourage geomorphologists around the world to engage with local Indigenous peoples to develop new approaches to geomorphic research. In Aotearoa New Zealand, we hope to inspire geomorphologists to embark on research journeys in genuine partnership with Māori that promote toitū te mātauranga – the enduring protection, promotion and respect of mātauranga Māori.

Description: Suspended sediment load in rivers is highly uncertain because sediment production and transport at catchment scale are strongly variable in space and time, and they are affected by catchment hydrology, topography, and land cover. Among the main sources of this variability are the spatially distributed nature of overland flow as an erosion driver and of surface erodibility given by soil type and vegetation cover distribution. Temporal variability mainly results from the time sequence of rainfall intensity during storms and snowmelt leading to soil saturation and overland flow. We present a new spatially distributed soil erosion and suspended sediment transport module integrated into the computationally efficient physically based hydrological model TOPKAPI-ETH, with which we investigate the effects of the two erosion drivers – precipitation and surface erodibility – on catchment sediment fluxes in a typical pre-Alpine mesoscale catchment. By conducting a series of numerical experiments, we quantify the impact of spatial variability in the two key erosion drivers on erosion–deposition patterns, sediment delivery ratio, and catchment sediment yields. Main findings are that the spatial variability in erosion drivers affects sediment yield by (i) increasing sediment production due to a spatially variable precipitation, while decreasing it due to a spatially variable surface erodibility, (ii) favouring the clustering of sediment source areas in space by surface runoff generation, and (iii) decreasing their connectivity to the river network by magnifying sediment buffers. The results highlight the importance of resolving spatial gradients controlling hydrology and sediment processes when modelling sediment dynamics at the mesoscale, in order to capture the key effects of sediment sources, buffers, and hillslope hydrological pathways in determining the sediment signal.

Description: Agricultural production in sub-Saharan Africa remains dependent on high inputs of human labor, a situation associated with direct exposure to daylight heat during critical periods of the agricultural calendar. We ask the question: how is the Wet-Bulb Globe Temperature (WBGT) going to be distributed in the future, and how will this affect the ability of smallholder farmers to perform agricultural activities? Data from general circulation models are used to estimate the distribution of WBGT in 2000, 2050 and 2100, and for high activity periods in the agricultural calendar. The distribution of WBGT is divided into recommended maximum WBGT exposure levels (°C) at different work intensities, and rest/work ratios for an average acclimatized worker wearing light clothing (ISO, 18). High WBGTs are observed during the two periods of the East African. In February to March, eastern and coastal regions of Kenya and Tanzania witness high WBGT values—some necessitating up to 75% rest/hour work intensities in 2050 and 2100. In August to September, eastern and northern Kenya and north and central Uganda are vulnerable to high WBGT values. Designing policies to address this key challenge is a critical element in adaptation methods to address the impact of climate change.

Description: Micrometeorites (MMs) are small particles that account for most of the extraterrestrial material deposited on Earth. Synchrotron X-ray fluorescence and diffraction allowed for chemical and mineral characterization to distinguish MM from atmospheric particulate. The relative components of iron, nickel, and other elements were considered in the identification of ferrous MM while high amounts of titanium were considered an indication that the particles were of atmospheric origin. Out of 100 samples collected by high school students and teachers, eight were taken to a synchrotron for analysis. Of those eight, three exhibited extraterrestrial compositions. X-ray absorption near-edge structure analysis revealed that the same three samples contained sulfide, the main sulfur form constituent in MM. X-ray microdiffraction analysis showed the presence of the minerals pentlandite and forsterite. Collectively, these results support the extraterrestrial nature of the three particles.

Description: The Yiliu tungsten polymetallic deposit, located in the south central portion of the Nanling nonferrous metal metallogenic province in South China, is an area with common Yanshanian tectonothermal events. Early Yanshanian magmatism leads to the emplacement of voluminous tungsten-bearing granite intrusions, such as the Baoshan, Benggangling and Junye plutons, which are considered temporally and spatially associated with W-polymetallic mineralization in the Yiliu region. Here, we investigate the basic geological and petrological characteristics of the Junye granites, and present major and trace element geochemical data and bulk-rock Sr-Nd-Pb-Hf isotopic data to gain insight into the petrogenesis and tectonic setting of granitic intrusions in the region. The Junye granites are high-K calc-alkaline and metaluminous to weakly peraluminous [A/CNK = molar ratios of Al2O3/(CaO + Na2O + K2O) = 0.97–1.02] with enrichment in SiO2 (75.68–76.44 wt.%), relatively high total alkalis (K2O + Na2O = 8.06–8.45 wt.%) with K2O/Na2O ratios ranging from 1.12 to 1.42, and moderate Al2O3 (12.62–13.00 wt.%), but low in P2O5 (

Description: We used observational data and the results from a high-resolution numerical simulation model to analyze the occurrence and development of an extreme precipitation event in the Ili Valley, Xinjiang, China on 26 June 2015. We analyzed the horizontal wavelength, period, speed, ducting, energy propagation and feedback mechanism of inertial gravity waves. A low-level convergence line was formed in the valley by the northerly and westerly winds as a result of Central Asian vortices and the trumpet-shaped topography of the Ili Valley. There was sufficient water vapor in the valley for the precipitation event to develop. A mesoscale vortex formed and developed on the low-level convergence line and the rainfall was distributed either near the convergence line or the mesoscale vortex. The low-level convergence line and the uplift caused by the terrain triggered convection, and then the convection triggered waves at lower levels. The combination of ascending motion induced by the lower level waves and the mesoscale vortex led to the development of convection, causing the precipitation to intensify. When the convection moved eastward to Gongliu County, it was coupled with the ascending phase of upper level waves, causing both the convection and precipitation to intensify again. We applied spectral analysis methods to verify that the waves were inertial gravity waves. The upper level inertial gravity waves propagated westward at a mean speed of −12 m s−1 with periods of 73–179 min and horizontal wavelengths of 50–55 km. The lower level inertial gravity waves propagated eastward at a mean speed of 8 m s−1 with periods of 73–200 min and a horizontal wavelength of 85 km. The more (less) favorable waveguide conditions determined whether the gravity waves persisted for a long (short) time and propagated for a longer (shorter) distance. Based on the mesoscale Eliassen–Palm flux theory, the wave energy of inertial gravity waves had an important effect on the maintenance and development of convection and precipitation by affecting wind strength and wind divergence. Feedback was mainly through the meridional and vertical transport of zonal momentum and the meridional transport of heat.

Description: In convergent zones, several parts of the geodynamic system (e.g., continental margins, back-arc regions) can be deformed, uplifted, and eroded through time, each of them potentially delivering clastic sediments to neighboring basins. Tectonically driven events are mostly recorded in syntectonic clastic systems accumulated into different kinds of basins: trench, fore-arc, and back-arc basins in subduction zones and foredeep, thrust-top, and episutural basins in collisional settings. The most widely used tools for provenance analysis of synorogenic sediments and for unraveling the tectonic evolution of convergent zones are sandstone petrography and U–Pb dating of detrital zircon. In this paper, we present a comparison of previously published data discussing how these techniques are used to constrain provenance reconstructions and contribute to a better understanding of the tectonic evolution of (i) the Cretaceous transition from extensional to compressional regimes in the back-arc region of the southern Andean system; and (ii) the involvement of the passive European continental margin in the Western Alps subduction system during impending Alpine collision. In both cases, sediments delivered from the down-bending continental block are significantly involved. Our findings highlight its role as a detrital source, which is generally underestimated or even ignored in current tectonic models.

Description: Urban air pollution from gaseous pollutants is a growing public health problem in many countries including South Africa. Examining the levels, trends and health risk of exposure to ambient gaseous pollutants will assist in understanding the effectiveness of existing control measures and plan for suitable management strategies. This study determined the concentration levels and non-cancer risk of CO, SO2, NO2, and O3 at an industrial area in Pretoria West, South Africa. We utilised a set of secondary data for CO, NO2, SO2, and O3 that was obtained from a monitoring station. Analysis of the hourly monitored data was done. Their non-cancer risk (HQ) was determined using the human health risk assessment model for different age categories. The annual levels of NO2 (39.442 µg/m3), SO2 (22.464 µg/m3), CO (722.003 µg/m3) and the 8-hour concentration of CO (649.902 µg/m3) and O3 (33.556 µg/m3) did not exceed the South African National Ambient Air Quality Standards for each pollutant. The HQ for each pollutant across exposed groups (except children) was less than 1. This indicates that the recorded levels could not pose non-cancer risk to susceptible individuals.

Description: We report on the application of the U-Th-He method for the direct dating of pyrite and provide an original methodological approach for measurement of U, Th and He in single grains without loss of parent nuclides during thermal extraction of He. The U-Th-He age of ten samples of high-crystalline stoichiometric pyrite from unoxidized massive ores of the Uzelga volcanogenic massive sulfide (VMS) deposit, South Urals, is 382 ± 12 Ma (2σ) (U concentrations ~1–5 ppm; 4He ~10−4 cm3 STP g−1). This age is consistent with independent (biostratigraphic) estimations of the age of ore formation (ca, 389–380 Ma) and is remarkably older than the probable age of the regional prehnite-pumpellyite facies metamorphism (~340–345 Ma). Our results indicate that the U-Th-He dating of ~1 mg weight pyrite sample is possible and open new perspectives for the dating of ore deposits. The relative simplicity of U-Th-He dating in comparison with other geochronological methods makes this approach interesting for further application.

Description: Human activity events are often recorded with their geographic locations and temporal stamps, which form spatial patterns of the events during individual time periods. Temporal attributes of these events help us understand the evolution of spatial processes over time. A challenge that researchers still face is that existing methods tend to treat all events as the same when evaluating the spatiotemporal pattern of events that have different properties. This article suggests a method for assessing the level of spatiotemporal clustering or spatiotemporal autocorrelation that may exist in a set of human activity events when they are associated with different categorical attributes. This method extends the Voronoi structure from 2D to 3D and integrates a sliding-window model as an approach to spatiotemporal tessellations of a space-time volume defined by a study area and time period. Furthermore, an index was developed to evaluate the partial spatiotemporal clustering level of one of the two event categories against the other category. The proposed method was applied to simulated data and a real-world dataset as a case study. Experimental results show that the method effectively measures the level of spatiotemporal clustering patterns among human activity events of multiple categories. The method can be applied to the analysis of large volumes of human activity events because of its computational efficiency.

Description: During February 2019, two severe storms affected the island of Crete, located in south Greece. Both storms produced excessive rainfall, provoking severe damages, especially in the western part of Crete. The role of the prevailing synoptic patterns and the interaction of the flow with the high mountains of Crete were investigated. For this purpose, a variety of observational and numerical model data were exploited, including data from a dense rain gauge network, satellite imagery, and model analysis of various parameters describing the stability of the impinging flow. The first storm was a long-lasting event, with convective outbreaks embedded in a more stratiform rainfall pattern. The second storm was brief but mostly convection dominated. The analysis of the available data underlined the role of the low-level convergence upstream of the mountains during both storms, highlighting similarities and differences, as well as the role of the stability of the impinging flow. High soil moisture content was also evidenced as a key ingredient for the severe flooding that occurred during the second storm. This work complements similar studies on the role of Mediterranean islands and their topography on the spatial and temporal distribution of extreme rainfall.