Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Climate gradient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

Median bed-material sediment particle size across rivers in the contiguous US

Abstract. Bed-material sediment particle size data, particularly the median sediment particle size (D50), are critical for understanding and modeling riverine sediment transport. However, sediment particle size observations are primarily available at individual sites. Large-scale modeling and assessment of riverine sediment transport are limited by the lack of continuous regional maps of bed-material sediment particle size. We hence present a map of D50 over the contiguous US in a vector format that corresponds to approximately 2.7 million river segments (i.e., flowlines) in the National Hydrography Dataset Plus (NHDPlus) dataset. We develop the map in four steps: (1) collect and process the observed D50 data from 2577 U.S. Geological Survey stations or U.S. Army Corps of Engineers sampling locations; (2) collocate these data with the NHDPlus flowlines based on their geographic locations, resulting in 1691 flowlines with collocated D50 values; (3) develop a predictive model using the eXtreme Gradient Boosting (XGBoost) machine learning method based on the observed D50 data and the corresponding climate, hydrology, geology, and other attributes retrieved from the NHDPlus dataset; and (4) estimate the D50 values for flowlines without observations using the XGBoost predictive model. We expect this map to be useful for various purposes, such as research in large-scale river sediment transport using model- and data-driven approaches, teaching environmental and earth system sciences, planning and managing floodplain zones, etc. The map is available at https://doi.org/10.5281/zenodo.4921987 (Li et al., 2021a).

54 ENVIRONMENTAL SCIENCES↗

Sea breezes and advective effects in southwest James Bay

Observations from a transect extending 100 km inland during the Northern Wetlands Study (NOWES) in 1990 show that the sea breeze develops on approximately 25% of days during summer and may penetrate up to 100 km inland on occasions. The sea breeze exhibits a marked diurnal clockwise rotation as a result of the Coriolis effect along the unobstructed coastline. The marine advective effect is shown to depend on gradient wind direction. With northwesterly upper level flow the sea breeze tends to be northeasterly in direction and is associated with decreased temperatures and vapor pressure deficits (VPD). With southwesterly upper level flow the sea breeze tends to have a southeasterly direction and less effect on temperatures and VPD. This is attributed to shorter residence times of air parcels over water. For two cases, Colorado State University mesoscale model simulations show good agreement with surface wind observations and suggest that under northwesterly gradient flow, Bowen ratios are increased in the onshore flow along western James Bay, while during southwesterly gradient flow these effects are negligible. These results have implications for the interpretation of local climate, ecology, and hydrology as well as land-based and airborne turbulent flux measurements made during NOWES.

Mckendry, Ian↗

Disturbances in North American boreal forest and Arctic tundra: impacts, interactions, and responses

Ecosystems in the North American Arctic-Boreal Zone (ABZ) experience a diverse set of disturbances associated with wildfire, permafrost dynamics, geomorphic processes, insect outbreaks and pathogens, extreme weather events, and human activity. Climate warming in the ABZ is occurring at over twice the rate of the global average, and as a result the extent, frequency, and severity of these disturbances are increasing rapidly. Disturbances in the ABZ span a wide gradient of spatiotemporal scales and have varying impacts on ecosystem properties and function. However, many ABZ disturbances are relatively understudied and have different sensitivities to climate and trajectories of recovery, resulting in considerable uncertainty in the impacts of climate warming and human land use on ABZ vegetation dynamics and in the interactions between disturbance types. Here we review the current knowledge of ABZ disturbances and their precursors, ecosystem impacts, temporal frequencies, spatial extents, and severity. We also summarize current knowledge of interactions and feedbacks among ABZ disturbances and characterize typical trajectories of vegetation loss and recovery in response to ecosystem disturbance using satellite time-series. We conclude with a summary of critical data and knowledge gaps and identify priorities for future study.

54 ENVIRONMENTAL SCIENCES↗

Disturbances in North American Boreal Forest and Arctic Tundra: Impacts, Interactions, and Responses

Ecosystems in the North American Arctic-Boreal Zone (ABZ) experience a diverse set of disturbances associated with wildfire, permafrost dynamics, geomorphic processes, insect outbreaks and pathogens, extreme weather events, and human activity. Climate warming in the ABZ is occurring at over twice the rate of the global average, and as a result the extent, frequency, and severity of these disturbances are increasing rapidly. Disturbances in the ABZ span a wide gradient of spatiotemporal scales and have varying impacts on ecosystem properties and function. However, many ABZ disturbances are relatively understudied and have different sensitivities to climate and trajectories of recovery, resulting in considerable uncertainty in the impacts of climate warming and human land use on ABZ vegetation dynamics and in the interactions between disturbance types. Here we review the current knowledge of ABZ disturbances and their precursors, ecosystem impacts, temporal frequencies, spatial extents, and severity. We also summarize current knowledge of interactions and feedbacks among ABZ disturbances and characterize typical trajectories of vegetation loss and recovery in response to ecosystem disturbance using satellite time-series. We conclude with a summary of critical data and knowledge gaps and identify priorities for future study.

high-latitude↗

Molecular Nature of Mineral-Organic Associations within Redox-Active Mountainous Floodplain Sediments

Floodplains are critical terrestrial-aquatic interfaces that act as hotspots of organic carbon (OC) cycling, regulating ecosystem carbon storage as well as export to riverine systems. Within floodplain sediments, regular flooding and textural gradients interact to create dynamic redox conditions. While anaerobic protection of OC upon burial is a well-recognized carbon storage mechanism in redox-active floodplain sediments, the impact of protective mineral-organic associations is relatively unknown. Here we determined the quantitative importance and chemical composition of mineral-organic associations along well-defined redox gradients emerging from textural variations and depth within meander sediments of the subalpine East River watershed (Gothic, CO). Here we characterized mineral-organic associations using a combination of sequential extractions, physical fractionation, and high-resolution mass spectrometry. Across the meander, we found that mineral-associated OC constitutes a significant fraction of total OC, and that extractable iron (Fe) and aluminum (Al) phases as well as high-density isolates were strongly correlated with total OC content, suggesting that mineral-organic associations are quantitatively important for floodplain sediment OC protection. Our mass spectrometry results showed OC associated with increasingly ordered Fe and Al phases are relatively enriched in low-molecular weight, oxidized, aromatic compounds. Surprisingly, however, total OC content showed weak or no correlation with indicators of anaerobic protection, such as relatively bioavailable OC pools (water-extractable and particulate OC) or the molecular weight and oxidation state of OC. Overall, this work highlights that protection of OC bound to reactive mineral phases-in addition to anaerobic protection-can play a quantitatively important role in controlling soil carbon storage in redox-active floodplain sediments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The effects of surface evaporation parameterizations on climate sensitivity to solar constant variations

The effects of two different evaporation parameterizations on the sensitivity of simulated climate to solar constant variations are investigated by using a zonally averaged climate model. One parameterization is a nonlinear formulation in which the evaporation is nonlinearly proportional to the sensible heat flux, with the Bowen ratio determined by the predicted vertical temperature and humidity gradients near the earth's surface (model A). The other is the formulation of Saltzman (1968) with the evaporation linearly proportional to the sensible heat flux (model B). The computed climates of models A and B are in good agreement except for the energy partition between sensible and latent heat at the earth's surface. The difference in evaporation parameterizations causes a difference in the response of temperature lapse rate to solar constant variations and a difference in the sensitivity of longwave radiation to surface temperature which leads to a smaller sensitivity of surface temperature to solar constant variations in model A than in model B. The results of model A are qualitatively in agreement with those of the general circulation model calculations of Wetherald and Manabe (1975).

Chou, S.-H.↗

Amplified Extreme Floods and Shifting Flood Mechanisms in the Delaware River Basin in Future Climates

Abstract Historical records in the Delaware River Basin reveal complex and spatially diverse flood generating mechanisms influenced by the region's mountains‐to‐plains gradients. This study focuses on predicting future flood hazards and understanding the underlying drivers of changes across the region. Using a process‐based hydrological model, we analyzed the hydrometeorological condition of each historical and future flood event. For each event, at the subbasin scale, we identified the dominant flood generating mechanism, including snowmelt, rain‐on‐snow, short‐duration rain, and long‐duration rain. The rain‐induced floods are further categorized based on the soil's Antecedent Moisture Condition (AMC) before the event, whether dry, normal, or wet. Our historical analysis suggests that rain‐on‐snow is the primary flood mechanism of the Upper Basin. Although most frequent, the magnitude of rain‐on‐snow floods is often less severe than short rain floods. In contrast, historical floods in the Lower Basin are primarily caused by short rain under normal AMC. Given the uncertainties in climate projections, we used an ensemble of future climate scenarios for flood projections. Despite variations in regional climate projections, coherent perspectives emerge: the region will shift toward a warmer, wetter climate, with a projected intensification of extreme floods. The Upper Basin is projected to experience a marked decrease in rain‐on‐snow floods, but a substantial increase in short rain floods with wet AMC. The largest increase in flood magnitude will be driven by short rains with wet AMC in the Upper Basin and by short rains with normal AMC in the Lower Basin.

54 ENVIRONMENTAL SCIENCES↗

Baroclinic adjustment

A detailed comparison is presented of the actual shear in the atmosphere with the critical shear given by the two-layer model of Phillips (1954), in which there is a critical temperature gradient separating stable conditions from baroclinically unstable ones. A very simple parameterization of the effect of eddy fluxes on atmospheric temperature is suggested, where the parameterization includes beta effects. The parameterization is illustrated by applying it in a one-dimensional heat-balance climate model. Enhancement of the eddy flux in a continuous atmosphere under supercritical conditions is stressed. This enhancement leads to a negative feedback between the meridional eddy flux of heat and the meridional temperature gradient. The feedback restricts gradients to values near the critical value, a process referred to as baroclinic adjustment. This should facilitate the development of simple climate models involving feedbacks associated with both the meridional vertical temperature structure.

Stone, P. H.↗

Climate, soil mineralogy and mycorrhizal fungi influence soil organic matter fractions in eastern US temperate forests

Identifying the primary controls of particulate (POM) and mineral-associated organic matter (MAOM) content in soils is critical for determining future stocks of soil carbon (C) and nitrogen (N) across the globe. However, drivers of these soil organic matter fractions are likely to vary among ecosystems in response to climate, soil type and the composition of local biological communities. We tested how soil factors, climate and plant–fungal associations influenced the distribution and concentrations of C and N in MAOM and POM in seven temperate forests in the National Ecological Observatory Network (NEON) across the eastern United States. Samples of upper mineral horizon soil within each forest were collected in plots representing a gradient of dominant tree–mycorrhizal association, allowing us to test how plant and microbial communities influenced POM and MAOM across sites differing in climate and soil conditions. We found that concentrations of C and N in soil organic matter were primarily driven by soil mineralogy, but the relative abundance of MAOM versus POM C was strongly linked to plot-level mycorrhizal dominance. Furthermore, the effect of dominant tree mycorrhizal type on the distribution of N among POM and MAOM fractions was sensitive to local climate: in cooler sites, an increasing proportion of ectomycorrhizal-associated trees was associated with lower proportions of N in MAOM, but in warmer sites, we found the reverse. As an indicator of soil carbon age, we measured radiocarbon in the MAOM fraction but found that within and across sites, Δ 14 C was unrelated to mycorrhizal dominance, climate, or soil factors, suggesting that additional site-specific factors may be primary determinants of long-term SOM persistence. Our results indicate that while soil mineralogy primarily controls SOM C and N concentrations, the distribution of SOM among density fractions depends on the composition of vegetation and microbial communities, with these effects varying across sites with distinct climates. Here we also suggest that within biomes, the age of mineral-associated soil carbon is not clearly linked to the factors that control concentrations of MAOM C and N.

54 ENVIRONMENTAL SCIENCES↗

Record-breaking 2020 summer marine heatwaves in the western North Pacific

Record-breaking marine heatwaves (MHWs) occurred in the western North Pacific during the summer of 2020. These unprecedented MHWs were consistent with favorable large-scale conditions that are linked to an anomalous western North Pacific subtropical high (WNPSH), resulting mainly from sea surface temperature (SST) anomalies across the tropical oceans. In addition, a moderate La Niña-like pattern was also conducive to transporting warm seawater to the western North Pacific. Mixed-layer heat budgets suggest that surface heat flux contributed to the SST anomaly in the western subtropical Pacific. In contrast, oceanic heat advection dominanted in the South China Sea and the western equatorial Pacific. Numerical model experiments indicated that the tropical Indian Ocean SST anomalies were responsible for the enhanced WNPSH. The increased zonal SST gradient across the tropical Pacific also played an important role. Inter-ocean interactions can modulate climate variability through ocean-atmospheric coupling and deserve more attention when predicting MHWs within the context of global warming. Additionally, it is critical to consider MHWs as a powerful tool in detecting acute, intense thermal stress events in the coral bleaching pre-warning system.

Coral bleaching↗

Oxygen Isotopes in Fresh Water Biogenic Opal: Northeastern US Alleroed-Younger Dryas Temperature Shift

The first oxygen isotope analysis of biogenic opal from lake sediments, from the Allerod/Younger Dryas transition in a core from Linsley Pond, Connecticut, gives an average estimate of a 6 C drop in temperature during the Younger Dryas. This shift represents temperatures during the bloom season, and may be less than the winter temperature drop. The sharp transition itself, with a duration of about 200 years, suggests that the temperature decrease may have been as large as 12 C. Previous estimates of the Allerod/Younger Dryas temperature shifts are controversial, and range from 3-20 C, suggesting that further interdisciplinary research on the same samples is warranted. One way that global climate change manifests itself is by redistributing energy throughout the globe. The Northern Hemisphere latitudinal temperature gradient during the late-glacial is at present a controversial topic. The magnitude of air temperature shifts during the Allerod/Younger Dryas (YD) oscillation are estimated from mid-latitude pollen records surrounding the North Atlantic to be 3-5 C in Europe [Lowe et al., 19941 and 3-4 C in the eastern US [Peteet et al., 1993]. In contrast, lake temperatures estimates derived from aquatic midge larvae in the Canadian eastern maritimes and Maine range from 6-20 C, with larger shifts at more southern sites [Levesque et al., 1997]. The magnitude of YD cooling in Greenland ice cores ranges from at least 7 C from the Bolling warming [Dansgaard et al., 1989] to 15 C - a more recent estimate from borehole temperatures [Cuffey et al., 1995]. The ice core geochemical records reveal that massive frequent and short-term (decadal or less) changes in atmospheric composition occurred throughout this event, suggesting a very dynamic circulation [Mayewski et al., 1993).

Shemesh, Aldo↗

On the Dome Effect of Flux Radiometers to Radiative Forcing

Since the introduction of thermopile, pyranometers (solar, e.g., 0.3-3.0 micrometers) and pyrgeometers (terrestrial, e.g., 4-50 micrometers) have become instruments commonly used for measuring the broadband hemispherical irradiances at the surface in a long-term, monitoring mode for decades. These commercially available radiometers have been manufactured in several countries such as from the United States, Asia, and Europe, and are generally reliable and economical. These worldwide distributions of surface measurements become even more important in the era of Earth remote sensing in studying climate forcing. However, recent studies from field campaigns have pointed out that erroneous factors (e.g., temperature gradients between the filter dome and detector, emissivity of the thermopile) are responsible for the unacceptable level of uncertainty (e.g., 10-20 W m (exp -2)). Using a newly developed instrument of Quantum Well Infrared Photodetector (QWIP), we have characterized the brightness temperature fields of pyranometers and pyrgeometers under various sky conditions. The QWIP is based on the superlattice (GaAs/AlGaAs) technology and has a noise equivalent temperature (NE delta T) less than 0.1 K. The quality of pyranometer and pyrgeometer measurements can be improved largely by applying proper knowledge of the thermal parameters affecting the operation of the thermopile systems. For example, we show a method to determine the "dome factor" (the longwave emission divided by the longwave transmission of a pyrgeometer dome) from field measurements. The results show, and are verified independently by the QWIP, that our dome factors of 0.59 and 0.90 are much smaller than the value of 4.0 assumed by the WMO (World Meteorological Organization). Data correction procedure and algorithm will be presented and discussed.

Tsay, S.-C.↗

On the Dome Effect of Flux Radiometers to Radiative Forcing

Since the introduction of thermopile, pyranometers (solar, e.g., 0.3-3.0 microns) and pyrgeometers (terrestrial, e.g., 4-50 microns) have become instruments commonly used for measuring the broadband hemispherical irradiances at the surface in a long-term, monitoring mode for decades. These commercially available radiometers have been manufactured in several countries such as from the United States, Asia, and Europe, and are generally reliable and economical. These worldwide distributions of surface measurements become even more important in the era of Earth remote sensing in studying climate forcing. However, recent studies from field campaigns have pointed out that erroneous factors (e.g., temperature gradients between the filter dome and detector, emissivity of the thermopile) are responsible for the unacceptable level of uncertainty (e.g., 10-20 W/square Meter). Using a newly developed instrument of Quantum Well Infrared Photodetector (QWIP), we have characterized the brightness temperature fields of pyranometers and pyrgeometers under various sky conditions. The QWIP is based on the superlattice (GaAs/AlGaAs) technology and has a noise equivalent temperature (NE delta T) less than 0.1 K. The quality of pyranometer and pyrgeometer measurements can be improved largely by applying proper knowledge of the thermal parameters affecting the operation of the thermopile systems. For example, we show a method to determine the "dome factor" (the longwave emission divided by the longwave transmission of a pyrgeometer dome) from field measurements. The results show, and are verified independently by the QWIP, that our dome factors of 0.59 and 0.90 are much smaller than the value of 4.0 assumed by the WMO. Data correction procedure and algorithm will be presented and discussed.

Tsay, S.-C.↗

Characterization of Thermal Parameters for Improving Pyranometer and Pyrgeometer Measurements

Since the introduction of thermopile, pyranometers (solar, e.g., 0.3-3.0 micrometers) and pyrgeometers (terrestrial, e.g., 4-50 micrometers) have become instruments commonly used for measuring the broadband hemispherical irradiances at the surface in a long-term, monitoring mode for decades. These commercially available radiometers have been manufactured in several countries such as from the United States, Asia, and Europe, and are generally reliable and economical. These worldwide distributions of surface measurements become even more important in the era of Earth remote sensing in studying climate change. However, recent studies from field campaigns have pointed out that erroneous factors (e.g., temperature gradients between the filter dome and detector, emissivity of the thermopile) are responsible for the unacceptable level of uncertainty (e.g., 20 W m(exp -2)). Using a newly developed instrument of Quantum Well Infrared Photodetector (QWTP), we have characterized the brightness temperature fields of pyranometers and pyrgeometers under various sky conditions. The QWIP is based on the superlattice (GaAs/AlGaAs) technology and has a noise equivalent temperature (NEAT) less than 0.1 K. The quality of pyranometer and pyrgeometer measure- ments can be improved largely by applying proper knowledge of the thermal parameters affecting the operation of the thermopile systems. Data correction procedure and algorithm will be presented and discussed.

Tsay, Si-Chee↗

Thermal Characteristics of Pyranometers and Pyrgeometers in Atmosphere-Surface Energetic Measurements

Since the introduction of thermopile, pyranometers (solar, e.g., 0.3 - 3.0 microns) and pyrgeometers (terrestrial, e.g., 4 - 50 microns) have become instruments commonly used for measuring the broadband hemispherical irradiances at the surface in a long-term, monitoring mode for decades. These commercially available radiometers have been manufactured in several countries such as from the United States, Asia, and Europe, and are generally reliable and economical. These worldwide distributions of surface measurements become even more important in the era of Earth remote sensing in studying climate change. However, recent studies from field campaigns have pointed out that erroneous factors (e.g., temperature gradients between the filter dome and detector, emissivity of the thermopile) are responsible for the unacceptable level of uncertainty (e.g., 20 W/sq m). It is best to utilize an energy balance equation to describe the thermal dome effect of pyranometers and pyrgeometers. Therefore, quality of pyranometer and pyrgeometer measurements can be improved largely by applying proper knowledge of the thermal parameters affecting the operation of the thermopile systems. Data correction procedure and algorithm will be presented and discussed.

Tsay, Si-Chee↗

Electrical Breakdown in a Martian Gas Mixture

The high probability for dust interactions during Martian dust storms and dust devils combined with the cold, dry climate of Mars most likely result in airborne dust that is highly charged. On Earth, potential gradients up to 5 kV/m have been recorded and in some cases resulted in lightning. Although the Martian atmosphere is not conducive to lightning generation, it is widely believed that electrical discharge in the form of a corona occurs. In order to understand the breakdown of gases, Paschen measurements are taken which relate the minimum potential required to spark across a gap between two electrodes. The minimum potential is plotted versus the pressure-distance value for electrodes of a given geometry. For most gases, the potential decreases as the pressure decreases. For CO2, the minimum in the curve happens to be at Mars atmospheric pressures (5-7 mm Hg) for many distances and geometries. However, a very small amount (<0.1%) of mixing gases radically changes the curve, as noted by Leach. Here, we present the first experimental results of a Paschen curve for a Mars gas mixture compared with 100% pure CO2.

Buhler, C. R.↗

Forest Age Rivals Climate to Explain Reproductive Allocation Patterns in Forest Ecosystems Globally

ABSTRACT Forest allocation of net primary productivity (NPP) to reproduction (carbon required for flowers, fruits, and seeds) is poorly quantified globally, despite its critical role in forest regeneration and a well‐supported trade‐off with allocation to growth. Here, we present the first global synthesis of a biometric proxy for forest reproductive allocation (RA) across environmental and stand age gradients from a compiled dataset of 824 observations across 393 sites. We find that ecosystem‐scale RA increases ~60% from boreal to tropical forests. Climate shows important non‐linear relationships with RA, but is not the sole predictor. Forest age effects are comparable to climate in magnitude (MAT: ß = 0.24,p = 0.021; old growth forest: ß = 0.22,p < 0.001), while metrics of soil fertility show small but significant relationships with RA (soil pH: ß = 0.07,p = 0.001; soil N: ß = −0.07,p = 0.001). These results provide strong evidence that ecosystem‐scale RA is mediated by climate, forest age, and soil conditions, and is not a globally fixed fraction of positive NPP as assumed by most vegetation and ecosystem models. Our dataset and findings can be used by modellers to improve predictions of forest regeneration and carbon cycling.

Environmental Sciences & Ecology↗

Escarpment evolution on high-elevation rifted margins: Insights derived from a surface processes model that combines diffusion, advection and reaction

Experiments with a surface processes model of large-scale (1-1000 km) long-term (1-100 m.y.) erosional denudation are used to establish the controls on the evolution of a model escarpment that is related to the rifting of a continent. The mdoel describes changes in topographic form as a result of sumultaneous short- and long-range mass transport representing hillslope (diffusive) processes and fluvial transport (advection), repsectively. Fluvial entrainment is modeled as a first-order kinetic reaction which reflects the erodibility of the substrate, and therefore the fluvial system is not necessarily carrying at capacity. One dimensional and planform models demonstrate that the principal controls on the evolution of an initially steep model escarpment are (1) antecedent topography/drainage; (2) the timesale (or equivalently a length scale) in the fluvial entrainment reaction; (3) the flexural response of the lithosphere to denudation; and (4) the relative efficiencies of the short- and long-range transport processes. When rainfall and substrate lithology are uniform, a significant amount of discharge draining over the escarpment top causes it to degrade. Only when the top of the model escarpment coincides with a drainage divide can escarpment retreat occur for these conditions. An additional requirement for retreat of a model escarpment without decline is a long reaction time scale for fluvial entrainment. This corresponds to a substrate that is hard to detach by flucial erosion, and therefore to fluvial erosion that is not transport limited. Coninuous backtilting of an escarpment due ot flexural isostatic uplift in response to denudational unloading helps maintain the scarp top as a divide. It is essntial if the escarpment gradient is to be preserved during retreat in a uniform lithology. Low flexural rigidieties propote steep and slowly retreating escarpments. For given rainfall and substrate conditions, the morphology of a retraeating model escarpment is determined by the ratio of hte short-range diffusive and long range advective transport efficiencies. A low ration (which is interpreted to correspond to a more humid, temperate climate) produces a convex upper slope, and concave lower slope morphology and only major escarpments are predicted to preserve a high scarp gradient. Lithological contrasts in hte model produce more complex morphologies and predict the formation of scarps crowned by an erosionaly resistant caprock. However, resistant caprocks are not an essential requirement for model scarps to retreat. We conclude that the inferred controls and model behavior are both consistent with the present-day morphology of rifted continental margins and with modern conceptual models of landscape evolution.

Kooi, Henk↗