Investigating the use of early satellite data to test historical reconstructions of sea surface temperature

Email: t.hall@pgr.reading.ac.uk

Observations of sea surface temperature (SST) form one of the key components of the climate record. There are a number of different in-situ based reconstructions of SST extending back over 150 years, but they are not truly independent of each other because the observations they are based on are largely the same (Berry et al., 2018). Datasets of SST retrieved from satellite radiometers exist for the 1980s onwards, providing an independent record of SST. Before this, SST reconstructions are based on sparse, ship-based measurements.

There were meteorological measurements being made from satellites in the 1960s and 70s, however, some of which can potentially be used to retrieve SST. My PhD focuses on investigating if we can retrieve SST from one of these early satellite instruments, to test the reliability of the SST reconstructions before the 1980s. This instrument is the Infrared Interferometer Spectrometer (IRIS), which made measurements of atmospheric emission spectra on-board the Nimbus 4 satellite from April 1970 to January 1971. IRIS had over 800 thermal infrared (IR) channels, covering the 400-1600 cm-1 spectral region. Figure 1 shows an example of two typical IRIS radiance spectra, with the main spectral absorption features labelled as well as the atmospheric window regions, which are the main spectral regions used for SST retrieval.

blog_fig1
Figure 1: Example of two typical IRIS radiance spectra; the main spectral absorption features are labelled as well as the atmospheric window regions.

Before using the IRIS data to retrieve SST, it was necessary to apply a series of quality assurance tests to filter out bad data. A few months into my PhD, work by Bantges et al. (2016) revealed evidence for a wavelength dependent cold bias of up to 2K in the data. A large part of my PhD was spent trying to quantify this bias. This was done by comparing clear-sky IRIS spectra with spectra simulated with a radiative transfer model. Unfortunately, this meant that the SSTs eventually retrieved from IRIS are not totally independent of the SST reconstructions as the simulations are based on reanalysis data forced by the HadISST2 reconstruction. Figure 2 compares our estimate of the IRIS spectral bias with globally averaged spectral differences between IRIS, the Interferometric Monitor for Greenhouse Gases (IMG) and the Infrared Atmospheric Sounding Instrument (IASI) from Bantges et al. (2016). This shows close agreement between our bias estimate and the IRIS-IMG and IRIS-IASI differences outside of the ozone spectral region, which is not relevant for SST retrieval.

It cannot just be assumed that the bias is the same for each IRIS measurement. Comparison of IRIS (bias-corrected using our initial bias estimate) with window channel data from the Temperature-Humidity Infrared Radiometer (THIR), also on-board Nimbus 4, reveals that the relative IRIS-THIR bias varies with window brightness temperature and orbit angle. The THIR, however, may have biases of its own, so these biases cannot be attributed to IRIS.

blog_fig2
Figure 2: Area-weighted global mean brightness temperature difference averaged over AMJ between IRIS (1970), IMG (1997) and IASI (2012) (black and blue lines) from Bantges et al. (2016), compared with our IRIS bias estimate, also area-weighted and averaged over AMJ (red line). The ozone absorption band is not used for SST retrieval, so is shaded grey.

The technique of optimal estimation was applied to retrieve SST from IRIS. This uses the observation-simulation differences together with information about the sensitivity of the simulations to the state of the atmosphere and ocean to estimate the SST. IR satellite retrievals of SST are usually performed in clear-sky conditions only. However, the low spatial resolution of IRIS means that very few cases are fully clear-sky. For this reason, we had to adapt the retrieval method to be tolerant of some cloud. This involves retrieving SST simultaneously with cloud fraction (CF). The retrieval method was then tested on partly cloudy (≤0.2 CF) IASI spectra made ‘IRIS-like’ by spatial averaging, spectral smoothing and simulating IRIS-like errors. The retrieved IRIS-like SSTs were validated against quality-controlled drifting buoy SSTs. This revealed latitudinal biases in the retrieved SSTs for the partly cloudy cases, not present in the SSTs for clear-sky cases.

SSTs were then retrieved for all IRIS cases with an expected CF ≤ 0.2. Figure 3 shows the difference between the gridded, monthly average IRIS SSTs and two of the SST reconstructions (HadSST3 and HadISST2) for July 1970. There are large, spatially correlated differences between the IRIS SSTs and reconstructions. We expect a latitudinal bias in the IRIS SSTs and some level of remaining bias in the IRIS spectra is likely, contributing to further SST bias. It is therefore likely that the differences in Figure 3 are mainly due to bias in the IRIS SSTs rather than the reconstructions.

iris_recon_07_70
Figure 3: Gridded IRIS-HadSST3 (left) and IRIS-HadISST2 (right) SST for July 1970. HadISST2 is a globally complete, interpolated dataset whereas HadSST3 is not globally complete.

Despite being unable to retrieve bias-free SST estimates from IRIS, my work has contributed towards better understanding the characteristics of IRIS. This ties in with a current project aiming to recover and assess the quality of data from a number of different historic satellite sensors, including IRIS, for assimilation in the next generation of climate reanalyses.

References

Bantges, R., H. Brindley, X. H. Chen, X. L. Huang, J. Harries, J. Murray (2016), On the detection of robust multi-decadal changes in the Earth’s Outgoing Longwave Radiation spectrum. J. Climate, 29, 4939-4947. https://doi.org/10.1175/JCLI-D-15-0713.1

Berry, D. I., G. K. Corlett, O. Embury, C. J. Merchant (2018), Stability assessment of the (A)ATSR Sea Surface Temperature climate dataset from the European Space Agency Climate Change Initiative. Remote Sens., 10, 126. https://doi.org/10.3390/rs10010126

Evaluating aerosol forecasts in London

Email: e.l.warren@pgr.reading.ac.uk

Aerosols in urban areas can greatly impact visibility, radiation budgets and our health (Chen et al., 2015). Aerosols make up the liquid and solid particles in the air that, alongside noxious gases like nitrogen dioxide, are the pollution in cities that we often hear about on the news – breaking safety limits in cities across the globe from London to Beijing. Air quality researchers try to monitor and predict aerosols, to inform local councils so they can plan and reduce local emissions.

Figure 1: Smog over London (Evening Standard, 2016).

Recently, large numbers of LiDARs (Light Detection and Ranging) have been deployed across Europe, and elsewhere – in part to observe aerosols. They effectively shoot beams of light into the atmosphere, which reflect off atmospheric constituents like aerosols. From each beam, many measurements of reflectance are taken very quickly over time – and as light travels further with more time, an entire profile of reflectance can be constructed. As the penetration of light into the atmosphere decreases with distance, the reflected light is usually commonly called attenuated backscatter (β). In urban areas, measurements away from the surface like these are sorely needed (Barlow, 2014), so these instruments could be extremely useful. When it comes to predicting aerosols, numerical weather prediction (NWP) models are increasingly being considered as an option. However, the models themselves are very computationally expensive to run so they tend to only have a simple representation of aerosol. For example, for explicitly resolved aerosol, the Met Office UKV model (1.5 km) just has a dry mass of aerosol [kg kg-1] (Clark et al., 2008). That’s all. It gets transported around by the model dynamics, but any other aerosol characteristics, from size to number, need to be parameterised from the mass, to limit computation costs. However, how do we know if the estimates of aerosol from the model are actually correct? A direct comparison between NWP aerosol and β is not possible because fundamentally, they are different variables – so to bridge the gap, a forward operator is needed.

In my PhD I helped develop such a forward operator (aerFO, Warren et al., 2018). It’s a model that takes aerosol mass (and relative humidity) from NWP model output, and estimates what the attenuated backscatter would be as a result (βm). From this, βm could be directly compared to βo and the NWP aerosol output evaluated (e.g. see if the aerosol is too high or low). The aerFO was also made to be computationally cheap and flexible, so if you had more information than just the mass, the aerFO would be able to use it!

Among the aerFO’s several uses (Warren et al., 2018, n.d.), was the evaluation of NWP model output. Figure 2 shows the aerFO in action with a comparison between βm and observed attenuated backscatter (βo) measured at 905 nm from a ceilometer (a type of LiDAR) on 14th April 2015 at Marylebone Road in London. βm was far too high in the morning on this day. We found that the original scheme the UKV used to parameterise the urban surface effects in London was leading to a persistent cold bias in the morning. The cold bias would lead to a high relative humidity, so consequently the aerFO condensed more water than necessary, onto the aerosol particles as a result, causing them to swell up too much. As a result, bigger particles mean bigger βm and an overestimation. Not only was the relative humidity too high, the boundary layer in the NWP model was developing too late in the day as well. Normally, when the surface warms up enough, convection starts, which acts to mix aerosol up in the boundary layer and dilute it near the surface. However, the cold bias delayed this boundary layer development, so the aerosol concentration near the surface remained high for too long. More mass led to the aerFO parameterising larger sizes and total numbers of particles, so overestimated βm. This cold bias effect was reflected across several cases using the old scheme but was notably smaller for cases using a newer urban surface scheme called MORUSES (Met Office – Reading Urban Surface Exchange Scheme). One of the main aims for MORUSES was to improve the representation of energy transfer in urban areas, and at least to us it seemed like it was doing a better job!

Figure 2: Vertical profiles of attenuated backscatter [m−1 sr−1] (log scale) that are (a, g) observed (βo) with estimated mixing layer height (red crosses, Kotthaus and Grimmond,2018) and (b, h) forward modelled (βm) using the aerFO (section 2).(c, i) Attenuated backscatter difference (βm – βo) calculated using the hourly βm vertical profile and the vertical profile of βo nearest in time; (d, j) aerosol mass mixing ratio (m) [μg kg−1]; (e, k) relative humidity (RH) [%] and (f, l) air temperature (T) [°C] at MR on 14th April 2015.

References

Barlow, J.F., 2014. Progress in observing and modelling the urban boundary layer. Urban Clim. 10, 216–240. https://doi.org/10.1016/j.uclim.2014.03.011

Chen, C.H., Chan, C.C., Chen, B.Y., Cheng, T.J., Leon Guo, Y., 2015. Effects of particulate air pollution and ozone on lung function in non-asthmatic children. Environ. Res. 137, 40–48. https://doi.org/10.1016/j.envres.2014.11.021

Clark, P.A., Harcourt, S.A., Macpherson, B., Mathison, C.T., Cusack, S., Naylor, M., 2008. Prediction of visibility and aerosol within the operational Met Office Unified Model. I: Model formulation and variational assimilation. Q. J. R. Meteorol. Soc. 134, 1801–1816. https://doi.org/10.1002/qj.318

Warren, E., Charlton-Perez, C., Kotthaus, S., Lean, H., Ballard, S., Hopkin, E., Grimmond, S., 2018. Evaluation of forward-modelled attenuated backscatter using an urban ceilometer network in London under clear-sky conditions. Atmos. Environ. 191, 532–547. https://doi.org/10.1016/j.atmosenv.2018.04.045

Warren, E., Charlton-Perez, C., Kotthaus, S., Marenco, F., Ryder, C., Johnson, B., Lean, H., Ballard, S., Grimmond, S., n.d. Observed aerosol characteristics to improve forward-modelled attenuated backscatter. Atmos. Environ. Submitted


Understanding the Processes Involved in Electrifying Convective Clouds

All clouds within our atmosphere are charged to some extent (Nicoll & Harrison, 2016), caused by the build-up of charge at the cloud edges caused by the charge travelling from the top of the atmosphere towards the surface. On the other hand, nearly all convective clouds are actively charged, caused by charge separation mechanisms that exchange charge between different sized hydrometeors within the cloud (Saunders, 1992). If these actively charged clouds can separate enough charge, then electrical breakdown can occur in the atmosphere and lightning can be initiated.

As lightning is a substantial hazard to both life and infrastructure, understanding the processes that cause a cloud to become charged are crucial for being able to forecast it. Even though most convective clouds within the UK are charged, most will never produce lightning. This puts us in a unique position of observing clouds which can and cannot produce lightning, providing a contrast of the cloud processes involved.

Cumliform_Electrification_Hypothesis
Figure 1: A conceptual image of a cloud showing the main processes that are thought to be involved in the electrification of a convective cloud. Based on current literature. See MacGorman & Rust (1998) for an overview.

Currently, lightning must be observed before the identification of a thunderstorm can be made, leaving zero lead time. As the first lightning strike can often be the most powerful, (especially in UK winter thunderstorms), improvements to understand the processes that charge a cloud are important to provide a thunderstorm lead time.

There are many processes involved in charging a cloud, a summary of the processes that are thought to have the greatest contribution can be found in Figure 1. The separation of charge through collisions is thought to be the dominant electrification mechanism which occurs in the ice phase of the cloud. The most successful charge separation occurs between growing ice hydrometeors of different sizes (Emersic & Saunders, 2010). The growth of the ice forms an outer shell of supercooled liquid water which contains negative charge. Collisions between different sized hydrometeors cause a net exchange of mass and charge, creating positively and negatively charged ice.

The liquid phase is crucial to maintaining a high moisture content in the ice phase of the cloud brought up by the updraught of the cloud. Once the electric field within the cloud is strong enough, the liquid drops can become polarised which can also help separate charge from collisions of different sized liquid hydrometeors in the later stages of convective cloud development.

Once the charge has been separated, the different polarities must be moved to different regions of the cloud to enhance the electric field. The most well-established method involves gravitational separation (Mason & Dash, 2000). Under the assumption that smaller hydrometeors will often contain negative charge and larger hydrometeors will contain positive charge, there is a distinct separation of hydrometeor sizes. After the formation of an updraught, the smaller hydrometeors can be lifted higher into the cloud overcoming the gravitational forces. Under the right conditions, the larger hydrometeors will be too heavy, and gravity will force the hydrometeors to remain lower in the cloud.

Field_Campaign_Instruments_Small
Figure 2: The Electrostatic Biral Thunderstorm Detector (BTD-300) (a), Electrostatic JCI 131 Field Mill (b), and 35 GHz “Copernicus” Radar installed at Chilbolton Observatory, UK. The Field Mill and BTD-300 were installed on 06/10/16 and 16/11/16 respectively.

To understand if any of the processes discussed in Figure 1 are evident within real-world convective clouds, a field campaign was set-up at Chilbolton Observatory (CO) to measure the properties of the cloud. Two electrical instruments were installed at CO, the Biral Thunderstorm Detector (BTD-300, Figure 2a) and the electrostatic Field Mill (FM, Figure 2b) which can measure the displacement current (jD) and the potential gradient (PG) respectively. The 35GHz “Copernicus” dopplerised radar was used to measure the properties of the cloud. Through a two-year field campaign, 653 convective clouds were identified over CO, with 524 clouds (80.2%) found to be charged and 129 clouds (19.8%) found to be uncharged.

To understand the importance of hydrometeor size for charging a cloud, the 95th percentile of the radar reflectivity (Z) was used. The Z is strongly related to the diameter (sixth power) and number concentration of hydrometeors (first power). Figure 3 shows a boxplot of all 653 clouds, classified by the cloud charge as measured at the surface. For each cloud, the liquid and ice regions were separated (purple and blue boxplots respectively) to highlight the dominant region for charge separation.

Z_Grouped_BoxPlot_LinearScale_v22
Figure 3: A boxplot of the 95th percentile reflectivity for 635 identified clouds (black box) classified into no charge, small charge and large charge groups. The cloud was also decoupled into the liquid (purple box) and ice (blue box) phases. The box-plot shows the mean (purple bar), median (red bar) and upper and lower quartiles (upper and lower limits of black box).

For all phases of the cloud, there was a substantial increase in the reflectivity for all regions of the cloud, especially the liquid phase. This suggests that the size of the hydrometeors in both the ice and liquid phase are indeed important for charging a cloud.

In the remainder of my PhD, the relative importance of each process discussed in Figure 1 will be addressed to try and decouple each process. Further to these observations made at the surface, ten radiosonde flights will be made (from the Reading University Atmospheric Observatory) inside convectively charged clouds. Measurements of the charge, optical thickness, amount of supercooled liquid water and turbulence will be used to increase the robustness of the results presented in the previous works.

Email: james.gilmore@pgr.reading.ac.uk

References

Bouniol, D., Illingworth, A. J. & Hogan, R. J., 2003. Deriving turbulent kinetic energy dissipation rate within clouds using ground-based 94 GHz radar. Seattle, 31st International Conference on Radar Meteorology.

Emersic, C. & Saunders, C. P., 2010. Further laboratory investigations into the Relative Diffusional Growth Rate theory of thunderstorm electrification. Atmos. Res., Volume 98, pp. 327-340.

MacGorman, D. R. & Rust, D. W., 1998. The Electrical Nature of Storms. 1st ed. New York: Oxford University Press.

Mason, B. L. & Dash, J. G., 2000. Charge and mass transfer in ice–ice collisions: Experimental observations of a mechanism in thunderstorm electrification. J. Geophys. Res., Volume 105, pp. 10185-10192.

Nicoll, K. A. & Harrison, R. G., 2016. Stratiform cloud electrification: comparison of theory with multiple in-cloud measurements. Q.J.R. Meteorol. Soc., Volume 142, pp. 2679-2691.

Renzo, M. D. & Urzay, J., 2018. Aerodynamic generation of electric fields in turbulence laden with charged inertial particles. Nat. Comms., 9(1676).

Saunders, C. P. R., 1992. A Review of Thunderstorm Electrification Processes. J. App. Meteo., Volume 32, pp. 642-655.

Baroclinic and Barotropic Annular Modes of Variability

Email: l.boljka@pgr.reading.ac.uk

Modes of variability are climatological features that have global effects on regional climate and weather. They are identified through spatial structures and the timeseries associated with them (so-called EOF/PC analysis, which finds the largest variability of a given atmospheric field). Examples of modes of variability include El Niño Southern Oscillation, Madden-Julian Oscillation, North Atlantic Oscillation, Annular modes, etc. The latter are named after the “annulus” (a region bounded by two concentric circles) as they occur in the Earth’s midlatitudes (a band of atmosphere bounded by the polar and tropical regions, Fig. 1), and are the most important modes of midlatitude variability, generally representing 20-30% of the variability in a field.

Southern_Hemi_Antarctica
Figure 1: Southern Hemisphere midlatitudes (red concentric circles) as annulus, region where annular modes have the largest impacts. Source.

We know two types of annular modes: baroclinic (based on eddy kinetic energy, a proxy for eddy activity and an indicator of storm-track intensity) and barotropic (based on zonal mean zonal wind, representing the north-south shifts of the jet stream) (Fig. 2). The latter are usually referred to as Southern (SAM or Antarctic Oscillation) or Northern (NAM or Arctic Oscillation) Annular Mode (depending on the hemisphere), have generally quasi-barotropic (uniform) vertical structure, and impact the temperature variations, sea-ice distribution, and storm paths in both hemispheres with timescales of about 10 days. The former are referred to as BAM (baroclinic annular mode) and exhibit strong vertical structure associated with strong vertical wind shear (baroclinicity), and their impacts are yet to be determined (e.g. Thompson and Barnes 2014, Marshall et al. 2017). These two modes of variability are linked to the key processes of the midlatitude tropospheric dynamics that are involved in the growth (baroclinic processes) and decay (barotropic processes) of midlatitude storms. The growth stage of the midlatitude storms is conventionally associated with increase in eddy kinetic energy (EKE) and the decay stage with decrease in EKE.

ThompsonWoodworth_Fig2a_SAM_2f_BAM(1)
Figure 2: Barotropic annular mode (right), based on zonal wind (contours), associated with eddy momentum flux (shading); Baroclinic annular mode (left), based on eddy kinetic energy (contours), associated with eddy heat flux (shading). Source: Thompson and Woodworth (2014).

However, recent observational studies (e.g. Thompson and Woodworth 2014) have suggested decoupling of baroclinic and barotropic components of atmospheric variability in the Southern Hemisphere (i.e. no correlation between the BAM and SAM) and a simpler formulation of the EKE budget that only depends on eddy heat fluxes and BAM (Thompson et al. 2017). Using cross-spectrum analysis, we empirically test the validity of the suggested relationship between EKE and heat flux at different timescales (Boljka et al. 2018). Two different relationships are identified in Fig. 3: 1) a regime where EKE and eddy heat flux relationship holds well (periods longer than 10 days; intermediate timescale); and 2) a regime where this relationship breaks down (periods shorter than 10 days; synoptic timescale). For the relationship to hold (by construction), the imaginary part of the cross-spectrum must follow the angular frequency line and the real part must be constant. This is only true at the intermediate timescales. Hence, the suggested decoupling of baroclinic and barotropic components found in Thompson and Woodworth (2014) only works at intermediate timescales. This is consistent with our theoretical model (Boljka and Shepherd 2018), which predicts decoupling under synoptic temporal and spatial averaging. At synoptic timescales, processes such as barotropic momentum fluxes (closely related to the latitudinal shifts in the jet stream) contribute to the variability in EKE. This is consistent with the dynamics of storms that occur on timescales shorter than 10 days (e.g. Simmons and Hoskins 1978). This is further discussed in Boljka et al. (2018).

EKE_hflux_cross_spectrum_blog
Figure 3: Imaginary (black solid line) and Real (grey solid line) parts of cross-spectrum between EKE and eddy heat flux. Black dashed line shows the angular frequency (if the tested relationship holds, the imaginary part of cross-spectrum follows this line), the red line distinguishes between the two frequency regimes discussed in text. Source: Boljka et al. (2018).

References

Boljka, L., and T. G. Shepherd, 2018: A multiscale asymptotic theory of extratropical wave, mean-flow interaction. J. Atmos. Sci., in press.

Boljka, L., T. G. Shepherd, and M. Blackburn, 2018: On the coupling between barotropic and baroclinic modes of extratropical atmospheric variability. J. Atmos. Sci., in review.

Marshall, G. J., D. W. J. Thompson, and M. R. van den Broeke, 2017: The signature of Southern Hemisphere atmospheric circulation patterns in Antarctic precipitation. Geophys. Res. Lett., 44, 11,580–11,589.

Simmons, A. J., and B. J. Hoskins, 1978: The life cycles of some nonlinear baroclinic waves. J. Atmos. Sci., 35, 414–432.

Thompson, D. W. J., and E. A. Barnes, 2014: Periodic variability in the large-scale Southern Hemisphere atmospheric circulation. Science, 343, 641–645.

Thompson, D. W. J., B. R. Crow, and E. A. Barnes, 2017: Intraseasonal periodicity in the Southern Hemisphere circulation on regional spatial scales. J. Atmos. Sci., 74, 865–877.

Thompson, D. W. J., and J. D. Woodworth, 2014: Barotropic and baroclinic annular variability in the Southern Hemisphere. J. Atmos. Sci., 71, 1480–1493.

Trouble in paradise: Climate change, extreme weather and wildlife conservation on a tropical island.

Joseph Taylor, NERC SCEARNIO DTP student. Zoological Society of London.

Email: J.Taylor5@pgr.reading.ac.uk

Projecting the impacts of climate change on biodiversity is important for informing

Mauritius Kestrel by Joe Taylor
Male Mauritius kestrel (Falco punctatus) in the Bambous Mountains, eastern Mauritius. Photo by Joe Taylor.

mitigation and adaptation strategies. There are many studies that project climate change impacts on biodiversity; however, changes in the occurrence of extreme weather events are often omitted, usually because of insufficient understanding of their ecological impacts. Yet, changes in the frequency and intensity of extreme weather events may pose a greater threat to ecosystems than changes in average weather regimes (Jentsch and Beierkuhnlein 2008). Island species are expected to be particularly vulnerable to climate change pressures, owing to their inherently limited distribution, population size and genetic diversity, and because of existing impacts from human activities, including habitat destruction and the introduction of non-native species (e.g. Fordham and Brook 2010).

Mauritius is an icon both of species extinction and the successful recovery of threatened species. However, the achievements made through dedicated conservation work and the investment of substantial resources may be jeopardised by future climate change. Conservation programmes in Mauritius have involved the collection of extensive data on individual animals, creating detailed longitudinal datasets. These provide the opportunity to conduct in-depth analyses into the factors that drive population trends.

My study focuses on the demographic impacts of weather conditions, including extreme events, on three globally threatened bird species that are endemic to Mauritius. I extended previous research into weather impacts on the Mauritius kestrel (Falco punctatus), and applied similar methods to the echo parakeet (Psittacula eques) and Mauritius fody (Foudia rubra). The kestrel and parakeet were both nearly lost entirely in the 1970s and 1980s respectively, having suffered severe population bottlenecks, but all three species have benefitted from successful recovery programmes. I analysed breeding success using generalised linear mixed models and analysed survival probability using capture-mark-recapture models. Established weather indices were adapted for use in this study, including indices to quantify extreme rainfall, droughts and tropical cyclone activity. Trends in weather indices at key conservation sites were also analysed.

The results for the Mauritius kestrel add to a body of evidence showing that precipitation is an important limiting factor in its demography and population dynamics. The focal population in the Bambous Mountains of eastern Mauritius occupies an area in which rainfall is increasing. This trend could have implications for the population, as my analyses provide evidence that heavy rainfall during the brood phase of nests reduces breeding success, and that prolonged spells of rain in the cyclone season negatively impact the survival of juveniles. This probably occurs through reductions in hunting efficiency, time available for hunting and prey availability, so that kestrels are unable to capture enough prey to sustain themselves and feed their young (Nicoll et al. 2003, Senapathi et al. 2011). Exposure to heavy and prolonged rainfall could also be a direct cause of mortality through hypothermia, especially for chicks if nests are flooded (Senapathi et al. 2011). Future management of this species may need to incorporate strategies to mitigate the impacts of increasing rainfall.

References:

Fordham, D. A. and Brook, B. W. (2010) Why tropical island endemics are acutely susceptible to global change. Biodiversity and Conservation 19(2): 329‒342.

Jentsch, A. and Beierkuhnlein, C. (2008) Research frontiers in climate change: Effects of extreme meteorological events on ecosystems. Comptes Rendus Geoscience 340: 621‒628.

Nicoll, M. A. C., Jones, C. G. and Norris, K. (2003) Declining survival rates in a reintroduced population of the Mauritius kestrel: evidence for non-linear density dependence and environmental stochasticity. Journal of Animal Ecology 72: 917‒926.

Senapathi, D., Nicoll, M. A. C., Teplitsky, C., Jones, C. G. and Norris, K. (2011) Climate change and the risks associated with delayed breeding in a tropical wild bird population. Proceedings of the Royal Society B 278: 3184‒3190.

Climate model systematic biases in the Maritime Continent

Email: y.y.toh@pgr.reading.ac.uk

The Maritime Continent commonly refers to the groups of islands of Indonesia, Borneo, New Guinea and the surrounding seas in the literature. My study area covers the Maritime Continent domain from 20°S to 20°N and 80°E to 160°E as shown in Figure 1. This includes Indonesia, Malaysia, Brunei, Singapore, Philippines, Papua New Guinea, Solomon islands, northern Australia and parts of mainland Southeast Asia including Thailand, Laos, Cambodia, Vietnam and Myanmar.

subsetF1
Figure 1: JJA precipitation (mm/day) and 850 hPa wind (m s−1) for (a) GPCP and ERA-interim, (b) MMM biases and (c)–(j) AMIP biases for 1979–2008 over the Maritime Continent region (20°S–20ºN, 80°E–160ºE). Third panel shows the Maritime Continent domain and land-sea mask

The ability of climate model to simulate the mean climate and climate variability over the Maritime Continent remains a modelling challenge (Jourdain et al. 2013). Our study examines the fidelity of Coupled Model Intercomparison Project phase 5 (CMIP5) models at simulating mean climate over the Maritime Continent. We find that there is a considerable spread in the performance of the Atmospheric Model Intercomparison Project (AMIP) models in reproducing the seasonal mean climate and annual cycle over the Maritime Continent region. The multi-model mean (MMM) (Figure 1b) JJA precipitation and 850hPa wind biases with respect to observations (Figure 1a) are small compared to individual model biases (Figure 1c-j) over the Maritime Continent. Figure 1 shows only a subset of Fig. 2 from Toh et al. (2017), for the full figure and paper please click here.

We also investigate the model characteristics that may be potential sources of bias. We find that AMIP model performance is largely unrelated to model horizontal resolution. Instead, a model’s local Maritime Continent biases are somewhat related to its biases in the local Hadley circulation and global monsoon.

cluster2
Figure 2: Latitude-time plot of precipitation zonally averaged between 80°E and 160°E for (a) GPCP, (b) Cluster I and (c) Cluster II. White dashed line shows the position of the maximum precipitation each month. Precipitation biases with respect to GPCP for (d) Cluster I and (e) Cluster II.

To characterize model systematic biases in the AMIP runs and determine if these biases are related to common factors elsewhere in the tropics, we performed cluster analysis on Maritime Continent annual cycle precipitation. Our analysis resulted in two distinct clusters. Cluster I (Figure 2b,d) is able to reproduce the observed seasonal migration of Maritime Continent precipitation, but it overestimates the precipitation, especially during the JJA and SON seasons. Cluster II (Figure 2c,e) simulate weaker seasonal migration of Intertropical Convergence Zone (ITCZ) than observed, and the maximum rainfall position stays closer to the equator throughout the year. Tropics-wide properties of clusters also demonstrate a connection between errors at regional scale of the Maritime Continent and errors at large scale circulation and global monsoon.

On the other hand, comparison with coupled models showed that air-sea coupling yielded complex impacts on Maritime Continent precipitation biases. One of the outstanding problems in the coupled CMIP5 models is the sea surface temperature (SST) biases in tropical ocean basins. Our study highlighted central Pacific and western Indian Oceans as the key regions which exhibit the most surface temperature correlation with Maritime Continent mean state precipitation in the coupled CMIP5 models. Future work will investigate the impact of SST perturbations in these two regions on Maritime Continent precipitation using Atmospheric General Circulation Model (AGCM) sensitivity experiments.

 

 

References:

Jourdain N.C., Gupta A.S., Taschetto A.S., Ummenhofer C.C., Moise A.F., Ashok K. (2013) The Indo-Australian monsoon and its relationship to ENSO and IOD in reanalysis data and the CMIP3/CMIP5 simulations. Climate Dynamics. 41(11–12):3073–3102

Toh, Y.Y., Turner, A.G., Johnson, S.J., & Holloway, C.E. (2017). Maritime Continent seasonal climate biases in AMIP experiments of the CMIP5 multimodel ensemble. Climate Dynamics. doi: 10.1007/s00382-017-3641-x

Should we be ‘Leaf’-ing out vegetation when parameterising the aerodynamic properties of urban areas?

Email: C.W.Kent@pgr.reading.ac.uk

When modelling urban areas, vegetation is often ignored in attempt to simplify an already complex problem. However, vegetation is present in all urban environments and it is not going anywhere… For reasons ranging from sustainability to improvements in human well-being, green spaces are increasingly becoming part of urban planning agendas. Incorporating vegetation is therefore a key part of modelling urban climates. Vegetation provides numerous (dis)services in the urban environment, each of which requires individual attention (Salmond et al. 2016). However, one of my research interests is how vegetation influences the aerodynamic properties of urban areas.

Two aerodynamic parameters can be used to represent the aerodynamic properties of a surface: the zero-plane displacement (zd) and aerodynamic roughness length (z0). The zero-plane displacement is the vertical displacement of the wind-speed profile due to the presence of surface roughness elements. The aerodynamic roughness length is a length scale which describes the magnitude of surface roughness. Together they help define the shape and form of the wind-speed profile which is expected above a surface (Fig. 1).

blogpostpic

Figure 1: Representation of the wind-speed profile above a group of roughness elements. The black dots represent an idealised logarithmic wind-speed profile which is determined using the zero-plane displacement (zd) and aerodynamic roughness length (z0) (lines) of the surface.

For an urban site, zd and z0 may be determined using three categories of methods: reference-based, morphometric and anemometric. Reference-based methods require a comparison of the site to previously published pictures or look up tables (e.g. Grimmond and Oke 1999); morphometric methods describe zd and z0 as a function of roughness-element geometry; and, anemometric methods use in-situ observations. The aerodynamic parameters of a site may vary considerably depending upon which of these methods are used, but efforts are being made to understand which parameters are most appropriate to use for accurate wind-speed estimations (Kent et al. 2017a).

Within the morphometric category (i.e. using roughness-element geometry) sophisticated methods have been developed for buildings or vegetation only. However, until recently no method existed to describe the effects of both buildings and vegetation in combination. A recent development overcomes this, whereby the heights of all roughness elements are considered alongside a porosity correction for vegetation (Kent et al. 2017b). Specifically, the porosity correction is applied to the space occupied and drag exerted by vegetation.

The development is assessed across several areas typical of a European city, ranging from a densely-built city centre to an urban park. The results demonstrate that where buildings are the dominant roughness elements (i.e. taller and occupying more space), vegetation does not obviously influence the calculated geometry of the surface, nor the aerodynamic parameters and the estimated wind speed. However, as vegetation begins to occupy a greater amount of space and becomes as tall as (or larger) than buildings, the influence of vegetation is obvious. Expectedly, the implications are greatest in an urban park, where overlooking vegetation means that wind speeds may be slowed by up to a factor of three.

Up to now, experiments such as those in the wind tunnel focus upon buildings or trees in isolation. Certainly, future experiments which consider both buildings and vegetation will be valuable to continue to understand the interaction within and between these roughness elements, in addition to assessing the parameterisation.

References

Grimmond CSB, Oke TR (1999) Aerodynamic properties of urban areas derived from analysis of surface form. J Appl Meteorol and Clim 38:1262-1292.

Kent CW, Grimmond CSB, Barlow J, Gatey D, Kotthaus S, Lindberg F, Halios CH (2017a) Evaluation of Urban Local-Scale Aerodynamic Parameters: Implications for the Vertical Profile of Wind Speed and for Source Areas. Boundary-Layer Meteorology 164: 183-213.

Kent CW, Grimmond CSB, Gatey D (2017b) Aerodynamic roughness parameters in cities: Inclusion of vegetation. Journal of Wind Engineering and Industrial Aerodynamics 169: 168-176.

Salmond JA, Tadaki M, Vardoulakis S, Arbuthnott K, Coutts A, Demuzere M, Dirks KN, Heaviside C, Lim S, Macintyre H (2016) Health and climate related ecosystem services provided by street trees in the urban environment. Environ Health 15:95.