Met Office Hadley Centre observations datasets
Home > HadISDH >

HadISDH Figures from Willett et al., (2014)

Please read the paper before using the material presented here.

All figures have been created using HadISDH.land*.2.0.0.2013p and are under Crown Copyright. Please acknowledge use of this material using the following reference:

Willett, K. M., Dunn, R. J. H., Thorne, P. W., Bell, S., de Podesta, M., Jones, P. D., Parker D. E. and Williams Jr., C. N.,submitted.: HadISDH land surface multi-variable humidity and temperature record for climate monitoring. Climate of the Past.
Open Access pdf and Supplementary material

Title Description Link Images
Figures from Willett et al. 2014
Fig. 1 HadISDH Processing Flow diagram of HadISDH.2.0.0 processing from hourly quality-controlled values to monthly mean homogenised station values. Fig1_HadISDH.2.0.0Process_v1.png Fig 1
Fig. 2 Station coverage Final station coverage maps with supersaturated and subzero error stations identified that are not included in the final HadISDH data sets. (a) Good stations common to all variables, (b) additional temperature stations, (c) additional dew point temperature stations, (d) additional specific humidity stations, and (e) additional relative humidity stations. See Fig. SM1 for other variables. Fig2_StationLocationsDiffs.landALL.2.0.0.2013p_SEP2014_v1.png Fig 2
Fig. 3 Station temporal coverage Station temporal coverage. Colours show missing station months by WMO ID (left y-axis) and the solid black line shows total station count (right y-axis). WMO IDs/colours are bounded as follows: Europe: 0-199999 (dark red), Russia/Eastern Europe: 200000-379999 (crimson), Central Asia/Middle East/India/Pakistan: 380000-439999 (dark orange), South East Asia/East Asia/China: 440000-599850 (orange), Africa: 600000-689999 (gold), USA/Canada: 690000-749 999 (dark green), Central America: 760000-799999 (olive green), South America: 800000-879999 (royal blue), Antarctica: 880000-899999 (sky blue), Pacific Islands (Inc. Hawaii): 911000-919999 (navy), Papua New Guinea/Australia/New Zealand: 930000-949999 (lilac), Indonesia/Philippines/Borneo: 960000-999999 (purple). (a) Temperature, (b) dew point temperature, (c) specific humidity and (d) relative humidity. See Fig. SM2 for other variables. Fig3a_TemporalCoverage.landT.2.0.0.2013p_SEP2014INV_v1.png, Fig3b_TemporalCoverage.landTd.2.0.0.2013p_SEP2014INV_v1.png, Fig3c_TemporalCoverage.landq.2.0.0.2013p_SEP2014INV_v1.png, Fig3d_TemporalCoverage.landRH.2.0.0.2013p_SEP2014INV_v1.png Fig 3a Fig 3b Fig 3c Fig 3d
Fig. 4 Pagri, Tibet Example of an inhomogeneity having a greater signal-to-noise ratio in DPD compared to q for station 557 730, Pagri, Tibet. (a) DPD annual mean series for raw data (red) and direct PHA homogenised data (blue) for 557730 and its raw neighbour series (black). (b) q annual mean series for raw data (red) and direct PHA homogenised data (blue) for 557730 and its raw neighbour series (black). (c) q annual mean series for raw data (red) and indirect PHA homogenised data (blue) for 557730 and its raw neighbour series (black). Fig4a_55773099999_trendcompPHAdpd_7313SEP2014abs_v1.png, Fig4b_55773099999_trendcompPHAq_7313SEP2014abs_v1.png, Fig4c_55773099999_trendcompIDPHAq_7313SEP2014abs_v1.png Fig 4a Fig 4b Fig 4c
Fig. 5 Homogeneity adjustment frequency and magnitude Adjustment sizes and locations applied by PHA (a, b, e, f) and ID PHA (c, d, g, h) to all monthly mean station time series combined. (a, c, e, g) Distribution of adjustment size (black line) with a Gaussian fit (grey line), a best-fit merged line taking the larger wings and Gaussian fit over the missing middle (red dashed line) and the difference between the best-fit line and actual distribution (blue dotted line) for temperature, dew point depression and relative humidity, respectively. The mean and standard deviation of the differences are shown. b,d,f,h) Distribution of adjustments over time for temperature (a, b, c, d), dew point depression (e, f) and relative humidity (g, h), respectively. Fig5ab_HadISDH.landT.2.0.0.2013p_adjspread_PHA_SEP2014_v1.png, Fig5cd_HadISDH.landT.2.0.0.2013p_adjspread_IDPHAMG_SEP2014_v1.png, Fig5ef_HadISDH.landDPD.2.0.0.2013p_adjspread_PHA_SEP2014_v1.png, Fig5gh_HadISDH.landRH.2.0.0.2013p_adjspread_IDPHA_SEP2014_v1.png Fig 5ab Fig 5cd Fig 5ef Fig 5gh
Fig. 6 Annual climatologies Annual climatologies (1976-2005) for (a) temperature, (b) dew point temperature, (c) specific humidity and (d) relative humidity. See Fig. SM5 for other variables. Fig6a_HadISDH.landT.2.0.0.2013p_FLATgridIDPHA5by5_SEP2014_climsANN7605_v1.png, Fig6b_HadISDH.landTd.2.0.0.2013p_FLATgridPHADPD5by5_SEP2014_climsANN7605_v1.png, Fig6c_HadISDH.landq.2.0.0.2013p_FLATgridIDPHA5by5_SEP2014_climsANN7605_v1.png, Fig6d_HadISDH.landRH.2.0.0.2013p_FLATgridIDPHA5by5_SEP2014_climsANN7605_v1.png Fig 6a Fig 6b Fig 6c Fig 6c
Fig. 7 Comparing raw and homogenised trends in temperature Comparison of decadal trends from 1973-2013 between raw and homogenised monthly mean temperature. (a) Ratio of raw to homogenised trends for each grid box (a value of > 1 means that the raw trend is greater than the homogenised trend). (b) Scatter plot of raw verses homogenised decadal trends for each grid box with percentage presence in each quadrant. (c) Trend size distribution for raw and homogenised grid boxes. (d) Large-scale area-average time series (weighted by the cosine of the latitude) and decadal trends calculated using the median of pairwise slopes with 5th to 95th percentile slopes as con?dence intervals. Fig7_HadISDH.landT.2.0.0.2013p_FLATgridIDPHA5by5_MAY2014_MPtrendsscat_19732013_v1.png Fig 7
Fig. 8 Comparing raw and homogenised trends in dew point temperature Comparison of decadal trends from 1973-2013 between raw and homogenised monthly mean dew point temperature. (a) Ratio of raw to homogenised trends for each grid box (a value of > 1 means that the raw trend is greater than the homogenised trend). (b) Scatter plot of raw verses homogenised decadal trends for each grid box with percentage presence in each quadrant. (c) Trend size distribution for raw and homogenised grid boxes. (d) Large-scale area-average time series (weighted by the cosine of the latitude) and decadal trends calculated using the median of pairwise slopes with 5th to 95th percentile slopes as con?dence intervals. Fig8_HadISDH.landTd.2.0.0.2013p_FLATgridPHADPD5by5_MAY2014_MPtrendsscat_19732013_v1.png Fig 8
Fig. 9 Comparing raw and homogenised trends in specific humidity Comparison of decadal trends from 1973-2013 between raw and homogenised monthly mean specific humidity. (a) Ratio of raw to homogenised trends for each grid box (a value of > 1 means that the raw trend is greater than the homogenised trend). (b) Scatter plot of raw verses homogenised decadal trends for each grid box with percentage presence in each quadrant. (c) Trend size distribution for raw and homogenised grid boxes. (d) Large-scale area-average time series (weighted by the cosine of the latitude) and decadal trends calculated using the median of pairwise slopes with 5th to 95th percentile slopes as con?dence intervals. Fig9_HadISDH.landq.2.0.0.2013p_FLATgridIDPHA5by5_MAY2014_MPtrendsscat_19732013_v1.png Fig 9
Fig. 10 Comparing raw and homogenised trends in relative humidity Comparison of decadal trends from 1973-2013 between raw and homogenised monthly mean relative humidity. (a) Ratio of raw to homogenised trends for each grid box (a value of > 1 means that the raw trend is greater than the homogenised trend). (b) Scatter plot of raw verses homogenised decadal trends for each grid box with percentage presence in each quadrant. (c) Trend size distribution for raw and homogenised grid boxes. (d) Large-scale area-average time series (weighted by the cosine of the latitude) and decadal trends calculated using the median of pairwise slopes with 5th to 95th percentile slopes as con?dence intervals. Fig10_HadISDH.landRH.2.0.0.2013p_FLATgridIDPHA5by5_MAY2014_MPtrendsscat_19732013_v1.png Fig 10
Fig. 11 Gridbox decadal trends 1973-2013 Decadal trends from 1973-2013 for monthly mean climate anomalies relative to 1976-2005. Trends are ?tted using the median of pairwise slopes. Black dots show trends that are considered to be significantly different from a zero trend where the 5th and 95th percentiles of the pairwise slopes are in the same direction. (a) Temperature, (b) dew point temperature, (c) specific humidity and (d) relative humidity. See Fig. SM6 for other variables. Fig11a_HadISDH.landT.2.0.0.2013p_FLATgridIDPHA5by5_SEP2014_MPdectrends_19732013_v1.png, Fig11b_HadISDH.landTd.2.0.0.2013p_FLATgridPHADPD5by5_SEP2014_MPdectrends_19732013_v1.png, Fig11c_HadISDH.landq.2.0.0.2013p_FLATgridIDPHA5by5_SEP2014_MPdectrends_19732013_v1.png, Fig11d_HadISDH.landRH.2.0.0.2013p_FLATgridIDPHA5by5_SEP2014_MPdectrends_19732013_v1.png, Fig 11a Fig 11b Fig 11c Fig 11d
Fig. 12 Global average time series Comparison global (70° N to 70° S) average time series of annual mean climate anomalies from HadISDH.2.0.0 with difference series for other monitoring products spatially matched to HadISDH grid box size and coverage. HadISDH uncertainty ranges for station (blue), sampling (red) and coverage (gold) uncertainty contributions combined are shown. All series have been given a zero mean over the common 1981-2010 period. Decadal trends with 5th to 95th percentile confidence ranges are shown for HadISDH and correlation coeficients of the standardised climate anomalies (not detrended) are shown. (a, h) Temperature, (b, i) wet bulb temperature, (c, j) dew point temperature, (d, k) specific humidity, (e, l) vapour pressure, (f, m) relative humidity, (g, n) dew point depression. Fig12_PlotCompsTimeSeries.2.0.0.2013p_annual_SEP2014_v1.png Fig 12
Fig. 13 Regional average time series HadISDH.2.0.0 regional average time series for all variables for with combined uncertainties broken down by coverage (gold), temporal sampling (red) and station (gold). Decadal trends are shown with 5th to 95th percentile confidence intervals. (a, h, o) Temperature, (b, i, p) wet bulb temperature, (c, j, q) dew point temperature, (d, k, r) specific humidity, (e, l, s) vapour pressure, (f, m ,t) relative humidity, (g, n, u) dew point depression. Fig13_PlotRegionsTimeSeries.2.0.0.2013p_annual_MAY2014_v1.png Fig 13

Further plots and text will appear here over time