The Laguna del Maule volcanic field is a large rhyolitic magmatic system in the Chilean Andes, which has exhibited frequent eruptions during the past 20 ka. Rapid surface uplift (>20 cm/year) has been observed since 2007 accompanied by localized earthquake swarms and microgravity changes, indicating the inflating magma reservoir may interact with a preexisting weak zone (i.e., Troncoso fault). In this investigation, we model the magma reservoir by data assimilation with Interferometric Synthetic Aperture Radar data. The reservoir geometry is comparable to the magma body inferred by seismic tomography, magnetotelluric, and gravity studies. The models also suggest that a weak zone, which has little effect on surface displacement, is important as a fluid transport channel to promote earthquakes and microgravity changes. In particular, concentrated dilatancy within the weak zone facilitates the microfracture formation during reservoir inflation. High-pressure fluid can inject into the weak zone from the magma reservoir to trigger earthquakes and further migrate upward to create positive gravity changes by occupying unsaturated storages. The pore pressure will then decrease, halting the seismicity swarm until the next cycle. This "hydrofracturing" process may release some accumulated stress along the magma reservoir delaying an eventual eruption in turn. Besides, the resultant models are propagated forward in time to evaluate potential stress trajectories for future unrest.

In volcano gravimetry, when analyzing residual spatiotemporal (time-lapse) gravity changes, the accurate deformation-induced topographic effect (DITE) should be used to account for the gravitational effect of surface deformation. Numerical realization of DITE requires the deformation field available in grid form. We compute the accurate DITE correction for gravity changes observed at the Laguna del Maule volcanic field in Chile over three nearly annual periods spanning 2013?2016 and compare it numerically with the previously used free-air effect (FAE) correction. We assess the impact of replacing the FAE by DITE on the model source parameters of analytic inversion solutions and apply a new inversion approach based on model exploration and growing source bodies. The new inversion results based on the DITE correction shift the position of the mass intrusion upwards by a few hundred meters and lower the total mass of the migrated fluids to roughly a half, compared to the inversion results based on the local-FAE correction. Our new Growth inversion results indicate that vertical dip slip faults beneath the lake, as well as the Troncoso fault play active roles in hosting migrating liquid. We also show that for the study period, the DITE at Laguna del Maule can be accurately evaluated by the planar Bouguer approximation, which only requires the availability of elevation changes at gravity network benchmarks. We hypothesize that this finding may be generalized to all volcanic areas with flatter or less rugged terrain and may modify interpretations based on the commonly used FAE corrections. ? 2021 Elsevier B.V. All rights reserved.

Improved understanding of the impact of crystal mush rheology on the response of magma chambers to magmatic events is critical for better understanding crustal igneous systems with abundant crystals. In this study, we extend an earlier model by Liao et al. (2018); which considers the mechanical response of a magma chamber with poroelastic crystal mush, by including poroviscoelastic rheology of crystal mush. We find that the coexistence of the two mechanisms of poroelastic diffusion and viscoelastic relaxation causes the magma chamber to react to a magma injection event with more complex time-dependent behaviors. Specifically, we find that the system's short-term evolution is dominated by the poroelastic diffusion process, while its long-term evolution is dominated by the viscoelastic relaxation process. We identify two post-injection timescales that represent these two stages and examine their relation to the material properties of the system. We find that better constraints on the poroelastic diffusion time are more important for the potential interpretation of surface deformation using the model.

The Laguna del Maule volcanic field in Chile has been exhibiting unrest since 2005. New GPS and InSAR data reveal a second episode of accelerated deformation beginning in late 2016 and continuing through May 2020, with an uplift rate > 290 mm/year between 2019 and 2020. To explain the spatial and temporal pattern of deformation, we apply a dynamic model of viscous magma flowing through a conduit into a fluid-filled reservoir surrounded by a heterogeneous, viscoelastic crust. A Monte Carlo procedure optimizes the ellipsoid reservoir geometry and the inlet pressure history. The two episodes of accelerating uplift are each modeled with a pressure increase rate of similar to 9 MPa/year. Since 2016, 0.10 km(3) of magma was injected into the system for a total of 0.37 km(3) since 2005.

Classical mechanisms of volcanic eruptions mostly involve pressure buildup and magma ascent towards the surface'. Such processes produce geophysical and geochemical signals that may be detected and interpreted as eruption precursors(1-3) . On 22 May 2021, Mount Nyiragongo (Democratic Republic of the Congo), an open-vent volcano with a persistent lava lake perched within its summit crater, shook up this interpretation by producing an approximately six-hour-long flank eruption without apparent precursors, followed-rather than preceded-by lateral magma motion into the crust. Here we show that this reversed sequence was most likely initiated by a rupture of the edifice, producing deadly lava flows and triggering a voluminous 25-km-long dyke intrusion. The dyke propagated southwards at very shallow depth (less than 500 m) underneath the cities of Goma (Democratic Republic of the Congo) and Gisenyi (Rwanda), as well as Lake Kivu. Thisvolcanic crisis raises new questions about the mechanisms controlling such eruptions and the possibility of facing substantially more hazardous events, such as effusions within densely urbanized areas, phreato-magmatism or a limnic eruption from the gas-rich Lake Kivu. It also more generally highlights the challenges faced with open-vent volcanoes for monitoring, early detection and risk management when a significant volume of magma is stored close to the surface.

The degree to which elevated CO2 concentrations (e[CO2]) increase the amount of carbon (C) assimilated by vegetation plays a key role in climate change. However, due to the short-term nature of CO2 enrichment experiments and the lack of reconciliation between different ecological scales, the effect of e[CO2] on plant biomass stocks remains a major uncertainty in future climate projections. Here, we review the effect of e[CO2] on plant biomass across multiple levels of ecological organization, scaling from physiological responses to changes in population-, community-, ecosystem-, and global-scale dynamics. We find that evidence for a sustained biomass response to e[CO2] varies across ecological scales, leading to diverging conclusions about the responses of individuals, populations, communities, and ecosystems. While the distinct focus of every scale reveals new mechanisms driving biomass accumulation under e[CO2], none of them provides a full picture of all relevant processes. For example, while physiological evidence suggests a possible long-term basis for increased biomass accumulation under e[CO2] through sustained photosynthetic stimulation, population-scale evidence indicates that a possible e[CO2]-induced increase in mortality rates might potentially outweigh the effect of increases in plant growth rates on biomass levels. Evidence at the global scale may indicate that e[CO2] has contributed to increased biomass cover over recent decades, but due to the difficulty to disentangle the effect of e[CO2] from a variety of climatic and land-use-related drivers of plant biomass stocks, it remains unclear whether nutrient limitations or other ecological mechanisms operating at finer scales will dampen the e[CO2] effect over time. By exploring these discrepancies, we identify key research gaps in our understanding of the effect of e[CO2] on plant biomass and highlight the need to integrate knowledge across scales of ecological organization so that large-scale modeling can represent the finer-scale mechanisms needed to constrain our understanding of future terrestrial C storage.

The origin of major volatiles nitrogen, carbon, hydrogen, and sulfur in planets is critical for understanding planetary accretion, differentiation, and habitability. However, the detailed process for the origin of Earth's major volatiles remains unresolved. Nitrogen shows large isotopic fractionations among geochemical and cosmochemical reservoirs, which could be used to place tight constraints on Earth's volatile accretion process. Here we experimentally determine N-partitioning and -isotopic fractionation between planetary cores and silicate mantles. We show that the core/mantle N-isotopic fractionation factors, ranging from -4 parts per thousand to +10 parts per thousand, are strongly controlled by oxygen fugacity, and the core/mantle N-partitioning is a multi-function of oxygen fugacity, temperature, pressure, and compositions of the core and mantle. After applying N-partitioning and -isotopic fractionation in a planetary accretion and core-mantle differentiation model, we find that the N-budget and -isotopic composition of Earth's crust plus atmosphere, silicate mantle, and the mantle source of oceanic island basalts are best explained by Earth's early accretion of enstatite chondrite-like impactors, followed by accretion of increasingly oxidized impactors and minimal CI chondrite-like materials before and during the Moon-forming giant impact. Such a heterogeneous accretion process can also explain the carbon-hydrogen-sulfur budget in the bulk silicate Earth. The Earth may thus have acquired its major volatile inventory heterogeneously during the main accretion phase.

The FourStar infrared camera is a 1.0-2.5 mu m (JHK(s)) near infrared camera for the Magellan Baade 6.5m telescope at Las Campanas Observatory (Chile). It is being built by Carnegie Observatories and the Instrument Development Group and is scheduled for completion in 2009. The instrument uses four Teledyne HAWAII-2RG arrays that produce a 10.9'x 10.9' field of view. The outstanding seeing at the Las Campanas site coupled with FourStar's high sensitivity and large field of view will enable many new survey and targeted science programs.

We report the detection of the first two planets from the N2K Doppler planet search program at the Magellan telescopes. The first planet has a mass of M sin i = 4.96 M-Jup and orbits the G3 IV star HD154672 with an orbital period of 163.9 days. The second planet orbits the F7 V star HD205739 with an orbital period of 279.8 days and has a mass of M sin i = 1.37 M-Jup. Both planets are in eccentric orbits, with eccentricities e = 0.61 and e = 0.27, respectively. Both stars are metal rich and appear to be chromospherically inactive, based on inspection of their Ca II H and K lines. Finally, the best Keplerian model fit to HD205739b shows a trend of 0.0649 m s(-1) day(-1), suggesting the presence of an additional outer body in that system.

Exploring the origin of Ly alpha nebulae ("blobs") at high redshift requires measurements of their gas kinematics that are impossible with only the resonant, optically thick Ly alpha line. To define gas motions relative to the systemic velocity of the blob, the Ly alpha line must be compared with an optically thin line like H alpha lambda 6563, which is not much altered by radiative transfer effects and is more concentrated about the galaxies embedded in the nebula's core. We obtain optical and near-IR (NIR) spectra of the two brightest Ly alpha blobs (CDFS-LAB01 and CDFS-LAB02) from the Yang et al. sample using the Magellan/Magellan Echellette Spectrograph optical and Very Large Telescope/SINFONI NIR spectrographs. Both the Ly alpha and H alpha lines confirm that these blobs lie at the survey redshift, z similar to 2.3. Within each blob, we detect several H alpha sources, which roughly correspond to galaxies seen in Hubble Space Telescope rest-frame UV images. The H alpha detections show that these galaxies have large internal velocity dispersions (sigma(upsilon) = 130-190 km s(-1)) and that, in the one system (LAB01), where we can reliably extract profiles for two H alpha sources, their velocity difference is Delta upsilon similar to 440 km s(-1). The presence of multiple galaxies within the blobs, and those galaxies' large velocity dispersions and large relative motion, is consistent with our previous finding that Ly alpha blobs inhabit massive dark matter halos that will evolve into those typical of present-day rich clusters and that the embedded galaxies may eventually become brightest cluster galaxies. To determine whether the gas near the embedded galaxies is predominantly infalling or outflowing, we compare the Ly alpha and H alpha line centers, finding that Ly alpha is not offset (Delta upsilon(Ly alpha) = +0 km s(-1)) in LAB01 and redshifted by only +230 km s(-1) in LAB02. These offsets are small compared to those of Lyman break galaxies, which average +450 km s(-1) and extend to about +700 km s(-1). In LAB02, we detect C II lambda 1334 and Si II lambda 1526 absorption lines, whose blueward shifts of similar to 200 km s(-1) are consistent with the small outflow implied by the redward shift of Ly alpha. We test and rule out the simplest infall models and those outflow models with super/hyperwinds, which require large outflow velocities. Because of the unknown geometry of the gas distribution and the possibility of multiple sources of Ly alpha emission embedded in the blobs, a larger sample and more sophisticated models are required to test more complex or a wider range of infall and outflow scenarios.