Design principles for origin-destination flow maps

Origin-destination flow maps are often difficult to read due to overlapping flows. Cartographers have developed design principles in manual cartography for origin-destination flow maps to reduce overlaps and increase readability. These design principles are identified and documented using a quantitative content analysis of 97 geographic origin-destination flow maps without branching or merging flows. The effectiveness of selected design principles is verified in a user study with 215 participants. Findings show that (a) curved flows are more effective than straight flows, (b) arrows indicate direction more effectively than tapered line widths, and (c) flows between nodes are more effective than flows between areas. These findings, combined with results from user studies in graph drawing, conclude that effective and efficient origin-destination flow maps should be designed according to the following design principles: overlaps between flows are minimized; symmetric flows are preferred to asymmetric flows; longer flows are curved more than shorter or peripheral flows; acute angles between crossing flows are avoided; sharp bends in flow lines are avoided; flows do not pass under unconnected nodes; flows are radially distributed around nodes; flow direction is indicated with arrowheads; and flow width is scaled with represented quantity.     

Climate change links fate of glaciers and an endemic alpine invertebrate

Climate warming in the mid- to high-latitudes and high-elevation mountainous regions is occurring more rapidly than anywhere else on Earth, causing extensive loss of glaciers and snowpack. However, little is known about the effects of climate change on alpine stream biota, especially invertebrates. Here, we show a strong linkage between regional climate change and the fundamental niche of a rare aquatic invertebrate—the meltwater stonefly Lednia tumana—endemic to Waterton-Glacier International Peace Park, Canada and USA. L. tumana has been petitioned for listing under the U.S. Endangered Species Act due to climate-change-induced glacier loss, yet little is known on specifically how climate impacts may threaten this rare species and many other enigmatic alpine aquatic species worldwide. During 14 years of research, we documented that L. tumana inhabits a narrow distribution, restricted to short sections (~500 m) of cold, alpine streams directly below glaciers, permanent snowfields, and springs. Our simulation models suggest that climate change threatens the potential future distribution of these sensitive habitats and the persistence of L. tumana through the loss of glaciers and snowfields. Mountaintop aquatic invertebrates are ideal early warning indicators of climate warming in mountain ecosystems. Research on alpine invertebrates is urgently needed to avoid extinctions and ecosystem change.     

Scaling Ecological Dynamics: Self-Organization, Hierarchical Structure, and Ecological Resilience

Assessing impacts of global change is complicated by the problems associated with translating models and data across spatial and temporal scales. One of the major problems of ecological scaling is the dynamic, self-organized nature of ecosystems. Ecological organization emerges from the interaction of structures and processes operating at different scales. The resilience of ecological organization to changes in key cross-scale processes can be used to assess the contexts within which scaling methods function well, need adjustment, and break down.     

U.S. National Forests adapt to climate change through Science–Management partnerships

Developing appropriate management options for adapting to climate change is a new challenge for land managers, and integration of climate change concepts into operational management and planning on United States national forests is just starting. We established science–management partnerships on the Olympic National Forest (Washington) and Tahoe National Forest (California) in the first effort to develop adaptation options for specific national forests. We employed a focus group process in order to establish the scientific context necessary for understanding climate change and its anticipated effects, and to develop specific options for adapting to a warmer climate. Climate change scientists provided the scientific knowledge base on which adaptations could be based, and resource managers developed adaptation options based on their understanding of ecosystem structure, function, and management. General adaptation strategies developed by national forest managers include: (1) reduce vulnerability to anticipated climate-induced stress by increasing resilience at large spatial scales, (2) consider tradeoffs and conflicts that may affect adaptation success, (3) manage for realistic outcomes and prioritize treatments that facilitate adaptation to a warmer climate, (4) manage dynamically and experimentally, and (5) manage for structure and composition. Specific adaptation options include: (1) increase landscape diversity, (2) maintain biological diversity, (3) implement early detection/rapid response for exotic species and undesirable resource conditions, (4) treat large-scale disturbance as a management opportunity and integrate it in planning, (5) implement treatments that confer resilience at large spatial scales, (6) match engineering of infrastructure to expected future conditions, (7) promote education and awareness about climate change among resource staff and local publics, and (8) collaborate with a variety of partners on adaptation strategies and to promote ecoregional management. The process described here can quickly elicit a large amount of information relevant for adaptation to climate change, and can be emulated for other national forests, groups of national forests with similar resources, and other public lands. As adaptation options are iteratively generated for additional administrative units on public lands, management options can be compared, tested, and integrated into adaptive management. Science-based adaptation is imperative because increasing certainty about climate impacts and management outcomes may take decades.     

A critique of the internal tracer method for estimating contaminant degradation rates

The internal tracer method for estimating contaminant degradation rates separates the attenuation effects not associated with degradation by using a codisposed recalcitrant internal tracer to normalize the degrading contaminant concentration. The remaining attenuation between the internal tracer and degrading contaminant is attributed to degradation and the degradation rate half-life is estimated from the first-order decay equation. An analytical solution of the advection- dispersion equation was used to evaluate flow-and-transport conditions that could result in incorrect estimates of contaminant degradation rate constants. Flow-and-transport characteristics that result in overestimating degradation rates were of particular interest because the internal tracer method often used to demonstrate natural attenuation can achieve remedial objectives. The analytical solution was also used to estimate the magnitude of error associated with using the internal tracer method at an example site and to explain different degradation rates estimated using tracers with different decay rate constants.     

Chemical equilibrium and mass balance relationships associated with the Long Valley hydrothermal system, California, U.S.A.

Recent drilling and sampling of hydrothermal fluids from Long Valley permit an accurate characterization of chemical concentrations and equilibrium conditions in the hydrothermal reservoir. Hydrothermal fluids are thermodynamically saturated with secondary quartz, calcite, and pyrite but are in disequilibrium with respect to aqueous sulfide-sulfate speciation. Hydrothermal fluids are enriched in ¹⁸O by approximately 1‰ relative to recharge waters. ¹⁸O and Cl concentrations in well cuttings and core from high-temperature zones of the reservoir are extensively depleted relative to fresh rhyolitic tuff compositions. Approximately 80% of the Li and 50% of the B are retained in the altered reservoir rock. Cl mass balance and open-system ¹⁸O fractionation models produce similar water-rock ratios of between 1.0 and 2.5 kg kg⁻¹. These water-rock ratios coupled with estimates of reservoir porosity and density produce a minimum fluid residence time of 1.3 ka. The low fluid Cl concentrations in Long Valley correlate with corresponding low rock concentrations. Mass balance calculations indicate that leaching of these reservoir rocks accounts for Cl losses during hydrothermal activity over the last 40 ka.     

Gravity variations along the Southeast Bohol Ophiolite Complex (SEBOC), Central Philippines: Implications on Ophiolite Emplacement

Abstract The basement complex of Bohol Island consists of the Southeast Bohol Ophiolite Complex (SEBOC), Cansiwang Melange and Alicia Schist. The SEBOC is a complete, but dismembered ophiolite with outcrops generally trending northeast– southwest and dipping north-west. The harzburgite units of the SEBOC are almost always observed to be thrusted onto the Cansiwang Melange, which in turn is thrusted onto the Alicia Schist. Bouguer gravity values on Bohol range from about +60 mGal in the west to +120 mGal in the east, in the region to the north-east of the SEBOC outcrops. Based on the present distribution of the SEBOC units and their thrust fault relationship with the Cansiwang Melange and Alicia Schist, it is proposed that the SEBOC was emplaced by onramping towards the south-eastward direction. However, the orientation of the Bouguer highs suggests that the thrusting direction of the ophiolite units is towards the south-west and not towards the south-east.     

The long-term salinity field in San Francisco Bay

Data are presented on long-term salinity behaviour in San Francisco Bay, California. A two-level, width averaged model of the tidally averaged salinity and circulation has been written in order to interpret the long-term (days to decades) salinity variability. The model has been used to simulate daily averaged salinity in the upper and lower levels of a 51 segment discretization of the Bay over the 22-yr period 1967–1988. Monthly averaged surface salinity from observations and monthly-averaged simulated salinity are in reasonable agreement. Good agreement is obtained from comparison with daily averaged salinity measured in the upper reaches of North Bay.The salinity variability is driven primarily by freshwater inflow with relatively minor oceanic influence. All stations exhibit a marked seasonal cycle in accordance with the Mediterranean climate, as well as a rich spectrum of variability due to extreme inflow events and extended periods of drought. Monthly averaged salinity intrusion positions have a pronounced seasonal variability and show an approximately linear response to the logarithm of monthly averaged Delta inflow. Although few observed data are available for studies of long-term salinity stratification, modelled stratification is found to be strongly dependent on freshwater inflow; the nature of that dependence varies throughout the Bay. Near the Golden Gate, stratification tends to increase up to very high inflows. In the central reaches of North Bay, modelled stratification maximizes as a function of inflow and further inflow reduces stratification. Near the head of North Bay, lowest summer inflows are associated with the greatest modelled stratification. Observations from the central reaches of North Bay show marked spring-neap variations in stratification and gravitational circulation, both being stronger at neap tides. This spring-neap variation is simulated by the model. A feature of the modelled stratification is a hysteresis in which, for a given spring-neap tidal range and fairly steady inflows, the stratification is higher progressing from neaps to springs than from springs to neaps.The simulated responses of the Bay to perturbations in coastal sea salinity and Delta inflow have been used to further delineate the time-scales of salinity variability. Simulations have been performed about low inflow, steady-state conditions for both salinity and Delta inflow perturbations. For salinity perturbations a small, sinusoidal salinity signal with a period of 1 yr has been applied at the coastal boundary as well as a pulse of salinity with a duration of one day. For Delta inflow perturbations a small, sinusoidally varying inflow signal with a period of 1 yr has been superimposed on an otherwise constant Delta inflow, as well as a pulse of inflow with a duration of one day. Perturbations in coastal salinity dissipate as they move through the Bay. Seasonal perturbations require about 40–45 days to propagate from the coastal ocean to the Delta and to the head of South Bay. The response times of the model to perturbations in freshwater inflow are faster than this in North Bay and comparable in South Bay. In North Bay, time-scales are consistent with advection due to lower level, up-estuary transport of coastal salinity perturbations; for inflow perturbations, faster response times arise from both upper level, down-estuary advection and much faster, down-estuary migration of isohalines in response to inflow volume continuity. In South Bay, the dominant time-scales are governed by tidal dispersion.     

Geochemistry and origin of the Proterozoic ultrapotassic rocks of the Churchill Province,Canada

Summary Early Proterozoic ultrapotassic dikes, lava flows, and pyroclastic rocks of the Christopher Island Formation (CIF) erupted throughout an area 600 × 300 km within the Churchill Province of the Canadian Shield at 1.84 Ga. The rocks range from mafic lamprophyres (mg # 60; SiO₂ 47–54%, mean K₂O/Na₂O > 4) with phenocrysts of phlogopite + diopside + apatite ± olivine ± magnetite, to phenocryst-poor felsic rocks and sanidine porphyries (SiO₂55–69%). All samples have high incompatible element contents and display large depletions of high field strength elements relative to K, Rb, Sr, Ba, and Th. The CIF has geochemical and petrographic characteristics of both minettes and lamproites, but overall most closely resembles young Mediterranean lamproites. Felsic rocks of the CIF were produced by crystal fractionation and crustal contamination of mafic ultrapotassic magma, and include both high-silica lamproites strongly enriched in Zr, U, and Th, and weakly potassic to sodic rocks of trachytic composition. Flows and feeder dikes have relatively homogeneous _{Nd, 1840 Ma} (–6 to –11) but highly variable ES., 1840 Ma (–40 to + 100); samples classified as lamproites have higher average _Sr. Dike samples have highly variable present-day Pb isotope compositions, ranging from moderately to strongly nonradiogenic. Geochemical and isotopic data are consistent with contributions from depleted Archean lithospheric mantle, and OIB-type convecting mantle, both metasomatized by subduction-related processes during the Early Proterozoic. The lithospheric mantle probably contained Archean enriched domains as well. Proterozoic enrichment may have accompanied shallow underplating of subducted oceanic lithosphere beneath the Churchill Province during amalgamation of the Laurentian supercontinent. There are strong analogies in isotopic composition, and interpreted source region history, between the CIF and lamproites and minettes of the Wyoming Province and western Greenland, which suggest the existence of a Laurentian ultrapotassic superprovince.

---

Geochemie und Entstehung der Proterozoischen ultrapotassischen Gesteine der Churchill Provinz, Kanada

Zusammenfassung Altproterozoische, ultrapotassische Gänge, Lavaströme und pyroklastische Gesteine der Christopher Island Formation (CIF), eruptierten in einem Gebiet von 600 × 300 km in der Churchill Provinz des Kanadischen Schildes vor 1.84 Ga. Die Zusammensetzung dieser Gesteine variiert von mafischen Lamprophyren (mg > 60; SiO₂ = 47–54%, durchschnittliches K₂O/Na₂O > 4) mit Phänokristallent von Phlogopit + Diopsid + Apatit + Olivin + Magnetit, bis zu phänokristallarmen felsischen Gesteinen und Sanidinporphyren (SiO₂ = 55–69%). Alle Proben zeigen hohe Gehalte an inkompatiblen Elementen und zeigen beträchtliche Verarmung an high field strength Elementen relativ zu K, Rb, Sr, Ba und Th. Die CIF hat geochemische und petrographische Eigenschaften sowohl von Minetten wie von Lamproiten, aber im allgemeinen ähnelt sie am stärksten jungen mediterranen Lamproiten. Felsische Gesteine der CIF wurden durch Fraktionierung und Krustenkontamination aus mafischen ultrapotassischen Magmen gebildet. Letztere umfassen sowohl siliziumreiche Lamproite, die deutlich an Zr, U und Th angereichert sind und schwach potassische bis sodische Gesteine von trachytischer Zusammensetzung. Lavenergüsse und zufuhrgänge zeigent ein relativ homogenes _{Nd, 1840 Ma} (–6 bis –11) aber ein sehr variables _{Sr, 1840 Ma} (-40 bis + 100); Proben die als Lamproite klassifiziert wurden, zeigent höhere durchschnittliche _Sr-Werte. Proben von Gängen haben sehr variable Bleiisotopen-Zusammensetzungen, die von mäßig bis stark nichtradiogen variieren. Geochemische und Isotopendaten weisen auf Beiträge aus verarmtem archaischen lithosphärischen Mantel und aus konvektierendem OIB-Typ Mantel hin, die beide während des Alproterozoikums durch Subduktions-Vorgänge metasomatisiert wurden. Der lithosphärische Mantel enthielt wahrscheinlich auch angereicherte archaische Domänen. Proterozoische Anreicherungsvorgänge dürften seichtes Underplating subduzierter ozeanischer Lithosphäre unter der Churchill Provinz während der Amalgamation des laurentischen Superkontinentes begleitet haben. Es gibt starke Analogien in der Isotopenzusammensetzung und in der interpretierten Geschichte der Ursprungsregion, zwischen den CIF und Lamproiten und Minetten der Wyoming Provinz, und des westlichen Grönland. Diese weisen auf die Existenz einer laurentischen ultrapotassischen Superprovinz hin.

---

---

With 7 Figures