Refining data–data and data–model vegetation comparisons using the Earth mover's distance (EMD)

Comparing temporal and spatial vegetation changes between reconstructions or between reconstructions and model simulations requires carefully selecting an appropriate evaluation metric. A common way of comparing reconstructed and simulated vegetation changes involves measuring the agreement between pollen- or model-derived unary vegetation estimates, such as the biome or plant functional type (PFT) with the highest affinity scores. While this approach based on summarising the vegetation signal into unary vegetation estimates performs well in general, it overlooks the details of the underlying vegetation structure. However, this underlying data structure can influence conclusions since minor variations in pollen percentages modify which biome or PFT has the highest affinity score (i.e. modify the unary vegetation estimate). To overcome this limitation, we propose using the Earth mover's distance (EMD) to quantify the mismatch between vegetation distributions such as biome or PFT affinity scores. The EMD circumvents the issue of summarising the data into unary biome or PFT estimates by considering the entire range of biome or PFT affinity scores to calculate a distance between the compared entities. In addition, each type of mismatch can be given a specific weight to account for case-specific ecological distances or, said differently, to account for the fact that reconstructing a temperate forest instead of a boreal forest is ecologically more coherent than reconstructing a temperate forest instead of a desert. We also introduce two EMD-based statistical tests that determine (1) if the similarity of two samples is significantly better than a random association given a particular context and (2) if the pairing between two datasets is better than might be expected by chance. To illustrate the potential and the advantages of the EMD as well as the tests in vegetation comparison studies, we reproduce different case studies based on previously published simulated and reconstructed biome changes for Europe and capitalise on the advantages of the EMD to refine the interpretations of past vegetation changes by highlighting that flickering unary estimates, which give an impression of high vegetation instability, can correspond to gradual vegetation changes with low EMD values between contiguous samples (case study 1). We also reproduce data–model comparisons for five specific time slices to identify those that are statistically more robust than a random agreement while accounting for the underlying vegetation structure of each pollen sample (case study 2). The EMD and the statistical tests are included in the paleotools R package (https://github.com/mchevalier2/paleotools, last access: 3 May 2023).</p

Biome changes and their inferred climatic drivers in northern and eastern continental Asia at selected times since 40 cal ka bp

© 2017 Springer-Verlag GmbH Germany, part of Springer Nature Recent global warming is pronounced in high-latitude regions (e.g. northern Asia), and will cause the vegetation to change. Future vegetation trends (e.g. the “arctic greening”) will feed back into atmospheric circulation and the global climate system. Understanding the nature and causes of past vegetation changes is important for predicting the composition and distribution of future vegetation communities. Fossil pollen records from 468 sites in northern and eastern Asia were biomised at selected times between 40 cal ka bp and today. Biomes were also simulated using a climate-driven biome model and results from the two approaches compared in order to help understand the mechanisms behind the observed vegetation changes. The consistent biome results inferred by both approaches reveal that long-term and broad-scale vegetation patterns reflect global- to hemispheric-scale climate changes. Forest biomes increase around the beginning of the late deglaciation, become more widespread during the early and middle Holocene, and decrease in the late Holocene in fringe areas of the Asian Summer Monsoon. At the southern and southwestern margins of the taiga, forest increases in the early Holocene and shows notable species succession, which m ay have been caused by winter warming at ca. 7 cal ka bp. At the northeastern taiga margin (central Yakutia and northeastern Siberia), shrub expansion during the last deglaciation appears to prevent the permafrost from thawing and hinders the northward expansion of evergreen needle-leaved species until ca. 7 cal ka bp. The vegetation-climate disequilibrium during the early Holocene in the taiga-tundra transition zone suggests that projected climate warming will not cause a northward expansion of evergreen needle-leaved species

LegacyClimate 1.0: a dataset of pollen-based climate reconstructions from 2594 Northern Hemisphere sites covering the last 30&thinsp;kyr and beyond

Here we describe LegacyClimate 1.0, a dataset of the reconstruction of the mean July temperature (TJuly), mean annual temperature (Tann), and annual precipitation (Pann) from 2594 fossil pollen records from the Northern Hemisphere, spanning the entire Holocene, with some records reaching back to the Last Glacial Period. Two reconstruction methods, the modern analog technique (MAT) and weighted averaging partial least squares regression (WA-PLS), reveal similar results regarding spatial and temporal patterns. To reduce the impact of precipitation on temperature reconstruction, and vice versa, we also provide reconstructions using tailored modern pollen data, limiting the range of the corresponding other climate variables. We assess the reliability of the reconstructions, using information from the spatial distributions of the root mean squared error in the prediction and reconstruction significance tests. The dataset is beneficial for synthesis studies of proxy-based reconstructions and to evaluate the output of climate models and thus help to improve the models themselves. We provide our compilation of reconstructed TJuly, Tann, and Pann as open-access datasets at PANGAEA (https://doi.org/10.1594/PANGAEA.930512; Herzschuh et al., 2023a). The R code for the reconstructions is provided at Zenodo (https://doi.org/10.5281/zenodo.7887565; Herzschuh et al., 2023b), including the harmonized open-access modern and fossil datasets used for the reconstructions, so that customized reconstructions can be easily established.</p

The Aguablanca Ni–(Cu) sulfide deposit, SW Spain: geologic and geochemical controls and the relationship with a midcrustal layered mafic complex

The Aguablanca Ni–(Cu) sulfide deposit is hosted by a breccia pipe within a gabbro–diorite pluton. The deposit probably formed due to the disruption of a partially crystallized layered mafic complex at about 12– 19 km depth and the subsequent emplacement of melts and breccias at shallow levels (<2 km). The ore-hosting breccias are interpreted as fragments of an ultramafic cumulate, which were transported to the near surface along with a molten sulfide melt. Phlogopite Ar–Ar ages are 341– 332 Ma in the breccia pipe, and 338–334 Ma in the layered mafic complex, and are similar to recently reported U–Pb ages of the host Aguablanca Stock and other nearby calcalkaline metaluminous intrusions (ca. 350–330 Ma). Ore deposition resulted from the combination of two critical factors, the emplacement of a layered mafic complex deep in the continental crust and the development of small dilational structures along transcrustal strike-slip faults that triggered the forceful intrusion of magmas to shallow levels. The emplacement of basaltic magmas in the lower middle crust was accompanied by major interaction with the host rocks, immiscibility of a sulfide melt, and the formation of a magma chamber with ultramafic cumulates and sulfide melt at the bottom and a vertically zoned mafic to intermediate magmas above. Dismembered bodies of mafic/ultramafic rocks thought to be parts of the complex crop out about 50 km southwest of the deposit in a tectonically uplifted block (Cortegana Igneous Complex, Aracena Massif). Reactivation of Variscan structures that merged at the depth of the mafic complex led to sequential extraction of melts, cumulates, and sulfide magma. Lithogeochemistry and Sr and Nd isotope data of the Aguablanca Stock reflect the mixing from two distinct reservoirs, i.e., an evolved siliciclastic middle-upper continental crust and a primitive tholeiitic melt. Crustal contamination in the deep magma chamber was so intense that orthopyroxene replaced olivine as the main mineral phase controlling the early fractional crystallization of the melt. Geochemical evidence includes enrichment in SiO2 and incompatible elements, and Sr and Nd isotope compositions (87Sr/86Sri 0.708–0.710; 143Nd/144Ndi 0.512–0.513). However, rocks of the Cortegana Igneous Complex have low initial 87Sr/86Sr and high initial 143Nd/144Nd values suggesting contamination by lower crustal rocks. Comparison of the geochemical and geological features of igneous rocks in the Aguablanca deposit and the Cortegana Igneous Complex indicates that, although probably part of the same magmatic system, they are rather different and the rocks of the Cortegana Igneous Complex were not the direct source of the Aguablanca deposit. Crust–magma interaction was a complex process, and the generation of orebodies was controlled by local but highly variable factors. The model for the formation of the Aguablanca deposit presented in this study implies that dense sulfide melts can effectively travel long distances through the continental crust and that dilational zones within compressional belts can effectively focus such melt transport into shallow environments

The granite‑hosted Variscan gold deposit from Santo António mine in the Iberian Massif (Penedono, NW Portugal): constraints from mineral chemistry, fuid inclusions, sulfur and noble gases isotopes

The study area is located in the Central Iberian Zone, a major tectonic unit of the Iberian Massif (Variscan belt). In this region the basement is composed of Cambrian-Ordovician sedimentary and minor volcanic rocks that underwent deformation and metamorphism during the Carboniferous. These metamorphic rocks host ca. 331–308 Ma granitic plutons emplaced during the D2 extensional and D3–D4 contractional deformation phases. The gold-bearing quartz veins from the Santo António mine (Penedono region) occur in granite formed at 310.1 ± 1.1 Ma and post-dated the peak of metamorphism. Gold–silver alloy is included in quartz, but mainly occurs in spaces between grains or micro-fractures within arsenopyrite of all three generations and less in pyrite. Late sulphides and sulphosalts were deposited along fractures mainly in arsenopyrite, and locally surrounding the gold–silver alloy grains. Ferberite, scheelite and stolzite replace arsenopyrite. The abundant aqueous carbonic fluids and the occurrence of a low-salinity fluid and their minimum possible entrapment temperature of 360–380 °C suggest that this gold-forming event began during the waning stages of the Variscan orogeny. The mean δ34S values of arsenopyrite and pyrite are − 4.7‰ and − 3.8‰, respectively. He–Ar–Ne isotopic data suggest a crustal origin. The ascent of the granite magma has provided the heat for remobilization of gold, other metals and metalloids from the metamorphic rocks. This gold-arsenopyrite deposit has thus similar characteristics as other selected gold-arsenopyrite deposits from the Iberian Massif, but it contains tungstates.El área de estudio está ubicada en la Zona Centroibérica, una importante unidad tectónica del Macizo Ibérico (cinturón varisco). En esta región el basamento está compuesto por rocas sedimentarias y volcánicas del Cámbrico-Ordovícico tectonizadas y metamorfzadas durante el Carbonífero. Estas rocas metamórfcas sirven como caja de los plutones graníticos datados en torno a 331–308 Ma y que fueron emplazados durante la fase de deformación extensional D2 y las fases de deformación contraccional D3 y D4. Las venas de cuarzo ricas en oro de la mina de Santo António (región de Penedono) que aparecen en un granito datado a los 310.1 ± 1.1 Ma son posteriores al pico metamórfco regional. La aleación de oro y plata se incluye en el cuarzo, pero se produce principalmente en los espacios entre granos o micro-fracturas dentro de arsenopirita de las tres generaciones y menos en pirita. Los sulfuros y sulfuros tardíos se depositaron a lo largo de las fracturas principalmente en arsenopirita, y alrededor de los granos de aleación de oro y plata. Ferberita, scheelita y la estolzita sustituyen a la arsenopirita. Los abundantes líquidos acuosos carbónicos y la presencia de un fuido de baja salinidad y su posible temperatura de atrapamiento mínima en torno de 360-380 ºC sugieren que este evento de formación de oro comenzó durante las etapas fnales de la orogenia varisca. Los valores medios de S de arsenopirita y pirita son − 4.7 ‰ y − 3.8 ‰, respectivamente. Los datos isotópicos de He–Ar–Ne sugieren que en el origen de los fuidos mineralizados participa la corteza continental. El ascenso del magma granítico ha provisto el calor para la movilización del oro, otros metales y metaloides desde las rocas metamórfcas. Este depósito de oroarsenopirita tiene así características similares a otros yaciamientos con arsenopirita y oro del Macizo Ibérico, pero sin embargo contienen tungstates.This research was financially supported by Fundação para a Ciência e Tecnologia through the projects GOLDGranites, Orogenesis, Long-term strain/stress and Deposition of ore metals—PTDC/GEO-GEO/2446/2012: COMPETE: FCOMP-01-0124-FEDER-029192 and UID/GEO/04035/2013

Semeadura de milho em solo sob preparo reduzido

Para avaliar a viabilidade da implantação da cultura do milho em solo sob preparo reduzido, em dezembro de 1988, instalou-se um experimento em um Latossolo Vermelho-Escuro distrófico, no município de Passo Fundo (RS). Na testemunha, o preparo do solo foi realizado por meio de uma escarificação mais uma gradagem com grade leve. Nos demais cinco tratamentos, o solo foi preparado por meio de uma operação conjugada (escarificador equipado com rolo destorroador). Na testemunha e em um dos tratamentos sob preparo conjugado, não foram realizadas adaptações na semeadora-adubadora. Nos demais tratamentos, foram adaptados discos de corte ondulados, ou secções de rolos destorroadores, na frente das linhas da semeadora. A percentagem de cobertura do solo, dois dias após a semeadura, foi significativamente menor (Tukey ao nível de 5%) na testemunha, enquanto nos demais tratamentos, caracterizados como conservacionistas, não se perceberam diferenças. O diâmetro médio geométrico (DMG) dos agregados do solo não apresentou diferenças significativas entre os tratamentos, porém os agregados coletados nas linhas apresentaram DMG significativamente menores do que aqueles encontrados nas entrelinhas. A velocidade de absorção de água pelas sementes, o comprimento médio ponderado das radículas, o índice de velocidade de emergência e o estande inicial de plantas de milho não variaram significativamente em função dos tratamentos. Os resultados evidenciam que o preparo conjugado do solo por meio de um escarificador equipado com um cilindro destorroador permitiu implantar a cultura do milho, sem a necessidade de operações de preparo secundário (gradagens) e sem requerer adaptações na semeadora-adubadora