Tests of sunspot number sequences: 3. Effects of regression procedures on the calibration of historic sunspot data

We use sunspot group observations from the Royal Greenwich Observatory (RGO) to investigate the effects of intercalibrating data from observers with different visual acuities. The tests are made by counting the number of groups RB above a variable cut-off threshold of observed total whole-spot area (uncorrected for foreshortening) to simulate what a lower acuity observer would have seen. The synthesised annual means of RB are then re-scaled to the full observed RGO group number RA using a variety of regression techniques. It is found that a very high correlation between RA and RB (rAB > 0.98) does not prevent large errors in the intercalibration (for example sunspot maximum values can be over 30 % too large even for such levels of rAB). In generating the backbone sunspot number (RBB), Svalgaard and Schatten (2015, this issue) force regression fits to pass through the scatter plot origin which generates unreliable fits (the residuals do not form a normal distribution) and causes sunspot cycle amplitudes to be exaggerated in the intercalibrated data. It is demonstrated that the use of Quantile-Quantile (“Q  Q”) plots to test for a normal distribution is a useful indicator of erroneous and misleading regression fits. Ordinary least squares linear fits, not forced to pass through the origin, are sometimes reliable (although the optimum method used is shown to be different when matching peak and average sunspot group numbers). However, other fits are only reliable if non-linear regression is used. From these results it is entirely possible that the inflation of solar cycle amplitudes in the backbone group sunspot number as one goes back in time, relative to related solar-terrestrial parameters, is entirely caused by the use of inappropriate and non-robust regression techniques to calibrate the sunspot data.

Robust future changes in temperature variability under greenhouse gas forcing and the relationship with thermal advection

Recent temperature extremes have highlighted the importance of assessing projected changes in the variability of temperature as well as the mean. A large fraction of present day temperature variance is associated with thermal advection, as anomalous winds blow across the land-sea temperature contrast for instance. Models project robust heterogeneity in the 21st century warming pattern under greenhouse gas forcing, resulting in land-sea temperature contrasts increasing in summer and decreasing in winter, and the pole-to-equator temperature gradient weakening in winter. In this study, future monthly variability changes in the 17 member ensemble ESSENCE are assessed. In winter, variability in midlatitudes decreases while in very high latitudes and the tropics it increases. In summer, variability increases over most land areas and in the tropics, with decreasing variability in high latitude oceans. Multiple regression analysis is used to determine the contributions to variability changes from changing temperature gradients and circulation patterns. Thermal advection is found to be of particular importance in the northern hemisphere winter midlatitudes, where the change in mean state temperature gradients alone could account for over half the projected changes. Changes in thermal advection are also found to be important in summer in Europe and coastal areas, although less so than in winter. Comparison with CMIP5 data shows that the midlatitude changes in variability are robust across large regions, particularly high northern latitudes in winter and mid northern latitudes in summer.

Analysis of the distribution of the number of bidders in construction contract auctions

The number of bidders, N, involved in a construction procurement auction is known to have an important effect on the value of the lowest bid and the mark-up applied by bidders. In practice, for example, it is important for a bidder to have a good estimate of N when bidding for a current contract. One approach, instigated by Friedman in 1956, is to make such an estimate by statistical analysis and modelling. Since then, however, finding a suitable model for N has been an enduring problem for researchers and, despite intensive research activity in the subsequent 30 years, little progress has been made, due principally to the absence of new ideas and perspectives. The debate is resumed by checking old assumptions, providing new evidence relating to concomitant variables and proposing a new model. In doing this and in order to ensure universality, a novel approach is developed and tested by using a unique set of 12 construction tender databases from four continents. This shows the new model provides a significant advancement on previous versions. Several new research questions are also posed and other approaches identified for future study.

Contrasting fast precipitation responses to tropospheric and stratospheric ozone forcing

The precipitation response to radiative forcing (RF) can be decomposed into a fast precipitation response (FPR), which depends on the atmospheric component of RF, and a slow response, which depends on surface temperature change. We present the first detailed climate model study of the FPR due to tropospheric and stratospheric ozone changes. The FPR depends strongly on the altitude of ozone change. Increases below about 3 km cause a positive FPR; increases above cause a negative FPR. The FPR due to stratospheric ozone change is, per unit RF, about 3 times larger than that due to tropospheric ozone. As historical ozone trends in the troposphere and stratosphere are opposite in sign, so too are the FPRs. Simple climate model calculations of the time-dependent total (fast and slow) precipitation change, indicate that ozone's contribution to precipitation change in 2011, compared to 1765, could exceed 50% of that due to CO2 change.

Projected changes in medicanes in the HadGEM3 N512 high-resolution global climate model

Medicanes or “Mediterranean hurricanes” represent a rare and physically unique type of Mediterranean mesoscale cyclone. There are similarities with tropical cyclones with regard to their development (based on the thermodynamical disequilibrium between the warm sea and the overlying troposphere) and their kinematic and thermodynamical properties (medicanes are intense vortices with a warm core and even a cloud-free eye). Although medicanes are smaller and their wind speeds are lower than in tropical cyclones, the severity of their winds can cause substantial damage to islands and coastal areas. Concern about how human-induced climate change will affect extreme events is increasing. This includes the future impacts on medicanes due to the warming of the Mediterranean waters and the projected changes in regional atmospheric circulation. However, most global climate models do not have high enough spatial resolution to adequately represent small features such as medicanes. In this study, a cyclone tracking algorithm is applied to high resolution global climate model data with a horizontal grid resolution of approximately 25 km over the Mediterranean region. After a validation of the climatology of general Mediterranean mesoscale cyclones, changes in medicanes are determined using climate model experiments with present and future forcing. The magnitude of the changes in the winds, frequency and location of medicanes is assessed. While no significant changes in the total number of Mediterranean mesoscale cyclones are found, medicanes tend to decrease in number but increase in intensity. The model simulation suggests that medicanes tend to form more frequently in the Gulf of Lion–Genoa and South of Sicily.

Southern Ocean isopycnal mixing and ventilation changes driven by winds

Observed and predicted changes in the strength of the westerly winds blowing over the Southern Ocean have motivated a number of studies of the response of the Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation (MOC) to wind perturbations and led to the discovery of the``eddy-compensation" regime, wherein the MOC becomes insensitive to wind changes. In addition to the MOC, tracer transport also depends on mixing processes. Here we show, in a high-resolution process model, that isopycnal mixing by mesoscale eddies is strongly dependent on the wind strength. This dependence can be explained by mixing-length theory and is driven by increases in eddy kinetic energy; the mixing length does not change strongly in our simulation. Simulation of a passive ventilation tracer (analogous to CFCs or anthropogenic CO$_2$) demonstrates that variations in tracer uptake across experiments are dominated by changes in isopycnal mixing, rather than changes in the MOC. We argue that, to properly understand tracer uptake under different wind-forcing scenarios, the sensitivity of isopycnal mixing to winds must be accounted for.

Impact of atmospheric forcing data on simulations of the Laptev Sea polynya dynamics using the sea-ice ocean model FESOM

The polynyas of the Laptev Sea are regions of particular interest due to the strong formation of Arctic sea-ice. In order to simulate the polynya dynamics and to quantify ice production, we apply the Finite Element Sea-Ice Ocean Model FESOM. In previous simulations FESOM has been forced with daily atmospheric NCEP (National Centers for Environmental Prediction) 1. For the periods 1 April to 9 May 2008 and 1 January to 8 February 2009 we examine the impact of different forcing data: daily and 6-hourly NCEP reanalyses 1 (1.875° x 1.875°), 6-hourly NCEP reanalyses 2 (1.875° x 1.875°), 6-hourly analyses from the GME (Global Model of the German Weather Service) (0.5° x 0.5°) and high-resolution hourly COSMO (Consortium for Small-Scale Modeling) data (5 km x 5 km). In all FESOM simulations, except for those with 6-hourly and daily NCEP 1 data, the openings and closings of polynyas are simulated in principle agreement with satellite products. Over the fast-ice area the wind fields of all atmospheric data are similar and close to in situ measurements. Over the polynya areas, however, there are strong differences between the forcing data with respect to air temperature and turbulent heat flux. These differences have a strong impact on sea-ice production rates. Depending on the forcing fields polynya ice production ranges from 1.4 km3 to 7.8 km3 during 1 April to 9 May 2011 and from 25.7 km3 to 66.2 km3 during 1 January to 8 February 2009. Therefore, atmospheric forcing data with high spatial and temporal resolution which account for the presence of the polynyas are needed to reduce the uncertainty in quantifying ice production in polynyas.

The role of local atmospheric forcing on the modulation of the ocean mixed layer depth in reanalysis and a coupled single column ocean model

The role of the local atmospheric forcing on the ocean mixed layer depth (MLD) over the global oceans is studied using ocean reanalysis data products and a single-column ocean model coupled to an atmospheric general circulation model. The focus of this study is on how the annual mean and the seasonal cycle of the MLD relate to various forcing characteristics in different parts of the world's ocean, and how anomalous variations in the monthly mean MLD relate to anomalous atmospheric forcings. By analysing both ocean reanalysis data and the single-column ocean model, regions with different dominant forcings and different mean and variability characteristics of the MLD can be identified. Many of the global oceans' MLD characteristics appear to be directly linked to different atmospheric forcing characteristics at different locations. Here, heating and wind-stress are identified as the main drivers; in some, mostly coastal, regions the atmospheric salinity forcing also contributes. The annual mean MLD is more closely related to the annual mean wind-stress and the MLD seasonality is more closely to the seasonality in heating. The single-column ocean model, however, also points out that the MLD characteristics over most global ocean regions, and in particular the tropics and subtropics, cannot be maintained by local atmospheric forcings only, but are also a result of ocean dynamics that are not simulated in a single-column ocean model. Thus, lateral ocean dynamics are essentially in correctly simulating observed MLD.

Small global-mean cooling due to volcanic radiative forcing

In both the observational record and atmosphere-ocean general circulation model (AOGCM) simulations of the last ∼∼ 150 years, short-lived negative radiative forcing due to volcanic aerosol, following explosive eruptions, causes sudden global-mean cooling of up to ∼∼ 0.3 K. This is about five times smaller than expected from the transient climate response parameter (TCRP, K of global-mean surface air temperature change per W m−2 of radiative forcing increase) evaluated under atmospheric CO2 concentration increasing at 1 % yr−1. Using the step model (Good et al. in Geophys Res Lett 38:L01703, 2011. doi:10.​1029/​2010GL045208), we confirm the previous finding (Held et al. in J Clim 23:2418–2427, 2010. doi:10.​1175/​2009JCLI3466.​1) that the main reason for the discrepancy is the damping of the response to short-lived forcing by the thermal inertia of the upper ocean. Although the step model includes this effect, it still overestimates the volcanic cooling simulated by AOGCMs by about 60 %. We show that this remaining discrepancy can be explained by the magnitude of the volcanic forcing, which may be smaller in AOGCMs (by 30 % for the HadCM3 AOGCM) than in off-line calculations that do not account for rapid cloud adjustment, and the climate sensitivity parameter, which may be smaller than for increasing CO2 (40 % smaller than for 4 × CO2 in HadCM3).

Tests of sunspot number sequences: 1. Using ionosonde data

More than 70 years ago it was recognised that ionospheric F2-layer critical frequencies [foF2] had a strong relationship to sunspot number. Using historic datasets from the Slough and Washington ionosondes, we evaluate the best statistical fits of foF2 to sunspot numbers (at each Universal Time [UT] separately) in order to search for drifts and abrupt changes in the fit residuals over Solar Cycles 17-21. This test is carried out for the original composite of the Wolf/Zürich/International sunspot number [R], the new “backbone” group sunspot number [RBB] and the proposed “corrected sunspot number” [RC]. Polynomial fits are made both with and without allowance for the white-light facular area, which has been reported as being associated with cycle-to-cycle changes in the sunspot number - foF2 relationship. Over the interval studied here, R, RBB, and RC largely differ in their allowance for the “Waldmeier discontinuity” around 1945 (the correction factor for which for R, RBB and RC is, respectively, zero, effectively over 20 %, and explicitly 11.6 %). It is shown that for Solar Cycles 18-21, all three sunspot data sequences perform well, but that the fit residuals are lowest and most uniform for RBB. We here use foF2 for those UTs for which R, RBB, and RC all give correlations exceeding 0.99 for intervals both before and after the Waldmeier discontinuity. The error introduced by the Waldmeier discontinuity causes R to underestimate the fitted values based on the foF2 data for 1932-1945 but RBB overestimates them by almost the same factor, implying that the correction for the Waldmeier discontinuity inherent in RBB is too large by a factor of two. Fit residuals are smallest and most uniform for RC and the ionospheric data support the optimum discontinuity multiplicative correction factor derived from the independent Royal Greenwich Observatory (RGO) sunspot group data for the same interval.

Constraining the aerosol influence on cloud fraction

Aerosol-cloud interactions have the potential to modify many different cloud properties. There is significant uncertainty in the strength of these aerosol-cloud interactions in analyses of observational data, partly due to the difficulty in separating aerosol effects on clouds from correlations generated by local meteorology. The relationship between aerosol and cloud fraction (CF) is particularly important to determine, due to the strong correlation of CF to other cloud properties and its large impact on radiation. It has also been one of the hardest to quantify from satellites due to the strong meteorological covariations involved. This work presents a new method to analyze the relationship between aerosol optical depth (AOD) and CF. By including information about the cloud droplet number concentration (CDNC), the impact of the meteorological covariations is significantly reduced. This method shows that much of the AOD-CF correlation is explained by relationships other than that mediated by CDNC. By accounting for these, the strength of the global mean AOD-CF relationship is reduced by around 80%. This suggests that the majority of the AOD-CF relationship is due to meteorological covariations, especially in the shallow cumulus regime. Requiring CDNC to mediate the AOD-CF relationship implies an effective anthropogenic radiative forcing from an aerosol influence on liquid CF of −0.48 W m−2 (−0.1 to −0.64 W m−2), although some uncertainty remains due to possible biases in the CDNC retrievals in broken cloud scenes.

Fast and slow precipitation responses to individual climate forcers: a PDRMIP multi-model study

Precipitation is expected to respond differently to various drivers of anthropogenic climate change. We present the first results from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), where nine global climate models have perturbed CO2, CH4, black carbon, sulfate, and solar insolation. We divide the resulting changes to global mean and regional precipitation into fast responses that scale with changes in atmospheric absorption and slow responses scaling with surface temperature change. While the overall features are broadly similar between models, we find significant regional intermodel variability, especially over land. Black carbon stands out as a component that may cause significant model diversity in predicted precipitation change. Processes linked to atmospheric absorption are less consistently modeled than those linked to top-of-atmosphere radiative forcing. We identify a number of land regions where the model ensemble consistently predicts that fast precipitation responses to climate perturbations dominate over the slow, temperature-driven responses.

What controls the vertical distribution of aerosol? Relationships between process sensitivity in HadGEM3–UKCA and inter-model variation from AeroCom Phase II

The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3–UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3–UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN  >  3 nm), while the profiles of larger particles (e.g. CN  >  100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AOD-normalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions.

ENSO related variation of equatorial MRG and Rossby waves and forcing from higher latitudes

The contrasting behaviour of westward-moving mixed Rossby-gravity (WMRG) and the first Rossby (R1) waves in El Niño (EN) and La Niña (LN) seasons is documented with a focus on the Northern Hemisphere winter. The eastward-moving variance in the upper troposphere is dominated by WMRG and R1 structures that appear to be Doppler-shifted by the flow and are referred to as WMRG-E and R1-E. In the East Pacific and Atlantic the years with stronger equatorial westerly winds have the stronger WMRG and WMRG- E. In the East Pacific, R1 is also a maximum in LN. However, R1-E exhibits an eastward-shift between LN and EN. The changes with ENSO phase provide a test-bed for the understanding of these waves. In the East Pacific and Atlantic, the stronger WMRG-E and WMRG with stronger westerlies are in accord with the dispersion relation with simple Doppler-shifting by the zonal flow. The possible existence of free waves can also explain stronger R1 in EN in the Eastern Hemisphere. 1-D free wave propagation theory based on wave activity conservation is also important for R1. However, this theory is unable to explain the amplitude maxima for other waves observed in the strong equatorial westerly regions in the Western Hemisphere, and certainly not their ENSO-related variation. The forcing of equatorial waves by higher latitude wave activity and its variation with ENSO phase is therefore examined. Propagation of extratropical eastward-moving Rossby wave activity through the westerly ducts into the equatorial region where it triggers WMRG-E is favoured in the stronger westerlies, in LN in the East Pacific and EN in the Atlantic. It is also found that WMRG is forced by Southern Hemisphere westward-moving wavetrains arching into the equatorial region where they are reflected. The most significant mechanism for both R1 and R1-E appear to be lateral forcing by subtropical wavetrains.