Assessing Global Water Storage Variability from GRACE: Trends,Seasonal Cycle,Subseasonal Anomalies and Extremes

Throughout the past decade, the Gravity Recovery and Climate Experiment (GRACE) has given an unprecedented view on global variations in terrestrial water storage. While an increasing number of case studies have provided a rich overview on regional analyses, a global assessment on the dominant features of GRACE variability is still lacking. To address this, we survey key features of temporal variability in the GRACE record by decomposing gridded time series of monthly equivalent water height into linear trends, inter-annual, seasonal, and subseasonal (intra-annual) components. We provide an overview of the relative importance and spatial distribution of these components globally. A correlation analysis with precipitation and temperature reveals that both the inter-annual and subseasonal anomalies are tightly related to fluctuations in the atmospheric forcing. As a novelty, we show that for large regions of the world high-frequency anomalies in the monthly GRACE signal, which have been partly interpreted as noise, can be statistically reconstructed from daily precipitation once an adequate averaging filter is applied. This filter integrates the temporally decaying contribution of precipitation to the storage changes in any given month, including earlier precipitation. Finally, we also survey extreme dry anomalies in the GRACE record and relate them to documented drought events. This global assessment sets regional studies in a broader context and reveals phenomena that had not been documented so far.     

Soil carbon losses by sheet erosion: a potentially critical contribution to the global carbon cycle

Despite soil erosion through water being a ubiquitous process and its environmental consequences being well understood, its effects upon the global carbon cycle still remain largely uncertain. How much soil organic carbon (SOC) is removed each year from soils by sheet wash, an important if not the most efficient mechanism of detachment and transport of surficial soil material? What are the main environnemental controls worldwide? These are important questions which largely remain unanswered. Empirical data from 240 runoff plots studied over entire rainy seasons from different regions of the world were analysed to estimate particulate organic carbon (POC) losses (POC_L), and POC enrichment in the sediments compared to the bulk soil (ER), which can be used as a proxy of the fate of the eroded POC. The median POC_L was 9.9 g C m^‐2 y^‐1 with highest values observed for semi‐arid soils (POC_L = 10.8 g C m^‐2 y^‐1), followed by tropical soils (POC_L = 6.4 g C m^‐2 y^‐1) and temperate soils (POC_L = 1.7 g C m^‐2 y^‐1). Considering the mean POC_L of 27.2 g C m^‐2 y^‐1, the total amount of SOC displaced annually by sheet erosion from its source would be 1.32 ± 0.20 Gt C, i.e. 14.6% of the net annual fossil fuel induced C emissions of 9 Gt C. Because of low sediment enrichment in POC, erosion‐induced CO₂ emissions are likely to be limited in clayey environments while POC burial within hillslopes is likely to constitute an important carbon sink. In contrast, most of the POC displaced from sandy soils is likely to be emitted to the atmosphere. These results underpin the major role sheet wash plays in the displacement of SOC from its source and in the fate of the eroded SOC, with large variations across the different pedo‐climatic regions of the world. Copyright © 2015 John Wiley & Sons, Ltd.     

Fluid flow reconstruction in hanging and footwall carbonates: Compartmentalization by Cenozoic reverse faulting in the Northern Oman Mountains (UAE)

In this paper, the diagenesis from either side of a major Cenozoic reverse fault in the Northern Oman Mountains is documented. Detailed petrographical and geochemical analysis of calcite-filled fractures in carbonate strata of Late Triassic and Early Cretaceous age in the hanging wall and footwall in Wadi Ghalilah reflect a different diagenetic history. In both hanging wall and footwall most of the fractures are pre-burial, extensional in origin, formed by a crack-seal mechanism, and the calcite vein infill has a host-rock buffered signature. In the hanging wall, the fluid responsible for calcite precipitation of these extensional fractures was a marine fluid at 60 °C. These veins predate deep burial and contractional tectonic deformation and consequently do not provide any information about syntectonic fluid flow. Neither do the pre-burial extension fractures in the footwall which are also host-rock buffered. The fractures post-dating the tectonic stylolitization in the footwall, by contrast, show evidence of syntectonic migration of saline formation waters at temperatures between 80 and 160 °C during contractional deformation. These fluids probably were sourced from the subsurface via the reverse fault, which acted as a fluid conduit. At the same time, however, this fault functioned as a permeability barrier towards the hanging wall, since no evidence of syntectonic fluid flow is present here. In this way compartmentalization of the hanging wall and footwall block was realized.     

Jason-1 global statistical evaluation and performance assessment: Calibration and cross-calibration results

Since Jason-1launch, extensive validation of Jason-1 data and cross-calibration relative to TOPEX/Poseidon (T/P) have been performed by the CLS validation team within the CNES Jason-1 project. These validation activities are routinely operated as part of the Jason-1 ground segment, and often lead to in-depth studies to understand all validation conclusions. This paper presents the main results in terms of Jason-1 data quality: verification of data availability and validity, monitoring of the most relevant altimeter and radiometer parameters, assessment of the Jason-1 altimeter system performances. From global statistical analysis of more than 2 years of Jason-1 GDR data, results for all components of the altimeter measurement are derived in terms of bias, trend and precision. This work also represents a contribution to the estimation of the Jason-1 error budget. Thorough studies have been more focused on specific issues in relation to data quality: this is the case for the analysis of the high frequency content of the Jason-1 data and its impact on the T/P to Jason-1 comparison. From the results presented in this paper, it is demonstrated that the Jason-1 mission fulfils the requirements of high precision altimetry. In particular, it allows continuing the observation of the Mean Sea Level (MSL) variations at the same accuracy as T/P, which was one of the challenges of the Jason-1 mission. Potential improvements and open issues are also identified, with the objective of still making progress in terms of altimeter data quality.     

Phytoplankton community adaptation to changing light levels in the southern Beaufort Sea,Canadian Arctic

The chlorophyll a specific absorption coefficient of phytoplankton, a_φ^∗(λ) is an important parameter to determine for primary production models and for the estimation of phytoplankton physiological condition. Knowledge of this parameter at high latitudes where nutrient rich cold water submitted to low incident light is a common environment is almost nonexistent. To address this issue, we investigated the light absorption properties of phytoplankton as a function of irradiance, temperature, and nutrients using a large data set in the southern Beaufort Sea during the open water to ice cover transition period. The a_φ^∗(λ) tended to increase from autumn when open water still existed to early winter when sea ice cover was formed, resulting from a biological selection of smaller-size phytoplankton more efficient to absorb light. There was no significant correlation between a_φ^∗(λ) and irradiance or temperature for both seasons. However, a_φ^∗(λ) showed a significant positive correlation with NO₃ + NO₂. Implications of the results for phytoplankton community adaptation to changing light levels are discussed.     

Input of particulate heavy metals from rivers and associated sedimentary deposits on the Gulf of Lion continental shelf

Fluxes of the heavy metals chromium (Cr), cobalt (Co), nickel (Ni), copper (Cu), lead (Pb), cadmium (Cd) and zinc (Zn) delivered by rivers to the Gulf of Lion (NW Mediterranean Sea) were estimated over a three year study of the River Rhone and its smaller tributaries. Most of the particulate metal fluxes (80–90%) delivered by these rivers occurred within a very short period of time (less than 12%), a typical trend for the Mediterranean environment, where highly contrasting hydrological regimes were observed over the year. Temporal and spatial variations in the fluxes of these particulate metals were driven by the fluxes in both water discharge and suspended particulate matter load. On the shelf, these particulate metal fluxes, largely arising from the Rhone watershed, were two to ten times more important than those resulting from atmospheric deposition. Co, Cr and Ni in the rivers and on the shelf surface sediments were mainly natural and associated with the finest particles. Cd and Phosphorus appeared to be associated with the silt fraction and to be enriched in the prodelta areas. Pb, Zn and Cu were more closely associated with the organic matter content and also showed enrichment in the organic rich prodeltaic sediments. Anthropogenic influences diminished offshore, except for Pb and Zn which could be supplied from the atmosphere by man-made aerosols. Although most of the metals tended to be enriched in the prodelta areas, these did not constitute a permanent sink due to resuspension processes affecting these shallow depths. A resuspension experiment conducted on sediment cores from the Rhone prodelta demonstrated that metal deposited on the surface layer, especially those associated with the organic matter, may be resuspended; this should be taken into account for a complete understanding of the biogeochemical cycle of these metals.     

Vertical structure variability in a seasonal simulation of a medium-resolution regional model of the Eastern South Pacific

A seasonal simulation from a medium-resolution ocean general circulation mode (OGCM) is used to investigate the vertical structure variability of the Southeast Pacific (SEP). The focus is on the extra-tropical Rossby wave (ETRW) variability and associated forcing mechanism. Some aspects of the model mean state are validated from available observations, which justifies a vertical mode decomposition of the model variability. The analysis of the baroclinic mode contributions to sea level indicates that the gravest mode is dominant over most of the domain at all frequencies. Annual variability is on average twice as large as the semi-annual variability which is confined near the coast for all the modes. The first baroclinic mode contribution to the annual cycle exhibits a clear westward propagation north of the critical latitude. The higher-order modes only contribute near the coast where they are associated with vertically propagating energy. The residual variability, which is the energy at all timescales other than annual and semi-annual periods peaks offshore between 20°S and 30°S for all baroclinic modes. The third baroclinic mode also exhibits a relative maximum variability off the coast of Peru south of the critical latitude of the annual cycle (13°S), where the Peru–Chile Undercurrent is the most intense. Sensitivity experiments to the atmospheric and boundary forcing suggest that the residual variability results from the non-linear interaction between annual Rossby waves and the mean flow, while the annual ETRWs in the model result from the summed-contribution from both the local wind stress and remote equatorial forcing. Overall the study extends the classical analysis of sea level variability in the SEP based on linear theory, and suggests that the peculiarities of the baroclinic modes need to be taken into account for interpreting the sea level variability and understanding its connection with the equatorial variability.     

Local versus regional-scale characteristics of monsoon onset and post-onset rainfall over Indonesia

The austral summer monsoon onset and post-onset rainfall and their associated low-level winds are analyzed during the August–February season over Indonesia from 1979 to 2006 using surface and satellite products as well as reanalyses and regional climate model simulations. Onset date is defined using a local agronomic definition. Its leading empirical orthogonal function is found to exhibit a regional-scale spatially-coherent signal across “monsoonal” Indonesia, i.e. mostly south of the equator, with an asymmetric temporal behavior, such that delayed onsets are more intense than early ones. Associated anomalies in rainfall tend to weaken quickly after mid-to-late November or early December, especially over islands, while they tend to persist over ocean. This weakening is shown to be associated with the evolution of distinct weather types revealed by a k-means cluster analysis. In particular, late onsets—usually related to warm El Niño Southern oscillation events—are found to be accompanied by an increased prevalence of a weather type characterized by weak low-level daily-averaged winds across monsoonal Indonesia and increased (decreased) rainfall over most of the island orography and southern and western coasts (seas). The regional model simulations provide evidence that this land–sea rainfall contrast could be associated with an enhanced diurnal sea-land breeze circulation.     

Influence of the warm pool and cold tongue El Niños on the following Caribbean rainy season rainfall

Caribbean rainfall and associated regional-scale ocean–atmosphere anomalies are analyzed during and after warm pool (WP) and cold tongue (CT) El Niño (EN) events (i.e. from the usual peak of EN events in boreal winter to next summer from 1950 to 2011). During and after a CT event, a north–south dipolar pattern with positive (negative) rainfall anomalies over the northern (southern) Caribbean during the boreal winter tends to reverse in spring, and then to vanish in summer. On the contrary, during and after a WP event, weak rainfall anomalies during the boreal winter intensify themselves from spring, with anomalous wet conditions over most of the Caribbean basin observed during summer, except over the eastern coast of Nicaragua and Costa Rica. The Caribbean rainfall anomalies associated with WP and CT events are shaped by competition between at least four different, but interrelated, mechanisms; (1) the near-equatorial large-scale subsidence anomaly over the equatorial Atlantic linked to the zonal adjustment of the Walker circulation; (2) the extra-tropical wave-like train combining positive phase of the Pacific/North American mode and negative phase of the North Atlantic Oscillation; (3) the wind-evaporation-sea surface temperature (SST) positive feedback coupling warmer-than-normal SST with weaker-than-normal low level easterlies over the tropical North Atlantic; and (4) the air-sea coupling between the speed of low level easterlies, including the Caribbean low level jet, and the SST anomaly (SSTA) gradient between the Caribbean basin and the eastern equatorial Pacific. It seems that Caribbean rainfall anomalies are shaped mostly by mechanisms (1–3) during CT events from the boreal winter to spring. These mechanisms seem less efficient during WP events when the atmospheric response seems driven mostly by mechanism (4), coupling positive west-east SSTA gradient with weaker-than-normal low level easterlies, and secondary by mechanism (3), from the boreal spring to summer.     

Temporal changes in seasonal precipitation over the Sahara Desert from 1979 to 2016

Rainfall variability dominates livelihoods in all countries of Saharan Africa. To better understand the processes involved in Sahara precipitation changes, we used the Global Precipitation Climatology Center(GPCC) dataset to examine dry and wet seasonal trends in the Sahara region from 1979 to 2016. We also used the European Centre for Medium-Range Weather Forecasts(ECMWF) to evaluate the general atmospheric circulation associated with seasonal change of Sahara precipitation. The Mann-Kendall test and Theil sens' slope estimator methods were adopted to test and estimate the significance and weight of precipitation trend, respectively. The results revealed that Sahara precipitation has increased significantly.The seasonal evaluation shows a positive trend of 0.42 mm/decade and 1.43 mm/decade in JAS(June, August, and September) seasons for the northern and southern Saharan Desert, respectively. Moreover, the JFMA(January, February,March, and April) period shows a negative trend but not statistically significant. An examination of the general circulation and moisture transport changes suggested an increase of rainfall in southern Sahara. The wet period is also driven by northward penetration of moisture originating from the Sahel region, African Easterly Jet(AEJ), and weakening in the upper tropospheric zonal wind. Summer rainfall has also been likely associated with positive anomalies of sea surface temperature(SST) in the North Tropical Atlantic(NTA) and the Mediterranean Sea.