USGS Science⌕ Search

SEARCH · USGS Science

Results for “Evolution and Development”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 667 records · Page 37Linked to original sources

Possible pingos and a periglacial landscape in northwest Utopia Planitia

Hydrostatic (closed-system) pingos are small, elongate to circular, ice-cored mounds that are perennial features of some periglacial landscapes. The growth and development of hydrostatic pingos is contingent upon the presence of surface water, freezing processes and of deep, continuous, ice-cemented permafrost. Other cold-climate landforms such as small-sized, polygonal patterned ground also may occur in the areas where pingos are found. On Mars, landscapes comprising small, elongate to circular mounds and other possible periglacial features have been identified in various areas, including Utopia Planitia, where water is thought to have played an important role in landscape evolution. Despite the importance of the martian mounds as possible markers of water, most accounts of them in the planetary science literature have been brief and/or based upon Viking imagery. We use a high-resolution Mars Orbiter Camera image (EO300299) and superposed Mars Orbiter Laser Altimeter data tracks to describe and characterise a crater-floor landscape in northwest Utopia Planitia (64.8?? N/292.7?? W). The landscape comprises an assemblage of landforms that is consistent with the past presence of water and of periglacial processes. This geomorphological assemblage may have formed as recently as the last episode of high obliquity. A similar assemblage of landforms is found in the Tuktoyaktuk peninsula of northern Canada and other terrestrial cold-climate landscapes. We point to the similarity of the two assemblages and suggest that the small, roughly circular mounds on the floor of the impact crater in northwest Utopia Planitia are hydrostatic pingos. Like the hydrostatic pingos of the Tuktoyaktuk peninsula, the origin of the crater-floor mounds could be tied to the loss of ponded, local water, permafrost aggradation and the evolution of a sub-surface ice core. ?? 2004 Elsevier Inc. All rights reserved.

Icarus↗

Green‐wave surfing increases fat gain in a migratory ungulate

Each spring, migratory herbivores around the world track or ‘surf’ green waves of newly emergent vegetation to distant summer or wet‐season ranges. This foraging tactic may help explain the great abundance of migratory herbivores on many seasonal landscapes. However, the underlying fitness benefits of this life‐history strategy remain poorly understood. A fundamental prediction of the green‐wave hypothesis is that migratory herbivores obtain fitness benefits from surfing waves of newly emergent vegetation more closely than their resident counterparts. Here we evaluate whether this behavior increases body‐fat levels – a critically important correlate of reproduction and survival for most ungulates – in elk Cervus elaphus of the Greater Yellowstone Ecosystem. Using satellite imagery and GPS tracking data, we found evidence that migrants (n = 23) indeed surfed the green wave, occupying sites 12.7 days closer to peak green‐up than residents (n = 16). Importantly, individual variation in surfing may help account for up to 6 kg of variation in autumn body‐fat levels. Our findings point to a pathway for anthropogenic changes to the green wave (e.g. climate change) or migrants’ ability to surf it (e.g. development) to impact migratory populations. To explore this possibility, we evaluated potential population‐level consequences of constrained surfing with a heuristic model. If green‐wave surfing deteriorates by 5–15 days from observed, our model predicts up to a 20% decrease in pregnancy rates, a 2.5% decrease in population growth, and a 30% decrease in abundance over 50 years. By linking green‐wave surfing to fitness and illustrating potential effects on population growth, our study provides new insights into the evolution of migratory behavior and the prospects for the persistence of migratory ungulate populations in a changing world.

Wyoming↗

VS2DRTI: Simulating heat and reactive solute transport in variably saturated porous media

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid-rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult-to-use models. To address the need for a simple and easy-to-use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two-dimensional, constant-density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature-dependent cation exchange. VS2DRTI is freely available public domain software.

Groundwater↗

Development of a process-based littoral sediment transport model for Dauphin Island, Alabama

Dauphin Island, Alabama, located in the Northern Gulf of Mexico just outside of Mobile Bay, is Alabama’s only barrier island and provides an array of historical, natural, and economic resources. The dynamic island shoreline of Dauphin Island evolved across time scales while constantly acted upon by waves and currents during both storms and calm periods. Reductions in the vulnerability and enhancements to the resiliency of Dauphin Island—through offshore sand placement, breach closure, berm construction, and other means—have been used to protect the island and its vital resources. Planning for a resilient Dauphin Island requires predicting the long-term evolution of the barrier island system and the dominant, temporally varying processes that influence it, including littoral alongshore sediment transport under typical wave conditions, beach and dune erosion, the island overwash and breaching that occur rapidly during storm events, and the recovery of primary sand dunes through Aeolian transport over decadal time scales. Littoral sediment transport within the Dauphin Island decadal-scale framework was simulated using the Delft-3D modeling software suite. The influences of wind, waves, water levels, and sediment transport are incorporated into the model. Model skill in the prediction of waves, water levels, currents, volumetric flow rates through inlets, and shoreline position was assessed by using a set of deterministic and statistical hindcast simulations. The Delft-3D modeling application described here can be coupled with validated models of storm-response and dune recovery to predict the evolution of Dauphin Island on decadal time scales.

Alabama↗

Segmentation of plate coupling, fate of subduction fluids, and modes of arc magmatism in Cascadia, inferred from magnetotelluric resistivity

Five magnetotelluric (MT) profiles have been acquired across the Cascadia subduction system and transformed using 2-D and 3-D nonlinear inversion to yield electrical resistivity cross sections to depths of ∼200 km. Distinct changes in plate coupling, subduction fluid evolution, and modes of arc magmatism along the length of Cascadia are clearly expressed in the resistivity structure. Relatively high resistivities under the coasts of northern and southern Cascadia correlate with elevated degrees of inferred plate locking, and suggest fluid- and sediment-deficient conditions. In contrast, the north-central Oregon coastal structure is quite conductive from the plate interface to shallow depths offshore, correlating with poor plate locking and the possible presence of subducted sediments. Low-resistivity fluidized zones develop at slab depths of 35–40 km starting ∼100 km west of the arc on all profiles, and are interpreted to represent prograde metamorphic fluid release from the subducting slab. The fluids rise to forearc Moho levels, and sometimes shallower, as the arc is approached. The zones begin close to clusters of low-frequency earthquakes, suggesting fluid controls on the transition to steady sliding. Under the northern and southern Cascadia arc segments, low upper mantle resistivities are consistent with flux melting above the slab plus possible deep convective backarc upwelling toward the arc. In central Cascadia, extensional deformation is interpreted to segregate upper mantle melts leading to underplating and low resistivities at Moho to lower crustal levels below the arc and nearby backarc. The low- to high-temperature mantle wedge transition lies slightly trenchward of the arc.

Geochemistry, Geophysics, Geosystems↗

Implications of the miocene(?) crooked ridge river of northern arizona for the evolution of the colorado river and grand canyon

The southwesterly course of the probably pre–early Miocene and possibly Oligocene Crooked Ridge River can be traced continuously for 48 km and discontinuously for 91 km in northern Arizona (United States). The course is visible today in inverted relief. Pebbles in the river gravel came from at least as far northeast as the San Juan Mountains (Colorado). The river valley was carved out of easily eroded Jurassic and Cretaceous rocks whose debris overloaded the river with abundant detritus, probably steepening the gradient. After the river became inactive, the regional drainage network was rearranged three times, and the nearby Four Corners region was lowered 1–2 km by erosion. The river provides constraints on the early evolution of the Colorado River and Grand Canyon. Continuation of this river into lakes in Arizona or Utah is unlikely, as is integration through Grand Canyon by lake spillover. The downstream course of the river probably was across the Kaibab arch in a valley roughly coincident with the present eastern Grand Canyon. Beyond this point, the course may have continued to the drainage basin of the Sacramento River, or to the proto–Snake River drainage. Crooked Ridge River was beheaded by the developing San Juan River, which pirated its waters and probably was tributary to a proto–Colorado River, flowing roughly along its present course west of the Monument upwarp.

Geosphere↗

Sources of variation in breeding-ground fidelity of mallards (Anas platyrhynchos)

Generalizations used to support hypotheses about the evolution of fidelity to breeding areas in birds include the tendency for fidelity to be greater in adult birds than in yearlings. In ducks, in contrast to most bird species, fidelity is thought to be greater among females than males. Researchers have suggested that fidelity in ducks is positively correlated with pond availability. However, most estimates of fidelity on which these inferences have been based represent functions of survival and recapture-resighting probabilities in addition to fidelity. We applied the modeling approach developed by Burnham to recapture and band recovery data of mallard ducks to test the above hypotheses about fidelity. We found little evidence of sex differences in adult philopatry, with females being slightly more philopatric than males in one study area, but not in a second study area. However, yearling females were more philopatric than yearling males in both study areas. We found that adults were generally more philopatric than yearlings. We could find no relationship between fidelity and pond availability. Our results, while partially supporting current theory concerning sex and age differences in philopatry, suggest that adult male mallards are more philopatric than once thought, and we recommend that other generalizations about philopatry be revisited with proper estimation techniques.

Behavioral Ecology↗

Simulation of daily streamflow for 12 river basins in western Iowa using the Precipitation-Runoff Modeling System

The U.S. Geological Survey, in cooperation with the Iowa Department of Natural Resources, constructed Precipitation-Runoff Modeling System models to estimate daily streamflow for 12 river basins in western Iowa that drain into the Missouri River. The Precipitation-Runoff Modeling System is a deterministic, distributed-parameter, physical-process-based modeling system developed to evaluate the response of streamflow and general drainage basin hydrology to various combinations of climate and land use. Calibration periods for each basin varied depending on the period of record available for daily mean streamflow measurements at U.S. Geological Survey streamflow-gaging stations. A geographic information system tool was used to delineate each basin and estimate initial values for model parameters based on basin physical and geographical features. A U.S. Geological Survey automatic calibration tool that uses a shuffled complex evolution algorithm was used for initial calibration, and then manual modifications were made to parameter values to complete the calibration of each basin model. The main objective of the calibration was to match daily discharge values of simulated streamflow to measured daily discharge values. The Precipitation-Runoff Modeling System model was calibrated at 42 sites located in the 12 river basins in western Iowa. The accuracy of the simulated daily streamflow values at the 42 calibration sites varied by river and by site. The models were satisfactory at 36 of the sites based on statistical results. Unsatisfactory performance at the six other sites can be attributed to several factors: (1) low flow, no flow, and flashy flow conditions in headwater subbasins having a small drainage area; (2) poor representation of the groundwater and storage components of flow within a basin; (3) lack of accounting for basin withdrawals and water use; and (4) limited availability and accuracy of meteorological input data. The Precipitation-Runoff Modeling System models of 12 river basins in western Iowa will provide water-resource managers with a consistent and documented method for estimating streamflow at ungaged sites and aid in environmental studies, hydraulic design, water management, and water-quality projects.

Iowa↗

Bedform response to flow variability

Laboratory observations and computational results for the response of bedform fields to rapid variations in discharge are compared and discussed. The simple case considered here begins with a relatively low discharge over a flat bed on which bedforms are initiated, followed by a short high-flow period with double the original discharge, during which the morphology of the bedforms adjusts, followed in turn by a relatively long period of the original low discharge. For the grain size and hydraulic conditions selected, the Froude number remains subcritical during the experiment, and sediment moves predominantly as bedload. Observations show rapid development of quasi-two-dimensional bedforms during the initial period of low flow with increasing wavelength and height over the initial low-flow period. When the flow increases, the bedforms rapidly increase in wavelength and height, as expected from other empirical results. When the flow decreases back to the original discharge, the height of the bedforms quickly decreases in response, but the wavelength decreases much more slowly. Computational results using an unsteady two-dimensional flow model coupled to a disequilibrium bedload transport model for the same conditions simulate the formation and initial growth of the bedforms fairly accurately and also predict an increase in dimensions during the high-flow period. However, the computational model predicts a much slower rate of wavelength increase, and also performs less accurately during the final low-flow period, where the wavelength remains essentially constant, rather than decreasing. In addition, the numerical results show less variability in bedform wavelength and height than the measured values; the bedform shape is also somewhat different. Based on observations, these discrepancies may result from the simplified model for sediment particle step lengths used in the computational approach. Experiments show that the particle step length varies spatially and temporally over the bedforms during the evolution process. Assuming a constant value for the step length neglects the role of flow alterations in the bedload sediment-transport process, which appears to result in predicted bedform wavelength changes smaller than those observed. However, observations also suggest that three-dimensional effects play at least some role in the decrease of bedform wavelength, so incorporating better models for particle hop lengths alone may not be sufficient to improve model predictions.

Earth Surface Processes and Landforms↗

Archean and Proterozoic crustal evolution: Evidence from crustal seismology

Seismic-velocity models for Archean and Proterozoic provinces throughout the world are analyzed. The thickness of the crust in Archean provinces is generally found to be about 35 km (except at collisional boundaries), whereas Proterozoic crust has a significantly greater thickness of about 45 km and has a substantially thicker high-velocity (>7.0 km/s) layer at the base. We consider two models that may explain these differences. The first model attributes the difference to a change in the composition of the upper mantle. The higher temperatures in the Archean mantle led to the eruption of komatiitic lavas, resulting in an ultradepleted lithosphere unable to produce significant volumes of basaltic melt. Proterozoic crust developed above fertile mantle, and subsequent partial melting resulted in basaltic underplating and crustal inflation. In the second model, convection in the hot Archean mantle is considered to have been too turbulent to sustain stable long-lived subduction zones. By the Proterozoic the mantle had cooled sufficiently for substantial island and continental arcs to be constructed, and the high-velocity basal layer was formed by basaltic underplating.

Geology↗

A hierarchical Bayesian approach for handling missing classification data

Ecologists use classifications of individuals in categories to understand composition of populations and communities. These categories might be defined by demographics, functional traits, or species. Assignment of categories is often imperfect, but frequently treated as observations without error. When individuals are observed but not classified, these “partial” observations must be modified to include the missing data mechanism to avoid spurious inference. We developed two hierarchical Bayesian models to overcome the assumption of perfect assignment to mutually exclusive categories in the multinomial distribution of categorical counts, when classifications are missing. These models incorporate auxiliary information to adjust the posterior distributions of the proportions of membership in categories. In one model, we use an empirical Bayes approach, where a subset of data from one year serves as a prior for the missing data the next. In the other approach, we use a small random sample of data within a year to inform the distribution of the missing data. We performed a simulation to show the bias that occurs when partial observations were ignored and demonstrated the altered inference for the estimation of demographic ratios. We applied our models to demographic classifications of elk ( Cervus elaphus nelsoni ) to demonstrate improved inference for the proportions of sex and stage classes. We developed multiple modeling approaches using a generalizable nested multinomial structure to account for partially observed data that were missing not at random for classification counts. Accounting for classification uncertainty is important to accurately understand the composition of populations and communities in ecological studies.

Ecology and Evolution↗

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey ( Petromyzon marinus ). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site‐associated DNA sequencing (RAD‐Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00–0.018) were concordant with previous microsatellite‐based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age‐1 and age‐2 families of full and half‐siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non‐native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.

Michigan, Ontario, Wisconsin↗

Divergent evolution in fluviokarst landscapes of central Kentucky

Central Kentucky is characterized by a mixture of karst and fluvial features, typically manifested as mosaic of karst-rich/ channel-poor (KRCP) and channel-rich/karst-poor (CRKP) environments. At the regional scale the location and distribution of KRCP and CRKP areas are not always systematically related to structural, lithological, topographic, or other controls. This study examines the relationship of KRCP and CRKP zones along the Kentucky River gorge area, where rapid incision in the last 1??5 million years has lowered local base levels and modified slopes on the edge of the inner bluegrass plateau. At the scale of detailed field mapping on foot within a 4 km2 area, the development of karst and fluvial features is controlled by highly localized structural and topographic constraints, and can be related to slope changes associated with retreat of the Kentucky River gorge escarpment. A conceptual model of karst/fluvial transitions is presented, which suggests that minor, localized variations are sufficient to trigger a karst-fluvial or fluvial-karst switch when critical slope thresholds are crossed. ?? 2004 John Wiley and Sons, Ltd.

Earth Surface Processes and Landforms↗

Trophic strategies, animal diversity and body size

A primary difference between predators and parasites is the number of victims that an individual attacks throughout a life-history stage. A key division within natural enemies is whether a successful attack eliminates the fitness of the prey or the host. A third distinctive axis for parasites is whether the host must die to further parasite development. The presence or absence of intensity-dependent pathology is a fourth factor that separates macroparasites from microparasites; this also distinguishes between social and solitary predators. Combining these four dichotomies defines seven types of parasitism, seven corresponding parasites, three forms of predation and, when one considers obligate and facultative combinations of these forms, four types of predator. Here, we argue that the energetics underlying the relative and absolute sizes of natural enemies and their victims is the primary selective factor responsible for the evolution of these different trophic strategies.

Trends in Ecology and Evolution↗

Debris-flow process controls on steepland morphology in the San Gabriel Mountains, California

Steep landscapes evolve largely by debris flows, in addition to fluvial and hillslope processes. Abundant field observations document that debris flows incise valley bottoms and transport substantial sediment volumes, yet their contributions to steepland morphology remain uncertain. This has, in turn, limited the development of debris-flow incision rate formulations that produce morphology consistent with natural landscapes. In many landscapes, including the San Gabriel Mountains (SGM), California, steady-state fluvial channel longitudinal profiles are concave-up and exhibit a power-law relationship between channel slope and drainage area. At low drainage areas, however, valley slopes become nearly constant. These topographic forms result in a characteristically curved slope-area signature in log-log space. Here, we use a one-dimensional landform evolution model that incorporates debris-flow erosion to reproduce the relationship between this curved slope-area signature and erosion rate in the SGM. Topographic analysis indicates that the drainage area at which steepland valleys transition to fluvial channels correlates with measured erosion rates in the SGM, and our model results reproduce these relationships. Further, the model only produces realistic valley profiles when parameters that dictate the relationship between debris-flow erosion, valley-bottom slope, and debris-flow depth are within a narrow range. This result helps place constraints on the mathematical form of a debris-flow incision law. Finally, modeled fluvial incision outpaces debris-flow erosion at drainage areas less than those at which valleys morphologically transition from near-invariant slopes to concave profiles. This result emphasizes the critical role of debris-flow incision for setting steepland form, even as fluvial incision becomes the dominant incisional process.

California↗

Improving ecological data science with workflow management software

Pressing environmental research questions demand the integration of increasingly diverse and large-scale ecological datasets as well as complex analytical methods, which require specialized tools and resources. Computational training for ecological and evolutionary sciences has become more abundant and accessible over the past decade, but tool development has outpaced the availability of specialized training. Most training for scripted analyses focuses on individual analysis steps in one script rather than creating a scripted pipeline, where modular functions comprise an ecosystem of interdependent steps. Although current computational training creates an excellent starting place, linear styles of scripting can risk becoming labor- and time-intensive and less reproducible by often requiring manual execution. Pipelines, however, can be easily automated or tracked by software to increase efficiency and reduce potential errors. Ecology and evolution would benefit from techniques that reduce these risks by managing analytical pipelines in a modular, readily parallelizable format with clear documentation of dependencies. Workflow management software (WMS) can aid in the reproducibility, intelligibility and computational efficiency of complex pipelines. To date, WMS adoption in ecology and evolutionary research has been slow. We discuss the benefits and challenges of implementing WMS and illustrate its use through a case study with the targets r package to further highlight WMS benefits through workflow automation, dependency tracking and improved clarity for reviewers. Although WMS requires familiarity with function-oriented programming and careful planning for more advanced applications and pipeline sharing, investment in training will enable access to the benefits of WMS and impart transferable computing skills that can facilitate ecological and evolutionary data science at large scales.

Methods in Ecology and Evolution↗

Character and evolution of the ground-water flow system in the central part of the western San Joaquin Valley, California

The occurrence of selenium in agricultural drain water derived from the western San Joaquin Valley, California, has focused concern on the ground-water flow system of the western valley. In this investigation, previous work and recently collected texture and water-level data are used to evaluate the character and evolution of the regional ground-water flow system in the central part of the western valley, with particular emphasis on the deposits overlying the Corcoran Clay Member of the Tulare Formation. The Corcoran Clay Member, where present, divides the flow system into an upper semiconfined zone and a lower con-fined zone. Above the Corcoran, three hydrogeologic units can be recognized: Coast Range alluvium, Sierran sand, and flood-basin deposits. These units differ in texture, hydrologic properties, and oxidation state. The development of irrigated agriculture in the central part of the western valley has significantly altered the flow system. Percolation of irrigation water past crop roots has caused a rise in the altitude of the water table in midfan and distal-fan areas. Pumpage of ground water from wells has caused a lowering of the water table beneath parts of the fanheads and a lowering of the potentiometric surface of the confined zone over much of the western valley. The combination of percolation and pumpage has resulted in development of a large downward hydraulic head gradient in the semi-confined zone and has created a groundwater divide along the western margin of the valley. Surface-water deliveries from the California Aqueduct have allowed a decrease in pumpage and a consequent recovery in hydraulic head throughout the system.

California↗

Leaders' perspectives in the Yellowstone to Yukon Conservation Initiative

The Yellowstone to Yukon Conservation Initiative (Y2Y) was created in 1993 to advance conservation in a 1.2 million km2 portion of the North American Rocky Mountains. We assembled 21 people with influence over Y2Y in a workshop to elucidate perspectives on challenges and solutions for this organization at a key point in its evolution, and used Q method to define four perspectives on challenges and three on solutions. Participants were differentiated by four models for effecting change-vision-based advocacy, practice-based learning, political engagement, and scientific management-with emphasis on the first three. Those with authority in Y2Y aligned with vision-based advocacy and expressed ambivalence about practice-based adaptive learning and rigorous appraisals of existing strategies. Workshop results were consistent with an apparent trend toward organizational maturation focused on stabilizing revenues, developing formal organizational arrangements, and focusing strategies. Consolidation of power in Y2Y around a long-standing formula does not bode well for the effectiveness of Y2Y. We recommend that leaders in Y2Y and similar organizations focused on large-scale conservation to create and maintain an open system-philosophically and operationally-that capitalizes on the diverse perspectives and skills of individuals who are attracted to such efforts. We also recommend that the Y2Y initiative be followed closely to harvest additional lessons for potential application to large-scale conservation efforts elsewhere. ?? Springer Science+Business Media, LLC(outside the USA).2011.

Policy Sciences↗