Search

Presents results from mechanistic experiments to clarify the origin and maintenance of an extratropical climate extreme of drought. How during April and May 1998, sea surface temperature anomalies combined with atmospheric circulation to establish the drought affecting Oklahoma and Texas; Maintenance of the drought by lower evaporation and precipitation; End of the drought with weather systems that were able to penetrate the region and overwhelm the soil-moisture feedback; Potential of numerical models to make predictions of regional climate.

Predicting the temporal and spatial patterns of South Asian monsoon rainfall within a season is of critical importance due to its impact on agriculture, water availability, and flooding. The monsoon intraseasonal oscillation (MISO) is a robust northwardpropagating mode that determines the active and break phases of the monsoon and much of the regional distribution of rainfall. However, dynamical atmospheric forecast models predict this mode poorly. Data-driven methods for MISO prediction have shown more skill, but only predict the portion of the rainfall corresponding to MISO rather than the full rainfall signal. Here, we combine state-of-the-art ensemble precipitation forecasts from a high-resolution atmospheric model with data-driven forecasts of MISO. The ensemble members of the detailed atmospheric model are projected onto a lower-dimensional subspace corresponding to the MISO dynamics and are then weighted according to their distance from the data-drivenMISOforecast in this subspace. We thereby achieve improvements in rainfall forecasts over India, as well as the broader monsoon region, at 10-to 30-d lead times, an interval that is generally considered to be a predictability gap. The temporal correlation of rainfall forecasts is improved by up to 0.28 in this time range. Our results demonstrate the potential of leveraging the predictability of intraseasonal oscillations to improve extended-range forecasts; more generally, they point toward a future of combining dynamical and data-driven forecasts for Earth system prediction. [ABSTRACT FROM AUTHOR]

Hybrid systems have become the state of the art among data assimilation methods. These systems combine the benefits of two other systems that are traditionally used in operational weather forecasting: an ensemble-based system and a variational system. One of the most recently proposed hybrid approaches is called hybrid gain (HG). It obtains the final analysis as a linear combination of two analyses, assuming that the innovations (i.e. the forecast and the set of observations used) between the two data assimilation methods are identical. A perfect model experiment was performed using the HG in the SPEEDY model to show a new methodology to assign different weights to the two analyses, LETKF and 3D-Var in the generation of the final analysis. Our new approach uses, in the assignment of the weights, the ensemble spread, considered to be a measure of uncertainty in the LETKF. Thus, it is possible to use the estimation of the uncertainty of the analysis that the LETKF provides, to determine where the system should give more weight to the LETKF or the 3D-Var analysis. For this purpose, we define a geographically varying weighting factor alpha, which multiplies the 3D-Var analysis, as the normalised spread for each variable at each level. Then, (1-alpha), which decreases with increasing spread, becomes the factor that multiplies the LETKF analysis. The underlying mechanism of the spread–error relationship is explained using a toy model experiment. The results are very encouraging: the original HG and the new weighted HG analyses have similar high quality and are better than both 3D-Var and LETKF. However, the dynamically weighted HG analyses are significantly more balanced than the original HG analyses are, which has probably contributed to the consistently improved performance observed in the weighted HG, which increases with time throughout the 5-day forecasts. [ABSTRACT FROM AUTHOR]

The Human System is within the Earth System. They should be modeled bidirectionally coupled, as they are in reality. The Human System is rapidly expanding, mostly due to consumption of fossil fuels (approximately one million times faster than Nature accumulated them) and fossil water. This threatens not only other planetary subsystems but also the Human System itself. Carrying Capacity is an important tool to measure sustainability, but there is a widespread view that Carrying Capacity is not applicable to humans. Carrying Capacity has generally been prescribed a priori, mostly using the logistic equation. However, the real dynamics of human population and consumption are not represented by this equation or its variants. We argue that Carrying Capacity should not be prescribed but should insteadbe dynamically derived a posteriori from the bidirectional coupling of Earth System submodels with the Human System model. We demonstrate this approach with a minimal model of Human–Nature interaction (HANDY). ▪ The Human System is a subsystem of the Earth System, with inputs (resources) from Earth System sources and outputs (waste, emissions) to Earth System sinks. ▪ The Human System is growing rapidly due to nonrenewable stocks of fossil fuels and water and threatens the sustainability of the Human System and to overwhelm the Earth System. ▪ Carrying Capacity has been prescribed a priori and using the logistic equation, which does not represent the dynamics of the Human System. ▪ Our new approach to human Carrying Capacity is derived from dynamically coupled Earth System–Human System models and can be used to estimate the sustainability of the Human System. [ABSTRACT FROM AUTHOR]

Satellite observations are indispensable for weather forecasting, climate change monitoring, and environmental studies. Understanding and quantifying errors and uncertainties associated with satellite observations are essential for hardware calibration, data assimilation, and developing environmental and climate data records. Satellite observation errors can be classified into four categories: measurement, observation operator, representativeness, and preprocessing errors. Current methods for diagnosing observation errors still yield large uncertainties due to these complex errors. When simulating satellite errors, empirical errors are usually used, which do not always accurately represent the truth. We address these challenges by developing an error inventory simulator, the Satellite Error Representation and Realization (SatERR). SatERR can simulate a wide range of observation errors, from instrument measurement errors to model assimilation errors. Most of these errors are based on physical models, including existing and newly developed algorithms. SatERR takes a bottom-up approach: errors are generated from root sources and forward propagate through radiance and science products. This is different from, but complementary to, the top-down approach of current diagnostics, which inversely solves unknown errors. The impact of different errors can be quantified and partitioned, and a ground-truth testbed can be produced to test and refine diagnostic methods. SatERR is a community error inventory, open-source on GitHub, which can be expanded and refined with input from engineers, scientists, and modelers. This debut version of SatERR is centered on microwave sensors, covering traditional large satellites and small satellites operated by NOAA, NASA, and EUMETSAT. [ABSTRACT FROM AUTHOR]

Assimilation methods, meant to constrain divergence of model trajectory from reality using observations, do not exactly satisfy the physical laws governing the model state variables. This allows mismatches in the analysis in the vicinity of observation locations where the effect of assimilation is most prominent. These mismatches are usually mitigated either by the model dynamics in between the analysis cycles and/or by assimilation at the next analysis cycle. However, if the observations coverage is limited in space, as it was in the ocean before the Argo era, these mechanisms may be insufficient to dampen the mismatches, which we call shocks, and they may remain and grow. Here we show through controlled experiments, using real and simulated observations in two different ocean models and assimilation systems, that such shocks are generated in the ocean at the lateral boundaries of the moored buoy network. They thrive and propagate westward as Rossby waves along these boundaries. However, these shocks are essentially eliminated by the assimilation of near-homogenous global Argo distribution. These findings question the fidelity of ocean reanalysis products in the pre-Argo era. For example, a reanalysis that ignores Argo floats and assimilates only moored buoys, wrongly represents 2008 as a negative Indian Ocean Dipole year. [ABSTRACT FROM AUTHOR]

Due to the physical coupling between atmosphere and ocean, information about the ocean helps to better predict the future of the atmosphere, and in turn, information about the atmosphere helps to better predict the ocean. Here, we investigate the spatial and temporal nature of this predictability: where, for how long, and at what frequencies does the ocean significantly improve prediction of the atmosphere, and vice versa? We apply Granger causality, a statistical test to measure whether a variable improves prediction of another, to local time series of sea surface temperature (SST) and low-level atmospheric variables. We calculate the detailed spatial structure of the atmosphere-to-ocean and ocean-to-atmosphere predictability. We find that the atmosphere improves prediction of the ocean most in the extratropics, especially in regions of large SST gradients. This atmosphere-to-ocean predictability is weaker but longer-lived in the tropics, where it can last for several months in some regions. On the other hand, the ocean improves prediction of the atmosphere most significantly in the tropics, where this predictability lasts for months to over a year. However, we find a robust signature of the ocean on the atmosphere almost everywhere in the extratropics, an influence that has been difficult to demonstrate with model studies. We find that both the atmosphere-to-ocean and ocean-to-atmosphere predictability are maximal at low frequencies, and both are larger in the summer hemisphere. The patterns we observe generally agree with dynamical understanding and the results of the Kalnay dynamical rule, which diagnoses the direction of forcing between the atmosphere and ocean by considering the local phase relationship between simultaneous sea surface temperature and vorticity anomaly signals. We discuss applications to coupled data assimilation. [ABSTRACT FROM AUTHOR]

This study extends recent ocean reanalysis comparisons to explore improvements to several next-generation products, the Simple Ocean Data Assimilation, version 3 (SODA3); the Estimating the Circulation and Climate of the Ocean, version 4, release 3 (ECCO4r3); and the Ocean Reanalysis System 5 (ORAS5), during their 23-yr period of overlap (1993–2015). The three reanalyses share similar historical hydrographic data, but the forcings, forward models, estimation algorithms, and bias correction methods are different. The study begins by comparing the reanalyses to independent analyses of historical SST, heat, and salt content, as well as examining the analysis-minus-observation misfits. While the misfits are generally small, they still reveal some systematic biases that are not present in the reference Hadley Center EN4 objective analysis. We next explore global trends in temperature averaged into three depth intervals: 0–300, 300–1000, and 1000–2000 m. We find considerable similarity in the spatial structure of the trends and their distribution among different ocean basins; however, the trends in global averages do differ by 30%–40%, which implies an equivalent level of disagreement in net surface heating rates. ECCO4r3 is distinct in having quite weak warming trends while ORAS5 has stronger trends that are noticeable in the deeper layers. To examine the performance of the reanalyses in the Arctic we explore representation of Atlantic Water variability on the Atlantic side of the Arctic and upper-halocline freshwater storage on the Pacific side of the Arctic. These comparisons are encouraging for the application of ocean reanalyses to track ocean climate variability and change at high northern latitudes. [ABSTRACT FROM AUTHOR]

These models range from stylized mathematical models [[4]], system dynamics models [[5]], and agent-based models [[6]], to complex integrated assessment models (IAMs) [[7]] and coupled component models [[8]] such as the Earth system models [[9]]. The Earth is a very large and complex system that consists of human and natural components interacting bidirectionally with each other, thus forming coupled human and natural systems (CHANS) [[1]]. B Advances in understanding CHANS. b The increasing research efforts on CHANS at multiple scales [[11]], including theoretical, empirical, and modeling studies, enhance understanding of CHANS components, system behaviors, and feedback mechanisms. [Extracted from the article]

Despite dramatic improvements over the last decades, operational NWP forecasts still occasionally suffer from abrupt drops in their forecast skill. Such forecast skill 'dropouts' may occur even in a perfect NWP system because of the stochastic nature of NWP but can also result from flaws in the NWP system. Recent studies have shown that dropouts occur due not to a model's deficiencies but to misspecified initial conditions, suggesting that they could be mitigated by improving the quality control (QC) system so that the observation-minus-background (O-B) innovations that would degrade a forecast can be detected and rejected. The ensemble forecast sensitivity to observations (EFSO) technique enables for the quantification of how much each observation has improved or degraded the forecast. A recent study has shown that 24-h EFSO can detect detrimental O-B innovations that caused regional forecast skill dropouts and that the forecast can be improved by not assimilating them. Inspired by that success, a new QC method is proposed, termed proactive QC (PQC), that detects detrimental innovations 6 h after the analysis using EFSO and then repeats the analysis and forecast without using them. PQC is implemented and tested on a lower-resolution version of NCEP's operational global NWP system. It is shown that EFSO is insensitive to the choice of verification and lead time (24 or 6 h) and that PQC likely improves the analysis, as attested to by forecast improvements of up to 5 days and beyond. Strategies for reducing the computational costs and further optimizing the observation rejection criteria are also discussed. [ABSTRACT FROM AUTHOR]

Atmospheric inversion of high spatiotemporal surface CO2 fllux without dynamic constraints and sufficient observations is an ill-posed problem, and a priori flux from a "bottom-up" estimation is commonly used in "top-down" inversion systems for regularization purposes. Ensemble Kalman filter-based inversion algorithms usually weigh a priori flux to the background or directly replace the background with the a priori flux. However, the "bottom-up" flux estimations, especially the simulated terrestrial-atmosphere CO2 exchange, are usually systematically biased at different spatiotemporal scales because of the deficiencies in understanding of some underlying processes. Here, we introduced a novel regularization algorithm into the Carbon in Ocean--Land--Atmosphere (COLA) data assimilation system, which assimilates a priori information as a unique observation (AAPO). The a priori information is not limited to "bottom-up" flux estimation. With the comprehensive assimilation regularization approach, COLA can apply the spatial gradient of the "bottom-up" flux estimation as a priori information to reduce the bias impact and enhance the dynamic information concerning the a priori "bottom-up" flux estimation. Benefiting from the enhanced signal-to-noise ratio in the spatial gradient, the global, regional, and grided flux estimations using the AAPO algorithm are significantly better than those obtained by the traditional regularization approach, especially over highly uncertain tropical regions in the context of observing simulation system experiments (OSSEs). We suggest that the AAPO algorithm can be applied to other greenhouse gas (e.g., CH4, NO2) and pollutant data assimilation studies. [ABSTRACT FROM AUTHOR]

Atmospheric inversion of carbon dioxide (CO 2) measurements to better understand carbon sources and sinks has made great progress over the last 2 decades. However, most of the studies, including a four-dimensional variational ensemble Kalman filter and Bayesian synthesis approaches, directly obtain only fluxes, while CO 2 concentration is derived with the forward model as part of a post-analysis. Kang et al. (2012) used the local ensemble transform Kalman filter (LETKF), which updates the CO 2 , surface carbon flux (SCF), and meteorology fields simultaneously. Following this track, a system with a short assimilation window and a long observation window was developed (Liu et al., 2019). However, this data assimilation system faces the challenge of maintaining carbon mass conservation. To overcome this shortcoming, here we apply a constrained ensemble Kalman filter (CEnKF) approach to ensure the conservation of global CO 2 mass. After a standard LETKF procedure, an additional assimilation is used to adjust CO 2 at each model grid point and to ensure the consistency between the analysis and the first guess of the global CO 2 mass. Compared to an observing system simulation experiment without mass conservation, the CEnKF significantly reduces the annual global SCF bias from ∼ 0.2 to less than 0.06 Gt and slightly improves the seasonal and annual performance over tropical and southern extratropical regions. We show that this system can accurately track the spatial distribution of annual mean SCF. And the system reduces the seasonal flux root mean square error from a priori to analysis by 48 %–90 %, depending on the continental region. Moreover, the 2015–2016 El Niño impact is well captured with anomalies mainly in the tropics. [ABSTRACT FROM AUTHOR]

There are widespread concerns that current trends in resource-use are unsustainable, but possibilities of overshoot/collapse remain controversial. Collapses have occurred frequently in history, often followed by centuries of economic, intellectual, and population decline. Many different natural and social phenomena have been invoked to explain specific collapses, but a general explanation remains elusive. In this paper, we build a human population dynamics model by adding accumulated wealth and economic inequality to a predator–prey model of humans and nature. The model structure, and simulated scenarios that offer significant implications, are explained. Four equations describe the evolution of Elites, Commoners, Nature, and Wealth. The model shows Economic Stratification or Ecological Strain can independently lead to collapse, in agreement with the historical record. The measure “Carrying Capacity” is developed and its estimation is shown to be a practical means for early detection of a collapse. Mechanisms leading to two types of collapses are discussed. The new dynamics of this model can also reproduce the irreversible collapses found in history. Collapse can be avoided, and population can reach a steady state at maximum carrying capacity if the rate of depletion of nature is reduced to a sustainable level and if resources are distributed equitably. [Copyright &y& Elsevier]

Atmospheric inversion of carbon dioxide (CO2) measurements to understand carbon sources and sinks has made great progress over the last two decades. However, most of the studies, including four-dimension variational (4D-Var), Ensemble Kalman filter (EnKF), and Bayesian synthesis approaches, obtains directly only fluxes while CO2 concentration is derived with the forward model as post-analysis. Kang et al. (2012) used the Local Ensemble Transform Kalman Filter (LETKF) that updates the CO2, surface carbon fluxes (SCF), and meteorology field simultaneously. Following this track, a system with a short assimilation window and a long observation window was developed (Liu et al., 2019). However, this system faces the challenge of maintaining global carbon mass. To overcome this shortcoming, here we introduce a Constrained Ensemble Kalman Filter (CEnKF) approach to ensure the conservation of global CO2 mass. After a standard LETKF procedure, an additional assimilation process is applied to adjust CO2 at each model grid point and to ensure the consistency between the analysis and the first guess of global CO2 mass. In the context of observing system simulation experiments (OSSEs), we show that the CEnKF can significantly reduce the annual global SCF bias from ~0.2 gigaton to less than 0.06 gigaton by comparing between experiments with and without it. Moreover, the annual bias over most continental regions is also reduced. At the seasonal scale, the improved system reduced the flux root-mean-square error from priori to analysis by 48-90 %, depending on the continental region. Moreover, the 2015-2016 El Nino impact is well captured with anomalies mainly in the tropics. [ABSTRACT FROM AUTHOR]

Oscillatory modes of the climate system are among its most predictable features, especially at intraseasonal time scales. These oscillations can be predicted well with data-driven methods, often with better skill than dynamical models. However, since the oscillations only represent a portion of the total variance, a method for beneficially combining oscillation forecasts with dynamical forecasts of the full system was not previously known. We introduce Ensemble Oscillation Correction (EnOC), a general method to correct oscillatory modes in ensemble forecasts from dynamical models. We compute the ensemble mean—or the ensemble probability distribution—with only the best ensemble members, as determined by their discrepancy from a data-driven forecast of the oscillatory modes. We also present an alternate method that uses ensemble data assimilation to combine the oscillation forecasts with an ensemble of dynamical forecasts of the system (EnOC-DA). The oscillatory modes are extracted with a time series analysis method called multichannel singular spectrum analysis (M-SSA), and forecast using an analog method. We test these two methods using chaotic toy models with significant oscillatory components and show that they robustly reduce error compared to the uncorrected ensemble. We discuss the applications of this method to improve prediction of monsoons as well as other parts of the climate system. We also discuss possible extensions of the method to other data-driven forecasts, including machine learning. [ABSTRACT FROM AUTHOR]

The NOAA Science Advisory Board appointed a task force to prepare a white paper on the use of observing system simulation experiments (OSSEs). Considering the importance and timeliness of this topic and based on this white paper, here we briefly review the use of OSSEs in the United States, discuss their values and limitations, and develop five recommendations for moving forward: national coordination of relevant research efforts, acceleration of OSSE development for Earth system models, consideration of the potential impact on OSSEs of deficiencies in the current data assimilation and prediction system, innovative and new applications of OSSEs, and extension of OSSEs to societal impacts. OSSEs can be complemented by calculations of forecast sensitivity to observations, which simultaneously evaluate the impact of different observation types in a forecast model system. [ABSTRACT FROM AUTHOR]

This paper describes the lessons learned from the implementation of a regional ensemble data assimilation and forecast system during the intensive observing period of the Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign (central Argentina, November–December 2018). This system is based on the coupling of the Weather Research and Forecasting (WRF) model and the Local Ensemble Transform Kalman Filter (LETKF). It combines multiple data sources both global and locally available like high-resolution surface networks, AMDAR data from local aircraft flights, soundings, AIRS retrievals, high-resolution GOES-16 wind estimates, and local radar data. Hourly analyses with grid spacing of 10 km are generated along with warm-start 36-h ensemble-forecasts, which are initialized from the rapid refresh analyses every three hours. A preliminary evaluation shows that a forecast error reduction is achieved due to the assimilated observations. However, cold-start forecasts initialized from the Global Forecasting System Analysis slightly outperform the ones initialized from the regional assimilation system discussed in this paper. The system uses a multi-physics approach, focused on the use of different cumulus and planetary boundary layer schemes allowing us to conduct an evaluation of different model configurations over central Argentina. We found that the best combinations for forecasting surface variables differ from the best ones for forecasting precipitation, and that differences among the schemes tend to dominate the forecast ensemble spread for variables like precipitation. Lessons learned from this experimental system are part of the legacy of the RELAMPAGO field campaign for the development of advanced operational data assimilation systems in South America. • A LETKF-WRF system was run successfully in real-time to support RELAMPAGO operations. • A reduction in forecast error was shown due to data assimilation cycles. • 60-member RRR analyses and forecasts are available for the research community. [ABSTRACT FROM AUTHOR]

The Ensemble Mars Atmosphere Reanalysis System (EMARS) dataset version 1.0 contains hourly gridded atmospheric variables for the planet Mars, spanning Mars Year (MY) 24 through 33 (1999 through 2017). A reanalysis represents the best estimate of the state of the atmosphere by combining observations that are sparse in space and time with a dynamical model and weighting them by their uncertainties. EMARS uses the Local Ensemble Transform Kalman Filter (LETKF) for data assimilation with the GFDL/NASA Mars Global Climate Model (MGCM). Observations that are assimilated include the Thermal Emission Spectrometer (TES) and Mars Climate Sounder (MCS) temperature retrievals. The dataset includes gridded fields of temperature, wind, surface pressure, as well as dust, water ice, CO 2 surface ice and other atmospheric quantities. Reanalyses are useful for both science and engineering studies, including investigations of transient eddies, the polar vortex, thermal tides and dust storms, and during spacecraft operations., (© 2019 The Authors. Geoscience Data Journal published by Royal Meteorological Society and John Wiley & Sons Ltd.)

Wind and solar farms offer a major pathway to clean, renewable energies. However, these farms would significantly change land surface properties, and, if sufficiently large, the farms may lead to unintended climate consequences. In this study, we used a climate model with dynamic vegetation to show that large-scale installations of wind and solar farms covering the Sahara lead to a local temperature increase and more than a twofold precipitation increase, especially in the Sahel, through increased surface friction and reduced albedo. The resulting increase in vegetation further enhances precipitation, creating a positive albedo-precipitation-vegetation feedback that contributes ~80% of the precipitation increase for wind farms. This local enhancement is scale dependent and is particular to the Sahara, with small impacts in other deserts., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

The Chinese National Forest Inventory (NFI) has reported increased forest coverage in China since 2000, however, the new satellite-based dataset Global Forest Change (GFC) finds decreased forest coverage. In this study, four satellite datasets are used to investigate this discrepancy in forest cover change estimates in China between 2000 and 2013: forest cover change estimated from MODIS Normalized Burn Ratio (NBR), existing MODIS Land Cover (LC) and Vegetation Continuous Fields (VCF) products, and the Landsat-based GFC. Among these satellite datasets, forest loss shows much better agreement in terms of total change area and spatial pattern than do forest gain. The net changes in forest cover as a proportion of China's land area varied widely from increases of 1.56% in NBR, 1.93% in VCF, and 3.40% in LC to a decline of -0.40% in GFC. The magnitude of net forest increase derived from MODIS datasets (1.56-3.40%) is lower than that reported in NFI (3.41%). Algorithm parameters, different spatial resolutions, and inconsistent forest definitions could be important sources of the discrepancies. Although several MODIS datasets support an overall forest increase in China, the direction and magnitude of net forest change is still unknown due to the large uncertainties in satellite-derived estimates.

The second West African Monsoon Modeling and Evaluation Project Experiment (WAMME II) is designed to improve understanding of the possible roles and feedbacks of sea surface temperature (SST), land use land cover change (LULCC), and aerosols forcings in the Sahel climate system at seasonal to decadal scales. The project's strategy is to apply prescribed observationally based anomaly forcing, i.e., "idealized but realistic" forcing, in simulations by climate models. The goal is to assess these forcings' effects in producing/amplifying seasonal and decadal climate variability in the Sahel between the 1950s and the 1980s, which is selected to characterize the great drought period of the last century. This is the first multi-model experiment specifically designed to simultaneously evaluate such relative contributions. The WAMME II models have consistently demonstrated that SST forcing is a major contributor to the 20 th century Sahel drought. Under the influence of the maximum possible SST forcing, the ensemble mean of WAMME II models can produce up to 60% of the precipitation difference during the period. The present paper also addresses the role of SSTs in triggering and maintaining the Sahel drought. In this regard, the consensus of WAMME II models is that both Indian and Pacific Ocean SSTs greatly contributed to the drought, with the former producing an anomalous displacement of the Intertropical Convergence Zone (ITCZ) before the WAM onset, and the latter mainly contributes to the summer WAM drought. The WAMME II models also show that the impact of LULCC forcing on the Sahel climate system is weaker than that of SST forcing, but still of first order magnitude. According to the results, under LULCC forcing the ensemble mean of WAMME II models can produces about 40% of the precipitation difference between the 1980s and the 1950s. The role of land surface processes in responding to and amplifying the drought is also identified. The results suggest that catastrophic consequences are likely to occur in the regional Sahel climate when SST anomalies in individual ocean basins and in land conditions combine synergistically to favor drought.

The atmospheric carbon dioxide (CO2) record displays a prominent seasonal cycle that arises mainly from changes in vegetation growth and the corresponding CO2 uptake during the boreal spring and summer growing seasons and CO2 release during the autumn and winter seasons. The CO2 seasonal amplitude has increased over the past five decades, suggesting an increase in Northern Hemisphere biospheric activity. It has been proposed that vegetation growth may have been stimulated by higher concentrations of CO2 as well as by warming in recent decades, but such mechanisms have been unable to explain the full range and magnitude of the observed increase in CO2 seasonal amplitude. Here we suggest that the intensification of agriculture (the Green Revolution, in which much greater crop yield per unit area was achieved by hybridization, irrigation and fertilization) during the past five decades is a driver of changes in the seasonal characteristics of the global carbon cycle. Our analysis of CO2 data and atmospheric inversions shows a robust 15 per cent long-term increase in CO2 seasonal amplitude from 1961 to 2010, punctuated by large decadal and interannual variations. Using a terrestrial carbon cycle model that takes into account high-yield cultivars, fertilizer use and irrigation, we find that the long-term increase in CO2 seasonal amplitude arises from two major regions: the mid-latitude cropland between 25° N and 60° N and the high-latitude natural vegetation between 50° N and 70° N. The long-term trend of seasonal amplitude increase is 0.311 ± 0.027 per cent per year, of which sensitivity experiments attribute 45, 29 and 26 per cent to land-use change, climate variability and change, and increased productivity due to CO2 fertilization, respectively. Vegetation growth was earlier by one to two weeks, as measured by the mid-point of vegetation carbon uptake, and took up 0.5 petagrams more carbon in July, the height of the growing season, during 2001-2010 than in 1961-1970, suggesting that human land use and management contribute to seasonal changes in the CO2 exchange between the biosphere and the atmosphere.