Toward Integrated Seasonal Predictions of Land and Ocean Carbon Flux: Lessons Learned From NASA’s Subseasonal-to-Seasonal Predictions
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Simulations by six Coupled Model Intercomparison Project Phase 6 (CMIP6) Earth system models indicate that the seasonal cycle of baseline tropospheric ozone at northern midlatitudes has been shifting since the mid-20th century. Beginning in ∼ 1940, the magnitude of the seasonal cycle increased by ∼10 ppb (measured from seasonal minimum to maximum), and the seasonal maximum shifted to later in the year by about 3 weeks. This shift maximized in the mid-1980s, followed by a reversal – the seasonal cycle decreased in amplitude and the maximum shifted back to earlier in the year. Similar changes are seen in measurements collected from the 1970s to the present. The timing of the seasonal cycle changes is generally concurrent with the rise and fall of anthropogenic emissions that followed industrialization and the subsequent implementation of air quality emission controls. A quantitative comparison of the temporal changes in the ozone seasonal cycle at sites in both Europe and North America with the temporal changes in ozone precursor emissions across the northern midlatitudes found a high degree of similarity between these two temporal patterns. We hypothesize that changing precursor emissions are responsible for the shift in the ozone seasonal cycle; this is supported by the absence of such seasonal shifts in southern midlatitudes where anthropogenic emissions are much smaller. We also suggest a mechanism by which changing emissions drive the changing seasonal cycle: increasing emissions of NOx allow summertime photochemical production of ozone to become more important than ozone transported from the stratosphere, and increasing volatile organic compounds (VOCs) lead to progressively greater photochemical ozone production in the summer months, thereby increasing the amplitude of the seasonal ozone cycle. Decreasing emissions of both precursor classes then reverse these changes. The quantitative parameter values that characterize the seasonal shifts provide useful benchmarks for evaluating model simulations, both against observations and between models.
Dry season length strongly influences tropical rainforest vegetation and is largely determined by precipitation patterns. Over the Amazon, the dry season length has increased since 1979 and severe short-term droughts have occurred. However, similar changes have not been investigated for the world’s second largest rainforest, the Congo Basin, where long-term drying and large-scale declines in forest greenness and canopy water content were reported. Here we present observational evidence for widespread increases in the boreal summer (June–August) dry season length over the Congo Basin since the 1980s, from both hydrological and ecological perspectives. We analyzed both dry season onset and dry season end via multiple independent precipitation and satellite derived vegetation datasets for the period 1979–2015. The dry season length increased by 6.4–10.4 days per decade in the period 1988–2013, primarily attributed to an earlier dry season onset and a delayed dry season end. The earlier dry season onset was caused by long-term droughts due to decreased rainfall in the pre-dry season (April–June). The delayed dry season end resulted from insufficiently replenished soil moisture, which postpones the start of the next wet season and hinders vegetation regrowth. If such changes continue, the enhanced water stress in a warming climate may affect the carbon cycle and alter the composition and structure of evergreen rainforest.
Local agricultural production is a key element of food security in many agricultural countries in Africa. Climate change and variability is likely to adversely affect these countries, particularly as they affect the ability of smallholder farmers to raise enough food to feed themselves. Seasonality influences farmers' decisions about when to sow and harvest, and ultimately the success or failure of their crops. At a 2009 conference in the United Kingdom hosted by the Institute of Development Studies, Jennings and Magrath (2009) described farmer reports from East Asia, South Asia, Southern Africa, East Africa and Latin America. Farmers indicate significant changes in the timing of rainy seasons and the pattern of rains within seasons, including: More erratic rainfall, coming at unexpected times in and out of season; Extreme storms and unusually intense rainfall are punctuated by longer dry spells within the rainy season; Increasing uncertainty as to the start of rainy seasons in many areas; Short or transitional second rainy seasons are becoming stronger than normal or are disappearing altogether. These farmer perceptions of change are striking in that they are geographically widespread and are remarkably consistent across diverse regions (Jennings and Magrath, 2009). The impact of these changes on farmers with small plots and few resources is large. Farming is becoming riskier because of heat stress, lack of water, pests and diseases that interact with ongoing pressures on natural resources. Lack of predictability in the start and length of the growing season affects the ability of farmers to invest in appropriate fertilizer levels or improved, high yielding varieties. These changes occur at the same time as the demand for food is rising and is projected to continue to rise for the next fifty years (IAASTD, 2008). Long-term data records derived from satellite remote sensing can be used to verify these reports, providing necessary analysis and documentation required to plan effective adaptation strategies. Remote sensing data can also provide some understanding of the spatial extent of these changes and whether they are likely to continue. Given the agricultural nature of most economies on the African continent, agricultural production continues to be a critical determinant of both food security and economic growth (Funk and Brown, 2009). Crop phenological parameters, such as the start and end of the growing season, the total length of the growing season, and the rate of greening and senescence are important for planning crop management, crop diversification, and intensification. The World Food Summit of 1996 defined food security as: "when all people at all times have access to sufficient, safe, nutritious food to maintain a healthy and active life". Food security roughly depends on three factors: 1) availability of food; 2) access to food and 3) appropriate use of food, as well as adequate water and sanitation. The first factor is dependent on growing conditions and weather and climate. In a previous paper we have investigated this factor by evaluating the effect of large scale climate oscillation on land surface phenology (Brown et al., 2010). We found that all areas in Africa are significantly affected by at least one type of large scale climate oscillations and concluded that these somewhat predictable oscillations could perhaps be used to forecast agricultural production. In addition, we have evaluated changes in agricultural land surface phenology over time (Brown et al., 2012). We found that land surface phenology models, which link large-scale vegetation indices with accumulated humidity, could successfully predict agricultural productivity in several countries around the world. In this chapter we are interested in the effect of variability in peak timing of the growing season, or phenology, on the second factor of food security, food access. In this chapter we want to determine if there is a link between market prices and land surface phenology and to determine which markets are vulnerable to land surface phenology changes and variability and which market prices are not correlated.
Vegetation growing season and maximum photosynthetic state determine spatiotemporal variability of seasonal total gross primary productivity of vegetation. Recent warming induced impacts accelerate shifts on growing season and physiological status over Northern vegetated land. Thus, understanding and quantifying these changes are very important. Here, we first investigate how vegetation growing season and maximum photosynthesis state are evolved and how such components contribute on inter-annual variation of seasonal total gross primary productivity. Furthermore, seasonally different response of northern vegetation to changing temperature and water availability is also investigated. We utilized both long-term remotely sensed data to extract larger scale growing season metrics (growing season start, end and duration) and productivity (i.e., growing season summed vegetation index, GSSVI) for answering these questions. We find that regionally diverged growing season shift and maximum photosynthetic state contribute differently characterized productivity inter-annual variability and trend. Also seasonally different response of vegetation gives different view of spatially varying interaction between vegetation and climate. These results highlight spatially and temporally varying vegetation dynamics and are reflective of biome-specific responses of northern vegetation to changing climate.
There are many pathways for an air parcel to travel from the troposphere to the stratosphere, each of which takes different time. The distribution of all the possible transient times, i.e. the stratospheric age spectrum, contains important information on transport characteristics. However, it is computationally very expensive to compute seasonally varying age spectra, and previous studies have focused mainly on the annual mean properties of the age spectra. To date our knowledge of the seasonality of the stratospheric age spectra is very limited. In this study we investigate the seasonal variations of the stratospheric age spectra in the Goddard Earth Observing System Chemistry Climate Model (GEOSCCM). We introduce a method to significantly reduce the computational cost for calculating seasonally dependent age spectra. Our simulations show that stratospheric age spectra in GEOSCCM have strong seasonal cycles and the seasonal cycles change with latitude and height. In the lower stratosphere extratropics, the average transit times and the most probable transit times in the winter/early spring spectra are more than twice as old as those in the summer/early fall spectra. But the seasonal cycle in the subtropical lower stratosphere is nearly out of phase with that in the extratropics. In the middle and upper stratosphere, significant seasonal variations occur in the sUbtropics. The spectral shapes also show dramatic seasonal change, especially at high latitudes. These seasonal variations reflect the seasonal evolution of the slow Brewer-Dobson circulation (with timescale of years) and the fast isentropic mixing (with timescale of days to months).
Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.
There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.
Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20°S to 20°N) the spring onset of elevated CH4 emissions starts 3 d earlier, and the CH4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.
The length of the sea ice season summarizes in one number the ice coverage conditions for an individual location for an entire year. It becomes a particularly valuable variable when mapped spatially over a large area and examined for regional and interannual differences, as is done here for the Southern Ocean over the years 1979-1986, using the satellite passive microwave data of the Nimbus 7 scanning multichannel microwave radiometer. Three prominent geographic anomalies in ice season lengths occur consistently in each year of the data set, countering the general tendency toward shorter ice seasons from south to north: (1) in the Weddell Sea the tendency is toward shorter ice seasons from southwest to northeast, reflective of the cyclonic ice/atmosphere/ocean circulations in the Weddell Sea region. (2) Directly north of the Ross Ice Shelf anomalously short ice seasons occur, lasting only 245-270 days, in contrast to the perennial ice coverage at comparable latitudes in the southern Bellingshausen and Amundsen Seas and in the western Weddell Sea. The short ice season off the Ross Ice Shelf reflects the consistently early opening of the ice cover each spring, under the influence of upwelling along the continental slope and shelf and atmospheric forcing from winds blowing off the Antarctic continent. (3) In the southern Amundsen Sea, anomalously short ice seasons occur adjacent to the coast, owing to the frequent existence of coastal polynyas off the many small ice shelves bordering the sea. Least squares trends in the ice season lengths over the 1979-1986 period are highly coherent spatially, with overall trends toward shorter ice seasons in the northern Weddell and Bellingshausen seas and toward longer ice seasons in the Ross Sea, around much of East Antarctica, and in a portion of the south central Weddell Sea.
As a representative site of the southern African biomass-burning region, sun-sky data from the 15 year Aerosol Robotic Network (AERONET) deployment at Mongu, Zambia, was analyzed. For the biomass-burning season months (July-November), we investigate seasonal trends in aerosol single scattering albedo (SSA), aerosol size distributions, and refractive indices from almucantar sky scan retrievals. The monthly mean single scattering albedo at 440 nm in Mongu was found to increase significantly from approx.. 0.84 in July to approx. 0.93 in November (from 0.78 to 0.90 at 675 nm in these same months). There was no significant change in particle size, in either the dominant accumulation or secondary coarse modes during these months, nor any significant trend in the Angstrom exponent (440-870 nm; r(exp 2) = 0.02). A significant downward seasonal trend in imaginary refractive index (r(exp 2) = 0.43) suggests a trend of decreasing black carbon content in the aerosol composition as the burning season progresses. Similarly, burning season SSA retrievals for the Etosha Pan, Namibia AERONET site also show very similar increasing single scattering albedo values and decreasing imaginary refractive index as the season progresses. Furthermore, retrievals of SSA at 388 nm from the Ozone Monitoring Instrument satellite sensor show similar seasonal trends as observed by AERONET and suggest that this seasonal shift is widespread throughout much of southern Africa. A seasonal shift in the satellite retrieval bias of aerosol optical depth from the Moderate Resolution Imaging Spectroradiometer collection 5 dark target algorithm is consistent with this seasonal SSA trend since the algorithm assumes a constant value of SSA. Multi-angle Imaging Spectroradiometer, however, appears less sensitive to the absorption-induced bias.
The quantitative assessment of cloud cover and atmospheric constituents improves our ability to exploit the climate feedback into the Amazon basin. In the 21st century, three droughts have already occurred in the Amazonia (e.g. 2005, 2010, 2015), inducing regional changes in the seasonal patterns of atmospheric constituents. In addition to climate, the atmospheric dynamic and attenuation properties are long-term challenges for satellite-based remote sensing of this ecosystem: high cloudiness, abundant water vapor content and biomass burning season. Therefore, while climatology analysis supports the understanding of atmospheric variability and trends, it also offers valuable insights for remote sensing applications. In this study, we evaluate the seasonal and interannual variability of cloud cover and atmospheric constituents (aerosol loading, water vapor and ozone content) over the Amazon basin, with focus on both climate analysis and remote sensing implications. We take the advantage of new atmosphere daily products at 1 km resolution derived from Multi-Angle Implementation for Atmospheric Correction (MAIAC) algorithm developed for Moderate Resolution Imaging Spectroradiometer (MODIS) data. An intercomparison of Aerosol Robotic Network (AERONET) and MAIAC aerosol optical depth (AOD) and columnar water vapor (CWV) showed quantitative information with a correlation coefficient higher than 0.81. Our results show distinct regional patterns of cloud cover across the Amazon basin: northwestern region presets a persistent cloud cover (> 80%) throughout the year, while low cloud cover (0-20%) occurs in the southern Amazon during the dry season. The cloud-free period in the southern Amazon is followed by an increase in the atmospheric burden due to fire emissions. Our results reveal that AOD records are changing in terms of area and intensity. During the 2005 and 2010 droughts, the positive AOD anomalies (δ > 0.1) occurred over 39.03% (240.3 million ha) and 27.14% (165.99 million ha) of total basin in the SON season, respectively. In contrast, the recent 2015 drought occurred towards the end of year (October through December) and these anomalies were observed over 23.72% (145 million ha) affecting areas in the central and eastern Amazon unlike previous droughts. The water vapor presents high concentration values (4.0-5.0 g/sq.cm) in the wet season (DJF), while we observed a strong spatial gradient from northwestern to southeastern of the basin during the dry season. In addition, we also found a positive trend of water vapor content ( ∼0.3 g/sq.cm) between 2000 and 2015. The total ozone typically varies between 220 and 270 DU, and it has a seasonal change of ∼25-35 DU from wet season to dry season caused by large emissions of ozone precursors and long-range transport. Finally, while this study contributes to climatological analysis of atmospheric constituents, the remote sensing users can also understand the regional constraints caused by atmospheric attenuation, such as high aerosol loading and cloud obstacles for surface observations.
There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.
Evidence-based responses to climate change by society require operational and sustained information including biophysical indicator systems that provide up-to-date measures of trends and patterns against historical baselines. Two key components linking anthropogenic climate change to impacts on socio-ecological systems are the periodic inter- and intra-annual variations in physical climate systems (seasonality) and in plant and animal life cycles (phenology). We describe a set of national indicators that reflect sub-seasonal to seasonal drivers and responses of terrestrial physical and biological systems to climate change and variability at the national scale. Proposed indicators and metrics include seasonality of surface climate conditions (e.g., frost and freeze dates and durations), seasonality of freeze/thaw in freshwater systems (e.g., timing of stream runoff and durations of lake/river ice), seasonality in ecosystem disturbances (e.g., wildfire season timing and duration), seasonality in vegetated land surfaces (e.g., green-up and brown-down of landscapes), and seasonality of organismal life-history stages (e.g., timings of bird migration). Recommended indicators have strong linkages to variable and changing climates, include abiotic and biotic responses and feedback mechanisms, and are sufficiently simple to facilitate communication to broad audiences and stakeholders interested in understanding and adapting to climate change.
Climate change will affect both the mean state and seasonality of marine physical and biogeochemical properties, with important implications for the oceanic sink of atmospheric CO2. Here, we investigate the seasonal cycle of the air-sea exchange of CO2 and pCO2,sw (surface seawater pCO2) and their long term changes using the CMIP6 submission of the NASA-GISS modelE (GISS-E2.1-G). In comparison to the CMIP5 submission (GISS-E2-R), we find that on the global scale, the seasonal cycles of the CO2 flux and NPP have improved, while the seasonal cycles of DIC, alkalinity, and macronutrients have deteriorated. Moreover, for all ocean biogeochemietry fields, changes in skill between E2.1-G and E2-R display large regional variability. For E2.1-G, we find similar modeled and observed CO2 flux seasonal cycles in the subtropical gyres, where seasonal anomalies of pCO2,sw and the flux are temperature-driven, and the Southern Ocean, where anomalies are DIC-driven. Biases in these seasonal cycles are largest in the subpolar and equatorial regions, driven by a combination of biases in temperature, DIC, alkalinity, and wind speed. When comparing the historical simulation to a simulation with an idealized increase in atmospheric pCO2, we find that the seasonal amplitudes of the CO2 flux and pCO2,sw generally increase. These changes are produced by increases in the sensitivity of pCO2,sw to its respective drivers. These findings are consistent with the notion that the sensitivity of pCO2,sw is expected to increase due to the increase of atmospheric pCO2, with changes in the seasonality of temperature, DIC, and alkalinity having secondary influences.
At most latitudes, the seasonal cycle of zonal-mean surface air temperature is notably asymmetric: the length of the warming season is not equal to the length of the cooling season. The asymmetry varies spatially, with the cooling season being ∼40 days shorter than the warming season in the subtropics and the warming season being ∼100 days shorter than the cooling season at the poles. Furthermore, the asymmetry differs between the Northern Hemisphere and the Southern Hemisphere. Here, we show that these observed features are broadly captured in a simple model for the evolution of temperature forced by realistic insolation. The model suggests that Earth's orbital eccentricity largely determines the hemispheric contrast, and obliquity broadly dictates the meridional structure. Clouds, atmospheric heat flux convergence, and time-invariant effective surface heat capacity have minimal impacts on seasonal asymmetry. This simple, first-order picture has been absent from previous discussions of the surface temperature seasonal cycle.
Surface pressure measurement performed by the Viking landers show substantial variations in pressure on seasonal timescales that are characterized by two local minima and two local maxima. These variations have widely been attributed to the seasonal condensation and sublimation of CO2 in the two polar regions. It has been somewhat of a surprise that the amplitude of the minimum and maximum that is dominated by the CO2 cycle in the north was much weaker than the corresponding amplitude of the south-dominated extrema. Another surprise was that the seasonal pressure cycle during years 2 and 3 of the Viking mission was so similar to that for year 1, despite the occurrence of two global dust storms during year 1 and none during years 2 and 3. An energy balance model that incorporates dynamical factors from general circulation model (GCM) runs in which the atmospheric dust opacity and seasonal date were systematically varied was used to model the observed seasonal pressure variations. The energy balance takes account of the following processes in determining the rates of CO2 condensation and sublimation at each longitudinal and latitudinal grid point: solar radiation, infrared radiation from the atmosphere and surface, subsurface heat conduction, and atmospheric heat advection. Condensation rates are calculated both at the surface and in the atmosphere. In addition, the energy balance model also incorporates information from the GCM runs on seasonal redistribution of surface pressure across the globe. Estimates of surface temperature of the seasonal CO2 caps were used to define the infrared radiative losses from the seasonal polar caps. The seasonal pressure variations measured at the Viking lander sites were closely reproduced.