text:
Research ArticleENVIRONMENTAL SCIENCE
Improved estimates of ocean heat content from 1960 to 2015
Lijing Cheng1,*,
Kevin E. Trenberth2,
John Fasullo2,
Tim Boyer3,
John Abraham4 and
Jiang Zhu1
1International Center for Climate and Environment Sciences, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China.
2National Center for Atmospheric Research, Boulder, CO 80307, USA.
3National Oceanic and Atmospheric Administration/National Centers for Environmental Information, Silver Spring, MD 20910, USA.
4University of St. Thomas, St. Paul, MN 55105, USA.
↵*Corresponding author. Email: chenglij{at}mail.iap.ac.cn
Science Advances 10 Mar 2017:
Vol. 3, no. 3, e1601545
Lijing Cheng
International Center for Climate and Environment Sciences, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China.
ORCID record for Lijing Cheng
For correspondence: chenglij@mail.iap.ac.cn
Kevin E. Trenberth
National Center for Atmospheric Research, Boulder, CO 80307, USA.
ORCID record for Kevin E. Trenberth
John Fasullo
ORCID record for John Fasullo
Tim Boyer
National Oceanic and Atmospheric Administration/National Centers for Environmental Information, Silver Spring, MD 20910, USA.
University of St. Thomas, St. Paul, MN 55105, USA.
Jiang Zhu
Earth's energy imbalance (EEI) drives the ongoing global warming and can best be assessed across the historical record (that is, since 1960) from ocean heat content (OHC) changes. An accurate assessment of OHC is a challenge, mainly because of insufficient and irregular data coverage. We provide updated OHC estimates with the goal of minimizing associated sampling error. We performed a subsample test, in which subsets of data during the data-rich Argo era are colocated with locations of earlier ocean observations, to quantify this error. Our results provide a new OHC estimate with an unbiased mean sampling error and with variability on decadal and multidecadal time scales (signal) that can be reliably distinguished from sampling error (noise) with signal-to-noise ratios higher than 3. The inferred integrated EEI is greater than that reported in previous assessments and is consistent with a reconstruction of the radiative imbalance at the top of atmosphere starting in 1985. We found that changes in OHC are relatively small before about 1980; since then, OHC has increased fairly steadily and, since 1990, has increasingly involved deeper layers of the ocean. In addition, OHC changes in six major oceans are reliable on decadal time scales. All ocean basins examined have experienced significant warming since 1998, with the greatest warming in the southern oceans, the tropical/subtropical Pacific Ocean, and the tropical/subtropical Atlantic Ocean. This new look at OHC and EEI changes over time provides greater confidence than previously possible, and the data sets produced are a valuable resource for further study.
Earth's energy imbalance
ocean heat content
subsample test
decadal variability
Global warming is driven by Earth's energy imbalance (EEI). The EEI is likely forced to first order by a combination of greenhouse gas and aerosol forcing, which shapes the timing and magnitude of global warming (1). It is also linked to the internal variations of the climate system and episodic volcanic eruptions; the latter may provide episodic strong radiative forcing to the Earth system (2, 3). By definition, radiative forcing is the change in the net radiative flux due to a change in an external driver of climate change, such as greenhouse gas concentrations. More than 90% of EEI is stored in the ocean, increasing ocean heat content (OHC), while the residual heat is manifest in melting of both land and sea ice, and in warming of the atmosphere and land surface (1, 4, 5). It is therefore essential to provide estimates of OHC changes over time with high confidence to improve our knowledge of EEI and its variability (4).
How much has Earth really warmed in recent decades? The magnitude and location of the ocean warming have become an area of active research, because of the large historical uncertainty in estimated OHC changes (3, 6–10). For instance, tracking Earth's heat and ocean heat is one of the key topics of the so-called "global warming hiatus" research surge (11–20). There are several sources of uncertainty in OHC estimates. A detailed discussion of uncertainty can be found in section S1, supplementing the brief overview here. A key source of uncertainty is instrumental in nature, relating to the measurement process; another originates from methodological choices, for example, in dealing with gaps in sampling (6, 8, 21–27). The community has made progress in detecting the systematic errors in expendable bathythermograph (XBT) data and has provided recommendations to correct the associated errors (25). These recommendations have markedly reduced the impact of XBT biases on multidecadal OHC estimation (26). Another major uncertainty arises from insufficient data coverage, mainly during the pre-Argo era (before 2005), that has led to spatial sampling errors in global and regional OHC estimation (28–30).
There have been many efforts to improve estimates of historical OHC changes including reduction of sampling errors using various methodologies (3, 6, 8, 28, 29, 31–37). One group is referred to as mapping methods, which statistically infill data gaps by using available observations, a prior guess, and an adjustment term based on nearby observations and covariance characteristics, because grid points within the ocean are often mutually correlated (38). For example, a widely referenced method used by the National Centers for Environmental Information (NCEI) (34, 35, 39) objectively analyzed gridded temperature anomalies on the basis of the distances between analysis grid boxes and nearby observations, with zero anomalies used as a first guess. Domingues et al. (36) presented the main ocean spatial variability by means of an empirical orthogonal function analysis. Recently, Cheng and Zhu (40) used an ensemble optimal interpolation method (CZ16) combined with covariances from Coupled Model Intercomparison Project phase 5 (CMIP5) multimodel simulations to provide an improved prior guess (section S2). The other group of gap-filling techniques is referred to as dynamic gap filling, which constrains an ocean/climate model by observations via a data assimilation method, for instance, Ocean Reanalysis System 4 (ORAS4) (3). In this case, the prior guess is produced by bringing forward past information using a model. The accuracy of this method can be limited by model error, which varies from model to model (9), although this limitation may be partly mitigated by model bias corrections (for example, based on assessed bias during the Argo era). This study will focus on the first group—that is, mapping methods—because these have been the main method used for historical OHC estimation.
The estimates obtained on the basis of all of these methods can reveal substantial differences in OHC changes on monthly, interannual, and interdecadal time scales (9, 10, 21, 41), although all of them still indicate a robust ocean warming in recent decades. Therefore, designing the best way to infill the data gaps in both space and time while minimizing sampling error is a key task in global OHC assessment that demands a comprehensive analysis. Quantification of the reliability of obtained OHC estimates across temporal and spatial scales is also essential.
The success of a mapping method can be judged by how accurately it reconstructs the full ocean temperature domain. When the global ocean is divided into a monthly 1°-by-1° grid, the monthly data coverage is <10% before 1960, <20% from 1960 to 2003, and <30% from 2004 to 2015 (see Materials and Methods for data information and Fig. 1). Coverage is still <30% during the Argo period for a 1°-by-1° grid because the original design specification of the Argo network was to achieve 3°-by-3° near-global coverage (42). However, we typically use temporal persistence from 1 month to the next as well as spatial covariance (as carried out using all the mapping methods cited above). Because the early observations were mainly carried out on commercial and scientific research vessels, their locations are limited to the regions around developed countries and along shipping routes (6, 28, 43). From 1990 to 1997, the World Ocean Circulation Experiment (44, 45) extended the in situ data to a mostly one-time global network. Since late 1992, altimetry from space has provided essential complementary information on sea surface height. This century (mainly since 2005), the deployment of the Argo floats into a global ocean network markedly increased data coverage (46, 47), with the exception of coastal, marginal sea, and ice-covered regions (4, 48). Now, the Argo community is extending the Argo float deployment under sea ice and into marginal seas (49).
Fig. 1 Fractional coverage of the mapping method used in this study.
(A) Averaged fraction (0 to 700 m) of the global ocean considered sampled for temperature data in each month for this study. (B) Same as (A) but for 700 to 2000 m. The mean fractional coverage of observations is shown in blue if the global ocean is divided into 1°-by-1° grids, and fractional coverage for NCEI yearly and the 5-year method is shown in green and orange, respectively. This figure starts from 1940. (C) Fractional coverage of the global ocean for layers within 0 to 700 m. (D) Same as (C) but for 700 to 2000 m.
In this paper, we extend and improve a recently proposed mapping strategy (CZ16) to provide a complete gridded temperature field for 0- to 2000-m depths from 1960 to 2015. We carefully evaluate the estimate by comparing the fully analyzed temperature signals with those based on a "subsample test," in which subsamples of data in the Argo era are colocated with historical ocean observations and used to assess sampling error. The time scales for which reliable estimates are obtained are quantified. Using this improved reconstruction, we derived an updated historical (1960–2015) ocean energy budget and contributions to Earth's total energy budget.
The efficiency of a mapping method
We first address the fractional coverage of monthly ocean temperature data achieved with reconstructed signals, which is defined as the fraction of total ocean area obtained by the mapping method (50). The large fractional coverage helps to ensure that a near-complete global reconstruction can be reached. A mapping method uses the data only near the analyzed grid within an assumed area to perform the reconstruction at individual grid cells, with the size of the area defined by the influencing radius or spatial correlation length scales. The choice of influencing radius mainly governs the extent of fractional coverage.
When the present method (CZ16; see Materials and Methods and section S2) is applied, the analyzed fractional coverage is larger than 90% from the late 1950s to 2015 from the sea surface to 2000 m (Fig. 1). For comparison, the NCEI method has a fractional coverage of 60 to 90% for its pentadal estimate (compositing 5 years of data) before 1990 and 90% in the upper ocean after 1990 (0 to 700 m). Below 700 m, the coverage for the NCEI method increases from 50 to 80% from 1960 to 1990 for its pentadal estimate and sharply increases from 30 to 90% during the 2000–2005 period for its yearly estimate. The present method has larger fractional coverage than NCEI and most of the other published methods (21) mainly owing to the difference in the assumed influencing radius. For example, the NCEI method (39) uses an 888-km (~8°) radius, and many other methods adopt a similar radius (31–33, 37). These methods take account of the spectrum of ocean spatial variability and are effective owing to the large spatial correlations within 8° (>0.6; Fig. 2). Instead, CZ16 uses 20° for an influencing radius within 0 to 700 m (25° for 700 to 2000 m) by using the relatively weak correlations (0.0 to 0.6; Fig. 2) for longer spatial distances (see section S3.1 for the choice of influencing radius in the CZ16 method). The application of a large influencing radius helps to ensure a high fractional coverage and also filters out the smaller-scale signals (that is, smaller than 20°), resulting in a smooth spatial pattern (an example is shown in fig. S2). However, the ocean shows spatial variability at multiple length scales (33, 38). To improve methods for assessing this variability, we adopted an iterative strategy similar to that adopted by NCEI (39). Three iterative scans were performed successively using 20° (25° below 700 m), 8°, and 4° to effectively include ocean variability across those spatial scales. Quantitative assessment of the iterative strategy (section S3) shows a significant reduction of analysis error. However, because this fractional coverage does not give information regarding how well a mapping method can reconstruct temperature signals, quantificative assessment of the accuracy of the reconstruction is also needed.
Fig. 2 Observational correlations with distance.
(A and C) Zonal-mean (Z) and (B and D) meridional-mean (M) correlation as a function of distance calculated using the ORAS4 reanalysis data. A linear trend and the mean seasonal cycle have been removed when calculating the correlation. The ORAS4 1° by 1° and monthly data were used.
Quantification of the global mean sampling error
The first step in quantifying the global mean sampling error is to investigate its contribution to both global temperature and OHC changes. The sampling error is assessed on the basis of a subsample test, which is defined as the difference between the reconstructed fields and "truth fields." The truth field is taken to be a set of the gridded averaged temperature anomalies during the Argo era, where the anomalies are relative to 1997–2005 climatology (see Materials and Methods). Each truth field is subsampled according to the locations of historical observations and mapped to obtain the reconstructed fields. The sampling error is estimated from the difference between the reconstructed and truth fields (see Materials and Methods).
Figure 3A presents the global averaged sampling errors since the late 1950s for different choices of truth fields, with the uncertainty estimated as 2 SDs (2σ). The global mean error due to the historical sampling is around 0°C from the late 1950s to 2014, which is significantly different from the temperature signals of the truth field (Fig. 3A, stars). This indicates that no significant sampling bias exists despite the marked evolution of the observation system over the past 56 years.
Fig. 3 Global and basin-averaged sampling error compared with reconstructed temperature change.
(A to G) The sampling errors for different truth fields subsampled by historical observation locations are shown as green dots, with the green solid line and the error bar for the mean and ±2σ, respectively. Blue stars denote the 16 different truth fields. The gray line is the reconstructed monthly temperature anomaly time series from 1960 to 2015 with ±2σ interval in gray shading, and the dark line is the time series after a 7-year low-pass filter (see Materials and Methods). The orange curve is the NCEI result, along with a 1-SE bar (dashed orange). (A) Global, (B) tropical/subtropical Pacific, (C) North Pacific, (D) Indian, (E) tropical/subtropical Atlantic, (F) North Atlantic, and (G) southern oceans. (H and I) S/N ratio of the temporal variability of our reconstruction on (H) decadal scales (solid lines) and (I) interannual scales (triangle lines) compared to the sampling error.
Comparison between the sampling errors and ocean temporal variability on different time scales provides an indication of whether the variability from sampling error is detectable. The temperature time series (Fig. 3A, gray lines) for decadal and multidecadal scales exhibit variability that is much larger than the uncertainties induced by sampling error (Fig. 3A, green dots and 2σ error bars). The decadal/multidecadal variability, which is defined as 2σ of the time series after applying a 7-year low-pass filter (see Materials and Methods), is 0.076°C, at least two times larger than the 2σ sampling error, with a signal-to-noise (S/N) ratio larger than 3 (Fig. 3H). This indicates a robust reconstruction of ocean decadal/multidecadal variability.
The temporal variability of ocean temperature change on interannual scales is 0.008°C after 1990, comparable with the sampling error (S/N ratio between 0.5:1 and 3:1; Fig. 3I). Errors in the signals in both the truth and reconstructed fields are unavoidable because of the imperfect removal of sampling error, imperfect quality control processes, the impact of mesoscale signals, and the residuals of instrumental biases. The year 1990 is chosen as a boundary here to better represent interannual variability, because associated errors in the reconstructed field are reduced (but not removed) in the most recent periods owing to improvements in data quality and density. On the other hand, the estimated sampling error also reveals uncertainty from other error sources in recent decades. Therefore, the comparable S/N ratio between sampling error and interannual variability implies that the reconstructions of year-to-year global temperature changes remain uncertain and can easily be affected by observational error. A recent study (10) based on a comparison of month-to-month variations of global OHC for the upper 2000 m with radiative observations at the top of atmosphere (TOA) noted that short-term fluctuations in OHC contain large noise.
An analysis of these temperature reconstructions across different layers (section S3.2) addresses uncertainty in interannual variability estimates for all ocean layers. In the upper ocean (that is, 20 m in fig. S7), the interannual variability is significant and larger than the sampling error (S/N ratio >2) because the upper ocean experiences strong interannual variability mostly dominated by El Niño–Southern Oscillation (ENSO) (3). However, the S/N ratio decreases markedly below 300 m (shown in figs. S8 to S10), and interannual variability is of comparable magnitude to sampling error in the deep ocean, simply because there is much weaker variability (signal) than in the upper 300 m.
Furthermore, the S/N ratio of both decadal and interannual variability increases with time, starting in the 1960s, and the interannual variability becomes slightly stronger than the sampling error during the Argo era (S/N ratio between 2:1 and 5:1), suggesting the need to maintain and extend the integrated ocean observation system. The latter includes the current Argo network in the open ocean, a ship-based observation system, glider and coastal moorings, the coastal regions, the moored array in the tropics, instrumented pinnipeds in polar regions, and modified Argo floats and ice-tethered profilers in ice-covered regions.
Quantification of the regional sampling error
It is also worthwhile examining the sampling error of basin-scale averages because different mechanisms may be responsible for OHC changes in different ocean basins. Moreover, ocean heat transport within and between the basins is of great interest in understanding Earth's key energy flows. Here, we divide the oceans into the North Pacific Ocean (30° to 65°N), the North Atlantic Ocean (30° to 65°N), the tropical and subtropical Pacific Ocean (30°S to 30°N), the tropical and subtropical Atlantic Ocean (30°S to 30°N), the Indian Ocean (30°S to 20°N), and the southern oceans (south of 30°S). There are continual observations for the North Atlantic and North Pacific oceans in middle latitudes from mechanical bathythermographs (1940s to 1960s) and XBTs (late 1960s to 2000), mainly distributed along cruise ship lines. In the tropics, the Tropical Atmosphere Ocean moorings started in 1979, although they were only fully deployed by about 1992. The least sampled area is south of 30°S, especially in the winter half-year. Even in the most recent decade, there are minimal data in the ice-covered regions.
In the six ocean basins, the mean sampling error is approximately zero (Fig. 3), indicating a robust reconstruction of temperature signals. Five ocean basins (not including the Indian Ocean) have experienced significant decadal-scale temperature changes, which are all larger than the 2σ sampling error (S/N ratio is larger than 2), confirming that long-term basin-averaged temperature changes are robust and insignificantly affected by sampling error. In the Indian Ocean (where there is evidently low signal), decadal variability is not significant before 1970, with an S/N ratio between 1.6:1 and 2:1. In particular, the southern oceans and the Atlantic Ocean have experienced the most marked decadal variations, with larger S/N ratios than the other basins. The magnitude of interannual variability in all six basins is comparable to the 2σ sampling error (S/N ratio between 0.3:1 and 5:1; Fig. 3I), as was found for the global ocean.
In addition, a geographical analysis of sampling error is presented in section S3.3. The sampling error for each 1°-by-1° grid has mean errors centered mostly about zero (fig. S11), suggesting that no significant regional bias in the reconstructed field exists over most ocean regions (see also fig. S12 for the 20- and 1600-m layers). However, in boundary current systems, Antarctic Circumpolar Current regions, and the Southern Hemisphere midlatitudes, there are larger 2σ sampling errors than in other regions (fig. S11B).
Global OHC change for the upper 2000 m
On the basis of reconstructed temperature fields and associated error bars, monthly OHC changes within 0 to 700 and 700 to 2000 m (Fig. 4A) show significant warming in the past 56 years. A stronger ocean warming trend has existed since the late 1980s for both 0- to 700- and 700- to 2000-m depths compared with the 1960s to the 1980s. The linear trend of OHC at 0- to 700-m depth is 0.15 ± 0.08 × 1022 J/year (0.09 ± 0.05 W/m2) during 1960–1991 and 0.61 ± 0.04 × 1022 J/year (0.38 ± 0.03 W/m2) during 1992–2015, a warming trend four times stronger than the 1960–1991 period. The linear trend of OHC at 700- to 2000-m depth is 0.04 ± 0.08 × 1022 J/year (0.02 ± 0.05 W/m2) in the period 1960–1991 and 0.37 ± 0.02 × 1022 J/year (0.23 ± 0.02 W/m2) from 1992 to 2015 (nine times stronger than that from 1960–1991) (table S1). This indicates an accelerating heat input into both the 0- to 700-m and 700- to 2000-m layers. The acceleration is most probably linked to the increasing EEI with time. In particular, a new study (51) shows that the recent acceleration is partly due to recovery from the volcanic eruption in 1992 (Mt. Pinatubo), which led to strong ocean cooling.
Fig. 4 Global OHC change time series.
(A) OHC from 0 to 700 m (blue), 700 to 2000 m (red), and 0 to 2000 m (dark gray) from 1955 to 2015 as obtained by this study, with the uncertainty of the ±2σ interval shown in shading. All time series of the new analysis are smoothed by a 12-month running mean filter, relative to the 1997–2005 base period. (B) The new estimate is compared with an independent estimate from NCEI with its SE as dashed lines. Both OHC 0 to 700 m and OHC 700 to 2000 m are shown from 1957 to 2014. The baseline of the time series from NCEI is adjusted to the values of the current analysis within 2005–2014.
The new estimates for OHC 700 to 2000 m show a stronger long-term trend before 2005 than NCEI estimates, with 10 to 20% differences, although the two estimates are not distinguishable within the margin of error (Fig. 4B). For both the upper 700-m and 700- to 2000-m layers, our current OHC estimates show larger long-term changes in the period 2005–2010 than NCEI estimates. To see where the "extra heat" is coming from, we include the NCEI estimates in Fig. 3 (A to G, orange lines) for comparison. The major differences appear in the southern oceans and the tropical/subtropical Pacific Ocean: The current OHC estimate shows quicker warming since 1960 than do NCEI estimates, albeit with sparse data coverage before the Argo era. This is probably due to the prior assumption of zero anomalies used in NCEI, which biases the reconstructed field toward zero in data-sparse regions. Rapid warming in the southern oceans has been observed in all Argo-based data sets since 2004 (10) and in the altimetry sea level (ASL) change since 1993 (fig. S18).
Observational estimates (both NCEI and this study) show a very weak (insignificant) deep (700 to 2000 m) ocean warming before 1980, in apparent contradiction to the strong radiative forcing since the 1960s and many model simulations (1). There are several hypotheses for this moderate deep ocean warming that require a more careful analysis in the future: (i) The warming signal is small during the 1960s and can be influenced by sampling error (note that the S/N ratios during the 1960s are small on global and regional scales in Fig. 3); (ii) there is a robust ocean warming in the upper 700 m, but the warming does not penetrate to the deep ocean below 700 m. By contrast, climate models may be too diffusive and thus may prematurely warm the deep ocean (51–53).
Regional OHC change in recent decades
Regional OHC changes are of great interest to the broad climate community because they provide a basis for understanding the mechanisms of ocean energy flow, especially since 1998 (the so-called "hiatus" period), as previous studies have sometimes shown contradictory results based on different data sets (14, 16, 17, 54, 55). Our studies show that there has been no slowdown in global OHC change since 1998 compared with the previous decade (Fig. 5): There is an acceleration of ocean warming (for both OHC 0 to 700 m and OHC 0 to 2000 m), consistent with a previous study based on model-based data assimilation (3). Moreover, the percentage of 700- to 2000-m ocean heating relative to 0- to 2000-m OHC is 32% (32%) for 1960–1998 (1960–2005), but 38% (49%) for 1998–2015 (2005–2015), indicating that the deep ocean has played an increasingly important role in the ocean energy budget since 1998. The total OHC increase from 1998 to 2015 is 15.2 × 1022 J in the upper 2000 m, with 17% stored in the Pacific Ocean, 24% in the Indian Ocean (30°S northward), 31% in the Atlantic Ocean, and 28% in the southern oceans (south of 30°S). Again, the total OHC change calculated here is not well characterized by a linear trend because of the relatively short time period considered and the presence of strong decadal variability. It is evident that all six ocean basins have experienced significant warming since 1998 but that heat was mainly stored in the southern oceans, the tropical/subtropical Pacific Ocean, and the tropical/subtropical Atlantic Ocean from 1960 to 1998 (Fig. 5). Understanding how this heat has been transported or redistributed in the ocean continues to be an important research topic.
Fig. 5 OHC changes from 1960 to 2015 for different ocean basins.
(A) For 0 to 2000 m, (B) 0 to 700 m, and (C) 700 to 2000 m. All the time series are relative to the 1997–1999 base period and smoothed by a 12-month running filter. The curves are additive, and the OHC changes in different ocean basins are shaded in different colors.
Atlantic Ocean warming shows the largest OHC 0 to 2000 m increase (about 3.5 times larger than the Pacific Ocean), despite having an area only 47% as large as that of the Pacific (Figs. 3 and 5). The OHC change in the North Atlantic Ocean shows strong decadal variability and is likely linked to the strengthening of Atlantic meridional overturning circulation up to the middle 1990s and a subsequent weakening during the 2000s (Fig. 3) (56, 57), although the reasons for these changes are still debated. Other regions with significant decadal variability are the North Pacific Ocean and the tropical/subtropical Pacific Ocean. However, these changes are mainly limited to the upper 700 m (with very weak changes below 700 m; Fig. 5C and figs. S7 to S10), showing the importance of the shallow subtropical circulation (11, 58). The southern oceans and the tropical/subtropical Atlantic Ocean have experienced continuous and monotonic long-term warming since the 1960s, revealing a robust footprint of global warming (59). The Atlantic Ocean and the southern oceans are the major new heat reservoirs (59%) even though their total area is just 48% that of the global ocean. This implies that the meridional overturning circulation is a key mechanism involved in the sequestering of heat in the deeper (below 700 m) ocean.
For the upper ocean (0 to 700 m), the OHC increase from 1998 to 2015 is 10.11 × 1022 J, with 21% occurring in the Pacific Ocean, 21% in the Indian Ocean, 28% in the Atlantic Ocean, and 30% in the southern oceans, generally consistent with the OHC 0 to 2000 m. Our data reveal significant heat input into each basin in the upper ocean in the recent decade (Figs. 3 and 5), including the Indian Ocean (10).
Implications for the planetary energy imbalance
This updated OHC estimate offers the opportunity to revisit the planetary energy imbalance over the historical period, since more than 90% of EEI is stored in the ocean. The revised ocean energy budget (Fig. 6) based on OHC in the upper 0 to 2000 m from this study is combined with OHC below 2000 m adapted from Purkey and Johnson (60) (0.11 ± 0.10 × 1022 J/year), which has been updated by Desbruyères et al. (61). It is extended to 1960 by assuming a zero heating rate before 1990 below 2000 m, because OHC changes are found to be very small for 700 to 2000 m before 1990 (with large uncertainty), consistent with the assumption adopted in previous sea level analyses (62). The OHC 700 to 2000 m in this study (NCEI) indicates a one-ninth (one-third) smaller warming rate for the 1960–1991 versus the 1992–2015 period. Therefore, the maximum temperature change below 2000 m before 1991 should be one-third of its change after 1991 (that is, 0.04 × 1022 J/year). To be conservative, the uncertainty of the OHC change below 2000 m is set to ±0.04 × 1022 J/year.
Fig. 6 Estimate of the ocean energy budget.
The three major volcanic eruptions are marked. The energy budgets are relative to the 1958–1962 base period. The integrated net radiative imbalance from Allan et al. (65) estimated from the TOA is included in yellow and is multiplied by 0.93 to be comparable with the ocean energy budget. The TOA radiation is adjusted to the value of OHC within 2013–2014. The dashed gray lines encompass the 95% confidence interval.
The new result (Fig. 6) suggests a total full-depth ocean warming of 33.5 ± 7.0 × 1022 J (equal to a net heating of 0.37 ± 0.08 W/m2 over the global surface and over the 56-year period) from 1960 to 2015, with 36.5, 20.4, 30.3, and 12.8% contributions from the 0- to 300-m, 300- to 700-m, 700- to 2000-m, and below 2000-m layers, respectively. Here, we prefer to use the total energy change in the budget analyses rather than the linear trend because the change is not linear owing to the general increase in radiative forcing with time. The new reconstruction confirms the previous finding that the upper ocean experiences the most statistically significant warming, while the 0- to 2000-m layer contributes to the vast majority of the ocean warming since 1960 (1). By comparison, the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (1) reported an ocean heat gain of 25.5 ± 6.1 × 1022 J from 1971 to 2010, whereas the results of the current study are somewhat greater at 28.8 ± 4.4 × 1022 J for the same period, at the upper end of their uncertainty range (table S1). Furthermore, the full-depth OHC for 1970–2005 (26.5 ± 4.8 × 1022 J) agrees with a recent analysis (63) based on adjusted observation-based estimates, a simple transparent gap-filling method, and ORAS4 reanalyses (data assimilation method) (28.3 ± 1.8 × 1022 J); it also agrees with the ensemble mean of CMIP5 simulations (26.6 ± 4.4 × 1022 J) (table S1) (63). The consistency among these independent studies indicates the strength of the assessment of the ocean energy budget since 1960.
In addition to warming the oceans, the energy imbalance of the Earth system has also melted ice and warmed the atmosphere and land (4, 5, 64). These events account for ~7% of the EEI (1, 5, 64) but may be slightly smaller in the earlier years (64). Therefore, the total planetary energy imbalance incorporating all appreciable thermal reservoirs since 1960 is 36.0 ± 7.5×1022 J (equal to a net heating of 0.40 ± 0.09 W/m2 over the global surface during the 56-year period).
Although the EEI can best be estimated by OHC changes, the radiation measurements from space provide an important complement (4, 5). Tracking changes in the energy inventory of the planet gives an estimate of the energy imbalance. Figure 6 incorporates the accumulated net downward radiative flux imbalance at the TOA reconstructed by Allan et al. (65) since 1985. Because 7% of EEI is absorbed into other energy reservoirs, we multiply the TOA net radiation by 0.93 to estimate the ocean component in Fig. 6. The EEI based on OHC and TOA radiation shows consistent long-term changes since the late 1980s, further confirming the robustness of this new assessment of Earth's energy inventory (tables S1 and S2). The inconsistency of the interannual variations of EEI based on TOA and OHC reveals the impacts of sampling errors in OHC, deeper OHC change below 2000 m, interannual variations in atmosphere/land/ice heat change, and errors in TOA radiation observations (66). Although there are errors in the various observation systems, comparison and reconciliation among independent observations (that is, OHC, TOA radiation, sea level, and also surface flux data) and models (that is, CMIP) will promote a more complete understanding of climate variability change and provide an opportunity for independent validation (1). Improving the existing integrated observation system, including the ocean and TOA observations, is therefore deemed essential for future progress.
Documentation of a detailed and reliable structure of ocean temperature changes is highly valuable in understanding the mechanisms of ocean heat variability driven by both internal processes and external forcing. This study provides a new reconstruction of the ocean temperature field (0 to 2000 m) and a quantification of its uncertainty on global and regional scales and interannual/decadal time scales based on subsample tests. These tests suggest a small sampling error related to temperature signals on decadal time scales. This new reconstruction provides a basis for accurately detecting and quantifying both regional and global variability across recent decades while minimizing error due to irregular ocean sampling by in situ observations.
This study aims to provide reconstructions on larger spatial and longer time scales (than mesoscale signals, eddies, and internal tides). Small-scale signals have been taken into account by smoothing in two ways: (i) collecting data within a time window for estimates in a given month and (ii) collecting data within the influencing radius, as carried out by most mapping methods (that is, the current method and NCEI), which acts as a spatial smoother. However, the small-scale signals could become errors in the reconstruction if they are spuriously mapped onto the large scale. This is where quality control becomes crucial. It is inherently subjective and based on knowledge of the ocean circulation and temperature/salinity patterns, but it remains the primary means to removing aliased mesoscale signals. Mesoscale eddies can be important in moving heat and fresh water around but do not show up in mean fields. The subsample test used in this study partly takes account of this error because gridded averaged data are used as truth without applying a spatial smoothing. One deficiency of this method is that the historical sampling error can be underestimated. We have tested for this effect and have shown that it does not affect the key conclusions of this study (section S4). It is difficult to fully address this issue because the Argo network does not fully capture mesoscale signals and lacks observations in coastal areas. High-resolution data, such as output from eddy-resolving ocean models, sea surface temperature, ocean color, and high-resolution ASL data (but only after 1993), will be helpful in further investigations of the impact of mesoscale signals, although the instrumental error of in situ profile observations may not be fully represented by these data.
Note that the new OHC analysis indicates only a weak ocean cooling after the Mt. Pinatubo eruption event in 1992 (Figs. 3 and 4), in contrast to the imbalance suggested by satellite data. Our subsample test suggests that interannual variations during the pre-Argo era are of comparable magnitude to sampling errors, but the current Argo system begins to allow for a resolution of interannual OHC variability. Therefore, if this volcanic eruption event were to occur with today's observing network, we expect to be able to resolve its impact on OHC with greater confidence. Further reduction in the error in OHC estimates on interannual scales likely requires international coordination, including extending the current integrated observation system and understanding the influence of quality control processes, instrumental bias corrections, and ocean sampling. In addition, the improvement of methods that synthesize all available observations, such as reanalyses, will continue to be essential.
In situ temperature data from 1960 to 2015 were from the World Ocean Database (WOD) from the NCEI (67) and were prepared following CZ16. The quality flags provided by WOD were used to remove erroneous profiles and measurements. There are 41 standard depths in the upper 2000 m: 1 m, 5 m, 10 to 100 m in 10-m intervals, 120 to 200 m in 20-m intervals, 250 to 900 m in 50-m intervals, and 1000 to 2000 m in 100-m intervals. In the upper 700 m, we used data collected over a 3-month window to calculate monthly means, which helped to increase the data coverage without losing interannual signals, such as ENSO variations. To achieve improved spatial coverage before 2005 below 700 m, the data within 4-, 8-, 12-, and 18-month bins around the selected month were averaged together for the 750- to 1000-m, 1100- to 1500-m, 1600- to 1900-m, and 2000-m ocean layers, respectively, using moving averages. This choice was based on the amount of data available before 2005 in these layers. For models, we assembled 40 outputs of historical runs from 1940 to 2005 and 31 outputs using the Representative Concentration Pathways 4.5 scenario from 2006 to 2015 (68). Here, we examined the temperature reconstruction for 70°S to 65°N. For both observations and models, we constructed the climatology by using the data for the 1997–2005 period (see section S1 for a discussion of climatology).
The mapping method adopted in this study was proposed by Cheng and Zhu (40) and is called the ensemble optimal interpolation method with dynamic ensemble (EnOI-DE) provided by CMIP5 historical simulations. We modified this method to better represent the ocean variability on various spatial scales by using iterative strategies. Three iterative runs were carried out, where the analysis fields of the previous scan provided the prior estimate for the next scan (a thorough description can be found in section S2). Different influencing radii were set in each iterative run, that is, 20° in scan 1 (25° for 700- to 2000-m layers), 8° in scan 2, and 4° in scan 3. Furthermore, the present study extended the method to ocean depths below 700 m, where the observations were sparser than the upper 700 m. The temporal variations of deep ocean temperatures were much smaller relative to the upper ocean (37, 39). Therefore, to obtain a reliable assessment of the decadal and multidecadal variations of ocean temperature below 700 m, it is desirable to use the data within adjacent months to increase the spatial data coverage.
The subsample test was designed as follows. The 1°-by-1° gridded averaged temperature anomalies during the Argo period were used for each selected month (January and August from 2007 to 2014) as truth values (for a total of 16 truth fields). Each truth field was then subsampled according to the history of the observational locations (every 5 years from 1955 to 2014 for January and 1957 to 2014 for August). Thus, 24 subsampled fields were generated and then mapped using the CZ16 method and compared with the truth fields. The data during the Argo era not only observed ocean variability in the real world in each ocean layer but also contained instrumental errors in observations and insufficient coverage in coastal areas. However, because one could never repeat ocean observations made before 2000, the data during the Argo era were the best for examining the performance of the reconstruction for the historical period. One study (50) has subsampled ASL data to examine the sampling error in the OHC estimate. Using ASL has benefits for its complete global ocean coverage, and provides an important insight into the sampling error. However, ASL is a combination of full-depth thermal expansion and ocean mass change and is not directly representative of OHC change. Further, ASL is a two-dimensional metric without depth information, but the ocean sampling changes markedly with depth.
The S/N ratio for the ocean temporal variability compared to the sampling error is defined as
where the signal is defined as twice the SD of the filtered temperature time series (Thistorical), which represents the magnitude of the OHC temporal variation. The noise is also twice the SD of the sampling errors. When S/N > 2, it indicates a significant signal given the sampling error; this is similar to a t test estimate of the >95% confidence level. When the ocean decadal/multidecadal variation was examined, a 7-year low-pass filter was applied to remove the ocean variability with frequencies higher than 1/7 per year. The method filtered a time series via the Lanczos method (cosine filter) (69). When the interannual variations were examined, the 7-year high-pass filter was applied to remove the ocean variability with frequencies lower than 1/7 per year.
Error bars of the OHC estimates (in Figs. 3, 4, and 6) are given as the ±2σ range, two SDs of the ensemble members obtained by the mapping method. ±2σ is equal to an approximate 95% confidence interval when assuming a Gaussian distribution, which means that there is only a 5% chance that the OHC falls outside the given interval. The distribution of ensemble anomalies at 20 and 1200 m is shown in fig. S19; although not strictly Gaussian, ±2σ captures more than 95% of this distribution.
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/3/3/e1601545/DC1
Supplementary Materials and Methods
fig. S1. Illustration of the iterative EnOI-DE/CMIP5 method used in the current study.
fig. S2. An example of the reconstructed fields after three iterative scans.
fig. S3. Reconstruction of temperature field at 1200 m in August 2011 for the historical sampling (in February).
fig. S4. Reconstruction of temperature field at 1200 m in August 2011 for the historical sampling (in September).
fig. S5. Mean temperature error as a function of different choices of the influencing radii between the reconstructed and truth fields.
fig. S6. Six major ocean basins defined in this study.
fig. S7. Global and basin-averaged sampling error compared with reconstructed temperature change at 20 m.
fig. S8. Global and basin-averaged sampling error compared with reconstructed temperature change at 300 m.
fig. S10. Global and basin-averaged sampling error compared with reconstructed temperature change at 1200 m.
fig. S11. Geographical distribution of mean and 2σ sampling error for 0- to 2000-m average.
fig. S12. Geographical distribution of mean and 2σ sampling error in 1°-by-1° grid at 20 and 1600 m.
fig. S13. 2σ sampling error for different scans at different depths.
fig. S14. Global OHC time series for the reconstructions after scan 1, scan 2, and scan 3.
fig. S15. Frequency distribution of temperature anomalies for the years 1986 and 2015.
fig. S16. Sampling error as calculated by two subsample methods.
fig. S17. S/N ratio analysis for two methods of subsample test.
fig. S18. Comparison between OHC and sea level change since 1993.
fig. S19. Distribution of the ensemble anomalies.
table S1. OHC trends obtained in this study for the 1960–1991 and 1992–2015 periods.
table S2. Net OHC and EEI changes obtained in the current study compared with some independent estimates.
M. Rhein, S. R. Rintoul, S. Aoki, E. Campos, D. Chambers, R. A. Feely, S. Gulev, G. C. Johnson, S. A. Josey, A. Kostianoy, C. Mauritzen, D. Roemmich, L. D. Talley, F. Wang, Observations: Ocean, in Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, T. F. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S. K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex, P. M. Midgley, Eds. (Cambridge Univ. Press, 2013).
J. A. Church,
N. J. White,
J. M. Arblaster
, Significant decadal-scale impact of volcanic eruptions on sea level and ocean heat content. Nature 438, 74–77 (2005).
M. A. Balmaseda,
K. E. Trenberth,
E. Källén
, Distinctive climate signals in reanalysis of global ocean heat content. Geophys. Res. Lett. 40, 1754–1759 (2013).
K. von Schuckmann,
M. D. Palmer,
K. Trenberth,
A. Cazenave,
N. Champollion,
J. Hasen,
S. A. Josey,
N. Loeb,
P.-P. Mathieu,
B. Meyssignac,
M. Wild
, An imperative to monitor Earth's energy imbalance. Nat. Clim. Change 6, 138–144 (2016).
J. T. Fasullo,
M. A. Balmaseda
, Earth's energy imbalance. J. Clim. 27, 3129–3144 (2014).
J. P. Abraham,
M. Baringer,
N. L. Bindoff,
T. Boyer,
L. J. Cheng,
J. L. Conroy,
C. M. Domingues,
J. Gilson,
G. Goni,
S. A. Good,
J. M. Gorman,
V. Gouretski,
M. Ishii,
G. C. Johnson,
S. Kizu,
J. M. Lyman,
A. M. Macdonald,
W. J. Minkowycz,
S. E. Moffitt,
A. R. Piola,
F. Reseghetti,
K. Schuckmann,
I. Velicogna,
J. K. Willis
, A review of global ocean temperature observations: Implications for ocean heat content estimates and climate change. Rev. Geophys. 51, 450–483 (2013).
L.-J. Cheng,
J. Zhu,
J. Abraham
, Global upper ocean heat content estimation: Recent progress and the remaining challenges. Atmos. Oceanic Sci. Lett. 8, 333–338 (2015).
V. V. Gouretski,
D. M. Smith,
, Robust warming of the global upper ocean. Nature 465, 334–337 (2010).
C. D. Roberts,
M. Balmaseda,
Y.-S. Chang,
G. Chepurin,
N. Ferry,
Y. Fujii,
S. Guinehut,
K. Haines,
F. Hernandez,
A. Köhl,
T. Lee,
M. J. Martin,
S. Masina,
S. Masuda,
K. A. Peterson,
A. Storto,
T. Toyoda,
M. Valdivieso,
G. Vernieres,
O. Wang,
Y. Xue
, Ocean heat content variability and change in an ensemble of ocean reanalyses. Clim. Dynam. 1–22 (2015).
L. Cheng
, Insights into Earth's energy imbalance from multiple sources. J. Clim. 29, 7495–7505 (2016).
H. M. England,
S. McGregor,
P. Spence,
G. A. Meehl,
A. Timmermann,
W. Cai,
A. S. Gupta,
M. J. McPhaden,
A. Purich,
A. Santoso
, Recent intensification of wind-driven circulation in the Pacific and the ongoing warming hiatus. Nat. Clim. Change 4, 222–227 (2014).
J. T. Fasullo
, An apparent hiatus in global warming? Earth's Future 1, 19–32 (2013).
G. Branstator,
A. S. Phillips
, Seasonal aspects of the recent pause in surface warming. Nat. Clim. Change 4, 911–916 (2014).
X. Chen,
K.-K. Tung
, Varying planetary heat sink led to global-warming slowdown and acceleration. Science 345, 897–903 (2014).
G. Foster, J. Abraham, Lack of evidence for a slowdown in global temperature. Variations (US Climate Variability and Predictability Program (CLIVAR), Washington DC, 2015), pp. 6–9.
L. Cheng,
F. Zheng,
J. Zhu
, Distinctive ocean interior changes during the recent warming slowdown. Sci. Rep. 5, 14346 (2015).
S.-K. Lee,
W. Park,
A. L. Gordon,
B. Huber,
Y. Liu
, Pacific origin of the abrupt increase in Indian Ocean heat content during the warming hiatus. Nat. Geosci. 8, 445–449 (2015).
K. E. Trenberth
, Has there been a hiatus? Science 349, 691–692 (2015).
S. Lewandowsky,
J. S. Risbey,
N. Oreskes
, The "pause" in global warming: Turning a routine fluctuation into a problem for science. Bull. Am. Meteorol. Soc. 97, 723–733 (2015).
T. R. Karl,
A. Arguez,
B. Huang,
J. H. Lawrimore,
J. R. McMahon,
M. J. Menne,
T. C. Peterson,
R. S. Vose,
H.-M. Zhang
, Possible artifacts of data biases in the recent global surface warming hiatus. Science 348, 1469–1472 (2015).
V. Gourestki,
J. K. Willis,
J. Antonov,
S. Wijffels,
R. Cowley,
N. Bindoff
, Sensitivity of global upper-ocean heat content estimates to mapping methods, XBT bias corrections, and baseline climatologies. J. Clim. 29, 4817–4842 (2016).
K. P. Koltermann
, How much is the ocean really warming? Geophys. Res. Lett. 34, L01610 (2007).doi:10.1029/2006GL027834
F. Reseghetti
, On depth and temperature biases in bathythermograph data: Development of a new correction scheme based on analysis of a global ocean database. Deep Sea Res. Pt. I 57, 812–833 (2010).
G. C. Johnson
, Estimating global ocean heat content changes in the upper 1800 m since 1950 and the influence of climatology choice. J. Climate 27, 1945–1957 (2013).
J. Abraham,
S. Dong,
F. Bringas,
M. Goes,
L. Houpert,
J. Sprintall,
, XBT Science: Assessment of instrumental biases and errors. Bull. Am. Meteorol. Soc. 97, 924–933 (2016).
S. Wijffels
, Time, probe type and temperature variable bias corrections to historical expendable bathythermograph observations. J. Atmos. Oceanic Tech. 31, 1793–1825 (2014).
, Influences of the choice of climatology on ocean heat content estimation. J. Atmos. Oceanic Tech. 32, 388–394 (2015).
, Artifacts in variations of ocean heat content induced by the observation system changes. Geophys. Res. Lett. 41, 7276–7283 (2014).
P. J. Durack,
P. J. Gleckler,
F. W. Landerer,
K. E. Taylor
, Quantifying underestimates of long-term upper-ocean warming. Nat. Clim. Change 4, 999–1005 (2014).
S. A. Good
, The impact of observational sampling on time series of global 0–700 m ocean average temperature: A case study. Int. J. Climatol. 10.1002/joc.4654 (2016).
M. Kimoto
, Reevaluation of historical ocean heat content variations with time-varying XBT and MBT depth bias corrections. J. Oceanogr. 65, 287–299 (2009).
J. M. Murphy
, An objective ocean temperature and salinity analysis using covariances from a global climate model. J. Geophys. Res. 112, C02022 (2007).
D. Roemmich,
B. Cornuelle
, Interannual variability in upper ocean heat content, temperature, and thermosteric expansion on global scales. J. Geophys. Res. 109, C12036 (2004).
S. Levitus,
T. Boyer
, Warming of the world ocean, 1955–2003. Geophys. Res. Lett. 32, L02604 (2005).
J. I. Antonov,
T. P. Boyer,
R. A. Locarnini,
H. E. Garcia,
A. V. Mishonov
, Global ocean heat content 1955–2008 in light of recently revealed instrumentation problems. Geophys. Res. Lett. 36, L07608 (2009).
S. E. Wijffels,
P. M. Barker,
J. R. Dunn
, Improved estimates of upper-ocean warming and multi-decadal sea-level rise. Nature 453, 1090–1093 (2008).
N. A. Rayner
, EN4: Quality controlled ocean temperature and salinity profiles and monthly objective analyses with uncertainty estimates. J. Geophys. Res. 118, 6704–6716 (2013).
X. Y. Zang,
C. Wunsch
, Spectral description of low-frequency oceanic variability. J. Phys. Oceanogr. 31, 3073–3095 (2001).
O. K. Baranova,
A. V. Mishonov,
J. R. Reagan,
D. Seidov,
E. S. Yarosh,
M. M. Zweng
, World ocean heat content and thermosteric sea level change (0–2000 m), 1955–2010. Geophys. Res. Lett. 39, L10603 (2012).
, Benefits of CMIP5 multimodel ensemble in reconstructing historical ocean subsurface temperature variations. J. Climate 29, 5393–5416 (2016).
Y. Xue,
S. Good,
I. Ishikawa,
A. Kumar,
M. Rienecker,
A. J. Rosati,
Y. Yin
, A comparative analysis of upper ocean heat content variability from an ensemble of operational ocean reanalyses. J. Climate 25, 6905–6929 (2012).
H. Freeland, D. Roemmich, L. Garzoli, P. Le Traon, M. Ravichandran, S. Riser, V. Thierry, S. Wijffels, M. Belbéoch, J. Gould, F. Grant, M. Ignazewski, B. King, B. Klein, K. A. Mork, B. Owens, S. Pouliquen, A. Sterl, T. Suga, M. Suk, P. Sutton, A. Troisi, P. Joaquin V. Belchi, J. Xu, Argo—A decade of progress, in Proceedings of OceanObs'09: Sustained Ocean Observations and Information for Society, J. Hall, D. E. Harrison, D. Stammer, Eds. (ESA Publication WPP-306, 2009).
K. M. AchutaRao,
B. D. Santer,
K. E. Taylor,
T. P. Barnett,
D. W. Pierce,
R. J. Stouffer,
T. M. L. Wigley
, Simulated and observed variability in ocean temperature and heat content. Proc. Natl. Acad. Sci. U.S.A. 104, 10768–10773 (2007).
A. Ganachaud,
, Improved estimates of global ocean circulation, heat transport and mixing from hydrographic data. Nature 408, 453–457 (2000).
C. Wunsch,
R. Ferrari
, Vertical mixing, energy and the general circulation of the oceans. Annu. Rev. Fluid Mech. 36, 281–314 (2004).
J. Church,
D. Monselesan,
P. Sutton,
, Unabated planetary warming and its ocean structure since 2006. Nat. Clim. Change 5, 240–245 (2015).
S. C. Riser,
H. J. Freeland,
A. Troisi,
M. Belbéoch,
D. Gilbert,
J. Xu,
S. Pouliquen,
A. Thresher,
P.-Y. Le Traon,
G. Maze,
B. Klein,
M. Ravichandran,
F. Grant,
P.-M. Poulain,
T. Suga,
B. Lim,
A. Sterl,
K.-A. Mork,
P. J. Velez-Belchi,
I. Ansorge,
B. King,
J. Turton,
S. R. Jayne
, Fifteen years of ocean observations with the global Argo array. Nat. Clim. Change 6, 145–153 (2016).
J.-B. Sallée,
C. Cabanes,
F. Gaillard,
S. Speich,
M. Hamon
, Consistency of the current global ocean observing systems from an Argo perspective. Ocean Sci. 10, 547–557 (2014).
S. G. Purkey
, Informing deep Argo array design using Argo and full-depth hydrographic section data. J. Atmos. Oceanic Tech. 32, 2187–2198 (2015).
, Estimating annual global upper-ocean heat content anomalies despite irregular in situ ocean sampling. J. Climate 21, 5629–5641 (2008).
R. S. Nerem,
B. Hamlington
, Is the detection of accelerated sea level rise imminent? Sci. Rep. 6, 31245 (2016).
C. E. Forest,
P. H. Stone,
A. P. Sokolov
, Constraining climate model parameters from observed 20th century changes. Tellus A 60, 911–920 (2008).
T. Kuhlbrodt,
J. M. Gregory
, Ocean heat uptake and its consequences for the magnitude of sea level rise and climate change. Geophys. Res. Lett. 39, L18608 (2012).
V. Nieves,
W. C. Patzert
, Recent hiatus caused by decadal shift in Indo-Pacific heating. Science 349, 532–535 (2015).
W. Liu,
S.-P. Xie,
J. Lu
, Tracking ocean heat uptake during the surface warming hiatus. Nat. Commun. 7, 10926 (2016).
J. Robson,
P. Ortega,
R. Sutton
, A reversal of climatic trends in the North Atlantic since 2005. Nat. Geosci. 9, 513–517 (2016).
L. C. Jackson,
R. A. Wood
, Recent slowing of Atlantic overturning circulation as a recovery from earlier strengthening. Nat. Geosci. 9, 518–522 (2016).
J. M. Arblaster,
A. Hu,
, Model-based evidence of deep-ocean heat uptake during surface-temperature hiatus periods. Nat. Clim. Change 1, 360–364 (2011).
S. T. Gille
, Decadal-scale temperature trends in the Southern Hemisphere ocean. J. Climate 21, 4749–4765 (2008).
S. Purkey,
G. Johnson
, Warming of global abyssal and deep Southern Ocean waters between the 1990s and 2000s: Contributions to global heat and sea level rise budgets. J. Climate 23, 6336–6351 (2010).
D. G. Desbruyères,
S. G. Purkey,
E. L. McDonagh,
B. A. King
, Deep and abyssal ocean warming from 35 years of repeat hydrography. Geophys. Res. Lett. 43, 10356–10365 (2016).
L. F. Konikow,
J. G. Cogley,
E. Rignot,
J. M. Gregory,
M. R. van den Broeke,
A. J. Monaghan,
I. Velicogna
, Revisiting the Earth's sea-level and energy budgets from 1961 to 2008. Geophys. Res. Lett. 38, L18601 (2011).
, Observed and simulated full-depth ocean heat content changes for 1970–2005. Ocean Sci. 12, 925–935 (2016).
J. Hansen,
M. Sato,
P. Kharecha,
K. von Schuckmann
, Earth's energy imbalance and implications. Atmos. Chem. Phys. 11, 13421–13449 (2011).
R. P. Allan,
C. Liu,
N. G. Loeb,
M. Roberts,
D. Smith,
P.-L. Vidale
, Changes in global net radiative imbalance 1985–2012. Geophys. Res. Lett. 41, 5588–5597 (2014).
, Changes in tropical clouds and radiation. Science 296, 2095 (2002).
T. P. Boyer, J. I. Antonov, O. K. Baranova, C. Coleman, H. E. Garcia, A. Grodsky, D. R. Johnson, R. A. Locarnini, A. V. Mishonov, T. D. O'Brien, C. R. Paver, J. R. Reagan, D. Seidov, I. V. Smolyar, M. M. Zweng, World Ocean Database 2013, S. Levitus, A. Mishonov, Eds. (NOAA Atlas NESDIS 72, 2013), 209 pp.
G. A. Meehl
, An overview of CMIP5 and the experiment design. Bull. Am. Meteorol. Soc. 93, 485–498 (2012).
W. J. Emery, R. E. Thomson, in Data Analysis Methods in Physical Oceanography (Elsevier, ed. 3, 2004), pp. 533–539.
Acknowledgments: We thank the three anonymous reviewers for their insightful comments on the previous version of this manuscript. Funding: L.C. and J.Z. were supported by the Chinese Academy Sciences' Project "Western Pacific Ocean System: Structure, Dynamics and Consequences" (XDA11010405), the Natural Science Foundation of China (41506029 and 41476016), and 2016YFC1401800. K.E.T. and J.F. were partially sponsored by the U.S. Department of Energy (grant DE-SC0012711). J.F. was sponsored additionally by NSF Award ID 1243125 and NASA Award Number NNH11ZDA001N. The National Center for Atmospheric Research is sponsored by the U.S. NSF. T.B. was funded by the National Oceanographic and Atmospheric Administration Climate Program Office, OAR-CPO-NES-NCEI SLA 17-1-001, from 2016 to 2018. Author contributions: L.C. and J.Z. initiated this study, designed the experiments, and prepared figures and the first manuscript. K.E.T., J.F., T.B., and J.A. contributed to the design of the experiments and analysis of the results and refined the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: The WOD13 data set is available at the National Oceanic and Atmospheric Administration/NCEI website www.nodc.noaa.gov/OC5/WOD13/. The reconstructed data set is available at http://159.226.119.60/cheng/. Argo data were collected and made freely available by the International Argo Program and the national programs that contribute to it (www.argo.ucsd.edu, http://argo.jcommops.org). The Argo Program is part of the Global Ocean Observing System. All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
Vol 3, No. 3
You are going to email the following Improved estimates of ocean heat content from 1960 to 2015
By Lijing Cheng, Kevin E. Trenberth, John Fasullo, Tim Boyer, John Abraham, Jiang Zhu
Science Advances 10 Mar 2017 : e1601545
A new assessment of how much heat Earth has accumulated since 1960 is made by examining ocean heat content changes.


← A little before and after….
I have lived in the Arizona desert since 1959 and have seen a lot.
Your photos are some of the very best I've ever seen.
They really show the beauty of the desert that a lot of people miss.


Box of 50 Tatuaje Verocu No. 5 (4 x 40) packaged and sealed with Humidipak for Humidified Shipping!
New for 2010, a stronger version of the Havana VI blend, this smaller version of the Verocu will be a favorite short smoke!


This has become an annual tradition. I hope you readers enjoy this as much as my family and friends do. It has been a whirlwind year.. took me hours of scratching my head and trying to remember what I was doing each particular month of the year.
I am so stoked to be living in the four-seasons climate again. I was so crushed to have missed the first snowfall. My colleagues didn't get it. As a cold weather girl, winter is a wonderful season of holidays, family and friends, and snow. Oh, snow.
This picture embodies my perfect winter.Taken on the top of one of the trails in Whistler, Canada. Snowboarder is my cousin, looking like he is saying, "Oh my god. I have to go down that??" He was actually also taking a photo of the view.
This is one of my favourite photos. I actually use this as a signature photo now. It was an earlier instance that boosted my confidence in photo taking.
Hanoi, Vietnam. The driver will go somewhere in between.. hopefully.
I knew in going to Myanmar that my access to media and Internet could potentially be non-existent. If I had any, my activities may even be monitored. But I was curious and wanted to test, a little, the limits. So I brought my laptop, which I could use to at least sort and edit my photos as the trip progresses.
facebook: *bing* How'd you pull it off, Zuckerberg??
postscript: I asked a guide about how he accesses his Internet when he mentioned he uses gmail. By using "https" vice "http," all the above sites, I could access.
I am especially proud of this slideshow.. The photos aren't necessarily all the ones I consider the best from the trip. But the selection well represents each area, each aspect of Myanmar that caught my interest and held my fascination. I hope you have a moment to sit back and enjoy.


Astonishing X-Men gets a face-lift in the form of a new creative team. Lead by writer Matthew Rosenberg and artist Greg Land, this new lineup and direction is a very promising one.
X-Men fan-favorite, Alex Summers, returns to his normal self. Havok is looking to get back to his heroic ways with a new X-Men team. After a public humiliation at the hands of the Avengers, his search for teammates gets off to a rocky start.
Right away, it feels so good to have Havok back to normal and in his classic costume. His characterization is a bit puzzling at first, he reads more like Jamie Madrox (who Rosenberg is also currently writing on Multiple Man) but there's only a brief adjustment period. His new attitude sets a light tone that is sure to smack a smile onto the face of every die-hard X-fan.
Matthew Rosenberg has absolutely been on fire recently with his Marvel contributions, this seems to be another great addition to his collection of X-Men work. There's an injection of joy and excitement that was missing towards the end of Charles Soule's 12-issue run of Astonishing X-Men. With the heavy dose of humor, there's also a star-studded cast that bounce perfectly off of Summers as he takes a tour through some of Marvel's mutants.
The little character moments are fantastic. Alex and Hank talking about superhero fatigue, Kitty phasing through the chalk-board, Havok getting the Stark treatment, etc. Astonishing X-Men #13 has so much charm spread throughout all of the featured characters. Havok's thirst for being an X-Men and what it means to him puts him right in the reader's shoes.
Greg Land may be the biggest source of joy in this issue, and not just for the fantastic cover art. Havok looks incredible in his classic costume, his abilities are illustrated in glorious fashion. Land delivers on every character, making them all distinctly different and visually pleasing. Colorist Frank D'Armata has an electric chemistry with the comic's artist, treating us to some of the finest work for both creators.
As X-Men: Gold and X-Men: Blue trend more and more into obscurity, irrelevance, and apparent cancellation, Astonishing X-Men emerges to keep X-Men fans satisfied. Between Astonishing X-Men and X-Men: Red, Marvel has their new flagship X-books. It's safe to say at this point, any time you see Matthew Rosenberg writing mutants you need to pick it up.


For convenience, nothing beats driving your own car. So, at Hertz, we've been try- ing to find a way of renting cars that comes as close to this convenience as possible. It took several years and several million dol- lars, but we finally di d it. It's called Hertz Number One Club. Once you're a mem- ber you'll never again have to bother with a rent a car form and all the aggravation that goes with it. Here's ho\v Number One Club works. Apply for a free membership card at any Hertz counter or call toll-free 800-654-3131. Then, when you re- serve a quiet Ford or other car anywhere in the country, just say you're a member of Hertz Number One Club, and your rental form will be filled out and waiting. So all you do when you get to the Hertz counter is show your license and charge card, sign your name and go. Join Hertz Number One Club. It's the begin- ning ofa whole new era in rent- Hertz . lng cars. ou ont fill nvey ur wne III them u . " '. . e u rms u on have nve urs. . Here's how Number One Cub works Just show your license and charge card. Jk .>>. j.*x. -* :. Po , . .. " -.t. Sign your name. . ", " '" "<- And go. OHE rz YS _M Hertz #1 Club It the little things that make us big. -.v v >0 v' 1.:\ .. --:;;:;;J ---J - . -........ . - '"


text:
Items Nymph and satyr
Nymph and satyr
The Miriam and Ira D. Wallach Division of Art, Prints and Photographs: Picture Collection
Lebègue, Léon-G., 1863-1944
Wonders: Images of the Ancient World
Mythology -- Nymphs
Lebègue, Léon-G., 1863-1944 (Artist)
Shelf locator: PC-WON OX MYT-Nym
Nymphs (Greek deities)
Satyrs (Greek mythology)
Borders (Ornament areas)
Thyrsus
Content: Includes additional text in French.
Source note: Les contes de Jacques Tournebroche. (Paris : Calmann-Levy, 1909) Lebègue, Léon (1863-), Illustrator. France, Anatole (1844-1924), Author.
Chromotypographs
Extent: Sheet 19.5 x 12.7 cm
Universal Unique Identifier (UUID): 110bbe70-c5bf-012f-b2b0-58d385a7bc34
The Miriam and Ira D. Wallach Division of Art, Prints and Photographs: Picture Collection, The New York Public Library. "Nymph and satyr" The New York Public Library Digital Collections. 1909. http://digitalcollections.nypl.org/items/510d47e4-0eea-a3d9-e040-e00a18064a99
The Miriam and Ira D. Wallach Division of Art, Prints and Photographs: Picture Collection, The New York Public Library. "Nymph and satyr" New York Public Library Digital Collections. Accessed January 19, 2020. http://digitalcollections.nypl.org/items/510d47e4-0eea-a3d9-e040-e00a18064a99
The Miriam and Ira D. Wallach Division of Art, Prints and Photographs: Picture Collection, The New York Public Library. (1909). Nymph and satyr Retrieved from http://digitalcollections.nypl.org/items/510d47e4-0eea-a3d9-e040-e00a18064a99
<ref name=NYPL>{{cite web | url=http://digitalcollections.nypl.org/items/510d47e4-0eea-a3d9-e040-e00a18064a99 | title= (still image) Nymph and satyr, (1909) |author=Digital Collections, The New York Public Library |accessdate=January 19, 2020 |publisher=The New York Public Library, Astor, Lenox, and Tilden Foundations}}</ref>

redpajama_set_name: RedPajamaCommonCrawl


During science week pre-school children at Leicester, Blaby took part in a group activity to create a volcano and to watch this erupt using vinegar, bicarbonate of soda and food colouring. The child were so engaged in this activity that staff asked parents if they would like to create a volcano with their child and bring in for them to explode!! Children and parents have been very busy and created some amazing volcanoes, some complete with dinosaurs.
Today the preschoolers have started to erupt their volcanoes. Awe and wonder is seen in abundance as the lava flows freely!! Some of the comments we have heard include "lava comes out the top and goes down there" "smoke comes out at first" "it's going out the top, whoa!!" "wow, its bubbling out, to coming to me" It is safe to say this is one activity that is going to be spoken about alot!


text:
El tirà ornat (Myiotriccus ornatus) és un ocell de la família dels tirànids (Tyrannidae) i única espècie del gènere Myiotriccus.

Hàbitat i distribució
Boscos clars de les vessants dels Andes, des de l'oest, nord i est de Colòmbia, cap al sud, a través de l'oest i est de l'Equador fins l'est i sud-est del Perú.

Taxonomia
Alguns autors han considerat que es tracta en realitat de dues espècies:
 tirà ornat meridional  (Myiotriccus phoenicurus Sclater, PL, 1855). Del sud-est de Colòmbia i est del Perú.
 tirà ornat septentrional  (Myiotriccus ornatus sensu stricto). De Colòmbia central i occidental i oest de l'Equador.

Referències

Tirànids

redpajama_set_name: RedPajamaWikipedia


text:
Ethnic discrimination
Article | Last updated: 07/09/2020 | Ministry of Culture
Everyone must have the same rights, obligations and opportunities, irrespective of ethnic background. The Equality and Anti-Discrimination Act provides protection against ethnic discrimination.
Discrimination based on ethnicity, religion or belief
Everyone must have the same rights regardless of ethnicity, religion and belief. The Equality and Anti-Discrimination Act provides protection against these types of discrimination.
The integration of efforts to promote equality and prevent discrimination is a goal in all policy areas and at all administrative levels. The Ministry of Culture plays a coordinating role in this work. However, all the ministries have an independent responsibility to promote equality in their areas of responsibility.
The Norwegian Directorate for Children, Youth and Family Affairs is responsible for following up directorate tasks relating to the prevention of discrimination based on ethnicity, religion and belief.
UN Racial Discrimination Convention
The UN Convention on the Elimination of All Forms of Racial Discrimination prohibits discrimination on the basis of race, skin colour, descent or national or ethnic origin.
The convention was adopted by the UN General Assembly on 21 December 1965 and came into force in January 1969. Norway ratified the convention without reservation on 6 August 1970. The convention has been incorporated into Norwegian law through the Equality and Anti-Discrimination Act. National implementation of the UN Racial Discrimination Convention is monitored by the UN Committee on the Elimination of Racial Discrimination, which is made up of independent experts.
The Equality and Anti-Discrimination Act
The Equality and Anti-Discrimination Act came into force on 1 January 2018. It prohibits direct and indirect discrimination based on factors such as ethnicity (including national origin, descent, skin colour and language), religion and belief. The purpose of the act is to promote equality, ensure equal opportunities and rights, and prevent discrimination. The act applies in all areas of society, but is not enforced in the context of family life and other personal relationships. Employers must make active, targeted and systematic efforts to promote equality and prevent discrimination in their undertakings. In addition, they must report on implemented and planned equality measures (the activity and reporting duties). Public authorities also have an activity duty in their areas of responsibility.
The Equality and Anti-Discrimination Ombud and the Equality and Anti-Discrimination Tribunal
The Equality and Anti-Discrimination Ombud is mandated to promote equality and combat discrimination on the basis of, among other things, gender, ethnicity, religion, belief, disability, sexual orientation and age.
The Ombud gives advice and is a driving force for equality and diversity. The Equality and Anti-Discrimination Tribunal enforces the statutory prohibitions against discrimination. Contact the Ombud and the Tribunal.
In June 2019, the Storting (the Norwegian parliament) decided to authorise the Equality and Anti-Discrimination Tribunal to enforce the prohibition against sexual harassment in the Equality and Anti-Discrimination Act, and to award compensation in cases concerning work-related sexual harassment and certain other instances.
Duty to promote equality actively
The Equality and Anti-Discrimination Act contains rules on the duty of employers, public authorities and organisations in working life to promote equality actively.
The purpose of the duty to promote equality actively is to raise awareness through ongoing activities within an organisation, and thereby to prevent discriminatory practices and organisational structures.
The activity duty of employers applies in areas such as recruitment, pay and working conditions, promotion, development opportunities and protection against harassment. While the duty does not prescribe specific measures, public-sector employers and private-sector employers with more than 50 employees are required to follow a four-step methodology. It is up to each undertaking to decide which measures are relevant in view of its own individual challenges. Employers have to report on implemented and planned measures in their annual report or annual budget.
In June 2019, the Storting decided to reinforce the activity and reporting duties of employers. Following the statutory amendment, the duties of employers are now specified in greater detail. The Equality and Anti-Discrimination Ombud and the three regional centres for equality and diversity play an important role in providing advice and guidance on the activity and reporting duties.
The Norwegian Government's Action Plan against Racism and Discrimination on the Grounds of Ethnicity and Religion 2020-2023 (Extracted Version)
Department for Equality, Non-discrimination and International Affairs
E-mail: postmottak@kud.dep.no
Address: Postboks 8030 Dep, 0030 Oslo
Visitor address: Grubbegata 1, Oslo,
The following persons are responsible for the Ministry of Culture's website:
Editor-in-chief: Marius Bakke
Web editor: Wenche Stadven Nybo
Email: postmottak@kud.dep.no
Organisation number: 972 417 866


*Charging up to 80% battery capacity in ideal conditions.
1- Based on NEDC testing standards with a base-level equipped vehicle.Range may differ depending on various factors, in particular:individual driving style, route characteristics, outside temperature, heating / air conditionaing, seated and standing passenger weight.
** With outward swing passenger door.


This month, saw a couple of Everything Dinosaur team members visit a trade show. We had lots of meetings to attend and whilst at the event, we took the opportunity to shoot a quick video featuring four of the new for 2017 Schleich dinosaur models. We appreciate the lighting at such venues does not always help when it comes to creating videos, but we thought that we could at least provide a quick peek at some of the larger Schleich models coming into our warehouse in just a few days.
The dinosaur models featured in the video are the Schleich Brachiosaurus (interesting skin texture on the back and tail of this model), the Allosaurus with an articulated jaw (please excuse the grammatical error in the narration). Plus, the Stegosaurus and the feathered Utahraptor, Schleich getting to grips with feathered Dromaeosaurs by bring out a shaggy coated Utahraptor to join the existing range of Schleich models that includes "scaly raptors". The video is quite short, it is under forty-five seconds long, but at least it provides a preview.
Schleich prehistoric animal brochures will be available to acquire from Everything Dinosaur's new website, which is due to be launched next month (February 2017), the new for 2017 Schleich models, including those featured in this video, will be in stock at our warehouse shortly.
One of the new introductions is the Schleich Brachiosaurus. The designers have chosen to give this dinosaur model a very textured skin, it's as if this super-sized dinosaur had super-sized scales on its neck, its back and on the tail, plus partially on the upper limbs. It is certainly quite an interesting effect, one that will make this dinosaur model very popular with those young dinosaur fans who like to see large scales on models.
The new for 2017 Schleich Brachiosaurus dinosaur model.
As with all the Schleich models the paintwork is very good, and there is plenty of detail to be found within the range. Schleich dinosaurs are suitable for children from three years and upwards.
The new for 2017 Schleich Allosaurus dinosaur model.
Also featured in the video, is the new Schleich Allosaurus figure.
Look out for more updates on new for 2017 models from Everything Dinosaur.


Do we need to teach the IR paradigms at all?
Over the weekend, fellow guest contributor Luke Perez had an interesting post on whether we need to include the grand paradigms of international relations (realism, liberalism, and constructivism) in foreign policy classes. He makes some good points on how to customize courses for foreign policy students; be sure to read it if you haven't. I'd like to go further and ask whether we need to teach these paradigms at all.
I'm coming at this from a different perspective than Perez. I teach at an undergraduate focused institution. So I'm preparing students for a broad array of potential political science careers. But the issue with the paradigms' importance transcends any single realm of higher education.
While in graduate school, my time as a teaching assistant for intro to IR courses and my own PhD seminars in international relations let me observe several different ways to organize an intro to IR course. This gave me a few potential models when I started putting together my own syllabus.
I had a lot to work out, but I decided on one thing: I would not organize the entire course around the paradigms. It's possible to have an "-ism-focused" course, presenting multiple sessions on each theory, complemented by related topics. For example, realism could include theory lectures ranging from classical realism to neoclassical realism, followed by lectures on the balance of power and other topics.
That didn't fit my own understanding of international relations, however. I don't think, for example, that power politics is necessarily a realist topic (for a systematic discussion on this, see this article by Goddard and Nexon). So I'd want power politics to be separate from realism. And as someone who studies religion and international relations, I wouldn't feel comfortable saying that's part of constructivism. That's where it would seem to belong except for the fact that constructivism never really dealt with religion. But an -ism-focused syllabus would force me to squeeze it in with constructivism. So I cover the –isms in one unit alongside the levels of analysis, but then proceed with narrower mid-range theories and policy issues.
The paradigms do not "own" any topic: "Balance of power" seems realist, but involves plenty of norms and symbolic politics (the English School has plenty to say on how states maintain order). Likewise, international institutions may be neoliberals' specialty, but they can also be a venue for power politics. If we make it seem like the paradigms are the foundation of all studies, however, students will attempt to make connections back to the paradigms on each topic. This will either flatten interesting elements of many international relations topics or lead to confusion. Minimizing the significance of the paradigms, however, can allow us to study each topic on its own merits.
The paradigms don't cover all topics (but we try to pretend they do): Students will often ask about a neoliberal explanation for war, or what constructivists say about terrorism. And some paradigm-focused textbooks try to have parallel debates on all issues across the paradigms ("on X, realism would say A, liberalism B and constructivism C"). In these situations we may try to make up paradigm-based arguments that don't exist. I could tell you a convincing constructivist explanation for al-Qaeda, but that's not something they focused on. At that point, we end up teaching students what the paradigms might have said or–in my case–explaining the massive gaps in what each school of thought found important. If the paradigms were not the centerpiece of a course, we would avoid these problems.
The paradigms are not relevant to today's discipline: This is the big one. In graduate school I came up with all these ideas about how to draw on elements of realism to study religion and power politics, then got to my first academic conference and realized…no one cared. We weren't debating the paradigms anymore. You didn't have to ground your arguments in deep theoretical foundations. We still have theoretical debates, but they're narrower and cross-paradigmatic. For example, debates on international hierarchy involve ideas and norms without relying on constructivism, and actually provide a forum to integrate this approach with work that would otherwise be neoliberal. So if learning about the paradigms doesn't necessarily help future professors, is it really useful for undergraduates who may pursue international affairs careers outside of academia?
The paradigms come up in foreign policy discussions: There are realist and liberal foreign policies, so understanding these theories may help students working in that realm. But the foreign policies don't necessarily line up well with their academic equivalents. And there is no constructivist foreign policy, a source of endless frustration for students. We would be better served discussing these foreign policy approaches separately, such as in a grand strategy lecture.
They are needed for graduate school: I've never sat on a graduate admissions committee, but I can imagine they want students with some familiarity with the paradigms. That is partly why I still teach them, to make sure students who want to get a PhD are prepared. But this is, as constructivists would say, a social construction. The paradigms are necessary because we think they're necessary. PhD programs could easily rethink their requirements.
I won't stop teaching the paradigms anytime soon, but this why I'm tempted to minimize their importance. I imagine there are counterarguments (and would love to hear them). And I didn't really get into broader debates about whether international relations needs grand theory; in fact, I sometimes we wish were still fighting the paradigm wars. But until international relations emphasizes paradigmatic theoretical debates again, I think we need to re-evaluate their role in our courses.
Syllabus Musing—can we teach foreign policy without paradigms?
IR is an obsolete academic discipline that needs to be retired.
Hope that in 20, 30 years, no one will be teaching and studying IR.
You don't need paradigms or "theories" to understand global politics or behaviour of states.
You need HISTORY and understanding of CULTURE of the society of that state. That is what needs to be taught and studied. But IR relegates culture to irrelevance, which in my view makes IR to be complete garbage.
Two good posts (including Luke's on FP), but as someone who teaches the IR 'paradigms' at both undergrad and postgrad level in the UK (as well as US Diplomatic History), I think both posts miss three important things.
First, neither posts seem to include Marxism as a theory / paradigm, which is interesting. Are you saying there is no need to help students understand the role of geopolitical economic structures in foreign policy or world politics? This seems a major gap. Assuming you don't mean this then surely a (basic!) introduction to Marxist theory is required, not least to illustrate the evolution from productive force Marxist models to Critical Theory. It seems hard to me to teach international relations (even without the capital I and R) and not at least give some pointers towards traditional Marxism, critical theory and World Systems Theory / dependency theory.
In my view there is a danger that ignoring 'paradigms' ignores the academic politics and power relations that exist *within* the fields in question. Who defines the 'paradigms' that we should cover, and why? I remain seriously unconvinced that IR Theory / foreign policy / teaching international politics, properly allows for true epistemological and ontological nuance and range. It's not the best way, but I think a good way of allowing for that room is by examining the epistemological and ontological stances of 'paradigms'. Peter, your statement that you would not feel comfortable being linked to Constructivism in your research on religion and world politics is a *very* interesting thing to explore and teach! My students would get a lot out of that 'theoretical' point! What does that argument tell us (both methodologically and philosophically)?
Putting this another way, and somewhat getting back to my first point, my own view is that IR is too concentrated around neo-positivism, so if we don't re-examine the 'paradigms' we take as read too many features accepted within mainstream arguments. The (excellent) Nexon and Goddard article still lacks a full engagement with the challenges posed by post-modernism (as demonstrated by so much being wrapped in 'constructivism'). That paper posits three broad competing positions in security studies, which is fine, but that's the *main* positions, not them all, and there is a benefit to letting students know this. You never know, one of our undergrad or postgrads may find something new and beneficial to us all in theoretical epitaphs.


text:
import type { Modifier, ModifierArguments, Padding, Rect } from '../types';
import type { Placement } from '../enums';
import getBasePlacement from '../utils/getBasePlacement';
import getLayoutRect from '../dom-utils/getLayoutRect';
import contains from '../dom-utils/contains';
import getOffsetParent from '../dom-utils/getOffsetParent';
import getMainAxisFromPlacement from '../utils/getMainAxisFromPlacement';
import within from '../utils/within';
import mergePaddingObject from '../utils/mergePaddingObject';
import expandToHashMap from '../utils/expandToHashMap';
import { left, right, basePlacements, top, bottom } from '../enums';
import { isHTMLElement } from '../dom-utils/instanceOf';

// eslint-disable-next-line import/no-unused-modules
export type Options = {
  element: HTMLElement | string | null,
  padding:
    | Padding
    | (({|
        popper: Rect,
        reference: Rect,
        placement: Placement,
      |}) => Padding),
};

const toPaddingObject = (padding, state) => {
  padding =
    typeof padding === 'function'
      ? padding({ ...state.rects, placement: state.placement })
      : padding;

  return mergePaddingObject(
    typeof padding !== 'number'
      ? padding
      : expandToHashMap(padding, basePlacements)
  );
};

function arrow({ state, name, options }: ModifierArguments<Options>) {
  const arrowElement = state.elements.arrow;
  const popperOffsets = state.modifiersData.popperOffsets;
  const basePlacement = getBasePlacement(state.placement);
  const axis = getMainAxisFromPlacement(basePlacement);
  const isVertical = [left, right].indexOf(basePlacement) >= 0;
  const len = isVertical ? 'height' : 'width';

  if (!arrowElement || !popperOffsets) {
    return;
  }

  const paddingObject = toPaddingObject(options.padding, state);
  const arrowRect = getLayoutRect(arrowElement);
  const minProp = axis === 'y' ? top : left;
  const maxProp = axis === 'y' ? bottom : right;

  const endDiff =
    state.rects.reference[len] +
    state.rects.reference[axis] -
    popperOffsets[axis] -
    state.rects.popper[len];
  const startDiff = popperOffsets[axis] - state.rects.reference[axis];

  const arrowOffsetParent = getOffsetParent(arrowElement);
  const clientSize = arrowOffsetParent
    ? axis === 'y'
      ? arrowOffsetParent.clientHeight || 0
      : arrowOffsetParent.clientWidth || 0
    : 0;

  const centerToReference = endDiff / 2 - startDiff / 2;

  // Make sure the arrow doesn't overflow the popper if the center point is
  // outside of the popper bounds
  const min = paddingObject[minProp];
  const max = clientSize - arrowRect[len] - paddingObject[maxProp];
  const center = clientSize / 2 - arrowRect[len] / 2 + centerToReference;
  const offset = within(min, center, max);

  // Prevents breaking syntax highlighting...
  const axisProp: string = axis;
  state.modifiersData[name] = {
    [axisProp]: offset,
    centerOffset: offset - center,
  };
}

function effect({ state, options }: ModifierArguments<Options>) {
  let { element: arrowElement = '[data-popper-arrow]' } = options;

  if (arrowElement == null) {
    return;
  }

  // CSS selector
  if (typeof arrowElement === 'string') {
    arrowElement = state.elements.popper.querySelector(arrowElement);

    if (!arrowElement) {
      return;
    }
  }

  if (__DEV__) {
    if (!isHTMLElement(arrowElement)) {
      console.error(
        [
          'Popper: "arrow" element must be an HTMLElement (not an SVGElement).',
          'To use an SVG arrow, wrap it in an HTMLElement that will be used as',
          'the arrow.',
        ].join(' ')
      );
    }
  }

  if (!contains(state.elements.popper, arrowElement)) {
    if (__DEV__) {
      console.error(
        [
          'Popper: "arrow" modifier\'s `element` must be a child of the popper',
          'element.',
        ].join(' ')
      );
    }

    return;
  }

  state.elements.arrow = arrowElement;
}

// eslint-disable-next-line import/no-unused-modules
export type ArrowModifier = Modifier<'arrow', Options>;
export default ({
  name: 'arrow',
  enabled: true,
  phase: 'main',
  fn: arrow,
  effect,
  requires: ['popperOffsets'],
  requiresIfExists: ['preventOverflow'],
}: ArrowModifier);

redpajama_set_name: RedPajamaGithub


text:
Assistant Douglas County manager might apply for top job
Local | October 25, 2018
Robert Macias
On Thursday, Assistant County Manager Jenifer Davidson announced she was planning to apply for the county's top job.
Davidson shared her plans during a discussion on the progress toward hiring a new county manager.
"I am considering putting my name in the hat for the position," she said.
Because Davidson could be a candidate for the office, she withdrew as the person who would work with Human Resources Director Wendy Lang on hiring a search firm.
"I was previously involved and provided her with direction," Davidson said. "Now I'm going to take a step back from the process and have her report directly to the board."
Because the update was not slated for action, commissioners could only talk about it, but that didn't stop Nancy McDermid from suggesting that the county might not need a search firm after all.
"We laid out the process prior to learning that the interim county manager was interested in putting in for the county manager position," she said. "I have been concerned that if we get people in our key positions who do not know Douglas County or Nevada with the Legislature coming up."
Commissioners instructed Lang in September to work with Davidson to hire a search firm to aid in finding a new county manager.
Lang said she has reviewed several firms and will have the job description revised and posted by the end of October.
Commissioner Barry Penzel said he wanted future discussions on the county manager's position agendized so commissioners can take action.
Chairman Steve Thaler agreed with Penzel.
"The question is whether we're going to spend $50,000 on a recruiting firm," he said.
Commissioner Dave Nelson said he was happy Davidson was planning to apply for the job.
"I'm very happy she put her hat in the ring, but I do believe it's just one more hat in the ring," he said. "I'd like to see us go forward with the current plan of action."
County commissioner-elect John Engel was enthusiastic at the prospect of Davidson becoming county manager.
"I strongly support her being out new county manager," he said. "I look for several qualities when hiring someone, but two of them are intangible. Those are integrity and honesty. She has been forthright in every encounter I've had with her."
He also felt Davidson's experience with water as Minden Town manager would be valuable.
"She has a fundamental understanding that I haven't seen in any of the other candidates," he said.
Davidson is a 2000 graduate of Douglas High School and holds a degree from Boise State University and completed the University of Nevada, Reno, supervisory management program.
She was administrative services manager at China Spring Youth Center before becoming town manager in 2011.
She was named assistant county manager in March.
Roadwork schedule for the Lake Tahoe area: Jan. 12-18
Call for artists issued for new mural on Lake Tahoe's East Shore Trail
Grateful snowboarder describes being rescued after a night out of bounds at Heavenly


Upgrading to the Double Diamond Package adds 2500 Plutos Credits (PC) to your account and pays 1.0% PC per day according to your settings. Also, if you pay with $US you will receive a bonus of 150 Ambassador Credits. Please have funds available in your selected balance prior to placing order.


You are here: Home | Exhibitor search | Kitagawa Inc.
My philosophy is that a meaning of a company's existence is always demanded.
We always try to develop desired products and become a company needed from society and cultivate human material.
It's a meaning of a company's existence for me.
The products which fill a satisfaction of either sellers or customers, but not both doesn't sell well or last long.
Only the products which satisfy the customers, the sellers and the manufacturers last long.
Even if it is difficult to fulfill a hundred percent, we always try to improve our products to satisfy and please customers all over the world.
We will make efforts to spread our wings around the world.


text:
Capital Yorkshire
By The Bridge
By the Bridge's unique therapeutic approach to fostering puts foster parents at the centre of their team and helps a child move forward.
Seeing a child in care discover their potential, giving them the best possible chance at happiness and success - it's what makes fostering with By the Bridge so rewarding.
By the Bridge was set up by people who have been foster parents themselves and have many members of staff who have fostered too, so always have someone on hand who really understands what it is like to foster a child in their own home.
With expert training you'll learn the skills to become a professional foster parent, ready to deal effectively with your foster child's needs.
Contact the fostering team on 0845 450 99 44
Visit the By the Bridge website
Find your nearest contact centre here
Find them on Twitter
Find and train the best foster parents
so that young people in care get a better chance in life


Whether your girlfriend is a professional zombie head or someone new to the undead scene, the best zombie ideas can make the best gifts. Keep your girlfriend alive (and hopefully her interest in the undead scene) with these wonderful and colorful zombie gifts.
Survival. Whether you're a zombie or a terrified civilian, survival is the name of the game. And there's no better time to prepare taking care of yourself than now. And what better way to make sure that your head remains intact, whether that's a decomposing head or a panicky human head, than to utilize all the best zombie survival gift ideas around.
No time to sit down and read with all of the undead clawing their way towards you? Well the Zombie Survival Guide Deck will provide you with that all important emergency crash course in zombie apocalypse survival. Flashcards, they're a great idea and one of the best zombie survival gift ideas around.
Need to run? The Zombie Apocalypse Zombie Apocalypse Survival Kit Sardine Can is your answer to getting out of those tight spots fast. No need to pack, this sardine contains a whistle, non-aspirin pain reliever, compass, matches, first aid kit, fish hook and line, and many more. This sardine can will definitely keep you alive during the coming zombie apocalypse.
It's the Halloween, you and your girlfriend are planning to enjoy the occasion (or at least, you're planning to enjoy the occasion as soon as she says yes to your idea). What better way to get her into the spirit of Treat or Tricking than by giving her a Walking Dead Halloween present. After all, no girl in her right TWD-watching mind can resists a Walking Dead Halloween present.
One of the best zombie ideas that you can have is to go as your favorite Walking Dead Couple during the Halloween (or even any event, for that matter. Zombie-inspired weddings? Check!). And with The Walking Dead costumes, props and decoration, you can definitely go in undead style.
Want a weapon? Why not go all Rick Grimes on your neighbors with the TWD Rick Grimes Foam Machete? Your girlfriend can then go as your trusty Michonne complete with her Michonne Wig and crazy Katana. You can even go as a pair of zombies if you just want to lumber around mumbling brains and eating all of your friend's cake and candy. Holding hands not required but recommended.
This is especially great if your girl is a fan of TWD. You don't have to convince her too much to wear a costume. Unless of course, the both of you want to go as Rick Grimes.


Step into the fall and winter seasons in style, wearing this casual women's striped blouse with its rounded hemline. This tunic length blouse is made from our finest cotton blend fabric. It has long, close fitting sleeves in classic black. The neckline is a simple, round one, with double reinforced stitching to help it retain its original shape. The neckline is trimmed with black piping. The torso of this blouse is loose fitting and comfortable. Its longer length provides a figure flattering design with the potential to hide any number of unwanted figure flaws. The hemline is rounded, and ends well below the hips. The torso portion of the top has narrow horizontal stripes in black, across a solid white background. Pair this cute top with your favorite skinny jeans for a great look for the cooler months of the year.


OSTA is excited to announce the OSTA State Tournament, Qualifier for USAT National Championships, will be held on April 27th at the Summit High School Gym in Reynoldsburg Ohio.
Event Poster for OSTA 2019 Ohio State Championship.
The 2018 OSTA State Championships will be held May 5th, 2018 at the Ohio Expo Center (Lausche Building) in Columbus.


The 2018 Trees of Lights Campaign is underway! Each year residents of Fairfield and the surrounding area make donations to the Chamber of Commerce for the holiday decorations on the Fairfield square. The Chamber's Christmas Committee, which meets several times throughout the year, delegates Christmas funds to the appropriate area of need. The Trees of Lights Campaign pays for new decorations on the square as well as renovation of the long-standing traditional decorations that so many people have come to love and enjoy.
If you would like to make a donation this year, please contact the Chamber at 641-472-2111 or email our Event Coordinator, Chelsey Haynes, at Events@FairfieldIowa.com.


Reborn doll pacifier with Magnet!
This authentic diaper is perfect for Reborn dolls!
7739 - Accessories : Elastic Beading String 8 mtr transp.
Use with our elastic bead string for easy braceleting and necklace making.
Plenty enough to make custom made name bracelets for your baby dolls!
Perfect to use with alphabet beads!