ALBERT

Description: Tree species are predicted to track future climate by shifting their geographic distributions, but climate-mediated migrations are not apparent in a recent continental-scale analysis. To better understand the mechanisms of a possible migration lag, we analyzed relative recruitment patterns by comparing juvenile and adult tree abundances in climate space. One would expect relative recruitment to be higher in cold and dry climates as a result of tree migration with juveniles located further poleward than adults. Alternatively, relative recruitment could be higher in warm and wet climates as a result of higher tree population turnover with increased temperature and precipitation. Using the USDA Forest Service's Forest Inventory and Analysis data at regional scales, we jointly modeled juvenile and adult abundance distributions for 65 tree species in climate space of the eastern United States. We directly compared the optimal climate conditions for juveniles and adults, identified the climates where each species has high relative recruitment, and synthesized relative recruitment patterns across species. Results suggest that for 77% and 83% of the tree species, juveniles have higher optimal temperature and optimal precipitation, respectively, than adults. Across species, the relative recruitment pattern is dominated by relatively more abundant juveniles than adults in warm and wet climates. These different abundance-climate responses through life history are consistent with faster population turnover and inconsistent with the geographic trend of large-scale tree migration. Taken together, this juvenile-adult analysis suggests that tree species might respond to climate change by having faster turnover as dynamics respond to longer growing seasons and higher temperatures, before there is evidence of poleward migration at biogeographic scales. This article is protected by copyright. All rights reserved.

Description: [1]  We investigated the evolution of seismicity and deformation in the unstable eastern flank of Etna volcano over a thirty-year period (from 1980 to 2012). A significant temporal correlation has been revealed between periods of flank acceleration and intensified seismic activity by comparing seismicity along the northern border (Pernicana fault system) of the sliding flank and the deformation of the eastern flank. Two marked phases have been observed in 1984-1986 and in the years following 2002. These two phases are separated by an intermediate phase from 1987 to 2001, in which the flank sliding slowed down and the seismicity dropped drastically. This common temporal evolution in the deformation rate and seismic release supports the hypothesis that the seismicity in the northern border can be viewed as a marker of the response to accommodate the stress exerted by the traction of the eastern flank sliding. This interplay has also been corroborated by Finite Element Method (FEM) numerical computations that highlight a good correlation between the seismicity pattern and areas of positive stress changes induced by the sliding surface. The two intense phases of flank acceleration are representative of two main different sources: volcano flank instability stretching the eastern sector in the first 1984-1986 phase and magmatic intrusions pushing the eastern flank seaward since the 2002-03 eruption. Establishing the relationship between flank acceleration and seismic activation, therefore, contributes to understanding Etna's mechanical behavior, and provides insights into the processes regulating the unstable flank response.

Description: [1]  Numerous studies in the Central Pyrenees have provided evidence for a rapid phase of exhumation of this mountain belt during the Late Eocene (37–30 Ma). Simultaneously, the closure of the Ebro foreland basin allowed the accumulation of sediments at the southern Piedmont, which partially covered the fold-and-thrust belt from Late Eocene ( e . g . when it was still actively deforming) to Miocene times. We aim here at understanding the consequences of such syn-tectonic sedimentation on the Southern Pyrenean fold-and-thrust belt by using a 2-D numerical model that reproduces the development of a thin-skinned wedge subject to different modes of sedimentation and erosion. The results show contrasting fold-and-thrust belt behavior when applying aggrading or prograding sedimentation, which we link to the critical state of the wedge. When the sediments are sourced from the hinterland (progradation), the thrusting propagates toward the foreland; whereas when the sediments aggrade from the basin, the thrusting sequence migrates backward. This latter mode shows patterns of deformation that compare favorably to the Pyrenean thrusting sequence observed during Eocene-Miocene times.

Description: [1]  In southern California, fault slip rate estimates along the San Andreas fault (SAF) and Garlock fault from geodetically-constrained kinematic models are systematically at the low end or lower than geologic slip rate estimates. The sum of geodetic model slip rates across the Eastern California Shear Zone is higher than the geologic sum. However, the ranges of reported model and geologic slip rate estimates in the literature are sufficiently large that it remains unclear whether these apparent discrepancies are real, or attributable to epistemic uncertainties in the two types of estimates. We further examine uncertainties in geodetically-derived slip rate estimates on major faults in southern California by conducting a suite of inversions with four kinematic models. Long-term-rigid elastic block models constrained by the geologic slip rates cannot fit the present-day GPS-derived velocity field. Deforming (permanent off-fault strain) elastic block models and viscoelastic earthquake cycle block models constrained by geologic slip rates can fit the present-day GPS-derived velocity field with 28-33% of the total geodetic moment rate occurring as distributed deformation off of the major faults. Models incorporating viscoelastic mantle flow predict systematically higher slip rates than purely elastic models on many of the the major southern California faults with ranges of (elastic/viscoelastic) 29-34/30-37 mm/yr for the Carrizo SAF segment, 20-24/20-32 mm/yr for the Mojave SAF segment, 14-17/18-22 mm/yr for the Coachella SAF segment, 13-19/14-22 mm/yr for the San Jacinto fault, and 5-11/5-11 mm/yr for the western Garlock fault.

Description: Climate warming threatens to increase mass coral bleaching events, and several studies have projected the demise of tropical coral reefs this century. However, recent evidence indicates corals may be able to respond to thermal stress though adaptive processes (e.g., genetic adaptation, acclimatization, and symbiont shuffling). How these mechanisms might influence warming induced bleaching is largely unknown. This study compared how different adaptive processes could affect coral bleaching projections. We used the latest bias-corrected global sea surface temperature (SST) output from the NOAA/GFDL Earth System Model 2 (ESM2M) for the pre-industrial period though 2100 to project coral bleaching trajectories. Initial results showed that, in the absence of adaptive processes, application of a pre-industrial climatology to the NOAA Coral Reef Watch bleaching prediction method over-predicts the present day bleaching frequency. This suggests that corals may have already responded adaptively to some warming over the industrial period. We then modified the prediction method so that the bleaching threshold either permanently increased in response to thermal history (e.g., simulating directional genetic selection) or temporarily increased for 2-10 years in response to a bleaching event (e.g., simulating symbiont shuffling). A bleaching threshold that changes relative to the preceding 60 years of thermal history reduced the frequency of mass bleaching events by 20-80% compared with the ‘no adaptive response’ prediction model by 2100, depending on the emissions scenario. When both types of adaptive responses were applied, up to 14% more reef cells avoided high frequency bleaching by 2100. However, temporary increases in bleaching thresholds alone only delayed the occurrence of high frequency bleaching by ~10 years in all but the lowest emissions scenario. Future research should test the rate and limit of different adaptive responses for coral species across latitudes and ocean basins to determine if and how much corals can respond to increasing thermal stress. This article is protected by copyright. All rights reserved.

Description: The snow-masking effect of vegetation exerts strong control on albedo in northern high latitude ecosystems. Large-scale changes in the distribution and stature of vegetation in this region will thus have important feedbacks to climate. The snow-albedo feedback is controlled largely by the contrast between snow-covered and snow-free albedo (Δα), which influences predictions of future warming in coupled climate models, despite being poorly constrained at seasonal and century time scales. Here we compare satellite observations and coupled climate model representations of albedo and tree cover for the boreal and Arctic region. Our analyses reveal consistent declines in albedo with increasing tree cover, occurring south of latitudinal tree line, that are poorly represented in coupled climate models. Observed relationships between albedo and tree cover differ substantially between snow-covered and snow-free periods, and among plant functional type (PFT). Tree cover in models varies widely but surprisingly does not correlate well with model albedo. Further, our results demonstrate a relationship between tree cover and snow-albedo feedback that may be used to accurately constrain high latitude albedo feedbacks in coupled climate models under current and future vegetation distributions. This article is protected by copyright. All rights reserved.

Description: Because of global land surface warming, extreme temperature events are expected to occur more often and more intensely, affecting the growth and development of the major cereal crops in several ways, thus affecting the production component of food security. In this paper, we have identified rice and maize crop responses to temperature in different, but consistent, phenological phases and development stages. A literature review and data compilation of around 140 scientific articles have determined the key temperature thresholds and response to extreme temperature effects for rice and maize, complementing an earlier study on wheat. Lethal temperatures and cardinal temperatures, together with error estimates, have been identified for phenological phases and development stages. Following the methodology of previous work, we have collected and statistically analysed temperature thresholds of the three crops for the key physiological processes such as leaf initiation, shoot growth and root growth and for the most susceptible phenological phases such as sowing to emergence, anthesis and grain filling. Our summary shows that cardinal temperatures are conservative between studies and are seemingly well-defined in all three crops. Anthesis and ripening are the most sensitive temperature stages in rice as well as in wheat and maize. We call for further experimental studies of the effects of transgressing threshold temperatures so such responses can be included into crop impact and adaptation models. This article is protected by copyright. All rights reserved.

Description: Will warming lead to an increased use of older soil organic carbon (SOC) by microbial communities, thereby inducing C losses from C-rich alpine soils? We studied soil microbial community composition, activity and substrate use after three and four years of soil warming (+4°C, 2007-2010) at the alpine treeline in Switzerland. The warming experiment was nested in a free air CO 2 enrichment experiment using depleted 13 CO 2 (δ 13 C = –30‰, 2001-2009). We traced this depleted 13 C label in phospholipid fatty acids (PLFA) of the organic layer (0-5 cm soil depth) and in C mineralized from root-free soils to distinguish substrate ages used by soil microorganisms: fixed before 2001 (“old”), from 2001 to 2009 (“new”) or in 2010 (“recent”). Warming induced a sustained stimulation of soil respiration (+38%) without decline in mineralizable SOC. PLFA concentrations did not reveal changes in microbial community composition due to soil warming, but soil microbial metabolic activity was stimulated (+66%). Warming decreased the amount of new and recent C in the fungal biomarker 18:2ω6,9 and the amount of new C mineralized from root-free soils, implying a shift in microbial substrate use towards a greater use of old SOC. This shift in substrate use could indicate an imbalance between C inputs and outputs, which could eventually decrease SOC storage in this alpine ecosystem. This article is protected by copyright. All rights reserved.

Description: Circumpolar expansion of tall shrubs and trees into Arctic tundra is widely thought to be occurring as a result of recent climate warming, but little quantitative evidence exists for northern Siberia, which encompasses the world's largest forest-tundra ecotonal belt. We quantified changes in tall shrub and tree canopy cover in eleven, widely-distributed Siberian ecotonal landscapes by comparing very-high-resolution photography from the Cold War-era “Gambit” and “Corona” satellite surveillance systems (1965-1969) with modern imagery. We also analyzed within-landscape patterns of vegetation change to evaluate the susceptibility of different landscape components to tall shrub and tree increase. The total cover of tall shrubs and trees increased in nine of eleven ecotones. In northwest Siberia, alder ( Alnus ) shrubland cover increased 5.3 – 25.9% in five ecotones. In Taymyr and Yakutia, larch ( Larix ) cover increased 3.0 – 6.7% within three ecotones, but declined 16.8% at a fourth ecotone due to thaw of ice-rich permafrost. In Chukotka, the total cover of alder and dwarf pine ( Pinus ) increased 6.1% within one ecotone and was little-changed at a second ecotone. Within most landscapes, shrub and tree increase was linked to specific geomorphic settings, especially those with active disturbance regimes such as permafrost patterned-ground, floodplains, and colluvial hillslopes. Mean summer temperatures increased at most ecotones since the mid-1960s, but rates of shrub and tree canopy cover expansion were not strongly correlated with temperature trends and were better correlated with mean annual precipitation. We conclude that shrub and tree cover is increasing in tundra ecotones across most of northern Siberia, but rates of increase vary widely regionally and at the landscape-scale. Our results indicate that extensive changes can occur within decades in moist, shrub-dominated ecotones, as in northwest Siberia, while changes are likely to occur much more slowly in the highly continental, larch-dominated ecotones of central and eastern Siberia. This article is protected by copyright. All rights reserved.

Description: The adaptation of different species to warming temperatures has been increasingly studied. Moose ( Alces alces ) is the largest of the ungulate species occupying the northern latitudes across the globe, and in Finland it is the most important game species. It is very well adapted to severe cold temperatures, but has a relatively low tolerance to warm temperatures. Previous studies have documented changes in habitat use by moose due to high temperatures. In many of these studies the used areas have been classified according to how much thermal cover they were assumed to offer based on satellite/aerial imagery data. Here, we identified the vegetation structure in the areas used by moose under different thermal conditions. For this purpose we used airborne laser scanning (ALS) data extracted from the locations of GPS-collared moose. This provided us with detailed information about the relationships between moose and the structure of forests it uses in different thermal conditions and we were therefore able to determine and differentiate between the canopy structures at locations occupied by moose during different thermal conditions. We also discovered a threshold beyond which moose behaviour began to change significantly: as day temperatures began to reach 20 ○ C and higher, the search for areas with higher and denser canopies during daytime became evident. The difference was clear when compared to habitat use at lower temperatures, and was so strong that it provides supporting evidence to previous studies, suggesting that moose are able to modify their behaviour to cope with high temperatures, but also that the species is likely to be affected by warming climate. This article is protected by copyright. All rights reserved.

Description: Mountain ecosystems are particularly susceptible to climate change. Characterizing intraspecific variation of alpine plants along elevational gradients is crucial for estimating their vulnerability to predicted changes. Environmental conditions vary with elevation, which might influence plastic responses and affect selection pressures that lead to local adaptation. Thus, local adaptation and phenotypic plasticity among low and high elevation plant populations in response to climate, soil and other factors associated with elevational gradients might underlie different responses of these populations to climate warming. Using a transplant experiment along an elevational gradient, we investigated reproductive phenology, growth and reproduction of the nutrient-poor grassland species Ranunculus bulbosus , Trifolium montanum , and Briza media . Seeds were collected from low and high elevation source populations across the Swiss Alps and grown in nine common gardens at three different elevations with two different soil depths. Despite genetic differentiation in some traits, the results revealed no indication of local adaptation to the elevation of population origin. Reproductive phenology was advanced at lower elevation in low and high elevation populations of all three species. Growth and reproduction of T. montanum and B. media were rarely affected by garden elevation and soil depth. In R. bulbosus , however, growth decreased and reproductive investment increased at higher elevation. Furthermore, soil depth influenced growth and reproduction of low elevation R. bulbosus populations. We found no evidence for local adaptation to elevation of origin and hardly any differences in the responses of low and high elevation populations. However, the consistent advanced reproductive phenology observed in all three species shows that they have the potential to plastically respond to environmental variation. We conclude that populations might not be forced to migrate to higher elevations as a consequence of climate warming, as plasticity will buffer the detrimental effects of climate change in the three investigated nutrient-poor grassland species. This article is protected by copyright. All rights reserved.

Description: The impact of climate change on the stability of soil organic carbon (SOC) remains a major source of uncertainty in predicting future changes in atmospheric CO 2 levels. One unsettled issue is whether the mineralization response to temperature depends on SOC mineralization rate. Long-term (〉25 years) bare fallow experiments (LTBF) in which the soil is kept free of any vegetation and organic inputs, and their associated archives of soil samples represent a unique research platform to examine this issue as with increasing duration of fallow, the lability of remaining total SOC decreases. We retrieved soils from LTBF experiments situated at Askov (Denmark), Grignon (France), Ultuna (Sweden) and Versailles (France) and sampled at the start of the experiments and after 25, 50, 52, and 79 years of bare fallow, respectively. Soils were incubated at 4, 12, 20 and 35 °C and the evolved CO 2 monitored. The apparent activation energy ( Ea ) of SOC was then calculated for similar loss of CO 2 at the different temperatures. The Ea was always higher for samples taken at the end of the bare-fallow period, implying a higher temperature sensitivity of stable C than of labile C. Our results provide strong evidence for a general relationship between temperature sensitivity and SOC stability upon which significant improvements in predictive models could be based. This article is protected by copyright. All rights reserved.

Description: The temperature dependence of aerobic scope has been suggested to be a major determinant of how marine animals will cope with future rises in environmental temperature. Here we present data suggesting that in some animals, the temperature dependence of anaerobic scope (i.e. the capacity for surviving severe hypoxia) may determine present-day latitudinal distributions and potential for persistence in a warmer future. As a model for investigating the role of anaerobic scope, we studied two sibling species of coral-dwelling gobies, Gobiodon histrio and G. erythrospilus , with different latitudinal distributions, but which overlap in equal abundance at Lizard Island (14°40'S) on the Great Barrier Reef. These species did not differ in the temperature dependence of resting oxygen consumption or critical oxygen concentration (the lowest oxygen level where resting oxygen consumption can be maintained). By contrast, the more equatorial species ( G. histrio ) had a better capacity to endure anaerobic conditions at oxygen levels below the critical oxygen concentration at the high temperatures (32 – 33 °C) more likely to occur near the equator, or in a warmer future. These results suggest that anaerobic scope, in addition to aerobic scope, could be important in determining the impacts of global warming on some marine animals. This article is protected by copyright. All rights reserved.

Description: Most North American forests are at some stage of post-disturbance regrowth, subject to a changing climate, and exhibit growth and mortality patterns that may not be closely coupled to annual environmental conditions. Distinguishing the possibly interacting effects of these processes is necessary to put short-term studies in a longer-term context, and particularly important for the carbon-dense, fire-prone boreal forest. The goals of this study were to combine dendrochronological sampling, inventory records, and machine-learning algorithms to understand how tree growth and death have changed at one highly studied site (Northern Old Black Spruce, NOBS) in the central Canadian boreal forest. Over the 1999-2012 inventory period, mean tree diameter increased even as stand density and basal area declined significantly. Tree mortality averaged 1.4±0.6% yr −1 , with most mortality occurring in medium-sized trees; new recruitment was minimal. There have been at least two, and probably three, significant influxes of new trees since stand initiation, but none in recent decades. A combined tree ring chronology constructed from sampling in 2001, 2004, and 2012 showed several periods of extreme growth depression, with increased mortality lagging depressed growth by ~5 years. Higher minimum and maximum air temperatures exerted a negative influence on tree growth, while precipitation and climate moisture index had a positive effect; both current- and previous-year data exerted significant effects. Models based on these variables explained 23-44% of the ring-width variability. We suggest that past climate extremes led to significant mortality still visible in the current forest structure, with decadal dynamics superimposed on slower patterns of fire and succession. These results have significant implications for our understanding of previous work at NOBS, the carbon sequestration capability of old-growth stands in a disturbance-prone landscape, and the sustainable management of regional forests in a changing climate. This article is protected by copyright. All rights reserved.

Description: The 20th century was a pivotal period at high northern latitudes as it marked the onset of a rapid climatic warming brought on by major anthropogenic changes in global atmospheric composition. In parallel, Arctic sea ice extent has been decreasing over the period of available satellite data record. Here we document how these changes influenced vegetation productivity in adjacent eastern boreal North America. To do this, we used normalized difference vegetation index (NDVI) data, model simulations of net primary productivity (NPP), and tree-ring width measurements covering the last 300 years. Climatic and proxy-climatic datasets were used to explore the relationships between vegetation productivity and Arctic sea ice concentration and extent, and temperatures. Results indicate that an unusually large amount of black spruce ( Picea mariana ) trees entered into a period of growth decline during the late 20th century (68% of sampled trees; n = 724 cross-sections of age 〉 70 years). This finding is coherent with evidence encoded in NDVI and simulated NPP data. Analyses of climatic and vegetation productivity relationships indicate that the influence of recent climatic changes in the studied forests has been via the enhanced moisture stress (i.e. greater water demands) and autotrophic respiration amplified by the declining sea ice concentration in the Hudson Bay and Hudson Strait. The recent decline strongly contrasts with other growth reduction events that occurred during the 19 th century, which were associated with cooling and high sea ice severity. The recent decline of vegetation productivity is the first one to occur under circumstances related to excess heat in a 300-year period, and further culminates with an intensifying wildfire regime in the region. Our results concur with observations from other forest ecosystems about intensifying temperature-driven drought stress and tree mortality with ongoing climatic changes. This article is protected by copyright. All rights reserved.

Description: Successful species interactions require that both partners share a similar cue. For many species, spring warming acts as a shared signal to synchronize mutualist behaviors. Spring flowering plants and the ants that disperse their seeds respond to warming temperatures so that ants forage when plants drop seeds. However, where warm-adapted ants replace cold-adapted ants, changes in this timing might leave early seeds stranded without a disperser. We investigate plant seed dispersal south and north of a distinct boundary between warm- and cold-adapted ants to determine if changes in the ant species influence local plant dispersal. The warm-adapted ants forage much later than the cold-adapted ants, and so we first assess natural populations of early and late blooming plants. We then transplant these plants south and north of the ant boundary to test whether distinct ant climate requirements disrupt the ant-plant mutualism. Whereas the early blooming plant's inability to synchronize with the warm-adapted ant leaves its populations clumped and patchy and its seedlings clustered around the parents in natural populations, when transplanted into the range of the cold-adapted ant, effective seed dispersal recovers. In contrast, the mutualism persists for the later blooming plant regardless of location because it sets seed later in spring when both warm- and cold-adapted ant species forage, resulting in effective seed dispersal. These results indicate that the climate response of species interactions, not just the species themselves, is integral in understanding ecological responses to a changing climate. Data linking phenological synchrony and dispersal are rare, and these results suggest a viable mechanism by which a species’ range is limited more by biotic than abiotic interactions – despite the general assumption that biotic influences are buried within larger climate drivers. These results show that biotic partner can be as fundamental a niche requirement as abiotic resources. This article is protected by copyright. All rights reserved.

Description: Increasing ocean temperatures and strengthening boundary currents have caused the poleward migration of many marine species. Cubozoan jellyfish known to cause Irukandji syndrome have historically been confined to tropical waters but may be expanding into sub-tropical regions. Here we examine the interactive effects of warming and acidification on the population dynamics of polyps of an Irukandji jellyfish, Alatina nr mordens, and the formation of statoliths in newly metamorphosed medusae, to determine if this jellyfish could tolerate future conditions predicted for southeast Queensland (SEQ), Australia. Two experiments, examining the orthogonal factors of temperature and pH were undertaken. Experiment 1 mimicked the current, ca. 2050 and ca. 2100 summer temperature and pH conditions predicted for SEQ using A1F1 scenarios (temperature: 25, 27, 29°C; pH: 7.9, 7.8, 7.6) and Experiment 2 mimicked current and future winter conditions (18 and 22°C, pH 7.9, 7.8, 7.6). All polyps in Experiment 1 survived and budded. Fewer polyps budded in the lower pH treatments but patterns varied slightly among temperature treatments. Statoliths at pH 7.6 were 24% narrower than those at pH 7.8 and 7.9. Most polyps survived the winter conditions mimicked by Experiment 2 but only polyps in the 22°C, pH 7.9 treatment increased significantly. The current absence of A . nr mordens medusae in SEQ, despite the polyps’ ability to tolerate the current temperature and pH conditions, suggests that ecological, rather than abiotic factors currently limit their distribution. Observations that budding was lower under low pH treatments suggest that rates of asexual reproduction will likely be much slower in the future. We consider that A . nr mordens polyps are likely to tolerate future conditions but are unlikely to thrive in the long term. However, if polyps can overcome potential ecological boundaries and acidification proceeds slowly A . nr mordens could expand polewards in the short-term. This article is protected by copyright. All rights reserved.

Description: Soil CO 2 efflux ( F soil ) is the largest source of carbon from forests and reflects primary productivity as well as how carbon is allocated within forest ecosystems. Through early stages of stand development, both elevated [CO 2 ] and availability of soil nitrogen (N; sum of mineralization, deposition, and fixation) have been shown to increase gross primary productivity, but the long-term effects of these factors on F soil are less clear. Expanding on previous studies at the Duke Free Air CO 2 Enrichment (FACE) site, we quantified the effects of elevated [CO 2 ] and N fertilization on F soil using daily measurements from automated chambers over 10 years. Consistent with previous results, compared to ambient-unfertilized plots, annual F soil increased under elevated [CO 2 ] (~17%) and decreased with N (~21%). N fertilization under elevated [CO 2 ] reduced F soil to values similar to untreated plots. Over the study period, base respiration rates increased with leaf productivity but declined after productivity saturated. Despite treatment-induced differences in aboveground biomass, soil temperature and water content were similar among treatments. Inter-annually, low soil water content decreased annual F soil from potential values – estimated based on temperature alone assuming non-limiting soil water content – by ~0.7% per 1.0% reduction in relative extractable water. This effect was only slightly ameliorated by elevated [CO 2 ]. Variability of soil N availability among plots accounted for the spatial variability of F soil , showing a decrease of ~114 g C m -2 y -1 per 1 g m -2 increase in soil N availability, with consistently higher F soil in elevated [CO 2 ] plots ~127 g C per 100 ppm [CO 2 ] over the +200 ppm enrichment. Altogether, reflecting increased belowground carbon partitioning in response to greater plant nutritional needs, the effects of elevated [CO 2 ] and N fertilization on F soil in this stand are sustained beyond the early stages of stand development and through stabilization of annual foliage production. This article is protected by copyright. All rights reserved.

Description: To meet growing global food demand with limited land and reduced environmental impact, agricultural greenhouse gas (GHG) emissions are increasingly evaluated with respect to crop productivity, i.e. on a yield-scaled as opposed to area basis. Here, we compiled available field data on CH 4 and N 2 O emissions from rice production systems to test the hypothesis that in response to fertilizer nitrogen (N) addition, yield-scaled global warming potential (GWP) will be minimized at N rates that maximize yields. Within each study, yield N surplus was calculated to estimate deficit or excess N application rates with respect to the optimal N rate (defined as the N rate at which maximum yield was achieved). Relationships between yield N surplus and GHG emissions were assessed using linear and nonlinear mixed-effects models. Results indicate that yields increased moving from deficit to optimal N rates. At N rates contributing to a yield N surplus, N 2 O and yield-scaled N 2 O emissions increased exponentially. In contrast, CH 4 emissions were not impacted by N inputs. Accordingly, yield-scaled CH 4 emissions decreased with N addition. Overall, yield-scaled GWP was minimized at optimal N rates, decreasing by 21% compared to treatments without N addition. These results are unique compared to aerobic cropping systems in which N 2 O emissions are the primary contributor to GWP, meaning yield-scaled GWP may not necessarily decrease for aerobic crops when yields are optimized by N fertilizer additions. Balancing gains in agricultural productivity with climate change concerns, this work supports the concept that high rice yields can be achieved with minimal yield-scaled GWP through optimal N application rates. Moreover, additional improvements in N use efficiency may further reduce yield-scaled GWP, thereby strengthening the economic and environmental sustainability of rice systems. This article is protected by copyright. All rights reserved.

Description: Climate change is projected to push the limits of cropping systems and has the potential to disrupt the agricultural sector from local to global scales. This article introduces the Coordinated Climate-Crop Modeling Project (C3MP), an initiative of the Agricultural Model Intercomparison and Improvement Project (AgMIP) to engage a global network of crop modelers to explore the impacts of climate change via an investigation of crop responses to changes in carbon dioxide concentration ([CO 2 ]), temperature, and water. As a demonstration of the C3MP protocols and enabled analyses, we apply the Decision Support System for Agrotechnology Transfer (DSSAT) CROPGRO-Peanut crop model for Henry County, Alabama, to evaluate responses to the range of plausible [CO 2 ], temperature changes, and precipitation changes projected by climate models out to the end of the 21 st century. These sensitivity tests are used to derive crop model emulators that estimate changes in mean yield and the coefficient of variation for seasonal yields across a broad range of climate conditions, reproducing mean yields from sensitivity test simulations with deviations of ~2% for rainfed conditions. We apply these statistical emulators to investigate how peanuts respond to projections from various global climate models, time periods, and emissions scenarios, finding a robust projection of modest (〈10%) median yield losses in the middle of the 21 st century accelerating to more severe (〉20%) losses and larger uncertainty at the end of the century under the more severe representative concentration pathway 8.5. This projection is not substantially altered by the selection of the AgMERRA global gridded climate dataset rather than the local historical observations, differences between the Third and Fifth Coupled Model Intercomparison Project (CMIP3 and CMIP5), or the use of the delta method of climate impacts analysis rather than the C3MP impacts response surface and emulator approach. This article is protected by copyright. All rights reserved.

Description: Large-scale, long-term FACE (Free Air CO 2 -enrichment) experiments indicate that increases in atmospheric CO 2 concentrations will influence forest C cycling in unpredictable ways. It has been recently suggested that responses of mycorrhizal fungi could determine whether forest NPP (net primary production) is increased by elevated CO 2 over long time periods and if forests soils will function as sources or sinks of C in the future. We studied the dynamic responses of ectomycorrhizae to N fertilization and atmospheric CO 2 -enrichment at the Duke FACE experiment using minirhizotrons over a six year period (2005-2010). Stimulation of mycorrhizal production by elevated CO 2 was observed during only one (2007) of six years. This increased the standing crop of mycorrhizal tips during 2007 and 2008; during 2008, significantly higher mortality returned standing crop to ambient levels for the remainder of the experiment. It is therefore unlikely that increased production of mycorrhizal tips can explain the lack of progressive nitrogen limitations and associated increases in N uptake observed in CO 2 -enriched plots at this site. Fertilization generally decreased tree reliance on mycorrhizae as tip production declined with the addition of nitrogen as has been shown in many other studies. Annual NPP of mycorrhizal tips was greatest during years with warm January temperatures and during years with cool spring temperatures. A 2° C increase in average late spring temperatures (May and June) decreased annual production of mycorrhizal root tip length by 50%. This has important implications for ecosystem function in a warmer world in addition to potential for forest soils to sequester atmospheric C. This article is protected by copyright. All rights reserved.

Description: Permafrost thaw in the Arctic driven by climate change is mobilizing ancient terrigenous organic carbon (OC) into fluvial networks. Understanding the controls on metabolism of this OC is imperative for assessing its role with respect to climate feedbacks. In this study we examined the effect of inorganic nutrient supply and dissolved organic matter (DOM) composition on aquatic extracellular enzyme activities (EEAs) in waters draining the Kolyma River Basin (Siberia), including permafrost derived OC. Reducing the phenolic content of the DOM pool resulted in dramatic increases in hydrolase EEAs (e.g. phosphatase activity increased 〉 28 fold) supporting the idea that high concentrations of polyphenolic compounds in DOM (e.g. plant structural tissues) inhibit enzyme synthesis or activity, limiting OC degradation. EEAs were significantly more responsive to inorganic nutrient additions only after phenolic inhibition was experimentally removed. In controlled mixtures of modern OC and thawed permafrost endmember OC sources, respiration rates per unit dissolved OC were 1.3 – 1.6 times higher in waters containing ancient carbon, suggesting that permafrost derived OC was more available for microbial mineralization. In addition, waters containing ancient permafrost derived OC supported elevated phosphatase and glucosidase activities. Based on these combined results, we propose that both composition and nutrient availability regulates DOM metabolism in Arctic aquatic ecosystems. Our empirical findings are incorporated into a mechanistic conceptual model highlighting two key enzymatic processes in the mineralization of riverine OM: 1) the role of phenol oxidase activity in reducing inhibitory phenolic compounds; and 2) the role of phosphatase in mobilizing organic P. Permafrost derived DOM degradation was less constrained by this initial “phenolic-OM” inhibition; thus, informing reports of high biological availability of ancient, permafrost derived DOM with clear ramifications for its metabolism in fluvial networks and feedbacks to climate. This article is protected by copyright. All rights reserved.

Description: Ecosystem functioning is simultaneously affected by changes in community composition and environmental change such as increasing atmospheric carbon dioxide (CO 2 ) and subsequent ocean acidification. However, it largely remains uncertain how the effects of these factors compare to each other. Addressing this question, we experimentally tested the hypothesis that initial community composition and elevated CO 2 are equally important to the regulation of phytoplankton biomass. We full-factorially exposed three compositionally different marine phytoplankton communities to two different CO 2 levels and examined the effects and relative importance (ω 2 ) of the two factors and their interaction on phytoplankton biomass at bloom peak. The results showed that initial community composition had a significantly greater impact than elevated CO 2 on phytoplankton biomass, which varied largely among communities. We suggest that the different initial ratios between cyanobacteria, diatoms, and dinoflagellates might be the key for the varying competitive and thus functional outcome among communities. Furthermore, the results showed that depending on initial community composition elevated CO 2 selected for larger sized diatoms, which led to increased total phytoplankton biomass. Our study highlights the relevance of initial community composition, which strongly drives the functional outcome, when assessing impacts of climate change on ecosystem functioning. In particular, the increase in phytoplankton biomass driven by the gain of larger sized diatoms in response to elevated CO 2 potentially has strong implications for nutrient cycling and carbon export in future oceans. This article is protected by copyright. All rights reserved.

Description: Less than half of anthropogenic carbon emissions are accumulating in the atmosphere, due to large net fluxes into both the oceans and the land (Le Queré et al., 2012). The land sink in particular has increased markedly, doubling in strength since the 1960's, to reach 26 petagrams of carbon in the latest decade. However, the location and drivers of this large terrestrial sink are still relatively poorly constrained by atmospheric measurements (Ciais et al. 2013). Pan et al. (2011) recently utilised 〉1 million forest inventory plots to provide summaries of forest carbon stocks, and the first global bottom-up estimates of carbon fluxes for the world's forest biomes for the period 1990-2007. One key result was that almost all the residual global terrestrial carbon sink (i.e. carbon uptake after accounting for land use change), some 2.4 ± 0.4 Pg of carbon per year, is located in the world's established forests (Pan et al., 2011). The sink is distributed worldwide, with globally significant net fluxes into boreal and temperate forests, and a large sink in intact tropical forest, albeit with large uncertainty. Furthermore, Pan et al. (2011) showed that this tropical intact forest sink - may have faded from an estimated annual 1.3 ± 0.4 Pg C in the 1990's to 1.0 ± 0.5 Pg C for 2000-2007. The tropical intact sink is offset by a net land-use emission (1.5 Pg C yr −1 [1990-1999]) declining to 1.1 Pg C yr −1 [2000-2007]), and as a consequence aircraft measurements and inverse modelling studies indicate the tropics to be close to neutral in terms of net carbon fluxes (reviewed by Ciais et al. 2013). While the intact tropical forest sink values represent updates from similar values published previously (e.g., Lewis et al., 2009a), the fact that almost the entire residual terrestrial carbon sink is accounted for by the forests of the world was a notable discovery. Evidence from the ground now points to established forests being a net sink across almost every major forest region, including all extra-tropical forest regions analysed. This article is protected by copyright. All rights reserved.

Description: Forecasting how global warming will affect onset of the growing season is essential for predicting terrestrial productivity, but suffers from conflicting evidence. We show that accurate estimates require ways to connect discrete observations of changing tree status (e.g., pre- vs. post-budbreak) with continuous responses to fluctuating temperatures. By coherently synthesizing discrete observations with continuous responses to temperature variation, we accurately quantify how increasing temperature variation accelerates onset of growth. Application to warming experiments at two latitudes demonstrates that maximum responses to warming are concentrated in late winter, weeks ahead of the main budbreak period. Given that warming will not occur uniformly over the year, knowledge of when temperature variation has the most impact can guide prediction. Responses are large and heterogeneous, yet predictable. The approach has immediate application to forecasting effects of warming on growing season length, requiring only information that is readily available from weather stations and generated in climate models. This article is protected by copyright. All rights reserved.

Description: Climate change is expected to cause geographic redistributions of species. To the extent that species within assemblages have different niche requirements, assemblages may no longer remain intact and dis- and reassemble at current or new geographic locations. We explored how climate change projected by 2100 may transform the world's avian assemblages (characterized at a 110 km spatial grain) by modelling environmental niche-based changes to their dietary guild structure under 0 km, 500 km, and 2000 km dispersal distances. We examined guild structure changes at coarse (primary, high-level, and mixed consumers) and fine (frugivores, nectarivores, insectivores, herbivores, granivores, scavengers, omnivores, and carnivores) ecological resolutions to determine whether or not geographic co-occurrence patterns among guilds were associated with the magnitude to which guilds are functionally resolved. Dietary guilds vary considerably in their global geographic prevalence, and under broad-scale niche-based redistribution of species, these are projected to change very heterogeneously. A non-dispersal assumption results in the smallest projected changes to guild assemblages, but with significant losses for some regions and guilds, such as South American insectivores. Longer dispersal distances are projected to cause greater degrees of disassembly, and lead to greater homogenization of guild composition, especially in northern Asia and Africa. This arises because projected range gains and losses result in geographically heterogeneous patterns of guild compensation. Projected decreases especially of primary and mixed consumers most often are compensated by increases in high-level consumers, with increasing uncertainty about these outcomes as dispersal distance and degree of guild functional resolution increases. Further exploration into the consequences of these significant broad-scale ecological functional changes at the community or ecosystem level should be increasingly on the agenda for conservation science. This article is protected by copyright. All rights reserved.

Description: Climate change scenarios predict increases in the frequency and duration of ENSO-related droughts for parts of South-East Asia until the end of this century exposing the remaining rainforests to increasing drought risk. A pan-tropical review of recorded drought-related tree mortalities in more than 100 monitoring plots before, during and after drought events suggested a higher drought-vulnerability of trees in South-East Asian than in Amazonian forests. Here, we present the results of a replicated (n=3 plots) throughfall exclusion experiment in a perhumid tropical rainforest in Sulawesi, Indonesia. In this first large-scale roof experiment outside semi-humid eastern Amazonia, 60% of the throughfall was displaced during the first 8 months and 80% during the subsequent 17 months, exposing the forest to severe soil desiccation for about 17 months. In the experiment's second year, wood production decreased on average by 40% with largely different responses of the tree families (ranging from -100 to +100% change). Most sensitive were trees with high radial growth rates under moist conditions. In contrast, tree height was only a secondary factor and wood specific gravity had no influence on growth sensitivity. Fine root biomass was reduced by 35% after 25 months of soil desiccation while fine root necromass increased by 250% indicating elevated fine root mortality. Cumulative aboveground litter production was not significantly reduced in this period. The trees from this Indonesian perhumid rainforest revealed similar responses of wood and litter production and root dynamics as those in two semi-humid Amazonian forests subjected to experimental drought. We conclude that trees from paleo- or neotropical forests growing in semi-humid or perhumid climates may not differ systematically in their growth sensitivity and vitality under sub-lethal drought stress. Drought vulnerability may depend more on stem cambial activity in moist periods than on tree height or wood specific gravity. This article is protected by copyright. All rights reserved.

Description: Soil microbial communities in Chihuahuan Desert grasslands generally experience highly variable spatiotemporal rainfall patterns. Changes in precipitation regimes can affect belowground ecosystem processes such as decomposition and nutrient cycling by altering soil microbial community structure and function. The objective of this study was to determine if increased seasonal precipitation frequency and magnitude over a seven-year period would generate a persistent shift in microbial community characteristics and soil nutrient availability. We supplemented natural rainfall with large events (one/winter and three/summer) to simulate increased precipitation based on climate model predictions for this region. We observed a two year delay in microbial responses to supplemental precipitation treatments. In Years 3-5, higher microbial biomass, arbuscular mycorrhizae abundance, and soil enzyme C and P acquisition activities were observed in the supplemental water plots even during extended drought periods. In Years 5-7, available soil P was consistently lower in the watered plots compared to control plots. Shifts in soil P corresponded to higher fungal abundances, microbial C utilization activity, and soil pH. This study demonstrated that 25% shifts in seasonal rainfall can significantly influence soil microbial and nutrient properties, which in turn may have long-term effects on nutrient cycling and plant P uptake in this desert grassland. This article is protected by copyright. All rights reserved.

Description: It is proposed that increases in anthropogenic reactive nitrogen (N r )-deposition may cause temperate and boreal forests to sequester a globally significant quantity of carbon (C); however, long-term data from boreal forests describing how C sequestration responds to realistic levels of chronic N r -deposition are scarce. Using a long term (14-year) stand scale (0.1 ha) N-addition experiment (three levels: 0, 12.5, and 50 kg N ha −1 yr −1 ) in the boreal zone of northern Sweden, we evaluated how chronic N additions altered N uptake and biomass of understory communities, and whether changes in understory communities explained N uptake and C sequestration by trees. We hypothesized that understory communities (i.e. mosses and shrubs) serve as important sinks for low-level N additions, with the strength of these sinks weakening as chronic N addition rates increase, due to shifts in species composition. We further hypothesized that trees would exhibit non-linear increases in N acquisition, and subsequent C sequestration as N addition rates increased, due to a weakening understory N sink. Our data showed that understory biomass was reduced by 50% in response to the high N addition treatment, mainly due to reduced moss biomass. A 15 N labelling experiment showed that feather mosses acquired the largest fraction of applied label, with this fraction decreasing as the chronic N addition level increased. Contrary to our hypothesis, the proportion of label taken up by trees was equal (~8%) across all three N addition treatments. The relationship between N addition and C sequestration in all vegetation pools combined was linear, and had a slope of 16 kg C kg −1 N. While canopy retention of N r deposition may cause C sequestration rates to be slightly different than this estimate, our data suggests that a minor quantity of annual anthropogenic CO 2 emissions are sequestered into boreal forests as a result of N r deposition. This article is protected by copyright. All rights reserved.

Description: Urbanization is one of the most extensive and ecologically significant changes happening to terrestrial environments, as it strongly affects the distribution of biodiversity. It is well established that native species richness is reduced in urban and suburban areas, but the species traits that predict tolerance to urbanization are yet little understood. In birds, one of the most studied groups in this respect, evidence is appearing that acoustic traits influence urban living, but it is unknown how this compares to the effects of more obvious ecological traits that facilitate urban living. Therefore, it remains unclear whether acoustic communication is an important predictor of urban tolerance among species. Here, with a comparative study across 140 European and North American passerines, I show that high song frequency, which is less masked by the low-frequency anthropogenic noise, is associated with urban tolerance, with an effect size over half that of the most important ecological trait studied: off-ground nesting. Other nesting and foraging traits accepted to facilitate urban living did not differ for species occurring in urban environments. Thus, the contribution of acoustic traits for passerine urban tolerance approximates that of more obvious ecological traits. Nonetheless, effect sizes of the biological predictors of urban tolerance were low and the phylogenetic signal for urban tolerance was null, both of which suggest that factors other than phenotypic traits have major effects on urban tolerance. A simple possibility is exposure to urbanization, since there was a higher proportion of urban-tolerant species in Europe, which is more urbanized than North America. This article is protected by copyright. All rights reserved.

Description: [1]  Determining the scale-length, magnitude, and distribution of heterogeneity in the lowermost mantle is crucial to understanding whole mantle dynamics, and yet it remains a much debated and ongoing challenge in geophysics. Common shortcomings of current seismically-derived lowermost mantle models are incomplete raypath coverage, arbitrary model parameterization, inaccurate uncertainty estimates, and an ad hoc definition of the misfit function in the optimization framework. In response, we present a new approach to global tomography. Apart from improving the existing raypath coverage using only high quality cross-correlated waveforms, the problem is addressed within a Bayesian framework where explicit regularization of model parameters is notrequired. We obtain high resolution images, complete with uncertainty estimates, of the lowermost mantle P-wave velocity structure using a hand-picked dataset of PKPab-df, PKPbc-df, and PcP-P differential traveltimes. Most importantly, our results demonstrate that the root mean square of the P-wave velocity variations in the lowermost mantle is approximately 0.87%, which is three times larger than previous global-scale estimates.

Description: There is evidence that climate change induced tree mortalities in boreal and temperate forests and increased forest turnover rates (both mortality and recruitment rates) in Amazon forests. However, no study has examined China's tropical and subtropical evergreen broadleaved forests (TEBF) that cover 〉26% of China's terrestrial land . The sustainability of this biome is vital to the maintenance of local ecosystem services (e.g., carbon sequestration, biodiversity conservation, climatic regulation etc.), many of which may influence patterns of atmospheric circulation and composition at regional to global scales. Here we analyze time-series data collected from thirteen permanent plots within China's unmanaged TEBF to study whether and how this biome has changed over recent decades. We find that the numbers of individuals and species for shrub and small tree have increased since 1978, whereas the numbers of individuals and species for tree have decreased over this same time period. The shift in species composition is accompanied by a decrease in the mean DBH (diameter at breast height) for all individuals combined. China's TEBF may thereby be transitioning from cohorts of fewer and larger individuals to ones of more and smaller individuals, which shows a unique change pattern differing from the documented. Regional-scale drying is likely responsible for the biome's reorganization. This biome-wide reconstitution would deeply impact the regimes of carbon sequestration and biodiversity conservation and have implications for the sustainability of economic development in the area. This article is protected by copyright. All rights reserved.

Description: Freshwater ecosystems provide vital resources for humans and support high levels of biodiversity, yet are severely threatened throughout the world. The expansion of human land uses, such as urban and crop cover, typically degrades water quality and reduces freshwater biodiversity, thereby jeopardizing both biodiversity and ecosystem services. Identifying and mitigating future threats to freshwater ecosystems requires forecasting where land use changes are most likely. Our goal was to evaluate the potential consequences of future land use on freshwater ecosystems in the coterminous United States by comparing alternative scenarios of land use change (2001-2051) with current patterns of freshwater biodiversity and water-quality risk. Using an econometric model, each of our land use scenarios projected greater changes in watersheds of the eastern half of the country, where freshwater ecosystems already experience higher stress from human activities. Future urban expansion emerged as a major threat in regions with high freshwater biodiversity (e.g., the Southeast) or severe water-quality problems (e.g., the Midwest). Our scenarios reflecting environmentally-oriented policies had some positive effects. Subsidizing afforestation for carbon sequestration reduced crop cover and increased natural vegetation in areas that are currently stressed by low water quality, while discouraging urban sprawl diminished urban expansion in areas of high biodiversity. On the other hand, we found that increases in crop commodity prices could lead to increased agricultural threats in areas of high freshwater biodiversity. Our analyses illustrate the potential for policy changes and market factors to influence future land use trends in certain regions of the country, with important consequences for freshwater ecosystems. Successful conservation of aquatic biodiversity and ecosystem services in the U.S. into the future will require attending to the potential threats and opportunities arising from policies and market changes affecting land use. This article is protected by copyright. All rights reserved.

Description: [1]  The Dzhungarian strike-slip fault of Central Asia is one of a series of long, NW-SE right-lateral strike-slip faults that are characteristic of the northern Tien Shan region, and extends over 300 km from the high mountains into the Kazakh Platform. Our field-based and satellite observations reveal that the Dzhungarian fault can be characterised by three 100 km long sections based on variation in strike direction. Through morphological analysis of offset streams and alluvial fans, and through OSLdating, we find that the Dzhungarian fault has a minimum average late Quaternary slip rate of 2.2 ± 0.8 mm/yr and accommodates N-S shortening related to the India-Eurasia collision. This shortening may also be partly accommodated by counter-clockwise rotation about a vertical axis. Evidence for a possible paleo-earthquake rupture indicates that earthquakes up to at least Mw 7 can be associated with just the partitioned component of reverse slip on segments of the central section of the fault up to 30 km long. An event rupturing longer sections of the Dzhungarian fault has the potential to generate greater magnitude earthquakes ( Mw 8), however long time periods (e.g. thousands of years) are expected in order to accumulate enough strain to generate such earthquakes.

Description: Animals living in tropical regions may be at increased risk from climate change because current temperatures at these locations already approach critical physiological thresholds. Relatively small temperature increases could cause animals to exceed these thresholds more often, resulting in substantial fitness costs or even death. Oviparous species could be especially vulnerable because the maximum thermal tolerances of incubating embryos is often lower than adult counterparts, and in many species mothers abandon the eggs after oviposition, rendering them immobile and thus unable to avoid extreme temperatures. As a consequence, the effects of climate change might become evident earlier and be more devastating for hatchling production in the tropics. Loggerhead sea turtles ( Caretta caretta ) have the widest nesting range of any living reptile, spanning temperate to tropical latitudes in both hemispheres. Currently, loggerhead sea turtle populations in the tropics produce nearly 30% fewer hatchlings per nest than temperate populations. Strong correlations between empirical hatching success and habitat quality allowed global predictions of the spatiotemporal impacts of climate change on this fitness trait. Under climate change, many sea turtle populations nesting in tropical environments are predicted to experience severe reductions in hatchling production, whereas hatching success in many temperate populations could remain unchanged or even increase with rising temperatures. Some populations could show very complex responses to climate change, with higher relative hatchling production as temperatures begin to increase, followed by declines as critical physiological thresholds are exceeded more frequently. Predicting when, where, and how climate change could impact the reproductive output of local populations is crucial for anticipating how a warming world will influence population size, growth, and stability. This article is protected by copyright. All rights reserved.

Description: It has recently been found that the frequency distribution of remotely sensed tree cover in the tropics has three distinct modes, which seem to correspond to forest, savanna and treeless states. This pattern has been suggested to imply that these states represent alternative attractors, and that the response of these systems to climate change would be characterized by critical transitions and hysteresis. Here, we show how this inference is contingent upon mechanisms at play. We present a simple dynamical model that can generate three alternative tree cover states (forest, savanna and a treeless state), based on known mechanisms, and use this model to simulate patterns of tree cover under different scenarios. We use these synthetic data to show that the hysteresis inferred from remotely sensed tree cover patterns will be inflated by spatial heterogeneity of environmental conditions. On the other hand, we show that the hysteresis inferred from satellite data may actually underestimate real hysteresis in response to climate change if there exists a positive feedback between regional tree cover and precipitation. Our results also indicate that such positive feedback between vegetation and climate should cause direct shifts between forest and a treeless state (rather than through an intermediate savanna-state) to become more likely. Lastly, we show how directionality of historical change in conditions may bias the observed relationship between tree cover and environmental conditions. This article is protected by copyright. All rights reserved.

Description: [1]  We performed shock compression experiments on preheated forsterite liquid (Mg 2 SiO 4 ) at an initial temperature of 2273 K and have revised the equation of state (EOS) that was previously determined by shock melting of initially solid Mg 2 SiO 4 (300 K). The linear Hugoniot, U S  = 2.674 ± 0.188 + 1.64 ± 0.06 u p km/s, constrains the bulk sound speed within a temperature and composition space as yet unexplored by 1-bar ultrasonic experiments. We have also revised the EOS for enstatite liquid (MgSiO 3 ) to exclude experiments that may have been only partially melted upon shock compression and also the EOS for anorthite liquid, which now excludes potentially un-relaxed experiments at low pressure. The revised fits and the previously determined EOS of fayalite and diopside were used to produce isentropes in the multicomponent CaO-MgO-Al 2 O 3 -SiO 2 -FeO system at elevated temperatures and pressures. Our results are similar to those previously presented for peridotite and simplified “chondrite” liquids such that regardless of where crystallization first occurs, the liquidus solid sinks upon formation. This process is not conducive to the formation of a basal magma ocean. We also examined the chemical and physical plausibility of the partial melt hypothesis to explain the occurrence and characteristics of ultralow velocity zones. We determined that the ambient mantle cannot produce an equilibrium partial melt and residue that is sufficiently dense to be a ULVZ mush. The partial melt would need to be segregated from its equilibrium residue and combined with a denser solid component to achieve a sufficiently large aggregate density.

Description: [1]  High resolution sparker and crustal-scale airgun seismic reflection data, coupled with repeat bathymetric surveys, document a region of repeated coseismic uplift on the portion of the Alaska subduction zone that ruptured in 1964. This area defines the western limit of Prince William Sound. Differencing of vintage and modern bathymetric surveys shows that the region of greatest uplift related to the 1964 Great Alaska earthquake was focused along a series of sub-parallel faults beneath Prince William Sound and the adjacent Gulf of Alaska shelf. Bathymetric differencing indicates that 12 m of coseismic uplift occurred along two faults that reached the sea floor as submarine terraces on the Cape Cleare bank southwest of Montague Island. Sparker seismic reflection data provide cumulative Holocene slip estimates as high as 9 mm/yr along a series of splay thrust faults within both the inner wedge and transition zone of the accretionary prism. Crustal seismic data show that these megathrust splay faults root separately into the subduction zone décollement. Splay fault divergence from this megathrust correlates with changes in mid-crustal seismic velocity and magnetic susceptibility values, best explained by duplexing of the subducted Yakutat terrane rocks above Pacific plate rocks along the trailing edge of the Yakutat terrane. Although each splay fault is capable of independent motion, we conclude that the identified splay faults rupture in a similar pattern during successive megathrust earthquakes and that the region of greatest seismic coupling has remained consistent throughout the Holocene.

Description: [1]  We perform a time-lapse analysis of Rayleigh and Love wave anisotropy above an underground gas storage facility in the Paris Basin. The data were acquired with a three-component seismic array deployed during several days in April and November 2010. Phase velocity and back azimuth of Rayleigh and Love waves are measured in the frequency range 0.2-1.1 Hz using a three-component beamforming algorithm. In both snapshots, higher surface wave modes start dominating the signal above 0.4 Hz with a concurrent increase in back azimuth ranges. We fit anisotropy parameters to the array detections above 0.4 Hz using a bootstrap approach which also provides estimation uncertainty and enables significance testing. The isotropic phase velocity dispersion for Love and Rayleigh waves match for both snapshots. We also observe a stable fast direction of NNW-SSE for Love and Rayleigh waves which is aligned with the preferred orientation of known shallow (〈300 m) and deeper (~1000 m) fault systems in the area, as well as the maximum horizontal stress orientation. At lower frequencies corresponding to deeper parts of the basin, the anisotropic parameters exhibit higher magnitude in the November data. This may perhaps be caused by the higher pore-pressure changes in the gas reservoir in that depth range.

Description: [1]  Eruptive activity at the summit of Kilauea Volcano, Hawaii beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with a 11-station broadband network to model the source mechanism of signals accompanying two large rockfalls on August 29, 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal-to-noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1-1000 s. The VLP signals associated with the rockfalls originate in a source region ~1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80 ∘ to the north, intersecting a north striking crack (another dike) dipping 65 ∘ to the east. Each rockfall is marked by a similar step-like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200 − 4500 m 3 ), length of magma column (120-210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3 – 6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400  −  15, 000 kg s − 1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on August 29, 2012.

Description: The magnitude and direction of phenological shifts from climate warming could be predictably variable across the planet depending upon the nature of physiological controls on phenology, the thermal sensitivity of the developmental processes and global patterns in the climate warming. We tested this with respect to the flight phenology of adult nocturnal moths (3.33 million captures of 334 species) that were sampled at sites in southern and northern Finland during 1993–2012 (with years 2005–2012 treated as an independent model validation data set). We compared eight competing models of physiological controls on flight phenology to each species and found strong support for thermal controls of phenology in 66% of the species generations. Among species with strong thermal control of phenology in both the south and north, the average development rate was higher in northern vs. southern populations at 10 °C, but about the same at 15 and 20 °C. With a 3 °C increase in temperature (approximating A2 scenario of IPPC for 2090–2099 relative to 1980–1999) these species were predicted to advance their phenology on average by 17 (SE ± 0.3) days in the south vs. 13 (±0.4) days in the north. The higher development rates at low temperatures of poleward populations makes them less sensitive to climate warming, which opposes the tendency for stronger phenological advances in the north from greater increases in temperature.

Description: Desert annuals are a critically important component of desert communities and may be particularly responsive to increasing atmospheric [CO 2 ] because of their high potential growth rates and flexible phenology. During the ten-year life of the Nevada Desert FACE (Free-air CO 2 enrichment) Facility, we evaluated the productivity, reproductive allocation, and community structure of annuals in response to long-term elevated [CO 2 ] exposure. The dominant forb and grass species exhibited accelerated phenology, increased size, and higher reproduction at elevated [CO 2 ] in a wet El Niño year near the beginning of the experiment. However, a multi-year dry cycle resulted in no increases in productivity or reproductive allocation for the remainder of the experiment. At the community level, early indications of increased dominance of the invasive Bromus rubens at elevated [CO 2 ] gave way to an absence of Bromus in the community during a drought cycle, with a resurgence late in the experiment in response to higher rainfall and a corresponding high density of Bromus in a final soil seed bank analysis, particularly at elevated [CO 2 ]. This long-term experiment resulted in two primary conclusions: (1) elevated [CO 2 ] does not increase productivity of annuals in most years; and (2) relative stimulation of invasive grasses will likely depend on future precipitation, with a wetter climate favoring invasive grasses but currently predicted greater aridity favoring native dicots. This article is protected by copyright. All rights reserved.

Description: [1]  On 5 September 2012, a large thrust earthquake (M w 7.6) ruptured a densely-instrumented seismic gap on the shallow-dipping plate boundary beneath the Nicoya Peninsula, Costa Rica. Ground motion recordings directly above the rupture zone provide a unique opportunity to study the detailed source process of a large shallow megathrust earthquake using very nearby land observations. Hypocenter relocation using local seismic network data indicates that the event initiated with small emergent seismic waves from a hypocenter ~10 km offshore, 13 km deep on the megathrust. A joint finite-fault inversion using high-rate GPS, strong-motion ground velocity recordings, GPS static offsets, and teleseismic P waves reveals that the primary slip zone (slip 〉 1 m) is located beneath the peninsula. The rupture propagated down-dip from the hypocenter with a rupture velocity of ~3.0 km/s. The primary slip zone extends ~70 km along strike and ~30 km along dip, with an average slip of ~2 m. The associated static stress drop is ~3 MPa. The seismic moment is 3.5 x 10 20  Nm, giving M w  = 7.6. The co-seismic large-slip patch directly overlaps an onshore inter-seismic locked region indicated by geodetic observations, and extends down-dip to the intersection with the upper plate Moho. At deeper depths, below the upper plate Moho, seismic tremor and low frequency earthquakes have been observed. Most tremor locates in adjacent areas of the megathrust that have little co-seismic slip; a region of prior slow slip deformation to the southeast also has no significant co-seismic slip or aftershocks. An offshore locked patch indicated by geodetic observations does not appear to have experienced co-seismic slip, and aftershocks do not overlap this region, allowing the potential for a comparable size rupture offshore in the future.

Description: [1]  Seismic shear waves emitted by earthquakes can be modeled as plane (transverse) waves. When entering an anisotropic medium they can be split into two orthogonal components moving at different speeds. This splitting occurs along an axis, the fast polarization, that is determined by geologic conditions. We present here a comprehensive analysis of the Silver and Chan (1991) method, used to obtain shear wave splitting parameters, comprising theoretical derivations and statistical tests of the assumptions used to construct the standard errors. We find discrepancies in the derivations of equations in their article, with the most important being a mistake in how the standard errors are calculated. Our simulations suggest that the degrees of freedom are being overestimated by this method and consequently the standard errors are too small. Using a set of S waveforms from very similar shallow earthquakes on Reunion Island, we perform a statistical analysis on the noise of these replicates and find that the assumption of Gaussian noise does not hold. Further, the properties of background noise differ substantially from the noise obtained from the shear wave splitting analysis. However, we find that the standard errors for the fast polarization are comparable to the spread in the fast polarization parameters. Delay time errors appear to be comparable to delay time estimates once cycle skipping is accounted for. Future work using synthetic seismograms with simulated noise should be conducted to confirm this is the case for earthquakes in general.

Description: [1]  The sliver strike-slip Great Sumatra Fault (GSF) traverses mainland Sumatra from the Sunda Strait in the southeast to Banda Aceh in the northwest, and defines the present day plate boundary between the Sunda Plate in the north and the Burmese Sliver Plate in the south. It has been well studied on mainland Sumatra but poorly north of Banda Aceh in the Andaman Sea. Here we present deep seismic reflection images along the northward extension of the GSF over 700 km until it joins the Andaman Sea Spreading Centre and we interpret these images in the light of earthquake, gravity, and bathymetry data. We find that the GSF has two strands between Banda Aceh and Nicobar Island: a transpression in the south and a deep narrow active rift system in the north dotted with volcanoes in the center, suggesting that the volcanic arc is coincident with rifting. Further north of Nicobar Island, an active strike-slip fault, the Andaman-Nicobar Fault, cuts through a rifted deep basin until its intersection with the Andaman Sea Spreading Centre. The volcanic arc lies just east of the rift basin. The western margin of this basin seems to be a rifted continental margin, tilted westward, and flooring the Andaman-Nicobar forearc basin. The Andaman-Nicobar forearc basin is bounded in the west by backthrusts similar to the West Andaman and Mentawai faults. The cluster of seismicity after the 2004 great Andaman-Sumatra earthquake just north of Nicobar Island coincides with the intersection of two strike-slip fault systems.

Description: [1]  The Japan Tohoku-Oki earthquake (9.0 Mw) of 11 March 2011 has left signatures in the Earth's gravity field that are detectable by data of the GRACE mission. Because ESA's satellite gravity mission GOCE – launched in 2009 – aims at high spatial resolution, its measurements could complement the GRACE information on coseismic gravity changes, although time-variable gravity was not foreseen as goal of the GOCE mission. We modeled the coseismic earthquake geoid signal and converted this signal to vertical gravity gradients at GOCE satellite altitude. We combined the single gradient observations in a novel way reducing the noise level, required to detect the coseismic gravity change, subtracted a global gravity model, and applied tailored outlier detection to the resulting gradient residuals. Furthermore, the measured gradients were along-track filtered using different gradient bandwidths where in the space domain Gaussian smoothing has been applied. One year periods before and after earthquake occurrence have been compared with the modeled gradients. The comparison reveals that the earthquake signal is well above the accuracy of the vertical gravity gradients at orbital height. Moreover, the obtained signal from GOCE shows a 1.3 times higher amplitude compared with the modeled signal. Besides the statistical significance of the obtained signal, it has a high spatial correlation of ~0.7 with the forward modeled signal. We conclude therefore that the coseismic gravity change of the Japan Tohoku-Oki earthquake left a statistically significant signal in the GOCE measured gravity gradients.

Description: [1]  The Ninetyeast Ridge (NER), one of the longest linear volcanic features on the Earth, offers an excellent opportunity of understanding the isostatic response to the interactions of mantle plume with the migrating mid-ocean ridge. Bathymetry, geoid and gravity (ship-borne and satellite) data along 72 closely spaced transects and 17 overlapping grids on the NER are analyzed and modeled to determine the effective elastic thickness ( Te ) beneath the entire ridge. The results of 2-D and 3-D flexural modeling of the NER show large spatial variations in Te values ranging from 4 to 35 km, suggesting that the ridge was compensated along its length by different isostatic mechanisms. The southern (south of 22°S latitude) and northern (north of 2°N latitude) parts of the NER have Te values of 〉10 and 〉23 km, respectively, revealing that the southern part was emplaced on a lithosphere of intermediate strength possibly on flank of the Indian plate, whereas the northern part was emplaced in an intraplate setting. In contrast, in the central part of the NER (between latitudes 22°S and 2°N), highly variable Te values (4–22 km) are estimated. The scattered Te values in the central NER suggest that this part may have evolved due to the occurrence of more frequent ridge jumps caused by the interaction of Kerguelen hot spot with rapid northward migration of the Wharton spreading ridge. Residual Mantle Bouguer Anomaly (RMBA) map of the NER and adjacent basins reveals that the entire length of the NER is associated with a significant negative anomaly up to 200 mGal, indicating the presence of thickened crust or less dense mantle beneath the ridge. 3-D crustal thickness map of the NER, generated by inversion of the RMBA data, shows a thick crust ranging from 15 to 19 km. The present study clearly shows that NER possesses a highly segmented isostatic pattern with the occurrence of sub-crustal underplating or sub-surface loading.

Description: [1]  We conducted deep-sea magnetic measurements using autonomous underwater vehicles in the Bayonnaise knoll caldera, the Izu-Ogasawara island arc, which hosts the large Hakurei hydrothermal field. We improved the conventional correction method applied for removing the effect of vehicle magnetization, thus greatly enhancing the precision of the resulting vector anomalies. The magnetization distribution obtained from the vector anomaly data shows an ∼ 2-km-wide belt of high magnetization, trending NNW–SSE going through the caldera, and a low magnetization zone ∼ 300 m by ∼ 500 m in area, extending over the Hakurei site. Comparison between the results obtained using the vector anomaly and the total intensity anomaly shows that the magnetic field is determined more accurately, especially in areas of sparse data distribution, when the vector anomaly rather than the total intensity anomaly is used. We suggest a geologically motivated model that basaltic volcanism associated with the backarc rifting occurred after the formation of the caldera, resulting in the formation of the high magnetization belt underneath the silicic caldera. The Hakurei hydrothermal field lies in the intersection of the basaltic volcanism belt and the caldera wall fault, suggesting a mechanism that hot water generated by the heat of the volcanic activity has been spouting out through the caldera wall fault. The deposit apparently extends beyond the low magnetization zone, climbing up the caldera wall. This may indicate that hot water rising from the deep through the alteration zone is transported laterally when it comes near the seafloor along fissures and fractures in the caldera wall.

Description: [1]  Using multiple ScS reverberations we examine mantle reflectivity structure beneath northeast China and the northwest Pacific. We find several upper mantle discontinuities, including a melt layer with a mean thickness of 64 km atop the 410-km discontinuity, present on both sides of the subducting slab near the Nankai trench. The transition zone contains a split 520-km discontinuity in several paths, and tomographic images show stagnant slabs at this depth. We believe this may be slab-related based on experimental work (Saikia, A., Frost, D. J., Rubie, D. C., 2008. Splitting of the 520-kilometer seismic discontinuity and chemical heterogeneity in the mantle. Science 319 (5869), 1515–1518). A negative reflector is found in one path beneath the northeast China craton at a depth of 598 km. Mid-mantle reflectors are found in all of our paths and are present throughout a wide depth range (~750 – 1600 km).

Description: [1]  Three-dimensional P - and S -wave velocity (V P , V S ) models and high-resolution earthquake relocations are determined for the New Madrid Seismic Zone using double-difference local earthquake tomography. The data set consists of arrival times and differential times recorded by the Cooperative New Madrid Seismic Network (CNMSN) from 2000-2007 and the 1989-1992 Portable Array Network and Data Acquisition deployment. Waveform cross-correlation derived differential times for the CNMSN data are also incorporated. The velocity solutions are compatible with previous solutions centered on the active arms of seismicity and cover a broader area including mafic intrusions along the margin of the Reelfoot Rift. Major features include elevated V P and V S associated with the mafic plutons and reduced V P and V S along and southeast of the Axial fault (AF), a major arm of seismicity trending along the rift axis. Low V P extends to a depth of at least 20 km along the portion of the AF that extends south of the Missouri bootheel. A locally high V P /V S anomaly imaged along the central portion of the Reelfoot fault is spatially correlated with a significant change in fault trend and is interpreted as a region containing high pore pressure and/or water-filled microcracks.

Description: [1]  We present a catalog of InSAR constraints on deformation that occurred during earthquake sequences in southern Iran between 1992-2011, and explore the implications on the accommodation of large-scale continental convergence between Saudi Arabia and Eurasia within the Zagros Mountains. The Zagros Mountains, a salt-laden fold-and-thrust-belt involving ~10 km of sedimentary rocks overlying Precambrian basement rocks, have formed as a result of ongoing continental collision since 10-20 Ma that is currently occurring at a rate of ~3 cm/yr. We first demonstrate that there is a biased misfit in earthquake locations in global catalogs that likely results from neglect of 3D velocity structure. Previous work involving two M ~ 6 earthquakes with well-recorded aftershocks has shown that the deformation observed with InSAR may represent triggered slip on faults much shallower than the primary earthquake, which likely occurred within the basement rocks (〉10 km depth). We explore the hypothesis that most of the deformation observed with InSAR spanning earthquake sequences is also due to shallow, triggered slip above a deeper earthquake, effectively doubling the moment release for each event. We quantify the effects that this extra moment release would have on the discrepancy between seismically and geodetically constrained moment rates in the region, finding that even with the extra triggered fault slip, significant aseismic deformation during the interseismic period is necessary to fully explain the convergence between Eurasia and Saudi Arabia.

Description: During the late Miocene, a dramatic global expansion of C 4 plant distribution occurred with broad spatial and temporal variations. Although the event is well documented, whether subsequent expansions were caused by a decreased atmospheric CO 2 concentration or climate change is a contentious issue. In the present study, we used an improved inverse vegetation modeling approach that accounts for the physiological responses of C 3 and C 4 plants to quantitatively reconstruct the paleoclimate in the Siwalik of Nepal based on pollen and carbon isotope data. We also studied the sensitivity of the C 3 and C 4 plants to changes in the climate and the atmospheric CO 2 concentration. We suggest that the expansion of the C 4 plant distribution during the late Miocene may have been primarily triggered by regional aridification and temperature increases. The expansion was unlikely caused by reduced CO 2 levels alone. Our findings suggest that this abrupt ecological shift mainly resulted from climate changes related to the decreased elevation of the Himalayan foreland. This article is protected by copyright. All rights reserved.

Description: [1]  The identification and evaluation of trigger mechanisms for volcano flank instabilities and/or collapse represent a key issue for risk assessment in densely populated volcanic areas, as well as in long distance settings, particularly in case of island or coastal volcanoes. Here, we address quantitatively the effects of external (seismic) and inner (magmatic) forcing on the stress-strain state associated to flank instabilities at Mt. Etna (Sicily, southern Italy) by means of a 2-D Finite-Difference-Method numerical modelling. Modelled seismic actions include strong near-field, strong far-field and low-magnitude near-field earthquakes. Magmatic actions consider the inner presssure changes induced by energetic lava fountains in the summit crater area, sub-vertical and oblique dyke ascent below the summit area. Model results are validated in light of available monitoring data and recent eruptive activity. Numerical results show that the main strain effects are produced by high-magnitude near-field earthquakes (expected return time of ~10 3  yrs), and by vertical uprise of a magma dyke below the volcano summit area. Maximum displacements in the order of tens of centimetres may involve the summit area, up to some 10 6  m 3 /m over some km laterally. Stress releases up to 10 7  Pa may affect a limited portion of the magmatic conduit, thus favouring major effusive flank eruptions. Major catastrophic events, such as volcano flank collapse, should not be expected by applying, either individually or combined, the aforementioned actions.

Description: [1]  Seismicity closely related to hydrological impacts has been observed in several locations worldwide; particularly in intraplate areas where tectonic stressing rates are small. The triggering mechanism is usually explained by a poroelastic response of the seismogenic crust to surface water flux, leading to pore pressure changes at depth. To explain the earthquake triggering in response of those small stress changes, however, the crust has to be near a critical state in which other transient processes might be significant. One of the prominent examples is the Mt. Hochstaufen in SW Germany, where seismicity is known to vary seasonally. A previous analysis showed that the seismicity in 2002 was highly correlated with model forecasts based on fluid diffusion and rate- and state-dependent frictional nucleation. Here we revisit this case by accounting additionally for poroelastic effects, as well as for thermoelastic and tidal stresses. We also test whether the model can explain the observations of the subsequent eight years between 2003 and 2010. Our analysis confirms that rainfall is the dominant driving force in this region. The model not only fits the year 2002 activity very well, but provides with the same parameters a reasonable fit to the subsequent period, with a probability gain of about 4 per event in comparison to a time-independent Poisson model.

Description: ABSTRACT [1]  The absolute magnitude of stress in the crust and the shear strength of faults are poorly known, yet fundamental quantities, in lithospheric dynamics. While stress magnitude cannot be measured directly, deviatoric stress state can be inferred indirectly from focal mechanism solutions collected before and after an earthquake. We extend a standard stress inversion for normalized stresses to invert for the 3D spatial distribution of absolute deviatoric stress and variation of fault strength with depth using focal mechanism solutions and coseismic stress changes produced by large earthquakes. We apply the method to the 2011 M9 Tohoku-oki, Japan earthquake. The northern Japan forearc crust between 5 and 15 km depth appears to be weak with fault strength of 40–90 MPa, consistent with a coefficient of friction of 0.2-0.5. The M9 Tohoku-oki coseismic stress change was large enough, relative to the ambient stress, to rotate the principal stress directions typically ~20° in the upper 20 km of the crust. The data from Japan require a heterogeneous ambient deviatoric stress field with short wavelength (~20-50 km) fluctions in principal stress orientations.

Description: [1]  We review marine heat flow data along the Nankai Trough and show that observations 〉 30 km seaward of the deformation front are 20% below conductive predictions (129–94 mW m -2 ) but consistent with the global heat flow average for oceanic crust of the same age (16-28 Ma). Heat flow values  〈  30 km seaward of the deformation front are generally 20% higher than conductive predictions. This heat flow pattern is consistent with the advection of heat by fluid flow in the subducting oceanic crust and explains both the high heat flux in the vicinity of the trench, 〉 200 and 〉 140 mW m -2 , and steep landward declines to values of approximately 60 mW m -2 over distances of 65 and 50 km along the Muroto and Kumano transects, respectively. Along the Ashizuri transect the lack of heat flow data preclude a definitive interpretation. We conclude that fluid flow in the subducting oceanic crust leads to temperatures that are generally 25 ° C higher near the toe of the margin wedge and 50 - 100 ° C lower near the downdip limit of the seismogenic zone than estimated by purely conductive models.

Description: Projections of future changes in land carbon (C) storage using biogeochemical models depend on accurately modeling the interactions between the C and nitrogen (N) cycles. Here, we present a framework for analyzing N limitation in global biogeochemical models to explore how C-N interactions of current models compare to field observations, identify the processes causing model divergence, and identify future observation and experiment needs. We used a set of N fertilization simulations from two global biogeochemical models (CLM-CN and O-CN) that use different approaches to modeling C-N interactions. On the global scale, net primary productivity (NPP) in the CLM-CN model was substantially more responsive to N fertilization than in the O-CN model. The most striking difference between the two models occurred for humid tropical forests, where the CLM-CN simulated a 62% increase in NPP at high N addition levels (30 g N m −2 yr −1 ), while the O-CN predicted a 2% decrease in NPP due to N fertilization increasing plant respiration more than photosynthesis. Across 35 temperate and boreal forest sites with field N fertilization experiments, we show that the CLM-CN simulated a 46% increase in aboveground NPP in response to N, which exceeded the observed increase of 25%. In contrast, the O-CN only simulated a 6% increase in aboveground NPP at the N fertilization sites. Despite the small response of NPP to N fertilization, the O-CN model accurately simulated ecosystem retention of N and the fate of added N to vegetation when compared to empirical 15 N tracer application studies. In contrast, the CLM-CN predicted lower total ecosystem N retention and partitioned more losses to volatilization than estimated based from observed N budgets of small catchments. These results point to the need for model improvements for both models to enhance the accuracy with which global C-N cycle feedbacks can be simulated. This article is protected by copyright. All rights reserved.

Description: This study reports the first well-replicated analysis of continuous coral growth records from warmer-water reefs (mean annual SST 〉28.5°C) around the Thai-Malay Peninsula in Southeast Asia. Based on analyses of 70 colonies sampled from 15 reefs within six locations, region-wide declines in coral calcification rate (~18.6%), linear extension rate (~15.4%) and skeletal bulk density (~3.9%) were observed over a 31-year period from 1980–2010. Decreases in calcification and linear extension rates were observed at five of the six locations and ranged from ~17.2-21.6% and ~11.4–19.6% respectively, while decline in skeletal bulk density was a consequence of significant reductions at only two locations (~6.9% and ~10.7%). A significant link between region-wide growth rates and average annual SST was found, and Porites spp. demonstrated a high thermal threshold of ~29.4°C before calcification rates declined. Responses at individual locations within the region were more variable with links between SST and calcification rates being significant at only four locations. Rates of sea temperature warming at locations in the Andaman Sea (Indian Ocean) (~1.3°C decade −1 ) were almost twice those in the South China Sea (Pacific Ocean) (~0.7°C decade −1 ), but this was not reflected in the magnitude of calcification declines at corresponding locations. Considering that massive Porites spp. are major reef-builders around Southeast Asia, this region-wide growth decline is a cause for concern for future reef accretion rates and resilience. However, this study suggests that the future rates and patterns of change within the region are unlikely to be uniform or dependent solely on the rates of change in the thermal environment. This article is protected by copyright. All rights reserved.

Description: Although striking changes have been documented in plant and animal phenology over the past century, less is known about how the fungal kingdom's phenology has been changing. A few recent studies have documented changes in fungal fruiting in Europe in the last few decades, but the geographic and taxonomic extent of these changes, the mechanisms behind these changes, and their relationships to climate, are not well understood. Here, we analyzed herbarium data of 274 species of fungi from Michigan to test the hypotheses that fruiting times of fungi depend on annual climate, and that responses depend on taxonomic and functional groups. We show that the fungal community overall fruits later in warmer and drier years, which has led to a shift toward later fruiting dates for autumn-fruiting species, consistent with existing evidence. However, we also show that these effects are highly variable among species and are partly explained by basic life history characteristics. Resulting differences in climate sensitivities are expected to affect community structure as climate changes. This study provides a unique picture of the climate-dependence of fungal phenology in North America and an approach for quantifying how individual species and broader fungal communities will respond to ongoing climate change. This article is protected by copyright. All rights reserved.

Description: A unique long-term phenological dataset of over 110,000 records of 1st cutting dates for haymaking across Germany, spanning the years 1951-2011 was examined. In addition, we analysed a long-term dataset on the beginning of flowering of meadow foxtail ( Alopecurus pratensis ) covering the last 20 years. We tested whether hay cutting dates (based on a human decision when to cut) showed trends, temperature relationships and spatial distribution similar to the development of this grassland species, and if these trends could be related to climate change. The timing of 1st hay cut was strongly influenced (p 〈 0.001) by altitude, latitude and longitude, revealing in particular an east-west gradient. Over the past 60 years there have been changes in the timing of hay cutting, with the majority of German federal states having significant (p 〈 0.05) advances of approximately 1 day per decade. Overall, the response to mean March- May temperature was highly significant (-2.87 days °C −1 ; p 〈 0.001). However, in the last 20 years no federal state experienced a significant advance and two were even significantly delayed. The temperature response in this post-1991 period became less or non significant for most of the federal states. We suggest that differences in agricultural land use and unequal uptakes of Agri-Environment Schemes (AES, which encourage later cutting) were likely to be responsible for the regional differences, while the general increase in AES appears to have confounded the overall trend in hay cutting in the last 20 years. Trends over time and responses to temperature were small relative to those associated with the phenology of meadow foxtail. The advance in phenology of this species is greater than the advance in hay cutting, implying that hay cutting may not be keeping pace with a changing climate, which may have a positive effect on grassland ecology. This article is protected by copyright. All rights reserved.

Description: [1]  We investigate whether predictions of mantle structure from tectonic reconstructions are in agreement with a detailed tomographic image of seismic P-wave velocity structure under the Caribbean region. In the upper mantle, positive seismic anomalies are imaged under the Lesser Antilles and Puerto Rico. These anomalies are interpreted as remnants of Atlantic lithosphere subduction and confirm tectonic reconstructions that suggest at least 1100 km of convergence at the Lesser Antilles island arc during the past ~45 Myr. The imaged Lesser-Antilles slab consists of a northern and southern anomaly, separated by a low velocity anomaly across most of the upper mantle, which we interpret as the subducted North America-South America plate boundary. The southern edge of the imaged Lesser Antilles slab agrees with vertical tearing of South America lithosphere. The northern Lesser Antilles slab is continuous with the Puerto Rico slab along the northeastern plate boundary. This results in an amphitheatre-shaped slab and it is interpreted as westward subducting North America lithosphere that remained attached to the surface along the northeastern boundary of the Caribbean plate. At the Muertos Trough, however, material is imaged until a depth of only 100 km, suggesting a small amount of subduction. The location and length of the imaged South Caribbean slab agrees with proposed subduction of Caribbean lithosphere under the northern South America plate. An anomaly related to proposed Oligocene subduction at the Nicaragua rise is absent in the tomographic model. Beneath Panama, a subduction window exists across the upper mantle, which is related to the cessation of subduction of the Nazca plate under Panama since 9.5 Ma and possibly the preceding subduction of the extinct Cocos-Nazca spreading center. In the lower mantle two large anomaly patterns are imaged. The westernmost anomaly agrees with the subduction of Farallon lithosphere. The second lower mantle anomaly is found east of the Farallon anomaly and is interpreted as a remnant of the late Mesozoic subduction of North and South America oceanic lithosphere at the Greater Antilles, Aves ridge and Leeward Antilles. The imaged mantle structure does not allow us to discriminate between an ‘Intra-Americas’ origin and a ‘Pacific origin’ of the Caribbean plate.

Description: [1]  We have measured interseismic deformation across the Ashkabad strike-slip fault using 13 Envisat interferograms covering a total effective timespan of ~30 years. Atmospheric contributions to phase delay are significant and variable due to the close proximity of the Caspian Sea. In order to retrieve the pattern of strain accumulation, we show it is necessary to use data from Envisat's Medium Resolution Imaging Spectrometer (MERIS) instrument, as well numerical weather model outputs from the European Centre for Medium-Range Weather Forecasting (ECMWF), to correct interferograms for differences in water vapour and atmospheric pressure respectively. This has enabled us to robustly estimate the slip rate and locking depth for the Ashkabad fault using a simple elastic dislocation model. Our data are consistent with a slip rate of 5–12 mm/yr below a locking depth of 5.5–17 km for the Ashkabad fault, and synthetic tests support the magnitude of the uncertainties on these estimates. Our estimate of slip rate is 1.25–6 times higher than some previous geodetic estimates, with implications for both seismic hazard and regional tectonics, in particular supporting fast relative motion between the South Caspian Block and Eurasia. This result reinforces the importance of correcting for atmospheric contributions to interferometric phase for small strain measurements. We also attempt to validate a recent method for atmospheric correction based on ECMWF ERA-Interim model outputs alone and find that this technique does not work satisfactorily for this region when compared to the independent MERIS estimates.

Description: [1]  The seafloor within the Perth Abyssal Plain (PAP), offshore Western Australia, is the only section of crust that directly records the early spreading history between India and Australia during the Mesozoic breakup of Gondwana. However, this early spreading has been poorly constrained due to an absence of data, including marine magnetic anomalies and data constraining the crustal nature of key tectonic features. Here, we present new magnetic anomaly data from the PAP that shows that the crust in the western part of the basin was part of the Indian Plate – the conjugate flank to the oceanic crust immediately offshore the Perth margin, Australia. We identify a sequence of M2 and older anomalies in the west PAP within crust that initially moved with the Indian Plate, formed at intermediate half-spreading rates (35 mm/yr) consistent with the conjugate sequence on the Australian Plate. More speculatively, we reinterpret the youngest anomalies in the east PAP, finding that the M0-age crust initially formed on the Indian Plate was transferred to the Australian Plate by a westward jump or propagation of the spreading ridge shortly after M0 time. Samples dredged from the Gulden Draak and Batavia Knolls (at the western edge of the PAP) reveal that these bathymetric features are continental fragments rather than igneous plateaus related to Broken Ridge. These microcontinents rifted away from Australia with Greater India during initial breakup at ~130 Ma, then rifted from India following the cessation of spreading in the PAP (~101-103 Ma).

Description: [1]  We develop a three-step Maximum-A-Posteriori probability (MAP) method for coseismic rupture inversion, which aims at maximizing the a posterior probability density function (PDF) of elastic deformation solutions of earthquake rupture. The method originates from the Fully Bayesian Inversion (FBI) and Mixed linear-nonlinear Bayesian inversion (MBI) methods, shares the same posterior PDF with them, while overcoming difficulties with convergence when large numbers of low-quality data are used and greatly improving the convergence rate using optimization procedures. A highly-efficient global optimization algorithm, Adaptive Simulated Annealing (ASA), is used to search for the maximum of a posterior PDF (" mode " in statistics) in the first step. The second step inversion approaches the " true" solution further using the Monte Carlo Inversion (MCI) technique with positivity constraints, with all parameters obtained from step one as the initial solution. Then slip artifacts are eliminated from slip models in the third step using the same procedure of the second step, with fixed fault geometry parameters. [2]  We first design a fault model with 45°-dip angle and oblique slip, and produce corresponding synthetic InSAR datasets to validate the reliability and efficiency of the new method. We then apply this method to InSAR data inversion for the coseismic slip-distribution of the April 14, 2010 Mw 6.9 Yushu, China earthquake. Our preferred slip model is composed of three segments with most of the slip occurring within 15 km depth and the maximum slip reaches 1.38 m at the surface. The seismic moment released is estimated to be 2.32e + 19 Nm, consistent with the seismic estimate of 2.50e + 19 Nm.

Description: The rate of vegetation recovery from boreal wildfire influences terrestrial carbon cycle processes and climate feedbacks by affecting the surface energy budget and land-atmosphere carbon exchange. Previous forest recovery assessments using satellite optical-infrared normalized difference vegetation index (NDVI) and tower CO 2 eddy covariance techniques indicate rapid vegetation recovery within 5 to 10 years, but these techniques are not directly sensitive to changes in vegetation biomass. Alternatively, the vegetation optical depth (VOD) parameter from satellite passive microwave remote sensing can detect changes in canopy biomass structure and may provide a useful metric of post-fire vegetation response to inform regional recovery assessments. We analyzed a multi-year (2003-2010) satellite VOD record from the NASA AMSR-E (Advanced Microwave Scanning Radiometer for EOS) sensor to estimate forest recovery trajectories for 14 large boreal fires from 2004 in Alaska and Canada. The VOD record indicated initial post-fire canopy biomass recovery within 3 to 7 years, lagging NDVI recovery by 1 to 5 years. The VOD lag was attributed to slower non-photosynthetic (woody) and photosynthetic (foliar) canopy biomass recovery, relative to the faster canopy greenness response indicated from the NDVI. The duration of VOD recovery to pre-burn conditions was also directly proportional (p〈0.01) to satellite (MODIS) estimated tree cover loss used as a metric of fire severity. Our results indicate that vegetation biomass recovery from boreal fire disturbance is generally slower than reported from previous assessments based solely on satellite optical-infrared remote sensing, while the VOD parameter enables more comprehensive assessments of boreal forest recovery. This article is protected by copyright. All rights reserved.

Description: Shifts in precipitation regimes are an inherent component of climate change, but in low energy systems are often assumed to be less important than changes in temperature. Because soil moisture is the hydrological variable most proximally linked to plant performance during the growing season in arctic-alpine habitats, it may offer the most useful perspective on the influence of changes in precipitation on vegetation. Here we quantify the influence of soil moisture for multiple vegetation properties at fine spatial scales, to determine the potential importance of soil moisture under changing climatic conditions. A fine-scale dataset, comprising vascular species cover and field-quantified ecologically-relevant environmental parameters, was analysed to determine the influence of soil moisture relative to other key abiotic predictors. Soil moisture was strongly related to community composition, species richness and the occurrence patterns of individual species, having a similar or greater influence than soil temperature, pH and solar radiation. Soil moisture varied considerably over short distances, and this fine-scale heterogeneity may contribute to offsetting the ecological impacts of changes in precipitation for species not limited to extreme soil moisture conditions. In conclusion, soil moisture is a key driver of vegetation properties, both at the species- and community-level, even in this low energy system. Soil moisture conditions represent an important mechanism through which changing climatic conditions impact vegetation, and advancing our predictive capability will therefore require a better understanding of how soil moisture mediates the effects of climate change on biota. This article is protected by copyright. All rights reserved.

Description: Global nitrogen (N) enrichment has resulted in increased nitrous oxide (N 2 O) emission that greatly contributes to climate change and stratospheric ozone destruction, but little is known about the N 2 O emissions from urban river networks receiving anthropogenic N inputs. We examined N 2 O saturation and emission in the Shanghai city river network, covering 6300 km 2 , over 27 months. The overall mean saturation and emission from 87 locations was 770% and 1.91 mg N 2 O-N•m −2 •d −1 , respectively. N 2 O saturation did not exhibit a clear seasonality, but the temporal pattern was co-regulated by both water temperature and N loadings. Rivers draining through urban and suburban areas receiving more sewage N inputs had higher N 2 O saturation and emission than those in rural areas. Regression analysis indicated that water ammonium (NH 4 + ) and dissolved oxygen (DO) level had great control on N 2 O production and were better predictors of N 2 O emission in urban watershed. About 0.29 Gg N 2 O-N•yr −1 N 2 O was emitted from the Shanghai river network annually, which was about 131% of IPCC's prediction using default emission values. Given the rapid progress of global urbanization, more study efforts, particularly on nitrification and its N 2 O yielding, are needed to better quantify the role of urban rivers in global riverine N 2 O emission. This article is protected by copyright. All rights reserved.

Description: [1]  We conducted total magnetic field and Bouguer gravity measurements to investigate the shallow structure beneath the summit caldera of Kīlauea Volcano, Hawai‘i. Two significant and distinctive magnetic anomalies were identified within the caldera. One is interpreted to be associated with a long-lived pre-historic eruptive centre, the Observatory vent, located ~1 km east of the Hawaiian Volcano Observatory. The second magnetic anomaly corresponds to a set of eruptive fissures that strike northeast from Halema‘uma‘u Crater, suggesting this is an important transport pathway for magma. The Bouguer gravity data were inverted to produce 3D models of density contrasts in the upper 2 km beneath Kīlauea. The models detect 3.0 km 3 of material, denser than 2800 kg m -3 , beneath the caldera that may represent an intrusive complex centred northeast of Halema‘uma‘u. Recent temporal gravity studies indicate continual addition of mass beneath the caldera during 1975–2008 centred west of Halema‘uma‘u and suggest this is due to filling of void space. The growth of a large intrusive complex, apparent cyclical caldera formation, and continual mass addition without inflation, however, can also be explained by extensional rifting caused by the continual southward movement of Kīlauea's unstable south flank.

Description: [1]  The 2006-2007 doublet of M W  〉 8 earthquakes in the Kuril subduction zone caused postseismic transient motion in the asthenosphere, which we observed on the Kuril GPS Array in 2007–2011. Here we show that the Maxwell asthenospheric viscosity that best fits the geodetic data increased by nearly an order of magnitude over the interval of four years, from 2 × 10 17 to 1 × 10 18  Pa s. These effective values of viscosity can be explained by a power-law rheology for which strain rate is proportional to stress raised to a power n  〉 1. The apparent change in viscosity can also be caused by other factors such as coupling between afterslip and viscoelastic flow. The open and intriguing question in connection with postseismic data after the Kuril earthquake doublet is the magnitude of the long-term asthenospheric viscosity, which shall be revealed by continued observations. An asthenosphere with viscosity of about 1 × 10 19  Pa s is favored by the postseismic deformation still observed several decades after the 1960 Chile and 1964 Alaska M W ~9 earthquakes. However, postseismic deformation associated with the 1952 southern Kamchatka M W ~9 earthquake currently is not observed in the northern Kurils, an indication that the long-term asthenospheric viscosity in the Kurils is lower than in Chile and Alaska.

Description: [1]  The region of central Chile offers a unique opportunity to study the links between the subducting Juan Fernandez Ridge, the flat slab, the Double Seismic Zone (DSZ) and the absence of modern volcanism. Here, we report the presence and characteristics of the first observed DSZ within the intermediate-depth Nazca slab using two temporary seismic catalogues (OVA99 and CHARSME). The lower plane of seismicity (LP) is located 20–25 km below the upper plane (UP), begins at 50 km depth and merges with the lower plane at 120 km depth, where the slab becomes horizontal. Focal mechanism analysis and stress tensor calculations indicate that the slab's state of stress is dominantly controlled by plate convergence and overriding crust thickness: Above 60–70 km depth, the slab is in horizontal compression, and below, it is in horizontal extension, parallel to plate convergence, which can be accounted for by vertical loading of the overriding lithosphere. Focal mechanisms below 60–70 km depth are strongly correlated with offshore outer rise bend faults, suggesting the reactivation of pre-existing faults below this depth. The large interplane distances for all Nazca DSZs can be related to the slab's unusually cold thermal structure with respect to its age. Since LPs globally seem to mimic mantle mineral dehydration paths, we suggest that fluid migration and dehydration embrittlement provide the mechanism necessary to weaken the rock and that the stress field determines the direction of rupture.

Description: Our ability to project the impact of global change on marine ecosystem is limited by our poor understanding on how to predict species sensitivity. For example, the impact of ocean acidification is highly species-specific, even in closely related taxa. The aim of this study was to test the hypothesis that the tolerance range of a given species to decreased pH corresponds to their natural range of exposure. Larvae of the green sea urchin Strongylocentrotus droebachiensis were cultured from fertilization to metamorphic competence (29 days) under a wide range of pH (from pH T =8.0/ p CO 2 ≈480μatm to pH T =6.5/ p CO 2 ≈20000μatm) covering present (from pH T 8.7 to 7.6), projected near-future's variability (from pH T 8.3 to 7.2) and beyond. Decreasing pH impacted all tested parameters (mortality, symmetry, growth, morphometry and respiration). Development of normal, although showing morphological plasticity, swimming larvae was possible as low as pH T ≥7.0. Within that range, decreasing pH increased mortality and asymmetry and decreased body length growth rate. Larvae raised at lowered pH and with similar body length had shorter arms and a wider body. Relative to a given body length, respiration rates and stomach volume both increased with decreasing pH suggesting changes in energy budget. At the lowest pHs (pH T ≤6.5), all the tested parameters were strongly negatively affected and no larva survived past 13 days post-fertilization. In conclusion, sea urchin larvae appeared to be highly plastic when exposed to decreased pH until a physiological tipping point at pH T =7.0. However, this plasticity was associated with direct (increased mortality) and indirect (decreased growth) consequences for fitness. This article is protected by copyright. All rights reserved.

Description: With a pace of about twice the observed rate of global warming, the temperature on the Qinghai-Tibetan Plateau (Earth's “third pole”) has increased by 0.2 °C per decade over the past 50 years, which results in significant permafrost thawing and glacier retreat. Our review suggested that warming enhanced net primary production (NPP) and soil respiration, decreased methane (CH 4 ) emissions from wetlands and increased CH 4 consumption of meadows, but might increase CH 4 emissions from lakes. Warming induced permafrost thawing and glaciers melting would also result in substantial emission of old carbon dioxide (CO 2 ) and CH 4 . Nitrous oxide (N 2 O) emission was not stimulated by warming itself, but might be slightly enhanced by wetting. However, there are many uncertainties in such biogeochemical cycles under climate change. Human activities (e.g., grazing, land cover changes) further modified the biogeochemical cycles and amplified such uncertainties on the plateau. If the projected warming and wetting continues, the future biogeochemical cycles will be more complicated. So facing research in this field is an ongoing challenge of integrating field observations with process-based ecosystem models to predict the impacts of future climate change and human activities at various temporal and spatial scales. To reduce the uncertainties and improve the precision of the predictions of the impacts of climate change and human activities on biogeochemical cycles, efforts should focus on conducting more field observation studies, integrating data within improved models, and developing new knowledge about coupling among carbon, nitrogen, and phosphorus biogeochemical cycles as well as about the role of microbes in these cycles. This article is protected by copyright. All rights reserved.

Description: Some species are adapting to changing environments by expanding their geographic ranges. Understanding whether range shifts will be accompanied by increased exposure to other threats is crucial to predicting when and where new populations could successfully establish. If species overlap to a greater extent with human development under climate change, this could form ecological traps which are attractive to dispersing individuals, but the use of which substantially reduces fitness. Until recently, the core nesting range for the Critically Endangered Kemp's ridley sea turtle ( Lepidochelys kempii ) was ~1,000km of sparsely populated coastline in Tamaulipas, Mexico. Over the past twenty-five years, this species has expanded its range into populated areas of coastal Florida (〉1,500km outside the historical range), where nesting now occurs annually. Suitable Kemp's ridley nesting habitat has persisted for at least 140,000 years in the western Gulf of Mexico, and climate change models predict further nesting range expansion into the eastern Gulf of Mexico and northern Atlantic Ocean. Range expansion is 6-12% more likely to occur along uninhabited stretches of coastline than are current nesting beaches, suggesting that novel nesting areas will not be associated with high levels of anthropogenic disturbance. Although the high breeding-site fidelity of some migratory species could limit adaptation to climate change, rapid population recovery following effective conservation measures may enhance opportunities for range expansion. Anticipating the interactive effects of past or contemporary conservation measures, climate change, and future human activities will help focus long-term conservation strategies. This article is protected by copyright. All rights reserved.

Description: We combine satellite and ground observations during 1950-2011 to study the long-term links between multiple climate (air temperature and cryospheric dynamics) and vegetation (greenness and atmospheric CO2 concentrations) indicators of the growing season of northern ecosystems (〉45oN) and their connection with the carbon cycle. During the last three decades, the thermal potential growing season has lengthened by about 10.5 days ( p 〈 0.01, 1982–2011), which is unprecedented in the context of the past 60 years. The overall lengthening has been stronger and more significant in Eurasia (12.6 days, p 〈 0.01) than North America (6.2 days, p 〉 0.05). The photosynthetic growing season has closely tracked the pace of warming and extension of the potential growing season in spring, but not in autumn when factors such as light and moisture limitation may constrain photosynthesis. The autumnal extension of the photosynthetic growing season since 1982 appears to be about half that of the thermal potential growing season, yielding a smaller lengthening of the photosynthetic growing season (6.7 days at circumpolar scale, p 〈 0.01). Nevertheless, when integrated over the growing season, photosynthetic activity has closely followed the interannual variations and warming trend in cumulative growing season temperatures. This lengthening and intensification of the photosynthetic growing season, manifested principally over Eurasia rather than North America, is associated with a long-term increase (22.2% since 1972, p 〈 0.01) in the amplitude of the CO2 annual cycle at northern latitudes. The springtime extension of the photosynthetic and potential growing seasons has apparently stimulated earlier and stronger net CO2 uptake by northern ecosystems, while the autumnal extension is associated with an earlier net release of CO2 to the atmosphere. These contrasting responses may be critical in determining the impact of continued warming on northern terrestrial ecosystems and the carbon cycle. This article is protected by copyright. All rights reserved.

Description: [1]  Among the different types of waves embedded in seismic noise, body waves present appealing properties but are still challenging to extract. Here we first validate recent improvements in numerical modeling of microseismic compressional ( P ) body waves and then show how this tool allows fast detection and location of their sources. We compute sources at ~ 0.2 Hz within typical P teleseismic distances (30-90 degrees) from the South California Seismic Network (SCSN) and analyze the most significant discrete sources. The locations and relative strengths of the computed sources are validated by the good agreement with beam-forming analysis. These ~75 noise sources exhibit a highly heterogeneous distribution, and cluster along the usual storm tracks in the Pacific and Atlantic oceans. They are mostly induced in the open ocean, at or near water depths of 2800 and 5600 km, most likely within storms or where ocean waves propagating as swell meet another swell or wind sea. We then emphasize two particularly strong storms to describe how they generate noise sources in their wake. We also use these two specific noise bursts to illustrate the differences between microseismic body- and surface-waves in terms of source distribution and resulting recordable ground motion. The different patterns between body- and surface-waves result from distinctive amplification of ocean wave-induced pressure perturbation and different seismic attenuation. Our study demonstrates the potential of numerical modeling to provide fast and accurate constraints on where and when to expect microseismic body waves, with implications for seismic imaging and climate studies.

Description: [1]  Earthquakes that rupture across steps between faults can be larger than those predicted from individual fault lengths, making understanding multifault events critical to assessing earthquake hazard. Empirical data from earthquake surface ruptures suggest that the distances between faults that rupture together can range from 〈1 km to 5 km. Dynamic and quasi-static models of planar faults determine similar distances. However, studies of interactions between realistic, 3D non-planar faults are few. A general comparison of quasi-static stress perturbations and triggering potentials with mechanical models incorporating either planar or non-planar faults highlight the sensitivity of planar fault models to model parameters and reveal no clear relationship between mean fault slip and triggering potential. More specifically, planar fault models predict triggering across a 3 km extensional step, while models incorporating non-planar faults indicate that a connecting fault is necessary to transfer slip through a 3 km step along the 1992 Landers, California earthquake rupture. The mechanical approach taken captures the stress changes as well as the total stress following fault slip, improving the criterion used to determine triggered failure potential. This underscores the need for additional constraint on fault strength and cohesion. The focus on complex fault geometry restricts analyses to the quasi-static realm, limiting the application of results to fault interactions over the short distances and slow rupture velocities for which the quasi-static stress field is relevant or approximates the dynamic stress field.

Description: [1]  Batch and flow-through experiments were performed on quartz-feldspar granular aggregates and sandstone samples to investigate time-dependent effects of fluid-rock interactions on fluid and rock conductivity, respectively. The experiments were conducted at temperatures up to 164, at confining and pore pressures up to 10 and 5 MPa, respectively, and for up to 136 days. It showed that changes in rock conductivity were unequivocally related to changes in pore fluid conductivity. It is inferred that these changes were dependent on kinetically controlled dissolution reactions between the mineral grains and the fluid. The time-dependent signature of rock conductivity implied a detectable transition from initial dissolution towards some state of equilibrium. The response of rock conductivity to temperature changes followed an Arrhenius-type behavior. An exploratory kinetic evaluation of the conductivity data for sandstone samples yielded an apparent activation energy of approximately 32 kJ/mol. A concurrent chemical fluid analysis showed that this is an integrated value over all reactions occurring in parallel within a sample. These reactions namely concern silica and silicate dissolutionbut also the dissolution of accessory salt minerals. It is concluded that measuring the evolution of rock conductivity in combination with chemical pore fluid analysis constitutes a powerful and quantitative tool for monitoring time-dependent changesin pore fluid chemistry and thus fluid-rock interactions in real time.

Description: [1]  We describe a multi-parameter experiment at Erebus volcano, Antarctica, employing Doppler radar, video, acoustic, and seismic observations to estimate the detailed energy budget of large (up to 40-m-diameter) bubble bursts from a persistent phonolite lava lake. These explosions are readily studied from the crater rim at ranges of less than 500 m, and present an ideal opportunity to constrain the dynamics and mechanism of magmatic bubble bursts that can drive Strombolian and Hawaiian eruptions. We estimate the energy budget of the first second of a typical Erebus explosion as a function of time and energy type, and constrain gas pressures and forces using an analytic model for the expansion of a gas bubble above a conduit that incorporates conduit geometry and magma and gas parameters. The model, consistent with video and radar observations, invokes a spherical bulging surface with a base diameter equal to that of the lava lake. The model has no ad hoc free parameters, and geometrical calculations predict zenith height, velocity and acceleration during shell expansion. During explosions, the energy contained in hot over-pressured gas bubbles is freed and partitioned into other energy types, where by far the greatest non-thermal energy component is the kinetic and gravitational potential energy of the accelerated magma shell (〉10 9 J). Seismic source energy created by explosions is estimated from radar measurements and is consistent with source energy determined from seismic observations. For the generation of the infrasonic signal, a dual mechanism incorporating a terminally disrupted slug is proposed, which clarifies previous models and provides good fits to observed infrasonic pressures. A new and straightforward method is presented for determining gas volumes from slug explosions at volcanoes from remote infrasound recordings.

Description: [1]  We explore the application of GPS data to earthquake early warning and investigate whether the co-seismic ground deformation can be used to provide fast and reliable magnitude estimations and ground shaking predictions. We use an algorithm to extract the permanent static offset from GPS displacement time series and invert for the slip distribution on the fault plane, which is discretized into a small number of rectangular patches. We developed a completely “self-adapting” strategy in which the initial fault plane model is built based on a quick, approximate magnitude estimation, and is then allowed to increase in size based on the evolutionary magnitude estimation resulting from the slip inversion. Two main early warning outputs are delivered in real-time: magnitude and the along-strike extent of the rupture area. These are finally used to predict the expected ground shaking due to the finite source. We tested the proposed strategy by simulating real-time environments for three earthquakes. For the Mw 9.0, 2011 Tohoku-Oki earthquake our algorithm provides the first magnitude estimate of 8.2 at 39 sec after the origin time, and then gradually increases to 8.9 at 120 sec. The estimated rupture length remains constant from the outset at ~360 km. For the Mw 8.3, 2003 Tokachi-Oki earthquake the initial magnitude estimate is 8.5 at 24 sec and drops to 8.2 at 40 sec with a rupture length of 290 km. Finally, for the Mw 7.2, 2010 El Mayor-Cucapah earthquake the magnitude estimate is 7.0 from the outset with a rupture length of 140 km. The accuracy of the ground shaking prediction using the GPS-based magnitude and finite extent is significantly better than existing seismology-based point source approaches. This approach would also facilitate more rapid tsunami warnings

Description: Coastal wetlands have the capacity to retain and denitrify large quantities of reactive nitrogen (N), making them important in attenuating increased anthropogenic N flux to coastal ecosystems. The ability of coastal wetlands to retain and transform N is being reduced by wetland losses resulting from land development. Nitrogen retention in coastal wetlands is further threatened by the increasing frequency and spatial extent of saltwater-inundation in historically freshwater ecosystems, due to the combined effects of dredging, declining river discharge to coastal areas due to human water use, increased drought frequency, and accelerating sea-level rise. Because saltwater incursion may affect N cycling through multiple mechanisms, the impacts of salinization on coastal freshwater wetland N retention and transformation are not well understood. Here, we show that repeated annual saltwater incursion during late summer droughts in the coastal plain of North Carolina changed N export from organic to inorganic forms and led to a doubling of annual NH 4 + export from a 440 hectare former agricultural field undergoing wetland restoration. Soil solution NH 4 + concentrations in two mature wetlands also increased with salinization, but the magnitude of increase was smaller than in the former agricultural field. Long-term saltwater exposure experiments with intact soil columns demonstrated that much of the increase in reactive N released could be explained by exchange of salt cations with sediment NH 4 + . Using these findings together with the predicted flooding of 1661 km 2 of wetlands along the NC coast by 2100, we estimate that saltwater incursion into these coastal areas could release up to 18,077 Mg N, or approximately half the annual NH 4 + flux of the Mississippi River. Our results suggest that that saltwater incursion into coastal freshwater wetlands globally could lead to increased N loading to sensitive coastal waters. This article is protected by copyright. All rights reserved.

Description: Recently there have been several studies using open top chambers (OTCs) or cloches to examine the response of Arctic plant communities to artificially elevated temperatures. Few, however, have investigated multi-trophic systems, or the effects of both temperature and vertebrate grazing treatments on invertebrates. This study investigated trophic interactions between an herbivorous insect ( Sitobion calvulum , Aphididae), a woody perennial host plant ( Salix polaris ) and a selective vertebrate grazer (barnacle geese, Branta leucopsis ). In a factorial experiment, the responses of the insect and its host to elevated temperatures using open top chambers (OTCs) and to three levels of goose grazing pressure were assessed over two summer growing seasons (2004 and 2005). OTCs significantly enhanced the leaf phenology of Salix in both years and there was a significant OTC by goose presence interaction in 2004. Salix leaf number was unaffected by treatments in both years, but OTCs increased leaf size and mass in 2005. Salix reproduction and the phenology of flowers were unaffected by both treatments. Aphid densities were increased by OTCs but unaffected by goose presence in both years. While goose presence had little effect on aphid density or host plant phenology in this system, the OTC effects provide interesting insights into the possibility of phenological synchrony disruption. The advanced phenology of Salix effectively lengthens the growing season for the plant, but despite a close association with leaf maturity, the population dynamics of the aphid appeared to lack a similar phenological response except for the increased population observed. This article is protected by copyright. All rights reserved.

Description: Evidence is accumulating that species’ responses to climate changes are best predicted by modelling the interaction of physiological limits, biotic processes and the effects of dispersal-limitation. Using commercially harvested blacklip ( Haliotis rubra ) and greenlip abalone ( H. laevigata ) as case studies, we determine the relative importance of accounting for interactions among physiology, metapopulation dynamics and exploitation in predictions of range (geographical occupancy) and abundance (spatially explicit density) under various climate change scenarios. Traditional correlative ecological niche models (ENM) predict that climate change will benefit the commercial exploitation of abalone by promoting increased abundances without any reduction in range size. However, models that account simultaneously for demographic processes and physiological responses to climate-related factors result in future (and present) estimates of area of occupancy and abundance that differ from those generated by ENMs alone. Range expansion and population growth are unlikely for blacklip abalone because of important interactions between climate-dependent mortality and metapopulation processes; in contrast, greenlip abalone should increase in abundance despite a contraction in area of occupancy. The strongly non-linear relationship between abalone population size and area of occupancy has important ramifications for the use of ENM predictions that rely on metrics describing change in habitat area as proxies for extinction risk. These results show that predicting species’ responses to climate change often require physiological information to understand climatic range determinants, and a metapopulation model that can make full use of this data to more realistically account for processes such as local extirpation, demographic rescue, source-sink dynamics and dispersal-limitation. This article is protected by copyright. All rights reserved.

Description: ABSTRACT [1]  Analysis of Lake Bonneville shorelines using LIDAR digital elevation data challenges accepted models of Wasatch fault deformation since the late Pleistocene. While footwall deformation of the Weber segment of the Wasatch fault is consistent with back-rotation of the footwall block and greatest displacement rate towards the center of the segment, shorelines along the footwall of the Salt Lake City segment decrease in elevation towards the interior and are highest at the segment boundaries, an opposite pattern of footwall deformation than predicted for boundaries arresting or strongly inhibiting displacement during earthquakes. The spatial pattern of footwall rebound implies that some of the proposed persistent fault segment boundaries do not stop earthquake ruptures that originate on adjacent fault segments, nor constrain ruptures initiated within the Salt Lake City segment. Net vertical fault displacement at the boundary between the Salt Lake and Provo segments is 16—20 m over the past 16.3—18.5 ka, corresponding to a vertical displacement rate of 0.8—1.2 mm/yr, a net fault slip rate of 2.0—2.8 mm/yr and horizontal extension rate of 1.8—2.6 mm/yr on the 25 o west-southwest dipping fault that forms the southern Salt Lake City segment boundary. Shoreline analysis suggests isostatic rebound caused by a drop in lake level was concentrated during a relatively short (~2000 yr) time period following the Bonneville flood at ~16 ka. LIDAR-derived topography in conjunction with robust geomorphic datums improves our ability to map deformation associated with lithospheric flexure and faulting while demonstrating the limitation of lacustrine shorelines in this type of analysis.

Description: [1]  We use recent results on statistical analysis of seismicity to present a robust method for comprehensive detection and analysis of earthquake clusters. The method is based on nearest-neighbor distances of events in space-time-energy domain. The method is applied to a 1981–2011 relocated seismicity catalog of southern California having 111,981 events with magnitudes m  ≥ 2, and corresponding synthetic catalogs produced by the Epidemic Type Aftershock Sequence (ETAS) model. Analysis of the ETAS model demonstrates that the cluster detection results are accurate and stable with respect to (i) three numerical parameters of the method, (ii) variations of the minimal reported magnitude, (iii) catalog incompleteness, and (iv) location errors. Application of the method to the observed catalog separates the 111,981 examined earthquakes into 41,393 statistically significant clusters comprised of foreshocks , mainshocks and aftershocks . The results reproduce the essential known statistical properties of earthquake clusters, which provide overall support for the proposed technique. In addition, systematic analysis with our method allows us to detect several new features of seismicity that include (i) existence of a significant population of single-event clusters ; (ii) existence of foreshock activity in natural seismicity that exceeds expectation based on the ETAS model; and (iii) dependence of all cluster properties, except area, on the magnitude difference of events from mainshocks but not on their absolute values. The classification of detected clusters into several major types, generally corresponding to singles, burst-like and swarm-like sequences, and correlations between different cluster types and geographic locations is addressed in a companion paper.

Description: [1]  For the first time, we report the amplitude variation with angle (AVA) pattern of bottom simulating reflectors (BSRs) beneath fracture-filled gas hydrate deposits when the effective medium is anisotropic. The common depth point (CDP) gathers of two mutually perpendicular multi-channel seismic profiles, located in the vicinity of Site NGHP-01-10, are appropriately processed such that they are fit for AVA analysis. AVA analysis of the BSR shows normal-incidence reflection coefficients of -0.04 to -0.11 with positive gradients of 0.04 to 0.31 indicating class IV pattern. The acoustic properties from isotropic rock physics model predict class III AVA pattern which cannot explain the observed class IV AVA pattern in Krishna-Godavari basin due to the anisotropic nature of fracture-filled gas hydrate deposits. [2]  We modeled the observed class IV AVA of the BSR by assuming that the gas hydrate bearing sediment can be represented by horizontally transversely isotropic (HTI) medium after accounting for anisotropic wave propagation effects on BSR amplitudes. The effective medium properties are estimated using Backus averaging technique and the AVA pattern of BSRs is modeled using the properties of overlying HTI and underlying isotropy/HTI media with or without free gas. Anisotropic AVA analysis of the BSR from the inline seismic profile shows 5–30 % gas hydrate concentration (equivalent to fracture density) and the azimuth of fracture system (fracture orientation) with respect to the seismic profile is close to 45°. Free gas below the base of gas hydrate stability zone is interpreted in the vicinity of fault system (F1).

Description: An elevated atmospheric CO 2 concentration ([CO 2 ]) can reduce stomatal conductance of leaves for most plant species, including rice ( Oryza sativa L.). However, few studies have quantified seasonal changes in the effects of elevated [CO 2 ] on canopy evapotranspiration, which integrates the response of stomatal conductance of individual leaves with other responses, such as leaf area expansion, changes in leaf surface temperature, and changes in developmental stages, in field conditions. We conducted a field experiment to measure seasonal changes in stomatal conductance of the uppermost leaves and in the evapotranspiration, transpiration, and evaporation rates using a lysimeter method. The study was conducted for flooded rice under open-air CO 2 elevation. Stomatal conductance decreased by 27% under elevated [CO 2 ], averaged throughout the growing season, and evapotranspiration decreased by an average of 5% during the same period. The decrease in daily evapotranspiration caused by elevated [CO 2 ] was more significantly correlated with air temperature and leaf area index rather than with other parameters of solar radiation, days after transplanting, vapor-pressure deficit and FAO reference evapotranspiration. This indicates that higher air temperatures, within the range from 16 to 27 °C, and a larger leaf area index, within the range from 0 to 4 m 2 m −2 , can increase the magnitude of the decrease in evapotranspiration rate caused by elevated [CO 2 ]. The crop coefficient (i.e., the evapotranspiration rate divided by the FAO reference evapotranspiration rate) was 1.24 at ambient [CO 2 ] and 1.17 at elevated [CO 2 ]. This study provides the first direct measurement of the effects of elevated [CO 2 ] on rice canopy evapotranspiration under open-air conditions using the lysimeter method, and the results will improve future predictions of water use in rice fields. © 2013 Blackwell Publishing Ltd

Description: Forests around the world are subject to risk of high rates of tree growth decline and increased tree mortality from combinations of climate warming and drought, notably in semi-arid settings. Here we assess how climate warming has affected tree growth in one of the world's most extensive zones of semi-arid forests, in Inner Asia, a region where lack of data limits our understanding of how climate change may impact forests. We show that pervasive tree growth declines since 1994 in Inner Asia have been confined to semi-arid forests where growing season water stress has been rising due to warming-induced increases in atmospheric moisture demand. A causal link between increasing drought and declining growth at semi-arid sites is corroborated by correlation analyses comparing annual climate data to records of tree-ring widths. These ring-width records tend to be substantially more sensitive to drought variability at semi-arid sites than at semi-humid sites. Fire occurrence and insect/pathogen attacks have increased in tandem with the most recent (2007-2009) documented episode of tree mortality. If warming in Inner Asia continues, further increases in forest stress and tree mortality could be expected, potentially driving the eventual regional loss of current semi-arid forests. © 2013 Blackwell Publishing Ltd

Description: [1]  Sequestration of large amounts of CO 2 within deep underground reservoirs has been proposed as a potential approach for reducing atmospheric emissions of greenhouse gases. A CO 2 sequestration project should address the associated environmental and safety issues and, in this respect, the importance of geomechanics has recently been widely recognized. Geomechanics is even more important when fluid injection is planned in faulted reservoirs. How much CO 2 can be safely injected into multi-compartment reservoirs? Are geomechanical constraints more restrictive than flow-dynamic constraints? These and other questions are addressed in the present study using a three-dimensional Finite Element - Interface Element geomechanical model. We simulate the possible mechanical failure in both the injected formation and the caprock, the fault/thrust reactivation, and the ground surface displacement in a faulted reservoir of the off-shore northern Italy, where seismic surveys provided an accurate characterization of the faulted geological structure. Based on reliable petrophysical/geomechanical properties from well-logs and pore overpressure as predicted by a fluid-dynamic model, the results show that the injection of 1 × 10 6  ton/a of CO 2 may be performed over a few years only. Thereafter part of the injected formation fails by shear stress. A number of parametric scenarios are investigated to address the major uncertainties on the geomechanical response to CO 2 injection. The modeling outcome suggests that shear failure and faults/thrusts reactivation can occur much before attaining the hydraulic fracturing pressure, hence representing two major constraints for a safe and permanent containment.

Description: [1]  We present numerical subduction models to investigate overriding plate deformation at subduction zones. All models show forearc shortening, resulting predominantly from shear stresses at the subduction zone interface and opposite-sense mantle shear stresses at the base of the forearc lithosphere. Models dominated by backarc extension show that it results from trench-normal positive velocity gradients in the mantle below the overriding plate. Such gradients result from toroidal mantle flow induced by slab rollback, with velocities below the leading part of the backarc faster than the overriding plate velocity. The velocity gradients induce basal shear stresses that increase trenchward and cause trenchward overriding plate motion at a velocity ( v OP⊥ ) whose spatial average is below the trench retreat velocity ( v T⊥ ). The combination of basal shear stresses and average v OP⊥ 〈 v T⊥ causes trench-normal deviatoric tension in the backarc and backarc extension. Models dominated by backarc shortening show that it results from a relatively immobile subduction hinge and trenchward overriding plate motion driven by poloidal mantle flow. The poloidal mantle flow is induced by downdip slab sinking and causes the average v OP⊥ 〉 v T⊥ . This results in trench-normal deviatoric compression and shortening in the leading part of the overriding plate as it collides with the subduction hinge. Ultimately, the geodynamic models demonstrate that backarc extension is favored for narrow slabs and near lateral slab edges, and is driven by rollback induced toroidal mantle flow, while backarc shortening is favored for the center of wide slabs, and is driven by poloidal mantle flow resulting from downdip slab motion.

Description: [1]  The Gamburtsev Subglacial Mountains (GSM), located near the center of East Antarctica, remain one of the most enigmatic mountain ranges on earth. A lack of direct geologic samples renders theirtectonic history almost totally unconstrained. We utilize teleseismic Rayleigh wave data from a two-year deployment of broadband seismic stations across the region to image shear velocity structure and analyze the lithospheric age of the GSM and surrounding regions. We solve for 2-D phase velocitiesand invert these results for 3-D shear velocity structure. We perform a Monte Carlo simulation to improve constraints of crustal thickness and shear velocity structure.Beneath the core of the GSM, we find crustal thickness in excess of 55 km.Mantle shear velocities remain faster than global average models to a depth of approximately 250 km, indicating a thick lithospheric root. Thinner crust and slower upper mantle velocities are observed beneath the Lambert Rift System and the Polar Subglacial Basin.When compared with phase velocity curves corresponding to specific tectonothermal ages elsewhere in the world, average phase velocity results for the GSM are consistent with regions of Archean – Paleoproterozoic origin. Combined with radiometric ages of detrital zircons found offshore, these results indicate a region of old crust that has undergone repeated periods of uplift and erosion, most recently during the Mesozoic breakup of Gondwana. Lower crustal seismic velocities imply a moderately dense lower crust beneath the core of the GSM, but with lower density than suggested by recent gravity models.

Description: [1]  This is a second paper in a study of statistical identification and classification of earthquake clusters using a relocated catalog of 1981–2011 seismicity in southern California and synthetic catalogs produced by the ETAS model. Here we focus on classification of event families – statistically significant clusters comprised of foreshocks , mainshocks and aftershocks – that are detected with the methodology discussed in part I of the study. The families are analyzed using their representation as time oriented tree graphs. The results (i) demonstrate that the clustering associated with the largest earthquakes, m  〉 7, is statistically different from that of small-to-medium earthquakes; (ii) establish the existence of two dominant types of small-to-medium magnitude earthquake families– burst-like and swarm-like sequences – and a variety of intermediate cluster forms obtained as a mixture of the two dominant types; (iii) suggest a simple new quantitative measure for identifying the cluster type based on its topological structure; (iv) demonstrate systematic spatial variability of the cluster characteristics on a scale of tens of kilometers in relation to heat flow and other properties governing the effective viscosity of a region; and (v) establish correlation between the family topological structure and a dozen of metric properties traditionally considered in the literature (number of aftershocks, duration, spatial properties, b -value, parameters of Omori-Utsu and Båth law, etc .). The burst-like clusters likely reflect highly-brittle failures in relatively cold regions, while the swarm-like clusters are likely associated with mixed brittle-ductile failures in regions with relatively high temperature and/or fluid content. The results of this and paper I may be used to develop improved region-specific hazard estimates and earthquake forecasts.

Description: The Humboldt Current System (HCS) sustains the world′s largest small pelagic fishery. While a cooling of this system has been observed during recent decades, there is debate about the potential impacts of rising atmospheric CO 2 concentrations on upwelling dynamics and productivity. Recent studies suggest that under increased atmospheric CO 2 scenarios the oceanic stratification may strongly increase and upwelling-favorable winds may remain nearly constant off Peru and increase off Chile. Here we investigate the impact of such climatic conditions on egg and larval dispersal phases, a key stage of small pelagic fish reproduction. We used larval retention rate in a predefined nursery area to provide a proxy for the recruitment level. Numerical experiments are based on hydrodynamics downscaled to the HCS from global simulations forced by pre-industrial (PI), 2xCO 2 and 4xCO 2 scenarios. A biogeochemical model is applied to the PI and 4xCO 2 scenarios in order to define a time-variable nursery area where larval survival is optimum. We test two distinct values of the oxycline depth that limits larval vertical distribution: one corresponding to the present-day situation and the other corresponding to a shallower oxycline potentially produced by climate change. It appeared that larval retention over the continental shelf increases with enhanced stratification due to regional warming. However, this increase in retention is largely compensated for by a decrease of the nursery area and the shoaling of the oxycline. The underlying dynamics are explained by a combination of stratification effects and mesoscale activity changes. Our results therefore show that future climate change may significantly reduce fish capacity in the HCS with strong ecological, economic and social consequences. © 2013 Blackwell Publishing Ltd

Description: Adoption of reduced-impact logging (RIL) methods could reduce CO 2 emissions by 30-50% across at least 20% of remaining tropical forests. We developed two cost effective and robust indices for comparing the climate benefits (reduced CO 2 emissions) due to RIL. The indices correct for variability in the volume of commercial timber among concessions. We determined that a correction for variability in terrain slope was not needed. We found that concessions certified by the Forest Stewardship Council (FSC, N=3), when compared with non-certified concessions (N=6), did not have lower overall CO 2 emissions from logging activity (felling, skidding, and hauling). On the other hand, FSC certified concessions did have lower emissions from one type of logging impact (skidding), and we found evidence of a range of improved practices using other field metrics. One explanation for these results may be that FSC criteria and indicators, and associated RIL practices, were not designed to achieve overall emissions reductions. Also, commonly used field metrics are not reliable proxies for overall logging emissions performance. Further, the simple distinction between certified and non-certified concessions does not fully represent the complex history of investments in improved logging practices. To clarify the relationship between RIL and emissions reductions, we propose the more explicit term “RIL-C” to refer to the sub-set of RIL practices that can be defined by quantified thresholds and that result in measurable emissions reductions. If tropical forest certification is to be linked with CO 2 emissions reductions, certification standards need to explicitly require RIL-C practices. This article is protected by copyright. All rights reserved.

Description: Sequestration of atmospheric carbon (C) in soils through improved management of forest and agricultural land is considered to have high potential for global CO 2 mitigation. However, the potential of soils to sequester soil organic carbon (SOC) in a stable form, which is limited by the stabilization of SOC against microbial mineralization, is largely unknown. In this study, we estimated the C sequestration potential of soils in southeast Germany by calculating the potential SOC saturation of silt and clay particles according to Hassink (1997) on the basis of 516 soil profiles. The determination of the current SOC content of silt and clay fractions for major soil units and land uses allowed an estimation of the C saturation deficit corresponding to the long-term C sequestration potential. The results showed that cropland soils have a low level of C saturation of around 50% and could store considerable amounts of additional SOC. A relatively high C sequestration potential was also determined for grassland soils. In contrast, forest soils had a low C sequestration potential as they were almost C saturated. A high proportion of sites with a high degree of apparent oversaturation revealed that in acidic, coarse-textured soils the relation to silt and clay is not suitable to estimate the stable C saturation. A strong correlation of the C saturation deficit with temperature and precipitation allowed a spatial estimation of the C sequestration potential for Bavaria. In total, about 395 Mt CO 2 -equivalents could theoretically be stored in A horizons of cultivated soils – four times the annual emission of greenhouse gases in Bavaria. Although achieving the entire estimated C storage capacity is unrealistic, improved management of cultivated land could contribute significantly to CO 2 mitigation. Moreover, increasing SOC stocks have additional benefits with respect to enhanced soil fertility and agricultural productivity. This article is protected by copyright. All rights reserved.

Description: [1]  Paleoseismologic data from the southern Panamint Valley fault (PVF) reveal evidence of at least four surface ruptures during late Holocene time (0.33-0.48 ka, 0.9-3.0 ka, 3.3-3.6 ka and 〉4.1 ka). These paleo-earthquake ages indicate that the southern PVF has ruptured at least once and possibly twice during the ongoing (≤1.5 ka) seismic cluster in the Mojave section of the eastern California shear zone (ECSZ). The most recent event (MRE) on the PVF is also similar in age to the 1872 Owens Valley earthquake and the geomorphically youthful MRE on the Death Valley fault. The timing of the three oldest events at our site shows that the PVF ruptured at least once and possibly thrice during the well-defined 2-5 ka seismic lull in the Mojave section of the ECSZ. Interestingly, the 3.3-3.6 ka age of Event 3 overlaps with the 3.3-3.8 ka age of the penultimate (i.e., pre-1872) rupture on the central Owens Valley fault. These new PVF data support the notion that earthquake occurrence in the ECSZ may be spatially and temporally complex, with earthquake clusters occurring in different regions at different times. Coulomb Failure Function modeling of the Panamint Valley and Garlock faults reveals significant stress interactions between these two faults that may influence future earthquake occurrence. Specifically, our models suggest a possible rupture sequence whereby an event on the southern Panamint Valley fault can lead to the potential triggering of an event on the Garlock fault, which in turn could trigger the Mojave section of the San Andreas Fault.

Description: [1]  The lithospheric structure of the Indian plate has been investigated using converted wave techniques (P and S receiver functions, RFs) and a novel stacking analysis technique (without using deconvolution) applied to a large seismological dataset from permanent and temporary broadband seismic stations. We observe coherent energy from at least two seismic discontinuities i.e. the crust-mantle (Moho) and lithosphere-asthenosphere boundary (LAB) in the uppermost mantle. Here, we provide a novel seismic image of the Indian lithosphere showing definitive evidence of its flexure, which is interpreted to be primarily caused by the hard collision at ~55My resulting in the world's highest mountain chain - the Himalayas and the Tibetan plateau. Results from geoidal and gravity studies do suggest post-collisional flexuring of the Indian plate; however, the flexure lacks observational constraints. The observed wavelength of the flex is ~1000 km with the thickness of the Indian plate varying from ~70 km to 140 km; such a low value for a continent implies that the Indian plate has been reworked in the past. The plate deepens in the Himalayan region to a depth of ~150 km. Further, the converted phases are interpreted to be resulting from the bottom of the lithosphere. We clearly demonstrate that these are distinct and different from the mid-lithospheric discontinuity. For a large number of stations, the MLD and LAB are clearly separated in depth. Our observations suggest that the Archaean lithosphere is no longer intact and is prone to deformation.

Description: Drought affects more people than any other natural disaster but there is little understanding of how ecosystems react to droughts. This study jointly analyzed spatio-temporal changes of drought patterns with vegetation phenology and productivity changes between 1999 and 2010 in major European bioclimatic zones. The Standardized Precipitation and Evapotranspiration Index (SPEI) was used as drought indicator whereas changes in growing season length and vegetation productivity were assessed using remote sensing time-series of NDVI (Normalized Difference Vegetation Index). Drought spatio-temporal variability was analyzed using a Principal Component Analysis, leading to the identification of four major drought events between 1999 and 2010 in Europe. Correspondence Analysis showed that at the continental scale the productivity and phenology reacted differently to the identified drought events depending on ecosystem and land cover. Northern and Mediterranean ecosystems proved to be more resilient to droughts in terms of vegetation phenology and productivity developments. Western Atlantic regions and Eastern Europe showed strong agglomerations of decreased productivity and shorter vegetation growing season length, indicating that these ecosystems did not buffer the effects of drought well. In a climate change perspective, increase in drought frequency or intensity may result in larger impacts over these ecosystems, thus management and adaptation strategies should be strengthened in these areas of concerns. This article is protected by copyright. All rights reserved.

Description: ‘Humans are now the most significant driver of global change, propelling the planet into a new geological epoch, the Anthropocene’. This landmark statement from the Stockholm Memorandum (2011) is supported by an overwhelming consensus in the scientific literature (Cook et al ., 2013). It is crucial to acknowledge, however, that several of Earth's ecosystems are still little affected by direct human activity, and appropriate conservation measures are fully feasible and should be enforced accordingly (Caro et al ., 2012). Arctic marine ecosystems belong to this category. This article is protected by copyright. All rights reserved.