Search

Abstract The generation of cosmic ray particles with ultra‐relativistic energies is a crucial issue in astrophysical and space physics. Shocks have been acknowledged as efficient accelerators in the universe which can produce energetic particles through a variety of mechanisms, such as diffusive shock acceleration (DSA), shock surfing acceleration (SSA) and shock drift acceleration (SDA). Using the data set obtained by the Magnetospheric Multiscale mission, we report an unusually sharp enhancement of energetic electron flux in the upstream region of Earth's bow shock, associated with a bipolar variation of the interplanetary magnetic field. To explain these observations, we propose a two‐step acceleration scenario where electrons trapped in a flux rope (FR) first undergo head‐on collisions and are further accelerated by FR contraction as the FR crosses the shock. Such a scenario can well explain the spacecraft observations. This work can improve our understanding of particle acceleration processes at shocks.

The 2015–2016 Amazon drought was characterized by below-average regional precipitation for an entire year, which distinguishes it from the dry-season-only droughts in 2005 and 2010. Studies of vegetation indices (VIs) derived from optical remote sensing over the Amazonian forests indicated three stages in canopy response during the 2015–2016 drought, with below-average greenness during the onset and end of the drought, and above-average greenness during the intervening months. To date, a satisfactory explanation for this broad temporal pattern has not been found. A better understanding of rainforest behaviors during this unusually long drought should help predict their response to future droughts. We hypothesized that negative VI anomalies could be caused by water and heat stress exceeding the tolerance ranges of the rainforest. To test our hypothesis, based on monthly observations of terrestrial water storage (TWS), land surface temperature (LST), and vapor pressure deficit (VPD) for January 2003 to December 2016, we proposed an approach to categorize regions into two groups: (1) those exceeding normal hydrological and thermal ranges and (2) those within normal ranges. Accordingly, regions exceeding normal ranges during different stages of the 2015–2016 event were delineated. The results showed a gradual southward shift in these regions: from the northeastern Amazon during August to October 2015 to the north–central part during November 2015 to February 2016 and finally to the southern Amazon in July 2016. Over these regions exceeding normal ranges during droughts, negative VI anomalies were expected, irrespective of radiation anomalies. Over the regions within normal ranges, VI anomalies were assumed to respond positively to radiation anomalies, as is expected under normal conditions. We found that our proposed approach can explain more than 70 % of the observed spatiotemporal patterns in VI anomalies during the 2015–2016 drought. These results suggest that our “exceeding normal ranges”-based approach combining (i) water storage, (ii) temperature, and (iii) atmospheric moisture demand drivers can reasonably identify the most likely drought-affected regions at monthly to seasonal timescales. Using observation-based hydrological and thermal condition thresholds can help with interpreting the response of the Amazon rainforest to future drought events.

Abstract Dipolarization fronts (DFs), characterized by sharp increases in the northward magnetic field and usually preceded by magnetic dips, are suggested to play an important role in energy conversion and transport in the magnetotail. It has been documented that strong energy conversion typically develops right at the fronts. Here we present spacecraft observations of electron‐scale energy conversion (EEC) developed inside the dip region ahead of a DF, by using high‐cadence data from the Magnetospheric Multiscale Mission. The EEC, with magnitude comparable to that at the front, is primarily driven by ion current and electron‐scale electric field. The electric field inside the dip is provided by electrostatic waves fed by lower hybrid drift instability, which experiences temporal decaying. Such decaying leads to nonhomogeneity of EEC along the dawn‐dusk direction. These results, uncovering a new channel for DF‐driven energy conversion, can provide important insights into understanding energy transport in the magnetotail.

A whistler wave is a fundamental electromagnetic radiation. In the universe, it transmits information through diverse mediums at audio frequencies. To know the properties of such a radio, one should determine how its frequency (energy) changes with its wavelength (momentum). The relation between them is called the dispersion relation, which encodes essential information on the kinetics. Although the dispersion relation governs the majority of the whistler properties, the experimental determination of their entire dispersion is still a challenge today. Especially in high-temperature mediums, the group velocity dispersion properties of whistlers are still unverified by spacecraft experiments despite their practical importance. Here, for the first time, we experimentally determine the dispersion relation of whistler waves in superhot (T_{e}>5keV) plasmas by the magnetospheric multiscale mission located in a magnetotail plasma sheet. Our result unmasks the characteristics of whistler waves at elevated temperatures and exhibits the evolution of their group velocity dispersion from positive to negative, which agrees well with kinetic theory, opening the door to controlled radio excitation and transport.

Abstract The Earth's magnetopause is an ion‐scale boundary that separates the magnetosphere from the shocked solar wind. At such boundary, energy‐conversion processes frequently occur. Previous studies suggested that such energy conversions are related to the magnetic reconnection process. Here, we report a new mechanism that the coaction between the twisted magnetic structure and lower‐hybrid waves can also drive energy conversion (J⋅ E, J is current density and E is electric field) at the magnetopause boundary layer, by using high‐resolution measurements of the four Magnetospheric Multiscale spacecraft. We find that such energy conversion is efficient (with magnitude up to 20 nW/m3) and is attributed to an intense current filament (j ≈ 3,800 nA/m2) and a fluctuating electric field driven by lower‐hybrid waves. With the help of the First‐Order Taylor Expansion method, we find that the intense current filament is driven by a twisted magnetic structure at electron scale. Our study improves the understanding of energy‐conversion processes at the Earth's magnetopause.

Abstract To improve the efficiency of the discrete unified gas kinetic scheme (DUGKS) in capturing cross-scale flow physics, an adaptive partitioning-based discrete unified gas kinetic scheme (ADUGKS) is developed in this work. The ADUGKS is designed from the discrete characteristic solution to the Boltzmann-BGK equation, which contains the initial distribution function and the local equilibrium state. The initial distribution function contributes to the calculation of free streaming fluxes and the local equilibrium state contributes to the calculation of equilibrium fluxes. When the contribution of the initial distribution function is negative, the local flow field can be regarded as the continuous flow and the Navier–Stokes (N-S) equations can be used to obtain the solution directly. Otherwise, the discrete distribution functions should be updated by the Boltzmann equation to capture the rarefaction effect. Given this, in the ADUGKS, the computational domain is divided into the DUGKS cell and the N-S cell based on the contribution of the initial distribution function to the calculation of free streaming fluxes. In the N-S cell, the local flow field is evolved by solving the N-S equations, while in the DUGKS cell, both the discrete velocity Boltzmann equation and the corresponding macroscopic governing equations are solved by a modified DUGKS. Since more and more cells turn into the N-S cell with the decrease of the Knudsen number, a significant acceleration can be achieved for the ADUGKS in the continuum flow regime as compared with the DUGKS.

Satellite radar backscatter contains unique information on land surface moisture, vegetation features, and surface roughness and has thus been used in a range of Earth science disciplines. However, there is no single global radar data set that has a relatively long wavelength and a decades-long time span. We here provide the first long-term (since 1992), high-resolution (∼8.9 km instead of the commonly used ∼25 km resolution) monthly satellite radar backscatter data set over global land areas, called the long-term, high-resolution scatterometer (LHScat) data set, by fusing signals from the European Remote Sensing satellite (ERS; 1992–2001; C-band; 5.3 GHz), Quick Scatterometer (QSCAT, 1999–2009; Ku-band; 13.4 GHz), and the Advanced SCATterometer (ASCAT; since 2007; C-band; 5.255 GHz). The 6-year data gap between C-band ERS and ASCAT was filled by modelling a substitute C-band signal during 1999–2009 from Ku-band QSCAT signals and climatic information. To this end, we first rescaled the signals from different sensors, pixel by pixel. We then corrected the monthly signal differences between the C-band and the scaled Ku-band signals by modelling the signal differences from climatic variables (i.e. monthly precipitation, skin temperature, and snow depth) using decision tree regression. The quality of the merged radar signal was assessed by computing the Pearson r, root mean square error (RMSE), and relative RMSE (rRMSE) between the C-band and the corrected Ku-band signals in the overlapping years (1999–2001 and 2007–2009). We obtained high Pearson r values and low RMSE values at both the regional (r≥0.92, RMSE ≤ 0.11 dB, and rRMSE ≤ 0.38) and pixel levels (median r across pixels ≥ 0.64, median RMSE ≤ 0.34 dB, and median rRMSE ≤ 0.88), suggesting high accuracy for the data-merging procedure. The merged radar signals were then validated against the European Space Agency (ESA) ERS-2 data, which provide observations for a subset of global pixels until 2011, even after the failure of on-board gyroscopes in 2001. We found highly concordant monthly dynamics between the merged radar signals and the ESA ERS-2 signals, with regional Pearson r values ranging from 0.79 to 0.98. These results showed that our merged radar data have a consistent C-band signal dynamic. The LHScat data set (https://doi.org/10.6084/m9.figshare.20407857; Tao et al., 2023) is expected to advance our understanding of the long-term changes in, e.g., global vegetation and soil moisture with a high spatial resolution. The data set will be updated on a regular basis to include the latest images acquired by ASCAT and to include even higher spatial and temporal resolutions.

Dipolarization fronts (DFs), earthward-propagating magnetic transients with a strong magnetic field, are important regions favorable for energetic electron acceleration in the magnetotail. The DF-driven electron acceleration usually generates coherent pitch angle distributions (PADs) inside flux pileup regions (FPRs), i.e., strong magnetic field regions behind the DFs, such as pancake, butterfly, and cigar distributions, which dominate at different tail regions and often occur separately. Here we present unique observations of electron PAD evolution inside the FPR, showing that electron PAD underwent local transition from cigar distribution, to butterfly distribution, then toward pancake distribution, forming a U-shaped distribution. During the local transition, electron perpendicular flux (relative to the local magnetic field) is anticorrelated with magnetic field strength, contrary to traditional expectation. The unexpected feature of the electron U-shaped distribution is associated with multiple physical processes at different scales, including local expansion of flux tubes and pitch angle variation near the neutral sheet. These atypical observations can advance our current understanding of electron acceleration and transport in the magnetosphere.

The molybdenum (Mo) non-point source pollution in the mining area has an irreversible impact on the surrounding water and soil ecosystems. Herein, three integrated vertical subsurface flow constructed wetlands (CWs) were constructed to assess the effects of combination substrates and plant on the removal of Mo(VI). Results showed that CW1 with combination substrates and cattail exhibited a favorable removal performance for Mo(VI) at 80.90%. Moreover, most Mo(VI) retained in the CWs was retained in the substrate (58.13–88.04%), and the largest fraction of Mo(VI) retained was the water-soluble fraction on the surface of the combination substrates. Mo(VI) removal was also influenced by the microbial community composition in substrate, especially their co-occurrence networks. The species that showed significant positive correlation with Mo(VI) removal were Planctomycetes, Latescibacteria, Armatimonadetes, and Gemmatimonadetes. Moreover, CWs added plants showed that more co-occurrences interaction between taxa occurs, which means that the wetlands efficiently select recruitment of potential microbial consortia and change the co-occurrences to remove pollution in the substrate. These results could be useful in providing an ecology-based solution for the treatment of Mo(VI) in wastewater, especially in adjusting the microbial communities for Mo(VI) removal at the genetic level. HIGHLIGHTS Complex substrates coupled cattail enhanced Mo(VI) removal in constructed wetlands.; Adsorption was the main Mo(VI) removal mechanism by DWTRs and MCPP.; Mo(VI) removal was also influenced by microbial community composition.; Cattails changed the co-occurrences of microbes to remove pollution in substrates.;

Abstract Using high‐resolution Solar Wind Electron Analyzer data from the Mars Atmosphere and Volatile EvolutioN spacecraft, we report two new types of electron pitch angle distributions (PADs) in the Martian plasma environment. The first type of electron PADs, showing pitch angle primally around 45°, is termed funnel distribution, which was observed near the terminator and at about 1,200 km altitude. The second type of electron PADs, showing pitch angle primally around 135°, is termed skirt distribution, which was detected on the nightside and at about 1,250 km altitude. The electrons of funnel and skirt distributions do not exhibit any photoelectron signatures and are shown to originate from the solar wind. Through the fitting analysis, we find that the electrons showing both funnel and skirt PADs are thermal electron populations and suprathermal electron populations. In addition, the possibilities that such two types of electron pitch angle distributions corresponding to magnetic field configurations are also discussed.

The initiation and reception of shields are major risk events for shield construction in water-rich and weak strata. Although the freezing method and the steel sleeve method receiving-shield tunneling technologies both have engineering applications, the environmental safety control effect cannot meet the construction requirements of water-rich soft soil. Considering the shield construction of a typical soft soil layer in Suzhou, China, as a research target, the applicability and safety of the freezing method combined with steel sleeve receiving-shield technology in water-rich soft soil were evaluated based on a field test system. The test results show that, during soil freezing, the temperature change trend of each measuring point in the temperature measuring hole is roughly the same. The freezing process can be divided into five typical stages. The closer the active freezing period of the water-rich soft soil is to the inside of the frozen-soil curtain, the faster the development rate of the frozen wall. The soil cooling gradient increased with an increase in the radial depth. After freezing the curtain circle, the soil frost heave significantly accelerated until the frost heave amount peaked. During the construction process, special attention should be paid to the change in the value of the soil settlement during each stage change to prevent sudden changes in soil displacement. The freezing method, combined with steel sleeve receiving-shield technology, can effectively reduce the environmental disturbance caused by shield construction in water-rich soft soil.

Interplanetary discontinuities are the most common intermittent structures in the solar winds, whose jump conditions conform to the Rankine–Hugoniot relations. The identification and analysis of such planar structures rely on an accurate determination of the normal vectors. In this study, we perform a benchmark test of a technique designed to estimate the dimensionality and orientation of the magnetic fields based on four spacecraft measurements. Such a technique operates on the fact that in a strictly planar magnetic field, three columns of ∇ B are all aligned and the characteristic direction corresponding to the maximum aligned directional derivative (MADD; i.e., the maximum of $| {\boldsymbol{\nabla }}{\boldsymbol{B}}\cdot \widehat{{\boldsymbol{r}}}| $ ) can be considered the normal vector of the structure. A benchmark test of this technique, using both the simulation data and numerical models of planar structures, shows that the MADD technique can identify planar magnetic field structures efficiently. For a discontinuity with translational and rotational motion, MADD is able to estimate its normal vector with second-order precision and the speed of motion with first-order precision. Two demonstrations of the MADD application to the tangential discontinuity and shock structures are presented to show its efficiency and superiority. In principle, MADD is suitable for all types of planar structures in space. Compared with the minimum variance analysis and timing techniques, MADD is superior in the study of intermediate shocks that exhibit magnetic coplanarity and rapid spatial evolution.

Betatron acceleration, a plasma process obtaining particle energy in the perpendicular direction but reserving energy in the field-aligned direction, is the consequence of magnetic strength enhancement when the first adiabatic invariant is conserved. Such process has been widely reported in the terrestrial magnetosphere but is barely reported in other planetary environments. Here, based on the in situ measurements from the Mars Atmosphere and Volatile Evolution mission, we report two events of betatron acceleration upstream of the Martian bow shock. In both events, betatron accelerations increase the fluxes of suprathermal electrons. The acceleration processes in these events are quantitatively reproduced with an analytical model. Gratifyingly, we find the acceleration factors derived from the analytical model are well consistent with the observations of magnetic strength enhancement. These results for the first time show that the betatron acceleration is an active upstream of the Martian bow shock and is very useful to help us understand the generation of energetic electrons in the Martian environment.

Magnetic discontinuities are fundamental structures in space and laboratory plasmas where the changes in magnetic and velocity fields are constrained by Rankine–Hugoniot relations. Due to the absence of precise measurements for particles, some issues therein are hardly investigated. The nature of discontinuities driven by the magnetohydrodynamics (MHD) turbulence, and the intermediate shock are two puzzles to be solved. The MHD turbulence generates numerous discontinuities with both small normal magnetic fields and nearly constant magnetic field magnitudes in statistics. By utilizing the data from the Parker Solar Probe, we identify among the turbulence-driven discontinuities two components that exhibit diverse statistical characteristics of the plasma density, and reveal that these discontinuities comprise 80.2% rotational and 19.8% tangential discontinuities. Then, we note a special class of discontinuities within 0.35 au that have jump conditions similar to that of the rotational discontinuity and the shock simultaneously, including (1) positively correlated jumps in the plasma density and temperature, (2) a small change in the magnetic field magnitude, and (3) opposite tangential magnetic fields on two sides. These features conform to the theoretical intermediate shock, which previous studies have found to not practically exist due to the breakdown of the evolutionary condition. By the conservation law of the mass flux across a boundary, we calculate their propagation speeds and find three intermediate shock candidates with super-Alfvénic upstream and sub-Alfvénic downstream flows. This work can improve our understanding of plasma intermittencies and suggests reassessing conclusions based on ideal MHD Rankine–Hugoniot relations.

Magnetic reconnection is a crucial process responsible for energy conversion and particle acceleration in space, astrophysical, and laboratory plasmas. Turbulence and magnetic reconnection can be mutually driven, but the underlying nature of energy dissipation, intrinsic turbulence waves, and magnetic field topologies in turbulent magnetic reconnection is still poorly understood. Here, using advanced multi-spacecraft mission and innovative methods, we provide a few new perspectives to investigate the properties of the turbulence and topology in a turbulent magnetic reconnection in the magnetotail. Our results reveal that in turbulent magnetic reconnection: (1) cyclotron resonance, an important mechanism of energy dissipation, is more effective in the core region of the reconnection than in the outflow regions; (2) energy is deposited in the form of kinetic Alfvén waves (KAWs) and fast/slow waves, with KAWs corresponding to low-frequency (ion cyclotron scale) and fast/slow waves corresponding to high-frequency (low-hybrid scale); and (3) the number of spiral nulls (O-lines) were about 3.6 times more than radial ones (X-lines), and three-dimensional structures were about 8 times more prevalent than two-dimensional ones. These findings should help us better unravel the dynamics of turbulent magnetic reconnection.

Minimum variance analysis of the magnetic field (MVAB), among various techniques of planar structure analysis, is most widely used for its numerical simplicity and loose requirements for data. Through a large number of studies based on MVAB, a global picture of the solar wind intermittency has been established. However, the huge discrepancy between the results from MVAB and other techniques like timing/triangulation implies that the uncertainty of MVAB is a crucial issue that is not fully understood. Utilizing Cluster data, we establish a data set comprised of 6752 discontinuities, whose orientations are precisely determined by timing, as a benchmark for testing MVAB. We find that the scatter of the MVAB normals around the timing normal can be reduced by elevating the threshold for the eigenvalue ratio λ _2 / λ _3 and narrowing the data window to which MVAB is applied. The misidentification of discontinuities with $\left|{B}_{N}\right|/| {\boldsymbol{B}}| $ < 0.4, Δ∣ B ∣/∣ B ∣ < 0.2 as rotational discontinuities (RDs, identified by $\left|{B}_{N}\right|/| {\boldsymbol{B}}| $ > 0.4, Δ∣ B ∣/∣ B ∣ < 0.2) is proved to be a major and inherent defect of MVAB, which can occur even when λ _2 / λ _3 is large. Such a misidentification process is revealed to be related to a special discontinuity geometry. It also explains the false RD predominance reported by previous studies based on MVAB. Finally, we provide advice for the application of MVAB and discuss the possibility of obtaining the real statistical properties of interplanetary discontinuities by using MVAB.

We perform a statistical investigation of the occurrence rates of energetic electron (100–500 eV) pancake pitch-angle distributions (PADs) in the Martian space environment by utilizing 6 yr of MAVEN data. In the Martian ionosphere, we find the following: (1) at the same altitude in the terminator and night regions, the occurrences rates in the center of the southern magnetic anomaly regions are very low, but at the edges of strong magnetic fields, they increase significantly; (2) the occurrence rates of energetic electron perpendicular anisotropies on the Martian dayside increase with altitude; and (3) some closed magnetic lines in the 10°S–55°S, 30°W–125°W region at 400–800 km altitude gradually become open during the rotation of Mars from duskside to dawnside, while more closed magnetic lines are produced in the 40°S–65°S, 35°E–90°E region. In the Martian induced magnetosphere, we find the following: (1) the high-energy electron perpendicular anisotropy in the magnetosheath is the most significant; (2) the occurrence rates in the southern ( Z _MSO ≤−1 R _M ) magnetosheath are higher than those in the northern ( Z _MSO ≥ 1 R _M ) magnetosheath; (3) in the region of ∣ Z _MSO ∣ < 0.5 R _M , these high-energy electron pancake PADs are mainly concentrated in the magnetosheath region with Y _MSO ∈ [−1.4R _M , 2 R _M ]; (4) the occurrence rates in the dawnside ( Y _MSO ≤−1 R _M ) magnetosheath are higher than those in the duskside ( Y _MSO ≥ 1 R _M ) magnetosheath; and (5) in the region of ∣ Y _MSO ∣ < 0.5 R _M , the occurrence rates throughout the magnetosheath are very high.

Magnetic flux ropes or magnetic islands are important structures responsible for electron acceleration and energy conversion during turbulent reconnection. However, the evolution of flux ropes and the corresponding electron acceleration process still remain open questions. In this paper, we present a comparative study of flux ropes observed by the Magnetospheric Multiscale mission in the outflow region during an example of turbulent reconnection in Earth's magnetotail. Interestingly, we find the farther the flux rope is away from the X-line, the bigger the size of the flux rope and the slower it moves. We estimate the power density converted at the observed flux ropes via the three fundamental electron acceleration mechanisms: Fermi, betatron, and parallel electric field. The dominant acceleration mechanism at all three flux ropes is the betatron mechanism. The flux rope that is closest to the X-line, having the smallest size and the fastest moving velocity, is the most efficient in accelerating electrons. Significant energy also returns from particles to fields around the flux ropes, which may facilitate the turbulence in the reconnection outflow region.

Energetic electrons have been frequently observed during magnetic reconnection in the magnetotail. The acceleration process of the energetic electrons is not fully understood. In this paper, we select for a detailed study a case of energetic electron acceleration from the earlier reported interval of turbulent magnetic reconnection in Earth’s magnetotail observed by the Magnetospheric Multiscale mission. We use the first-order Taylor expansion method to reconstruct the magnetic topology of electron acceleration sites from the data. We find that the energetic electron fluxes increase inside the flux rope forming in front of the magnetic pileup region. We show that the energetic electrons are produced by a two-step process where two different acceleration mechanisms are successively operating outside and inside the flux rope. First, the thermal electrons are energized in the field-aligned direction inside the magnetic pileup region owing to the Fermi mechanism forming a cigar-like distribution. Second, those energized electrons are further accelerated predominately antiparallel to the magnetic field direction by a parallel electric field inside the flux rope. Our findings provide information for a better understanding of the generation of energetic electrons during turbulent reconnection process.

Wildfire can become a catastrophic natural hazard, especially during dry summer seasons in Australia. Severity is influenced by various meteorological, geographical, and fuel characteristics. Modified Mark 4 McArthur's Grassland Fire Danger Index (GFDI) is a commonly used approach to determine the fire danger level in grassland ecosystems. The degree of curing (DOC, i.e. proportion of dead material) of the grass is one key ingredient in determining the fire danger. It is difficult to collect accurate DOC information in the field, and therefore ground-observed measurements are rather limited. In this study, we explore the possibility of whether adding satellite-observed data responding to vegetation water content (vegetation optical depth, VOD) will improve DOC prediction when compared with the existing satellite-observed data responding to DOC prediction models based on vegetation greenness (normalised difference vegetation index, NDVI). First, statistically significant relationships are established between selected ground-observed DOC and satellite-observed vegetation datasets (NDVI and VOD) with an r2 up to 0.67. DOC levels estimated using satellite observations were then evaluated using field measurements with an r2 of 0.44 to 0.55. Results suggest that VOD-based DOC estimation can reasonably reproduce ground-based observations in space and time and is comparable to the existing NDVI-based DOC estimation models.

The ancient Xi’an, China, has been suffering severe land subsidence and ground fissure hazards since the 1960s, which has affected the safety of Subways. Multi-sensor SAR data are conducted to monitor the latest complex ground deformation and its influence on subway line No.3 over Xi’an. Annual deformation rates have been retrieved to reveal the spatiotemporal evolution of ground subsidence in Xi’an city from 2013 to 2017. Meanwhile, the correlation between land subsidence and ground fissures are analyzed by retrieving the deformation differences in both sides of the fissures. Besides, the deformation along subway line No. 3 is analyzed, and the fast deformation section is quantitatively studied. Finally, a flat lying sill model with distributed contractions is implemented to model the InSAR deformation over YHZ subsidence center, which manifests that the ground deformation is mainly caused by groundwater withdrawal.

The use of dynamic global vegetation models (DGVMs) to estimate CO2 emissions from land-use and land-cover change (LULCC) offers a new window to account for spatial and temporal details of emissions and for ecosystem processes affected by LULCC. One drawback of LULCC emissions from DGVMs, however, is lack of observation constraint. Here, we propose a new method of using satellite- and inventory-based biomass observations to constrain historical cumulative LULCC emissions (ELUCc) from an ensemble of nine DGVMs based on emerging relationships between simulated vegetation biomass and ELUCc. This method is applicable on the global and regional scale. The original DGVM estimates of ELUCc range from 94 to 273 PgC during 1901–2012. After constraining by current biomass observations, we derive a best estimate of 155 ± 50 PgC (1σ Gaussian error). The constrained LULCC emissions are higher than prior DGVM values in tropical regions but significantly lower in North America. Our emergent constraint approach independently verifies the median model estimate by biomass observations, giving support to the use of this estimate in carbon budget assessments. The uncertainty in the constrained ELUCc is still relatively large because of the uncertainty in the biomass observations, and thus reduced uncertainty in addition to increased accuracy in biomass observations in the future will help improve the constraint. This constraint method can also be applied to evaluate the impact of land-based mitigation activities.

Consistent forest loss estimates are important to understand the role of forest loss and deforestation in the global carbon cycle, for biodiversity studies, and to estimate the mitigation potential of reducing deforestation. To date, most studies have relied on optical satellite data and new efforts have greatly improved our quantitative knowledge on forest dynamics. However, most of these studies yield results for only a relatively short time period or are limited to certain countries. We have quantified large-scale forest loss over a 21-year period (1990–2010) in the tropical biomes of South America using remotely sensed vegetation optical depth (VOD). This passive microwave satellite-based indicator of vegetation water content and vegetation density has a much coarser spatial resolution than optical data but its temporal resolution is higher and VOD is not impacted by aerosols and cloud cover. We used the merged VOD product of the Advanced Microwave Scanning Radiometer (AMSR-E) and Special Sensor Microwave Imager (SSM/I) observations, and developed a change detection algorithm to quantify spatial and temporal variations in forest loss dynamics. Our results compared reasonably well with the newly developed Landsat-based Global Forest Change (GFC) maps, available for the 2001 onwards period (r2 = 0.90 when comparing annual country-level estimates). This allowed us to convert our identified changes in VOD to forest loss area and compute these from 1990 onwards. We also compared these calibrated results to PRODES (r2 = 0.60 when comparing annual state-level estimates). We found that South American forest exhibited substantial interannual variability without a clear trend during the 1990s, but increased from 2000 until 2004. After 2004, forest loss decreased again, except for two smaller peaks in 2007 and 2010. For a large part, these trends were driven by changes in Brazil, which was responsible for 56 % of the total South American forest loss area over our study period according to our results. One of the key findings of our study is that while forest loss decreased in Brazil after 2005, increases in other countries partly offset this trend suggesting that South American forest loss as a whole decreased much less than that in Brazil.

Water vapor tomography is a promising technique for reconstructing the 4-D moisture field, which is important to the weather forecasting and nowcasting as well as to the numerical weather prediction. A near real-time 4-D water vapor tomographic system is developed in this study. GPS slant water vapor (SWV) observations are derived by a sliding time window strategy using double-difference model and predicted orbits. Besides GPS SWV, surface water vapor measurements are also assimilated as real time observations into the tomographic system in order to improve the distribution of observations in the lowest layers of tomographic grid. A 1-year term experiment in Hong Kong was carried out. The feasibility of the GPS data processing strategy is demonstrated by the good agreement between the time series of GPS-derived Precipitable Water Vapor (PWV) and radio-sounding-derived PWV with a bias of 0.04 mm and a root-mean-square error (RMSE) of 1.75 mm. Using surface humidity observations in the tomographic system, the bias and RMSE between tomography and radiosonde data are decreased by half in the ground level, but such improved effects weaken gradually with the rise of altitude until becoming adverse above 4000 m. The overall bias is decreased from 0.17 to 0.13 g m−3 and RMSE is reduced from 1.43 to 1.28 g m−3. By taking the correlation coefficient and RMSE between tomography and radiosonde individual profile as the statistical measures, quality of individual profile is also improved as the success rate of tomographic solution is increased from 44.44 to 63.82% while the failure rate is reduced from 55.56 to 36.18%.

Previous work has failed to address fully the response of (autotrophic and heterotrophic) respiration to grazing in different ecosystems, particularly in alpine grasslands outside the growing season. From 2010 to 2011 a field experiment combined two methods (static closed chambers and a closed dynamic soil CO2 flux system) in alpine grasslands located in the Tianshan Mountains. We examined the effects of grazing regime on ecosystem respiration (Re) both outside (NGS) and during (GS) the growing season and determined the pattern of Re in relation to climate change. There was no significant change in CO2 emissions under grazing. Heterotrophic respiration (Rh) accounted for 78.5% of Re with short-term grazing exclusion and 93.2% of Re with long-term grazing exclusion. Re, Rh and autotrophic respiration (Ra) fluxes outside the growing season were equivalent to 12.9%, 14.1% and 11.4% of the respective CO2 fluxes during the growing season. In addition, our results indicate that soil water content played a critical role in Ra in the cold and arid environment. Both Rh and Re were sensitive to soil temperature. Moreover, our results suggest that grazing exerted no significant effect on CO2 emissions in these alpine grasslands.

Widespread nitrogen (N) enrichment resulting from anthropogenic activities has led to great changes in carbon exchange between the terrestrial biosphere and the atmosphere. Grassland is one of the most sensitive ecosystems to N deposition. However, the effect of N deposition on ecosystem respiration (Re) in grasslands has been conducted mainly in temperate grasslands, which are limited largely by water availability, with few studies focused on alpine grasslands that are primarily constrained by low temperatures. Failure to assess the magnitude of the response in Re outside the growing season (NGS) in previous studies also limits our understanding of carbon exchange under N deposition conditions. To address these knowledge gaps we used a combination of static closed chambers and gas chromatography in an alpine grassland from 2010 to 2011 to test the effects of N application on ecosystem respiration (Re) both inside and outside the growing season. There was no significant change in CO2 emissions under N application. Re outside the growing season was at least equivalent to 9.4% of the CO2 fluxes during the growing season (GS). Annual Re was calculated to be 279.0–403.9 g CO2 m−2 yr−1 in Bayinbuluk alpine grasslands. In addition, our results indicate that soil temperature was the dominant abiotic factor regulating variation in Re in the cold and arid environment. Our results suggest that short-term N additions exert no significant effect on CO2 emissions in alpine grassland.

The Discrete Global Grid System (DGGS) is a new type of global spatial data model and is the extension of the plane grid on a sphere. The hexagon is usually used in the construction of DGGS for its advantageous geometric structure. The paper principally focuses on the issue of modeling and expression of vector data in the hexagon DGGS. The precision of vector data is the basis of data recording and data expression, and data with different precision fall into the grid cells of corresponding sizes, making the gridding data themselves contain the precision and scale information. The present method of data recording is reserved, as far as possible, in the data recording process, and only the geometric information of vectors is substituted by the one-dimension coding of grids. This approach is more simple and effective than the digital coordinate recording method. The gridding expression of vector data differs from the traditional technique, mainly due to the subdivision of the durative space by grids as well as the obedience of the subdivision special rules, among which the point expression should activate the corresponding grid cells in the light of the point coordinates. Linear expression should activate the corresponding grid cells of every coordinate as well as the connected grids between every two node cells, and area expression should express both the boundary and internal regions by virtue of grid cells. For spherical expression, vector data have to solve the cell filling problem, but also the extension from planum to sphere. This paper puts forward a reasonable sphere extension approach, in which the vector data expression on the spherical grids was accomplished by the dismantling of vector data on different extended areas and the multi-times transformation. Besides, the algorithm in connection with the vector data was verified through experiments for its effect and efficiency. Moreover, the distance and direction of vector data on the grids would change in the mapping process from planum to sphere girds, leading to an inaccurate spherical gridding expression. So, the effects on the rectilinear direction in grids of the hexagon from the planum-sphere mapping process was investigated, and accuracy control of the spherical expression was processed to make sure that the drawing error of the spherical grids for vector data should be limited within one cell.

Drylands, covering nearly 30% of the global land surface, are characterized by high climate variability and sensitivity to land management. Here, two satellite-observed vegetation products were used to study the long-term (1988–2008) vegetation changes of global drylands: the widely used reflective-based Normalized Difference Vegetation Index (NDVI) and the recently developed passive-microwave-based Vegetation Optical Depth (VOD). The NDVI is sensitive to the chlorophyll concentrations in the canopy and the canopy cover fraction, while the VOD is sensitive to vegetation water content of both leafy and woody components. Therefore it can be expected that using both products helps to better characterize vegetation dynamics, particularly over regions with mixed herbaceous and woody vegetation. Linear regression analysis was performed between antecedent precipitation and observed NDVI and VOD independently to distinguish the contribution of climatic and non-climatic drivers in vegetation variations. Where possible, the contributions of fire, grazing, agriculture and CO2 level to vegetation trends were assessed. The results suggest that NDVI is more sensitive to fluctuations in herbaceous vegetation, which primarily uses shallow soil water, whereas VOD is more sensitive to woody vegetation, which additionally can exploit deeper water stores. Globally, evidence is found for woody encroachment over drylands. In the arid drylands, woody encroachment appears to be at the expense of herbaceous vegetation and a global driver is interpreted. Trends in semi-arid drylands vary widely between regions, suggesting that local rather than global drivers caused most of the vegetation response. In savannas, besides precipitation, fire regime plays an important role in shaping trends. Our results demonstrate that NDVI and VOD provide complementary information and allow new insights into dryland vegetation dynamics.

Combining information derived from satellite-based passive and active microwave sensors has the potential to offer improved estimates of surface soil moisture at global scale. We develop and evaluate a methodology that takes advantage of the retrieval characteristics of passive (AMSR-E) and active (ASCAT) microwave satellite estimates to produce an improved soil moisture product. First, volumetric soil water content (m3 m−3) from AMSR-E and degree of saturation (%) from ASCAT are rescaled against a reference land surface model data set using a cumulative distribution function matching approach. While this imposes any bias of the reference on the rescaled satellite products, it adjusts them to the same range and preserves the dynamics of original satellite-based products. Comparison with in situ measurements demonstrates that where the correlation coefficient between rescaled AMSR-E and ASCAT is greater than 0.65 ("transitional regions"), merging the different satellite products increases the number of observations while minimally changing the accuracy of soil moisture retrievals. These transitional regions also delineate the boundary between sparsely and moderately vegetated regions where rescaled AMSR-E and ASCAT, respectively, are used for the merged product. Therefore the merged product carries the advantages of better spatial coverage overall and increased number of observations, particularly for the transitional regions. The combination method developed has the potential to be applied to existing microwave satellites as well as to new missions. Accordingly, a long-term global soil moisture dataset can be developed and extended, enhancing basic understanding of the role of soil moisture in the water, energy and carbon cycles.

Understanding the error structures of remotely sensed soil moisture observations is essential for correctly interpreting observed variations and trends in the data or assimilating them in hydrological or numerical weather prediction models. Nevertheless, a spatially coherent assessment of the quality of the various globally available datasets is often hampered by the limited availability over space and time of reliable in-situ measurements. As an alternative, this study explores the triple collocation error estimation technique for assessing the relative quality of several globally available soil moisture products from active (ASCAT) and passive (AMSR-E and SSM/I) microwave sensors. The triple collocation is a powerful statistical tool to estimate the root mean square error while simultaneously solving for systematic differences in the climatologies of a set of three linearly related data sources with independent error structures. Prerequisite for this technique is the availability of a sufficiently large number of timely corresponding observations. In addition to the active and passive satellite-based datasets, we used the ERA-Interim and GLDAS-NOAH reanalysis soil moisture datasets as a third, independent reference. The prime objective is to reveal trends in uncertainty related to different observation principles (passive versus active), the use of different frequencies (C-, X-, and Ku-band) for passive microwave observations, and the choice of the independent reference dataset (ERA-Interim versus GLDAS-NOAH). The results suggest that the triple collocation method provides realistic error estimates. Observed spatial trends agree well with the existing theory and studies on the performance of different observation principles and frequencies with respect to land cover and vegetation density. In addition, if all theoretical prerequisites are fulfilled (e.g. a sufficiently large number of common observations is available and errors of the different datasets are uncorrelated) the errors estimated for the remote sensing products are hardly influenced by the choice of the third independent dataset. The results obtained in this study can help us in developing adequate strategies for the combined use of various scatterometer and radiometer-based soil moisture datasets, e.g. for improved flood forecast modelling or the generation of superior multi-mission long-term soil moisture datasets.

Vertisols are clay soils that are common in the monsoonal and dry warm regions of the world. One of the characteristics of these soil types is to form deep cracks during periods of extended dry, resulting in significant variation of the soil and hydrologic properties. Understanding the influence of these varying soil properties on the hydrological behavior of the system is of considerable interest, particularly in the retrieval or simulation of soil moisture. In this study we compare surface soil moisture (θ in m3 m−3) retrievals from AMSR-E using the VUA-NASA (Vrije Universiteit Amsterdam in collaboration with NASA) algorithm with simulations from the Community Land Model (CLM) over vertisol regions of mainland Australia. For the three-year period examined here (2003–2005), both products display reasonable agreement during wet periods. During dry periods however, AMSR-E retrieved near surface soil moisture falls below values for surrounding non-clay soils, while CLM simulations are higher. CLM θ are also higher than AMSR-E and their difference keeps increasing throughout these dry periods. To identify the possible causes for these discrepancies, the impacts of land use, topography, soil properties and surface temperature used in the AMSR-E algorithm, together with vegetation density and rainfall patterns, were investigated. However these do not explain the observed θ responses. Qualitative analysis of the retrieval model suggests that the most likely reason for the low AMSR-E θ is the increase in soil porosity and surface roughness resulting from cracking of the soil. To quantitatively identify the role of each factor, more in situ measurements of soil properties that can represent different stages of cracking need to be collected. CLM does not simulate the behavior of cracking soils, including the additional loss of moisture from the soil continuum during drying and the infiltration into cracks during rainfall events, which results in overestimated θ when cracks are present. The hydrological influence of soil physical changes are expected to propagate through the modeled system, such that modeled infiltration, evaporation, surface temperature, surface runoff and groundwater recharge should be interpreted with caution over these soil types when cracks might be present. Introducing temporally dynamic roughness and soil porosity into retrieval algorithms and adding a "cracking clay" module into models are expected to improve the representation of vertisol hydrology.