Evaluating the role of deep granitic rocks in generating anomalous temperatures in SW-Queensland

Across central Australia and south-west Queensland, a large (~800,000km2) subsurface temperature anomaly occurs (Figure 1). Temperatures are interpreted to be greater than 235°C at 5km depth, ca. 85°C higher than the average geothermal gradient for the upper continental crust (Chopra & Holgate, 2005; Holgate & Gerner, 2011). This anomaly has driven the development of Engineered Geothermal Systems (EGS) at Innamincka, where high temperatures have been related to the radiogenic heat production of High Heat Producing Granites (HHPG) at depth, below thermally insulative sedimentary cover (Chopra & Holgate, 2005; Draper & D’Arcy, 2006; Meixner & Holgate, 2009). To evaluate the role of granitic rocks at depth in generating the broader temperature anomaly in SW-Queensland, we sampled 25 granitic rocks from basement intervals of petroleum drill cores below thermal insulative cover along two transects (WNW–ESE and NNE–SSW — Figure 1) and performed a multidisciplinary study involving petrography, whole-rock chemistry, zircon dating and thermal conductivity measurements.

Two modes of dipole events in tropical Indian Ocean

By analyzing the distributions of subsurface temperature and the surface wind stress anomalies in the tropical Pacific and Indian Oceans during the Indian Ocean Dipole (IOD) events, two major modes of the IOD and their formation mechanisms are revealed. (1) The subsurface temperature anomaly (STA) in the tropical Indian Ocean during the IOD events can be described as a "<" -shaped and west-east-oriented dipole pattern; in the east side of the "<" pattern, a notable tongue-like STA extends westward along the equator in the tropical eastern Indian Ocean; while in the west side of the "<" pattern, the STA has opposite sign with two centers (the southern one is stronger than the northern one in intensity) being of rough symmetry about the equator in the tropical mid-western Indian Ocean. (2) The IOD events are composed of two modes, which have similar spatial pattern but different temporal variabilities due to the large scale air-sea interactions within two independent systems. The first mode of the IOD event originates from the air-sea interaction on a scale of the tropical Pacific-Indian Ocean and coexists with ENSO. The second mode originates from the air-sea interaction on a scale of the tropical Indian Ocean and is closely associated with changes in the position and intensity of the Mascarene high pressure. The strong IOD event occurs when the two modes are in phase, and the IOD event weakens or disappears when the two modes are out of phase. Besides, the IOD events are normally strong when either of the two modes is strong. (3) The IOD event is caused by the abnormal wind stress forcing over the tropical Indian Ocean, which results in vertical transports, leading to the upwelling and pileup of seawater. This is the main dynamic processes resulting in the STA. When the anomalous easterly exists over the equatorial Indian Ocean, the cold waters upwell in the tropical eastern Indian Ocean while the warm waters pileup in the tropical western Indian Ocean, hence the thermocline in the tropical Indian Ocean is shallowed in the east and deepened in the west. The off-equator component due to the Coriolis force in the equatorial area causes the upwelling of cold waters and the shallowing of the equatorial India Ocean thermocline. On the other hand, the anomalous anticyclonic circulations and their curl fields located on both sides of the equator, cause the pileup of warm waters in the central area of their curl fields and the deepening of the equatorial Indian Ocean thermocline off the equator. The above three factors lead to the occurrence of positive phase IOD events. When anomalous westerly dominates over the tropical Indian Ocean, the dynamic processes are reversed, and the negative-phase IOD event occurs.

Subsurface ocean temperature signals on west coast of North America [abstract]

EXTRACT (SEE PDF FOR FULL ABSTRACT): Data were extracted from the U.S. Navy Fleet Numerical Oceanographic Center Master Oceanographic Observation Data Set for a 200 km to 300 km wide coastal strip on the west coast of the United States. These data were averaged for the September through February (winter) and March through August (summer) intervals. The resulting winter temperature anomaly values show the El Nino signal in the CCC [Coastal California Current] as positive temperature anomalies from the surface to at least 300 m.

Heat-producing crust regulation of subsurface temperatures : a stochastic model re-evaluation of the geothermal potential in western Queensland, Australia

A large subsurface, elevated temperature anomaly is well documented in Central Australia. High Heat Producing Granites (HHPGs) intersected by drilling at Innamincka are often assumed to be the dominant cause of the elevated subsurface temperatures, although their presence in other parts of the temperature anomaly has not been confirmed. Geological controls on the temperature anomaly remain poorly understood. Additionally, methods previously used to predict temperature at 5 km depth in this area are simplistic and possibly do not give an accurate representation of the true distribution and magnitude of the temperature anomaly. Here we re-evaluate the geological controls on geothermal potential in the Queensland part of the temperature anomaly using a stochastic thermal model. The results illustrate that the temperature distribution is most sensitive to the thermal conductivity structure of the top 5 km. Furthermore, the results indicate the presence of silicic crust enriched in heat producing elements between and 40 km.

The response of a uniform horizontal temperature gradient to heating

The response of a uniform horizontal temperature gradient to prescribed fixed heating is calculated in the context of an extended version of surface quasigeostrophic dynamics. It is found that for zero mean surface flow and weak cross-gradient structure the prescribed heating induces a mean temperature anomaly proportional to the spatial Hilbert transform of the heating. The interior potential vorticity generated by the heating enhances this surface response. The time-varying part is independent of the heating and satisfies the usual linearized surface quasigeostrophic dynamics. It is shown that the surface temperature tendency is a spatial Hilbert transform of the temperature anomaly itself. It then follows that the temperature anomaly is periodically modulated with a frequency proportional to the vertical wind shear. A strong local bound on wave energy is also found. Reanalysis diagnostics are presented that indicate consistency with key findings from this theory.

Recent changes in solar outputs and the global mean surface temperature. III. Analysis of contributions to global mean air surface temperature rise

A multivariate fit to the variation in global mean surface air temperature anomaly over the past half century is presented. The fit procedure allows for the effect of response time on the waveform, amplitude and lag of each radiative forcing input, and each is allowed to have its own time constant. It is shown that the contribution of solar variability to the temperature trend since 1987 is small and downward; the best estimate is -1.3% and the 2sigma confidence level sets the uncertainty range of -0.7 to -1.9%. The result is the same if one quantifies the solar variation using galactic cosmic ray fluxes (for which the analysis can be extended back to 1953) or the most accurate total solar irradiance data composite. The rise in the global mean air surface temperatures is predominantly associated with a linear increase that represents the combined effects of changes in anthropogenic well-mixed greenhouse gases and aerosols, although, in recent decades, there is also a considerable contribution by a relative lack of major volcanic eruptions. The best estimate is that the anthropogenic factors contribute 75% of the rise since 1987, with an uncertainty range (set by the 2sigma confidence level using an AR(1) noise model) of 49–160%; thus, the uncertainty is large, but we can state that at least half of the temperature trend comes from the linear term and that this term could explain the entire rise. The results are consistent with the intergovernmental panel on climate change (IPCC) estimates of the changes in radiative forcing (given for 1961–1995) and are here combined with those estimates to find the response times, equilibrium climate sensitivities and pertinent heat capacities (i.e. the depth into the oceans to which a given radiative forcing variation penetrates) of the quasi-periodic (decadal-scale) input forcing variations. As shown by previous studies, the decadal-scale variations do not penetrate as deeply into the oceans as the longer term drifts and have shorter response times. Hence, conclusions about the response to century-scale forcing changes (and hence the associated equilibrium climate sensitivity and the temperature rise commitment) cannot be made from studies of the response to shorter period forcing changes.

Instability of surface-temperature filaments in strain and shear

The effects of uniform straining and shearing on the stability of a surface quasi-geostrophic temperature filament are investigated. Straining is shown to stabilize perturbations for wide filaments but only for a finite time until the filament thins to a critical width, after which some perturbations can grow. No filament can be stabilized in practice, since there are perturbations that can grow large for any strain rate. The optimally growing perturbations, defined as solutions that reach a certain threshold amplitude first, are found numerically for a wide range of parameter values. The radii of the vortices formed through nonlinear roll-up are found to be proportional to θ/s, where θ is the temperature anomaly of the filament and s the strain rate, and are not dependent on the initial size of the filament. Shearing is shown to reduce the normal-mode growth rates, but it cannot stabilize them completely when there are temperature discontinuities in the basic state; smooth filaments can be stabilized completely by shearing and a simple scaling argument provides the shear rate required. Copyright © 2010 Royal Meteorological Society

A simple relationship between volcanic sulfate aerosol optical depth and surface temperature change simulated in an atmosphere-ocean general circulation model.

In this study we quantify the relationship between the aerosol optical depth increase from a volcanic eruption and the severity of the subsequent surface temperature decrease. This investigation is made by simulating 10 different sizes of eruption in a global circulation model (GCM) by changing stratospheric sulfate aerosol optical depth at each time step. The sizes of the simulated eruptions range from Pinatubo‐sized up to the magnitude of supervolcanic eruptions around 100 times the size of Pinatubo. From these simulations we find that there is a smooth monotonic relationship between the global mean maximum aerosol optical depth anomaly and the global mean temperature anomaly and we derive a simple mathematical expression which fits this relationship well. We also construct similar relationships between global mean aerosol optical depth and the temperature anomaly at every individual model grid box to produce global maps of best‐fit coefficients and fit residuals. These maps are used with caution to find the eruption size at which a local temperature anomaly is clearly distinct from the local natural variability and to approximate the temperature anomalies which the model may simulate following a Tambora‐sized eruption. To our knowledge, this is the first study which quantifies the relationship between aerosol optical depth and resulting temperature anomalies in a simple way, using the wealth of data that is available from GCM simulations.

An investigation into the impact of using various techniques to estimate Arctic surface air temperature anomalies

Time series of global and regional mean Surface Air Temperature (SAT) anomalies are a common metric used to estimate recent climate change. Various techniques can be used to create these time series from meteorological station data. The degree of difference arising from using five different techniques, based on existing temperature anomaly dataset techniques, to estimate Arctic SAT anomalies over land and sea ice were investigated using reanalysis data as a testbed. Techniques which interpolated anomalies were found to result in smaller errors than non-interpolating techniques relative to the reanalysis reference. Kriging techniques provided the smallest errors in estimates of Arctic anomalies and Simple Kriging was often the best kriging method in this study, especially over sea ice. A linear interpolation technique had, on average, Root Mean Square Errors (RMSEs) up to 0.55 K larger than the two kriging techniques tested. Non-interpolating techniques provided the least representative anomaly estimates. Nonetheless, they serve as useful checks for confirming whether estimates from interpolating techniques are reasonable. The interaction of meteorological station coverage with estimation techniques between 1850 and 2011 was simulated using an ensemble dataset comprising repeated individual years (1979-2011). All techniques were found to have larger RMSEs for earlier station coverages. This supports calls for increased data sharing and data rescue, especially in sparsely observed regions such as the Arctic.

Subsurface urban heat island and its effects on horizontal ground-source heat pump potential under climate change

Recent urban air temperature increase is attributable to the climate change and heat island effects due to urbanization. This combined effects of urbanization and global warming can penetrate into the underground and elevate the subsurface temperature. In the present study, over-100 years measurements of subsurface temperature at a remote rural site were analysed, and an increasing rate of 0.17⁰C per decade at soil depth of 30cm due to climate change was identified in the UK, but the subsurface warming in an urban site showed a much higher rate of 0.85⁰C per decade at a 30cm depth and 1.18⁰C per decade at 100cm. The subsurface urban heat island (SUHI) intensity obtained at the paired urban-rural stations in London showed an unique 'U-shape', i.e. lowest in summer and highest during winter. The maximum SUHII is 3.5⁰C at 6:00 AM in December, and the minimum UHII is 0.2⁰C at 18:00PM in July. Finally, the effects of SUHI on the energy efficiency of the horizontal ground source heat pump (GSHP) were determined. Provided the same heat pump used, the installation at an urban site will maintain an overall higher COP compared with that at a rural site in all seasons, but the highest COP improvement can be achieved in winter.

Representing the Arctic in global surface temperature time series of recent climate change

The Arctic is an important region in the study of climate change, but monitoring surface temperatures in this region is challenging, particularly in areas covered by sea ice. Here in situ, satellite and reanalysis data were utilised to investigate whether global warming over recent decades could be better estimated by changing the way the Arctic is treated in calculating global mean temperature. The degree of difference arising from using five different techniques, based on existing temperature anomaly dataset techniques, to estimate Arctic SAT anomalies over land and sea ice were investigated using reanalysis data as a testbed. Techniques which interpolated anomalies were found to result in smaller errors than non-interpolating techniques. Kriging techniques provided the smallest errors in anomaly estimates. Similar accuracies were found for anomalies estimated from in situ meteorological station SAT records using a kriging technique. Whether additional data sources, which are not currently utilised in temperature anomaly datasets, would improve estimates of Arctic surface air temperature anomalies was investigated within the reanalysis testbed and using in situ data. For the reanalysis study, the additional input anomalies were reanalysis data sampled at certain supplementary data source locations over Arctic land and sea ice areas. For the in situ data study, the additional input anomalies over sea ice were surface temperature anomalies derived from the Advanced Very High Resolution Radiometer satellite instruments. The use of additional data sources, particularly those located in the Arctic Ocean over sea ice or on islands in sparsely observed regions, can lead to substantial improvements in the accuracy of estimated anomalies. Decreases in Root Mean Square Error can be up to 0.2K for Arctic-average anomalies and more than 1K for spatially resolved anomalies. Further improvements in accuracy may be accomplished through the use of other data sources.

Organic temperature proxies in the subpolat region around Iceland

Subpolar regions are key areas to study natural climate variability, due to their high sensitivity to rapid environmental changes, particularly through sea surface temperature (SST) variations. Here, we have tested three independent organic temperature proxies (UK'37, TEX86 and LDI) on their potential applicability for SST reconstruction in the subpolar region around Iceland. UK'37, TEX86 and TEXL86 temperature estimates from suspended particulate matter showed a substantial discrepancy with instrumental data, while long chain alkyl diols were below detection limit in most of the stations. In the northern Iceland Basin, sedimenting particles revealed a seasonality in lipid fluxes i.e. high fluxes of alkenones and GDGTs were measured during late spring-summer, and high fluxes of long chain alkyl diols during late summer. The flux-weighted average temperature estimates had a significant negative (ca. 2.3°C for UK'37) and positive (up to 5°C for TEX86) offset with satellite-derived SSTs and temperature estimates derived from the underlying surface sediment. UK'37 temperature estimates from surface sediments around Iceland correlate well with summer mean sea surface temperatures, while TEX86 derived temperatures correspond with both annual and winter mean 0-200 m temperatures, suggesting a subsurface temperature signal. Anomalous LDI-SST values in surface sediments, and low mass flux of 1,13- and 1,15-diols compared to 1,14-diols, suggest that Proboscia diatoms are the major sources of long chain alkyl diols in this area rather than eustigmatophyte algae, and therefore the LDI cannot be applied in this region.