Elasticity of MgO to 130 GPa: Implications for lower mantle mineralogy

The aggregate shear wave velocities of MgO (periclase) have been determined throughout Earth's lower mantle pressure regime approaching 130 GPa using Brillouin spectroscopy in conjunction with synchrotron X-ray diffraction technique in a diamond anvil cell apparatus. We found that the extrapolations of the high-pressure shear wave velocities and shear moduli to ambient pressure are highly consistent with earlier studies. However, the measurements over a wide pressure range revealed that the pressure derivative of the shear modulus (dG/dP = G₀′) of MgO is 1.92(2), which is distinctly lower than that of previous lower-pressure experiments. Compared with the previous results on (Mg,Fe)O ferropericlase, there is no clear correlation between iron content and G₀′. We calculate that the shear wave velocity profile of lower mantle along the adiabatic geotherm applied by the lower G₀′ value of periclase can remarkably well reproduce the global seismological 1-D velocity profile model with uniform composition model. The best-fitting result indicates the possibility of a lower mantle mineralogy with ~ 92 vol.% silicate perovskite phase, implying that the bulk composition of lower mantle is likely not to be pyrolitic but more chondritic. The present acoustic measurements performed over the large pressure range have thus led us to a better understanding of compositional model of the Earth's lower mantle.     

Sound velocity measurements in dhcp-FeH up to 70 GPa with inelastic X-ray scattering: Implications for the composition of the Earth's core

We have determined the density evolution of the sound velocity of dhcp-FeH_x (x ≈ 1) up to 70 GPa at room temperature, by inelastic X-ray scattering and by X-ray diffraction. We find that the variation of V_P with density is different for the ferromagnetic and nonmagnetic dhcp-FeH_x, and that only nonmagnetic dhcp-FeH_x follows Birch's law. Combining our results with Birch's law for iron and assuming an ideal two-component mixing model, we obtain an upper bound of the hydrogen content in the Earth's inner core, 0.23(6) wt.% H, corresponding to FeH_0.13(3). The iron alloy with 0.23(6) wt.% H can satisfy the density, and compressional and shear sound velocities of the PREM inner core, assuming that there are no other light elements in the inner core.     

Viscosity of peridotite liquid up to 13 GPa: Implications for magma ocean viscosities

The viscosity of synthetic peridotite liquid has been investigated at high pressures using in-situ falling sphere viscometry by combining a multi-anvil technique with synchrotron radiation. We used a newly designed capsule containing a small recessed reservoir outside of the hot spot of the heater, in which a viscosity marker sphere is embedded in a forsterite + enstatite mixture having a higher solidus temperature than the peridotite. This experimental setup prevents spheres from falling before a stable temperature above the liquidus is established and thus avoids difficulties in evaluating viscosities from velocities of spheres falling through a partially molten sample.

Experiments have been performed between 2.8 and 13 GPa at temperatures ranging from 2043 to 2523 K. Measured viscosities range from 0.019 (± 0.004) to 0.13 (± 0.02) Pa s. At constant temperature, viscosity increases with increasing pressure up to  8.5 GPa but then decreases between  8.5 and 13 GPa. The change in the pressure dependence of viscosity is likely associated with structural changes of the liquid that occur upon compression. By combining our results with recently published 0.1 MPa peridotite liquid viscosities [D.B. Dingwell, C. Courtial, D. Giordano, A. Nichols, Viscosity of peridotite liquid, Earth Planet. Sci. Lett. 226 (2004) 127–138.], the experimental data can be described by a non-Arrhenian, empirical Vogel-Fulcher-Tamman equation, which has been modified by adding a term to account for the observed pressure dependence of viscosity. This equation reproduces measured viscosities to within 0.08 log₁₀-units on average. We use this model to calculate viscosities of a peridotitic magma ocean along a liquid adiabat to a depth of  400 km and discuss possible effects on viscosity at greater pressures and temperatures than experimentally investigated.     

Thermal equation of state of magnesite to 32 GPa and 2073 K

Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg_0.975Fe_0.015Mn_0.006Ca_0.004CO₃) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K₀ = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V₀ = 279.55 ± 0.02 Å³. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a₀ + a₁T with a₀ = 4.03 (7) × 10⁻⁵ K⁻¹ and a₁ = 0.49 (10) × 10⁻⁸ K⁻². The Anderson–Grüneisen parameter is estimated to be δ_T = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (K_Tc = 47 ± 1 GPa) than the a-axis (K_Tc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a₀ = 6.8 (2) × 10⁻⁵ K⁻¹ and a₁ = 2.2 (4) × 10⁻⁸ K⁻²) is greater than that of the a-axis (a₀ = 2.7 (4) × 10⁻⁵ K⁻¹ and a₁ = −0.2 (2) × 10⁻⁸ K⁻²). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.     

Phase relations in hydrous MORB at 18-28 GPa: implications for heterogeneity of the lower mantle

The phase relations in hydrous and anhydrous mid-ocean ridge basalt were determined at pressures of 18-28 GPa. Liquidus phase relations in hydrous and anhydrous MORB are different. Garnet is the liquidus phase at pressures below 21 GPa, Ca-Al (CAS) phase and stishovite are the liquidus phases at pressures of 22-27 GPa, and stishovite and Ca-perovskite are the liquidus phases above 27 GPa, whereas Ca-perovskite is a liquidus phase of anhydrous MORB at pressures above 23 GPa. Under subsolidus conditions, we have found that in the hydrous MORB system the stability fields of Al-bearing perovskite and Na-Al (NAL) phase might shift to lower pressure by about 1.5 GPa compared to the dry MORB system. This shift could be explained by oxidation of a garnet-bearing assemblage by hydrous fluid and formation of Fe³⁺-bearing aluminous perovskite at lower pressures relative to the anhydrous system and/or differences in water solubility of the phases existing in perovskite-bearing assemblages. Our data indicate that hydrous basaltic crust remains denser than peridotite along the geotherm of a subducting slab, i.e. there is no density crossover between peridotite and basalt. Therefore, in slabs going through the 660 km discontinuity, basalt would gravitationally sink into the lower mantle under relatively hydrous conditions. The delamination of former basaltic crust near the 660 km discontinuity might be possible under relatively dry conditions of subduction. There are no stable highly hydrous phases in MORB above 10 GPa even at lower temperatures corresponding to subducting slabs. Therefore, MORB cannot be an important carrier of water to the deep Earth interior. However, it can be constantly supplied by water-bearing fluid from the underlying peridotite part of the descending slab. Thus, it is plausible that water can control subduction of the oceanic crust into the lower mantle.     

Phase relations of iron and iron–nickel alloys up to 300 GPa: Implications for composition and structure of the Earth's inner core

We have investigated the phase relations of iron and iron–nickel alloys with 18 to 50 wt.% Ni up to over 300 GPa using a laser-heated diamond-anvil cell. The synchrotron X-ray diffraction measurements show the wide stability of hcp-iron up to 301 GPa and 2000 K and 319 GPa and 300 K without phase transition to dhcp, orthorhombic, or bcc phases. On the other hand, the incorporation of nickel has a remarkable effect on expanding the stability field of fcc phase. The geometry of the temperature–composition phase diagram of iron–nickel alloys suggests that the hcp–fcc–liquid triple point is located at 10 to 20 wt.% Ni at the pressure of the inner core boundary. The fcc phase could crystallize depending on the nickel and silicon contents in the Earth's core, both of which are fcc stabilizer.     

Volatile and trace elements in basaltic glasses from Samoa: Implications for water distribution in the mantle

We report volatile (H₂O, CO₂, F, S, Cl) and trace element data for submarine alkalic basalt glasses from the three youngest Samoan volcanoes, Ta'u, Malumalu and Vailulu'u. Most samples are visibly sulfide saturated, so have likely lost some S during fractionation. Cl / K ratios (0.04–0.15) extend to higher values than pristine MORBs, but are suspected to be partly due to source differences since Cl / K roughly varies as a function of ⁸⁷Sr / ⁸⁶Sr. There are no resolvable differences in the relative enrichment of F among sources, and compatibility of F during mantle melting is established to be nearly identical to Nd. Shallow degassing has affected CO₂ in all samples, and H₂O only in the most shallowly erupted samples from Vailulu'u. Absolute water contents are high for Samoa (0.63–1.50 wt.%), but relative enrichment of water compared to equally incompatible trace elements (Ce, La) is low and falls entirely below normal MORB values. H₂O / Ce (58–157) and H₂O / La (120–350) correlate inversely with ⁸⁷Sr / ⁸⁶Sr compositions (0.7045–0.7089). This leads us to believe that, because of very fast diffusion of hydrogen in mantle minerals, recycled lithospheric material with high initial water and trace element content will lose water to the drier ambient mantle during storage within the inner Earth. The net result is the counter-intuitive appearance of greater dehydration with greater mantle enrichment. We expect that subducted slabs will experience a two-stage dehydration history, first within subduction zones and then in the ambient mantle during long-term convective mixing.     

Phase transformations in the earth's mantle and subducting slabs: Implications for their compositions, seismic velocity and density structures and dynamics

Abstract Phase transformations in model mantle compositions and those in subducting slabs have been reviewed to a depth of 800 km on the basis of recent high-pressure experimental data. Seismic velocity and density profiles in these compositions have also been calculated using these and other mineral physics data. The nature of the seismic velocity and density profiles calculated for a pyrolite composition was found to generally agree with those determined by seismic observations (e.g. PREM). The locations of the seismic discontinuities at 400 and 670 km correspond almost exactly to the depths where the transformations of the olivine component to denser phases take place. Moreover, the steep gradients in the seismic velocity/density profiles observed between these depths are qualitatively consistent with those expected from the successive transformations in the complementary pyroxene-garnet component in the pyrolite composition. Further, the calculated seismic velocity and density values agree well with those observed in the upper mantle and mantle transition region within the uncertainties attached to these calculations and observations. Pyrolite or peridotite compositions are thus most likely to represent the composition of the mantle above 670 km depth, although some degrees of chemical heterogeneity may exist in the transition region. The observed sharp discontinuous increases of seismic velocities and density at this depth may be attributed either to the phase transformation to a perovskite-bearing assemblage in pyrolite or to chemical composition changes. Density profiles in subducted slabs have been calculated along adequate geotherms assuming that the slabs are composed of the former oceanic crust underlain by a thicker harzburgitic layer. It is shown that the former oceanic crust is substantially less dense than the surrounding pyrolite mantle at depths below 670 km, while it is denser than pyrolite in the upper mantle and the transition region. The subducted former oceanic crust may be trapped in this region, forming a geochemically enriched layer at the upper mantle-lower mantle boundary. Thick and cool slabs may penetrate into the lower mantle, but the chemically derived buoyancy may result in strong deformation and formation of megalith structures around the 670 km seismic discontinuity. These structures are consistent with those detected by recent seismic tomography studies for subduction zones.     

Experimental investigation of phase transformations of olivine and enstatite at the lower part of the mantle transition zone: Implications for structure of the 660 km seismic discontinuity

High-pressure polymorphs of olivine and enstatite are major constituent minerals in the mantle transition zone(MTZ).The phase transformations of olivine and enstatite at pressure and temperature conditions corresponding to the lower part of the MTZ are import for understanding the nature of the 660 km seismic discontinuity.In this study,we determine phase transformations of olivine(MgSi2O4) and enstatite(MgSiO3) systematiclly at pressures between 21.3 and 24.4 GPa and at a constant temperature of 1600℃.The most profound discrepancy between olivine and enstatite phase transformation is the occurency of perovskite.In the olivine system,the post-spinel transformation occures at 23.8 GPa,corresponding to a depth of 660 km.In contrast,perovskite appears at 23 GPa(640 km) in the enstatite system.The ~1 GPa gap could explain the uplifting and/or splitting of the 660 km seismic discountinuity under eastern China.     

Tetrahedral occupancy of ferric iron in (Mg,Fe)O: Implications for point defects in the Earth's lower mantle

We investigated the concentration and site occupation of ferric iron (Fe³⁺) in (Mg,Fe)O to understand the influence of point defects on transport properties such as atomic diffusion, electrical conductivity and viscosity. We conducted Mössbauer spectroscopy of (Mg_0.8Fe_0.2)O single crystals synthesized at temperatures from 1673 to 2273 K and pressures from 5 to 15 GPa with Re–ReO₂ and Mo–MoO₂ oxygen fugacity buffers. The isomer shift of the Mössbauer spectra suggests that Fe³⁺ occupies mostly the tetrahedral site at reduced conditions and both the octahedral and tetrahedral sites at oxidized conditions. We formulate a thermodynamic model of point defect dissolution in (Mg,Fe)O which suggests that unassociated tetrahedral Fe³⁺ is more stable than unassociated octahedral Fe³⁺ at high-pressure and low oxygen fugacity due to the effect of configurational entropy. The pressure dependence of Fe³⁺ concentration indicates a change in the dominant site occupancy of Fe³⁺: (1) Fe³⁺ in the tetrahedral site, (2) Fe³⁺ in the octahedral site, and (3) defect clusters of Fe³⁺ and cation vacancy, in the order of increasing oxygen fugacity and decreasing pressure. This is in reasonable agreement with previously reported experiments on Fe³⁺ concentration, Mg–Fe interdiffusivity and electrical conductivity. We consider it plausible that (Mg,Fe)O accommodates Fe³⁺ in the tetrahedral site down to the lower mantle. Based on our results and available experimental data, we discuss the solubility competition between Fe³⁺ and protons (H⁺), and its implications for transport properties in the lower mantle.     

Apatite and phosphorus in mantle source regions: An experimental study of apatite/melt equilibria at pressures to 25 kbar

The solubility of fluorapatite in a wide variety of basic magmatic liquids was experimentally determined over a range of upper mantle P-T conditions (8–25 kbar, 1275–1350°C). Fluorapatite is stable over the entire range of conditions investigated, but its solubility in melts is variable, depending negatively on SiO₂ content of the melt and positively upon temperature, with relatively little sensitivity to pressure above 8 kbar. At upper mantle pressures and a temperature of 1250°C, molten basalt (50% SiO₂) will dissolve 3–4 wt.% P₂O₅ before saturation in apatite is reached. For a magma 100°C cooler or containing 10 wt.% more SiO₂, apatite saturation occurs at less than 2 wt.% dissolved P₂O₅. The observed high solubility of apatite in basic magmas at their normal near-liquidus temperatures virtually precludes the occurrence of residual apatite in mantle source regions. If relatively low-temperature melting conditions prevail (e.g., 1100°C), as might be possible in H₂O-bearing regions of the upper mantle, apatite could remain in the residue, but only in amounts too small to have significant effects on the rare earth patterns of the liquids.Because of the high solubility of apatite in basic magmas, phosphorus can be confidently treated as an incompatible element in peridotite melting models. Such models, in combination with observed characteristics of basic lavas, indicate that the upper mantle contains ～200 ppm of phosphorus, much less than the chondritic abundance of ～900 ppm.     

The last 5000 years of activity at Sete Cidades volcano (São Miguel Island,Azores): Implications for hazard assessment

Sete Cidades is a central volcano with a summit caldera at the western end of São Miguel Island, Azores. Its stratigraphy comprises two main geological groups: the Inferior Group, the units of which date from more than 200 000 years ago through to 36 000 years before present, consisting of thick lava flows and subaerial volcaniclastic deposits that built the base of the central volcano; and the Superior Group which comprises all the activity from the last 36 000 years, including pumice and scoria fallout and PDC deposits with minor lava flows. The volcanostratigraphy is divided into six main formations — Risco, Ajuda, Bretanha, Lombas, Santa Bárbara and Lagoas, each defined by different activity phases in the volcano's evolution.     

Homologous temperature of olivine: Implications for creep of the upper mantle and fabric transitions in olivine

The homologues temperature of a crystalline material is defined as T/T_m, where T is temperature and T_m is the melting(solidus) temperature in Kelvin. It has been widely used to compare the creep strength of crystalline materials. The melting temperature of olivine system,(Mg,Fe)_2SiO_4, decreases with increasing iron content and water content, and increases with confining pressure. At high pressure, phase transition will lead to a sharp change in the melting curve of olivine. After calibrating previous melting experiments on fayalite(Fe_2SiO_4), the triple point of fayalite-Fe_2SiO_4 spinel-liquid is determined to be at 6.4 GPa and 1793 K. Using the generalized means, the solidus and liquidus of dry olivine are described as a function of iron content and pressure up to 6.4 GPa. The change of T/T_m of olivine with depth allows us to compare the strength of the upper mantle with different thermal states and olivine composition. The transition from semi-brittle to ductile deformation in the upper mantle occurs at a depth where T/T_m of olivine equals 0.5. The lithospheric mantle beneath cratons shows much smaller T/T_m of olivine than orogens and extensional basins until the lithosphere-asthenosphere boundary where T/T_m 0.66, suggesting a stronger lithosphere beneath cratons. In addition, T/T_m is used to analyze deformation experiments on olivine. The results indicate that the effect of water on fabric transitions in olivine is closely related with pressure. The hydrogen-weakening effect and its relationship with T/T_m of olivine need further investigation. Below 6.4 GPa(200 km), T/T_m of olivine controls the transition of dislocation glide from [100] slip to [001] slip. Under the strain rate of 10~(-12)–10~(-15) s~(-1) and low stress in the upper mantle, the [100](010) slip system(A-type fabric) becomes dominant when T/T_m 0.55–0.60. When T/T_m 0.55–0.60, [001] slip is easier and low T/T_m favors the operation of [001](100) slip system(C-type fabric). This is consistent with the widely observed A-type olivine fabric in naturally deformed peridotites, and the C-type olivine fabric in peridotites that experienced deep subduction in ultrahigh-pressure metamorphic terranes. However, the B-type fabric will develop under high stress and relatively low T/T_m. Therefore, the homologues temperature of olivine established a bridge to extrapolate deformation experiments to rheology of the upper mantle. Seismic anisotropy of the upper mantle beneath cratons should be simulated using a four-layer model with the relic A-type fabric in the upper lithospheric mantle, the B-type fabric in the middle layer, the newly formed A- or B-type fabric near the lithosphere-asthenosphere boundary, and the asthenosphere dominated by diffusion creep below the Lehmann discontinuity. Knowledge about transition mechanisms of olivine fabrics is critical for tracing the water distribution and mantle flow from seismic anisotropy.     

Simulating the effects of geologic heterogeneity and transient boundary conditions on streambed temperatures — Implications for temperature-based water flux calculations

Analytical solutions to the one-dimensional heat transport equation for steady-state conditions can provide simple means to quantify groundwater surface water exchange. The errors in exchange flux calculations that are introduced when the underlying assumptions of homogeneous sediments and constant temperature boundary conditions are violated were systematically evaluated in a simulation study. Temperatures in heterogeneous sediments were simulated using a numerical model. Heterogeneity in the sediments was represented by discrete, binary geologic units. High contrasts between the hydraulic conductivities (K) of the geologic units were found to lead to large errors, while the influence of the structural arrangement of the units was smaller. The effects of transient temperature boundary conditions were investigated using an analytical equation. Errors introduced by transient boundary conditions were small for Darcy-velocities > 0.1 m d^− 1 in the period near maximum and minimum annual surface water temperatures. For smaller fluxes, however, errors can be large. Assuming steady-state conditions and vertical flow in homogeneous sediments is acceptable at certain times of the year and for medium to high exchange fluxes, but pronounced geologic heterogeneity can lead to large errors.     

Twenty-five million years of silicic volcanism in the southern central volcanic zone of the Andes: Geochemistry and magma genesis of ignimbrites from 25 to 27 °S, 67 to 72 °W

Silicic volcanism in the Andean Central Volcanic Zone (CVZ) produced one of the world's largest Neogene ignimbrite provinces. The largest and best-known CVZ ignimbrites are located on the Altiplano-Puna plateau north of 24 °S. Their compositions and huge erupted volumes suggest an origin by large-scale crustal melting, and present-day geophysical anomalies in this region suggest still active zones of partial melting in the middle crust. Farther south in the CVZ, the Cerro Galán complex erupted ignimbrites in the late Miocene and Pliocene that are quite similar in volume and composition to those from north of 24 °S and they have a similar origin. However, there are a great many other, smaller ignimbrites in the southern CVZ whose compositions and geodynamic significance are poorly known. These are the subject of this paper.     

Phase relations in the system LiFMgF2 at elevated pressures: Implications for the proposed mixed-oxide zone of the earth's mantle

Solvi and liquidi for various LiFMgF₂ mixtures have been determined at pressures up to 40 kbar by differential-thermal-analysis in a piston-cylinder high-pressure device. The melting curves of pure LiF and MgF₂ were also studied and the initial slopes (dT_m/dP)_{P = 0} were found to be 11.2 and 8.3°C/kbar, respectively. The eutectic composition (LiF)_0.64(MgF₂)_0.36 is independent of pressure to 35 kbar and the eutectic temperature rises approximately 6.3°C per kbar. Initial slopes of 11°C/kbar and 35°C/kbar are inferred for the melting curves of MgO and SiO₂ (stishovite) respectively, on the basis of data for their structural analogue compounds. The observed solid solution of LiF in MgF₂ and other evidence suggest the possibility of solid solution in the system (Mg,Fe)OSiO₂ (stishovite) under mantle conditions which may have important consequences for the elastic properties of a “mixed-oxide” zone of the earth's mantle.     

Deuterium‐excess determination of evaporation to inflow ratios of an alpine lake: Implications for water balance and modeling

The numerous lakes on the Tibetan Plateau play an important role in the regional hydrological cycle and water resources, but systematic observations of the lake water balance are scarce on the Tibetan Plateau. Here, we present a detailed study on the water cycle of Cona Lake, at the headwaters of the Nujiang‐Salween River, based on 3 years (2011–2013) of observations of δ¹⁸O and δ²H, including samples from precipitation, lake water, and outlet surface water. Short‐term atmospheric water vapor was also sampled for isotope analyses. The δ²H–δ¹⁸O relationship in lake water (δ²H = 6.67δ¹⁸O ? 20.37) differed from that of local precipitation (δ²H = 8.29δ¹⁸O + 12.50), and the deuterium excess (d‐excess) in the lake water (?7.5‰) was significantly lower than in local precipitation (10.7‰), indicating an evaporative isotope enrichment in lake water. The ratio of evaporation to inflow (E /I ) of the lake water was calculated using both d‐excess and δ¹⁸O. The E /I ratios of Cona lake ranged from 0.24 to 0.27 during the 3 years. Observations of atmospheric water vapor isotopic composition (δ _A ) improved the accuracy in E /I ratio estimate over a simple precipitation equilibrium model, though a correction factor method provided nearly identical estimates of E /I ratio. The work demonstrates the feasibility of d‐excess in the study of the water cycle for lakes in other regions of the world and provides recommendations on sampling strategies for accurate calculations of E /I ratio.     

Wet surface resistance of forest canopy in monsoon Asia: Implications for eddy‐covariance measurement of evapotranspiration

During the Asian monsoon period, intense precipitation commonly occurs for an extended period in accompaniment with a reduction in solar radiation. This suggests that wet surface evapotranspiration is an important contributor to the total evapotranspiration. Therefore, investigating evapotranspiration over a wet canopy surface is critical to achieve a better understanding of water and energy cycles in Asia. In this study, we estimated surface resistances under wet conditions in a mixed forest influenced by the East Asian monsoon system. We showed that the surface resistance had a non‐negligible magnitude of about 30 sm^?1 even under wet conditions. We also found that the ratio between the actual and potential evapotranspiration depended on the friction velocity regardless of the time of day. Our analyses suggest that this dependency is tightly related to the underestimation of turbulent fluxes by the eddy‐covariance system under wet surface conditions. Together, our findings suggest that the wet surface resistance, although small, should be considered in simulating evapotranspiration because the forest ecosystem is strongly coupled to the overlying atmosphere. This could significantly improve the shortcomings of evapotranspiration measurement and modeling in Asian forest canopies influenced by the monsoon system. Copyright © 2012 John Wiley & Sons, Ltd.