Refine
Has Fulltext
- no (237)
Year of publication
- 2018 (237) (remove)
Document Type
- Article (237) (remove)
Language
- English (237)
Is part of the Bibliography
- yes (237)
Keywords
- Himalaya (5)
- Germany (4)
- climate change (4)
- permafrost (4)
- Palynology (3)
- Quaternary (3)
- landslides (3)
- Arctic (2)
- Australia (2)
- Chile (2)
Institute
- Institut für Geowissenschaften (237) (remove)
High-latitude treeless ecosystems represent spatially highly heterogeneous landscapes with small net carbon fluxes and a short growing season. Reliable observations and process understanding are critical for projections of the carbon balance of the climate-sensitive tundra. Space-borne remote sensing is the only tool to obtain spatially continuous and temporally resolved information on vegetation greenness and activity in remote circumpolar areas. However, confounding effects from persistent clouds, low sun elevation angles, numerous lakes, widespread surface inundation, and the sparseness of the vegetation render it highly challenging. Here, we conduct an extensive analysis of the timing of peak vegetation productivity as shown by satellite observations of complementary indicators of plant greenness and photosynthesis. We choose to focus on productivity during the peak of the growing season, as it importantly affects the total annual carbon uptake. The suite of indicators are as follows: (1) MODIS-based vegetation indices (VIs) as proxies for the fraction of incident photosynthetically active radiation (PAR) that is absorbed (fPAR), (2) VIs combined with estimates of PAR as a proxy of the total absorbed radiation (APAR), (3) sun-induced chlorophyll fluorescence (SIF) serving as a proxy for photosynthesis, (4) vegetation optical depth (VOD), indicative of total water content and (5) empirically upscaled modelled gross primary productivity (GPP). Averaged over the pan-Arctic we find a clear order of the annual peak as APAR <= GPP < SIF < VIs/VOD. SIF as an indicator of photosynthesis is maximised around the time of highest annual temperatures. The modelled GPP peaks at a similar time to APAR. The time lag of the annual peak between APAR and instantaneous SIF fluxes indicates that the SIF data do contain information on light-use efficiency of tundra vegetation, but further detailed studies are necessary to verify this. Delayed peak greenness compared to peak photosynthesis is consistently found across years and land-cover classes. A particularly late peak of the normalised difference vegetation index (NDVI) in regions with very small seasonality in greenness and a high amount of lakes probably originates from artefacts. Given the very short growing season in circumpolar areas, the average time difference in maximum annual photosynthetic activity and greenness or growth of 3 to 25 days (depending on the data sets chosen) is important and needs to be considered when using satellite observations as drivers in vegetation models.
The Gutenberg-Richter relation for earthquake magnitudes is the most famous empirical law in seismology. It states that the frequency of earthquake magnitudes follows an exponential distribution; this has been found to be a robust feature of seismicity above the completeness magnitude, and it is independent of whether global, regional, or local seismicity is analyzed. However, the exponent b of the distribution varies significantly in space and time, which is important for process understanding and seismic hazard assessment; this is particularly true because of the fact that the Gutenberg-Richter b-value acts as a proxy for the stress state and quantifies the ratio of large-to-small earthquakes. In our work, we focus on the automatic detection of statistically significant temporal changes of the b-value in seismicity data. In our approach, we use Bayes factors for model selection and estimate multiple change-points of the frequency-magnitude distribution in time. The method is first applied to synthetic data, showing its capability to detect change-points as function of the size of the sample and the b-value contrast. Finally, we apply this approach to examples of observational data sets for which b-value changes have previously been stated. Our analysis of foreshock and after-shock sequences related to mainshocks, as well as earthquake swarms, shows that only a portion of the b-value changes is statistically significant.
Tundra be dammed
(2018)
Increasing air temperatures are changing the arctic tundra biome. Permafrost is thawing, snow duration is decreasing, shrub vegetation is proliferating, and boreal wildlife is encroaching. Here we present evidence of the recent range expansion of North American beaver (Castor canadensis) into the Arctic, and consider how this ecosystem engineer might reshape the landscape, biodiversity, and ecosystem processes. We developed a remote sensing approach that maps formation and disappearance of ponds associated with beaver activity. Since 1999, 56 new beaver pond complexes were identified, indicating that beavers are colonizing a predominantly tundra region (18,293km(2)) of northwest Alaska. It is unclear how improved tundra stream habitat, population rebound following overtrapping for furs, or other factors are contributing to beaver range expansion. We discuss rates and likely routes of tundra beaver colonization, as well as effects on permafrost, stream ice regimes, and freshwater and riparian habitat. Beaver ponds and associated hydrologic changes are thawing permafrost. Pond formation increases winter water temperatures in the pond and downstream, likely creating new and more varied aquatic habitat, but specific biological implications are unknown. Beavers create dynamic wetlands and are agents of disturbance that may enhance ecosystem responses to warming in the Arctic.
Many regions central to our understanding of tectonics and landscape evolution are active or ancient magmatic terranes, and robust interpretation of low-temperature thermochronologic ages in these settings requires careful attention to the drivers of rock heating and cooling, including magmatism. However, we currently lack a quantitative framework for evaluating the potential role of magmatic coolingthat is, post-magmatic thermal relaxationin shaping cooling age patterns in regions with a history of intrusive magmatism. Here we use analytical approximations and numerical models to characterize how low-temperature thermochronometers document cooling inside and around plutons in steadily exhuming environments. Our models predict that the thermal field a pluton intrudes into, specifically the ambient temperatures relative to the closure temperature of a given thermochronometer, is as important as the pluton size and temperature in controlling the pattern and extent of thermochronometer resetting in the country rocks around a pluton. We identify one advective and several conductive timescales that govern the relationship between the crystallization and cooling ages inside a pluton. In synthetic vertical age-elevation relationships (AERs), resetting next to plutons results in changes in AER slope that could be misinterpreted as past changes in exhumation rate if the history of magmatism is not accounted for. Finally, we find that large midcrustal plutons, such as those emplaced at similar to 10-15-km depth, can reset the low-temperature thermochronometers far above them in the upper crusta result with considerable consequences for thermochronology in arcs and regions with a history of magmatic activity that may not have a surface expression.
The variation of Rayleigh ellipticity versus frequency is gaining popularity in site characterization. It becomes a necessary observable to complement dispersion curves when inverting shear wave velocity profiles. Various methods have been proposed so far to extract polarization from ambient vibrations recorded on a single three-component station or with an array of three-component sensors. If only absolute values were recovered 10 yr ago, new array-based techniques were recently proposed with enhanced efficiencies providing also the ellipticity sign. With array processing, higher-order modes are often detected even in the ellipticity domain. We suggest to explore the properties of a high-resolution beamforming where radial and vertical components are explicitly included. If N is the number of three-component sensors, 2N x 2N cross-spectral density matrices are calculated for all presumed directions of propagation. They are built with N radial and N vertical channels. As a first approach, steering vectors are designed to fit with Rayleigh wave properties: the phase shift between radial and vertical components is either -Pi/2 or Pi/2. We show that neglecting the ellipticity tilt due to attenuation has only minor effects on the results. Additionally, we prove analytically that it is possible to retrieve the ellipticity value from the usual maximization of the high-resolution beam power. The method is tested on synthetic data sets and on experimental data. Both are reference sites already analysed by several authors. A detailed comparison with previous results on these cases is provided.
Woodlands and steppes
(2018)
Based on fossil organism remains including plant macrofossils, charcoal, pollen, and invertebrates preserved in syngenetic deposits of the Batagay permafrost sequence in the Siberian Yana Highlands, we reconstructed the environmental history during marine isotope stages (MIS) 6 to 2. Two fossil assemblages, exceptionally rich in plant remains, allowed for a detailed description of the palaeo-vegetation during two climate extremes of the Late Pleistocene, the onset of the last glacial maximum (LGM) and the last interglacial. In addition, altogether 41 assemblages were used to outline the vegetation history since the penultimate cold stage of MIS 6. Accordingly, meadow steppes analogue to modern communities of the phytosociological order Festucetalia lenensis formed the primary vegetation during the Saalian and Weichselian cold stages. Cold-resistant tundra-steppe communities (Carici rupestris-Kobresietea bellardii) as they occur above the treeline today were, in contrast to more northern locations, mostly lacking. During the last interglacial, open coniferous woodland similar to modern larch taiga was the primary vegetation at the site. Abundant charcoal indicates wildfire events during the last interglacial. Zoogenic disturbances of the local vegetation were indicated by the presence of ruderal plants, especially by abundant Urtica dioica, suggesting that the area was an interglacial refugium for large herbivores. Meadow steppes, which formed the primary vegetation during cold stages and provided potentially suitable pastures for herbivores, were a significant constituent of the plant cover in the Yana Highlands also under the full warm stage conditions of the last interglacial. Consequently, meadow steppes occurred in the Yana Highlands during the entire investigated timespan from MIS 6 to MIS 2 documenting a remarkable environmental stability. Thus, the proportion of meadow steppe vegetation merely shifted in response to the respectively prevailing climatic conditions. Their persistence indicates low precipitation and a relatively warm growing season throughout and beyond the late Pleistocene. The studied fossil record also proves that modern steppe occurrences in the Yana Highlands did not establish as late as in the Holocene but instead are relicts of a formerly continuous steppe belt extending from Central Siberia to Northeast Yakutia during the Pleistocene. The persistence of plants and invertebrates characteristic of meadow steppe vegetation in interior Yakutia throughout the late Quaternary indicates climatic continuity and documents the suitability of this region as a refugium also for other organisms of the Pleistocene mammoth steppe including the iconic large herbivores. (C)2018 Elsevier Ltd. All rights reserved.
Both enstatite spectra are very bright in the VIS and NIR and show almost neutral to slightly bluish spectral slopes with a steep absorption in the UV. Very low iron in the enstatites (below similar to 0.04 wt% FeO) already results in weak albeit noticeable absorptions in the VNIR between 0.4 and 0.9 mu m. Orthoenstatite and clinoenstatite are not distinguishable based only on their spectra in the VIS and NIR. At the Reststrahlen bands in the MIR a systematic difference in the number and exact position of local minima at similar to 10 mu m between clinoenstatite and orthoenstatite is evident. This can be used to discern between the polymorphs in this wavelength range. Additionally, we can distinguish between Fe-free low- and high-Ca pyroxenes in the MIR.
We present the first high-resolution trace element (Mg/Ca, Sr/Ca, Ba/Ca) record from a stalagmite in southwestern Romania covering the last 3.6 ka, which provides the potential for quantitative climate reconstruction. Precise age control is based on three independent dating methods, in particular for the last 250 yr, where chemical lamina counting is combined with the identification of the 20th century radiocarbon bomb peak and Th-230/U dating. Long-term cave monitoring and model simulations of drip water and speleothem elemental variability indicate that precipitation-related processes are the main drivers of speleothem Mg/Ca ratios. Calibration against instrumental climate data shows a significant anti-correlation of speleothem Mg/Ca ratios with autumn/winter (October to March) precipitation (r = -0.61, p < 0.01), which is statistically robust when considering age uncertainties and auto-correlation. This relationship is used to develop a quantitative reconstruction of autumn/winter precipitation. During the late Holocene, our data suggest a heterogeneous pattern of past regional winter hydroclimate in the Carpathian/Balkan realm, along with intermittent weakening of the dominant influence of North Atlantic forcing. In agreement with other regional paleo-hydrological reconstructions, the observed variability reveals periodically occurring strong NW-SE hydro-climate gradients. We hypothesize, that this pattern is caused by shifts of the eastern edge of the area of influence of the NAO across central eastern Europe due to the confluence of North Atlantic forcing, and other climatic features such as the East Atlantic/Western Russia (EAWR) pattern. (C) 2018 Elsevier B.V. All rights reserved.
Nearly 13,000 years ago, the warming trend into the Holocene was sharply interrupted by a reversal to near glacial conditions. Climatic causes and ecological consequences of the Younger Dryas (YD) have been extensively studied, however proxy archives from the Mediterranean basin capturing this period are scarce and do not provide annual resolution. Here, we report a hydroclimatic reconstruction from stable isotopes (delta O-18, delta C-13) in subfossil pines from southern France. Growing before and during the transition period into the YD (12 900-12 600 cal BP), the trees provide an annually resolved, continuous sequence of atmospheric change. Isotopic signature of tree sourcewater (delta O-18(sw)) and estimates of relative air humidity were reconstructed as a proxy for variations in air mass origin and precipitation regime. We find a distinct increase in inter-annual variability of sourcewater isotopes (delta O-18(sw)), with three major downturn phases of increasing magnitude beginning at 12 740 cal BP. The observed variation most likely results from an amplified intensity of North Atlantic (low delta O-18(sw)) versus Mediterranean (high delta O-18(sw)) precipitation. This marked pattern of climate variability is not seen in records from higher latitudes and is likely a consequence of atmospheric circulation oscillations at the margin of the southward moving polar front.
The Sierra Nevada de Santa Marta in NW Colombia is an isolated massif at the northernmost end of the Andes chain near the boundary with the Caribbean plate. Previous geologic mapping and K-Ar dating have shown that Jurassic plutonic and volcanic units make up a large part of the Santa Marta Massif (SMM). These rocks have been considered to be part of a Jurassic magmatic arc extending from NW Colombia to northern Chile, but without any geochemical basis for comparison. This paper reports on a geochemical and Sr-Nd-Pb isotope study of the Jurassic rocks in the SMM and provides 12 new U-Pb zircon ages from in-situ laser ICP-MS dating. The plutonic and volcanic units span a range from 45 to 78 wt.% SiO2, with a dominance of intermediate to felsic compositions with SiO2 > 57 wt.%. They classify as calc-alkaline, medium to high-K, metaluminous rocks with trace-element features typical for arc-derived magma series. In terms of their major and trace-element compositions, the SMM Jurassic units overlap with contemporary plutonic and volcanic rocks from other regions of the Central and Eastern Cordilleras of Colombia, and confirm an arc affinity. The new U-Pb ages range from 176 +/- 1 Ma to 192 +/- 2 Ma (n = 12), with most between 180 and 188 Ma (n = 7). The initial Sr isotope ratios (at 180 Ma) are between 0.7012 and 0.7071 (n = 29), with 3 outliers attributed to mobilization of Rb and/or Sr, Nearly all samples have negative( )epsilon Nd-(180) values between - 10.3 and 0.0 (n = 30), the two exceptions being only slightly positive (1.1 and 1.9). Measured Pb isotope ratios fall in a narrow range, with Pb-206/Pb-204 from 18.02 to 19.95, (207) Pb/(204) Pb from 15.56 to 15.67 and Pb-208/Pb-204 from 37.76 to 39.04 (n = 28). In the regional context of previous studies, these results confirm early Jurassic ages and an arc affinity for the widespread magmatism exposed in the eastern and northeastern Colombian Andes. We also note patterns in the distribution and composition of magmas. The magmatic activity in the Central Cordillera tends to be younger than in the Eastern Cordillera and is spatially more restricted to the vicinity of regional fault systems. In terms of composition, Jurassic igneous rocks in the Eastern Cordillera have systematically lower epsilon Nd-(180) values than those from the Central Cordillera, whereas the Pb isotope ratios overlap. We ascribe the Nd isotope variations to heterogeneity in the mantle source and/or degree of crustal contamination, whereas the Pb isotope ratios are crust-dominated and similar throughout the region. The spatio-temporal and compositional evolution of Jurassic magmatic rocks in the Northern Andes reflect the major plate kinematic readjustment between the Triassic and the Early Jurassic in the proto-Andean margin.
The article provides a contribution to the glacial geomorphology of tropical high-altitude mountains. It focuses on the diversity of glacial landforms in the Cordillera Blanca (Peru). The landforms are classified as individual landform types as well as in their hypsometric sequence. Their spatial arrangement is depicted in vertical series and summarised in a glacial-morphological altitudinal belt model. The glacial landforms are analysed in a genetic succession from the moraine-delimited glacier forefields to the excessive pedestal moraines and taking account of the topographical conditions in the verticality. The Pleistocene glacial relief formation is considered as one of the dominating factors in the development of the moraine types. The sediments are on the one hand in situ deposition forms and on the other depositions dislocated by mass movements and represent paraglacial landforms. Using the example of the Paron Valley in the Artesonraju-Huandoy Massif and neighbouring valleys, the article provides a typology of the range of glacial and paraglacial landforms based on a holistic glacial landsystem approach with special reference to pedestal moraines. (C) 2018 Elsevier B.V. All rights reserved.
Risk-based insurance is a commonly proposed and discussed flood risk adaptation mechanism in policy debates across the world such as in the United Kingdom and the United States of America. However, both risk-based premiums and growing risk pose increasing difficulties for insurance to remain affordable. An empirical concept of affordability is required as the affordability of adaption strategies is an important concern for policymakers, yet such a concept is not often examined. Therefore, a robust metric with a commonly acceptable affordability threshold is required. A robust metric allows for a previously normative concept to be quantified in monetary terms, and in this way, the metric is rendered more suitable for integration into public policy debates. This paper investigates the degree to which risk-based flood insurance premiums are unaffordable in Europe. In addition, this paper compares the outcomes generated by three different definitions of unaffordability in order to investigate the most robust definition. In doing so, the residual income definition was found to be the least sensitive to changes in the threshold. While this paper focuses on Europe, the selected definition can be employed elsewhere in the world and across adaption measures in order to develop a common metric for indicating the potential unaffordability problem.
Tomato (Solanum lycopersicum L.) being an important vegetable is cultivated and used throughout the world. It not only contributes in fulfilling the basic nutritional requirements of the human body but also has many health benefits due to its rich biochemical composition. However, its production at large scale is hampered by many limiting factors such as biotic and abiotic stresses. Among the different abiotic stresses, drought poses drastic impact on tomato yield. Drought stress is genetically regulated by many transcription factors that not only regulate the stress responsive mechanism but also facilitate the growth and development of tomato plants. NAC is an important stress related transcription factor genes family, and the ATAF1 gene, a member of this family, is involved in ABA signaling and stress response. In this study, tomato variety Rio Drande was transformed with drought resistant ATAF1 gene via Agrobacterium mediated gene transformation method. The ATAF1 gene was first cloned in the pK7WFG2 vector having kanamycin selectable marker and then it was introduced in the Agrobacterium tumefaciens strain GV3101 through heat shock method. The tomato cotyledon and hypocotyl ex-plants of variety "Rio Ggrande" were cultured on callus induction medium (MS + 2.5 mg/L IAA + 2 mg/L BAP). The calli were then infected with Agrobacterium tumefaciens strain GV3101 containing ATAF1 gene and selection was carried out on the kanamycin selectable medium (MS + 100 mg/L Kan), and were regenerated on MS medium with 1 mg/L IAA + 1 mg/L BAP. Out of 216 putative transformed calli, 13 calli were able to regenerate on the selection medium. Of the 13 calli, three transgenic tomato plantlets were recovered, and these were confirmed through PCR analysis for the presence of 432 bp fragment of ATAF1 gene. The transformation protocol reported here can be used to generate drought resistant tomato plants in future.
It is of major interest to estimate the feedback of arctic ecosystems to the global warming we expect in upcoming decades. The speed of this response is driven by the potential of species to migrate, tracking their climate optimum. For this, sessile plants have to produce and disperse seeds to newly available habitats, and pollination of ovules is needed for the seeds to be viable. These two processes are also the vectors that pass genetic information through a population. A restricted exchange among subpopulations might lead to a maladapted population due to diversity losses. Hence, a realistic implementation of these dispersal processes into a simulation model would allow an assessment of the importance of diversity for the migration of plant species in various environments worldwide. To date, dynamic global vegetation models have been optimized for a global application and overestimate the migration of biome shifts in currently warming temperatures. We hypothesize that this is caused by neglecting important fine-scale processes, which are necessary to estimate realistic vegetation trajectories. Recently, we built and parameterized a simulation model LAVESI for larches that dominate the latitudinal treelines in the northernmost areas of Siberia. In this study, we updated the vegetation model by including seed and pollen dispersal driven by wind speed and direction. The seed dispersal is modelled as a ballistic flight, and for the pollination of ovules of seeds produced, we implemented a wind-determined and distance-dependent probability distribution function using a von Mises distribution to select the pollen donor. A local sensitivity analysis of both processes supported the robustness of the model's results to the parameterization, although it highlighted the importance of recruitment and seed dispersal traits for migration rates. This individual-based and spatially explicit implementation of both dispersal processes makes it easily feasible to inherit plant traits and genetic information to assess the impact of migration processes on the genetics. Finally, we suggest how the final model can be applied to substantially help in unveiling the important drivers of migration dynamics and, with this, guide the improvement of recent global vegetation models.
The equatorial electrojet occasionally reverses during morning and afternoon hours, leading to periods of westward current in the ionospheric E region that are known as counter electrojet (CEJ) events. We present the first analysis of CEJ climatology and CEJ dependence on solar flux and lunar phase for the Brazilian sector, based on an extensive ground-based data set for the years 2008 to 2017 from the geomagnetic observatory Tatuoca (1.2 degrees S, 48.5 degrees W), and we compare it to the results found for Huancayo (12.0 degrees S, 75.3 degrees W) observatory in the Peruvian sector. We found a predominance of morning CEJ events for both sectors. The afternoon CEJ occurrence rate in the Brazilian sector is twice as high as in the Peruvian sector. The afternoon CEJ occurrence rate strongly depends on season, with maximum rates occurring during the northern-hemisphere summer for the Brazilian sector and during the northern-hemisphere winter for the Peruvian sector. Significant discrepancies between the two sectors are also found for morning CEJ rates during the northern-hemisphere summer. These longitudinal differences are in agreement with a CEJ climatology derived from contemporary Swarm satellite data and can be attributed in part to the well-known longitudinal wave-4 structure in the background equatorial electrojet strength that results from nonmigrating solar tides and stationary planetary waves. Simulations with the Thermosphere-Ionosphere-Electrodynamics General Circulation Model show that the remaining longitudinal variability in CEJ during northern summer can be explained by the effect of migrating tides in the presence of the varying geomagnetic field in the South Atlantic Anomaly.
TDS-1 GNSS Reflectometry
(2018)
This study presents the development and a systematic evaluation study of GNSS reflectometry wind speeds. After establishing a wind speed retrieval algorithm, UK TechDemoSat-1 (TDS-1) derived winds, from May 2015 to July 2017, are compared to the Advanced Scatterometer (ASCAT). ERA-Interim wind fields of the European Centre for Medium-range Weather Forecasts (ECMWF) and in situ observation from Tropical Atmosphere Ocean buoy array in the Pacific are taken as reference. One-year averaged TDS-1 global winds demonstrate small differences with ECMWF in a majority of areas as well as discuss under- and overestimations. The pioneering TDS-1 winds demonstrate a root-mean-squared error (RMSE) and bias of 2.77 and -0.33 m/s, which are comparable to the RMSE and bias derived by ASCAT winds, as large as 2.31 and 0.25 m/s, respectively. Using buoys measurements as reference, RMSE and bias of 2.23 and -0.03 m/s for TDS-1 as well as 1.40 and -0.68 m/s for ASCAT are obtained. Utilizing rain microwave-infrared estimates of the Tropical Rainfall Measuring Mission, rain-affected observation of both ASCAT and TDS-1 are collected and evaluated. Although ASCAT winds show a significant performance degradation resulting in an RMSE and bias of 3.16 and 1.03 m/s, respectively, during rain condition, TDS-1 shows a more reliable performance with an RMSE and bias of 2.94 and -0.21 m/s, respectively, which indicates the promising capability of GNSS forward scattering for wind retrievals during rain. A decrease in TDS-1-derived bistatic radar cross sections during rain events, at weak winds, is also demonstrated.
Recent studies indicate that the tropical belt has been expanding during recent decades, which can significantly influence precipitation in subtropical climates. Often the location of the tropical border is identified using the Hadley cell edge (HCE) or the subtropical jet stream (STJ), but most studies concentrated on the zonal-mean state, thereby missing regional impacts. Here we detect longitudinal-resolved trends in STJ cores and HCEs over 1979-2016 in both hemispheres at a higher spatial and temporal resolution than previous studies. Besides pronounced regional trend differences in both sign and magnitude, we show that winter HCE and STJ variability increased in the Mediterranean region and decreased over the American and Asian continents. Rainfall variability in these regions changed likewise, and most of those changes can be explained by changes in HCE/STJ variability. This highlights the importance of understanding future tropical belt changes both regionally and in terms of variability. Plain Language Summary We applied a new network-based method to detect motion of the tropical climate border with longitudinal resolution. Depending on the longitudinal position, there are differences in both direction and magnitude of the border motion. In addition, we demonstrate that the rainfall variability is increasing in the Mediterranean region and decreasing over the American and Asian continents, which can be explained by the variability of the tropical belt location.
Coupling land-use change and hydrologic models for quantification of catchment ecosystem services
(2018)
Representation of land-use and hydrologic interactions in respective models has traditionally been problematic. The use of static land-use in most hydrologic models or that of the use of simple hydrologic proxies in land-use change models call for more integrated approaches. The objective of this study is to assess whether dynamic feedback between land-use change and hydrology can (1) improve model performances, and/or (2) produce a more realistic quantification of ecosystem services. To test this, we coupled a land-use change model and a hydrologic mode. First, the land-use change and the hydrologic models were separately developed and calibrated. Then, the two models were dynamically coupled to exchange data at yearly time-steps. The approach is applied to a catchment in South Africa. Performance of coupled models when compared to the uncoupled models were marginal, but the coupled models excelled at the quantification of catchment ecosystem services more robustly.
The ability of water to transport and transform soil materials is one of the main drivers of soil and landscape development. In turn, soil and landscape properties determine how water is distributed in soil landscapes. Understanding the complex dynamics of this co-evolution of soils, landscapes and the hydrological system is fundamental in adapting land management to changes in climate. Soil-Landscape Evolution Models (SLEMs) are used to simulate the development and evolution of soils and landscapes. However, many hydrologic processes, such as preferential flow and subsurface lateral flow, are currently absent in these models. This limits the applicability of SLEMs to improve our understanding of feedbacks in the hydro-pedo-geomorphological system. Implementation of these hydrologic processes in SLEMs faces several complications related to calculation demands, limited methods for linking pedogenic and hydrologic processes, and limited data on quantification of changes in the hydrological system over time. In this contribution, we first briefly review processes and feedbacks in soil-landscape-hydrological systems. Next, we elaborate on the development required to include these processes in SLEMs. We discuss the state-of-the-art knowledge, identify complications, give partial solutions and suggest important future development. The main requirements for incorporating hydrologic processes in SLEMs are: (1) designing a model framework that can deal with varying timescales for different sets of processes, (2) developing and implementing methods for simulating pedogenesis as a function of water flow, (3) improving and implementing knowledge on the evolution and dynamics of soil hydraulic properties over different timescales, and (4) improving the database on temporal changes and dynamics of flow paths.
Thermal conductivity (lambda) is an essential physical property of minerals and rocks and fundamental in constraining the thermal field of the lithosphere. In case that adequate samples to measure lambda are not available, it could be indirectly inferred from calculation. One of the most widely applied indirect methods for rocks involve modal mineralogy and porosity as parameters that are incorporated into mathematical mean or mixing models. Robust inferences from these approaches for crystalline rocks were impeded by a small number of studied samples or restriction to certain rock types. We employ this method and examine its applicability to low-porosity plutonic rocks by calculating bulk thermal conductivity lambda(b) for 45 samples covering the entire range from gabbro/diorite to granite. We show that the use of the harmonic-mean model for both rock matrix and porosity provided a good match between lambda(b.meas) and lambda(b.calc) of <10% deviation (2 sigma), with relative and absolute errors amounting to 1.49.7% and 4.44.9%, respectively. The results of our study constitute a big step forward to a robust conclusion on the overall applicability of the harmonic-mean model for inferring lambda(b) of isotropic, low-porosity, mafic to silicic plutonic and metamorphic rocks with an acceptable magnitude of error. Drill cuttings and enclaves form particularly interesting objects for application of this method, as they are poorly suited for direct measurement. Well-derived lambda values for those rocks would permit to calculate heat flow and to model more profoundly the thermal state of the deeper lithosphere.