Regionally coupled atmosphere-ocean-sea ice-marine biogeochemistry model ROM: 1. Description and validation

The general circulation models used to simulate global climate typically feature resolution too coarse to reproduce many smaller-scale processes, which are crucial to determining the regional responses to climate change. A novel approach to downscale climate change scenarios is presented which includes the interactions between the North Atlantic Ocean and the European shelves as well as their impact on the North Atlantic and European climate. The goal of this paper is to introduce the global ocean-regional atmosphere coupling concept and to show the potential benefits of this model system to simulate present-day climate. A global ocean-sea ice-marine biogeochemistry model (MPIOM/HAMOCC) with regionally high horizontal resolution is coupled to an atmospheric regional model (REMO) and global terrestrial hydrology model (HD) via the OASIS coupler. Moreover, results obtained with ROM using NCEP/NCAR reanalysis and ECHAM5/MPIOM CMIP3 historical simulations as boundary conditions are presented and discussed for the North Atlantic and North European region. The validation of all the model components, i.e., ocean, atmosphere, terrestrial hydrology, and ocean biogeochemistry is performed and discussed. The careful and detailed validation of ROM provides evidence that the proposed model system improves the simulation of many aspects of the regional climate, remarkably the ocean, even though some biases persist in other model components, thus leaving potential for future improvement. We conclude that ROM is a powerful tool to estimate possible impacts of climate change on the regional scale.

Precipitation and floodiness

There are a number of factors that lead to non-linearity between precipitation anomalies and flood hazard; this non-linearity is a pertinent issue for applications that use a precipitation forecast as a proxy for imminent flood hazard. We assessed the degree of this non-linearity for the first time using a recently developed global-scale hydrological model driven by the ERA-Interim Land precipitation reanalysis (1980–2010). We introduced new indices to assess large-scale flood hazard, or floodiness, and quantified the link between monthly precipitation, river discharge and floodiness anomalies at the global and regional scales. The results show that monthly floodiness is not well correlated with precipitation, therefore demonstrating the value of hydrometeorological systems for providing floodiness forecasts for decision-makers. A method is described for forecasting floodiness using the Global Flood Awareness System, building a climatology of regional floodiness from which to forecast floodiness anomalies out to two weeks.

Recent observed and simulated changes in precipitation over Africa

Multiple observational data sets and atmosphere-only simulations from the Coupled Model Intercomparison Project Phase 5 are analyzed to characterize recent rainfall variability and trends over Africa focusing on 1983–2010. Data sets exhibiting spurious variability, linked in part to a reduction in rain gauge density, were identified. The remaining observations display coherent increases in annual Sahel rainfall (29 to 43 mm yr−1 per decade), decreases in March–May East African rainfall (−14 to −65 mm yr−1 per decade), and increases in annual Southern Africa rainfall (32 to 41 mm yr−1 per decade). However, Central Africa annual rainfall trends vary in sign (−10 to +39 mm yr−1 per decade). For Southern Africa, observed and sea surface temperature (SST)-forced model simulated rainfall variability are significantly correlated (r~0.5) and linked to SST patterns associated with recent strengthening of the Pacific Walker circulation.

Robust comparison of climate models with observations using blended land air and ocean sea surface temperatures

The level of agreement between climate model simulations and observed surface temperature change is a topic of scientific and policy concern. While the Earth system continues to accumulate energy due to anthropogenic and other radiative forcings, estimates of recent surface temperature evolution fall at the lower end of climate model projections. Global mean temperatures from climate model simulations are typically calculated using surface air temperatures, while the corresponding observations are based on a blend of air and sea surface temperatures. This work quantifies a systematic bias in model-observation comparisons arising from differential warming rates between sea surface temperatures and surface air temperatures over oceans. A further bias arises from the treatment of temperatures in regions where the sea ice boundary has changed. Applying the methodology of the HadCRUT4 record to climate model temperature fields accounts for 38% of the discrepancy in trend between models and observations over the period 1975–2014.

SPARE-ICE: synergistic Ice Water Path from passive operational sensors

This article presents SPARE-ICE, the Synergistic Passive Atmospheric Retrieval Experiment-ICE. SPARE-ICE is the first Ice Water Path (IWP) product combining infrared and microwave radiances. By using only passive operational sensors, the SPARE-ICE retrieval can be used to process data from at least the NOAA 15 to 19 and MetOp satellites, obtaining time series from 1998 onward. The retrieval is developed using collocations between passive operational sensors (solar, terrestrial infrared, microwave), the CloudSat radar, and the CALIPSO lidar. The collocations form a retrieval database matching measurements from passive sensors against the existing active combined radar-lidar product 2C-ICE. With this retrieval database, we train a pair of artificial neural networks to detect clouds and retrieve IWP. When considering solar, terrestrial infrared, and microwave-based measurements, we show that any combination of two techniques performs better than either single-technique retrieval. We choose not to include solar reflectances in SPARE-ICE, because the improvement is small, and so that SPARE-ICE can be retrieved both daytime and nighttime. The median fractional error between SPARE-ICE and 2C-ICE is around a factor 2, a figure similar to the random error between 2C-ICE ice water content (IWC) and in situ measurements. A comparison of SPARE-ICE with Moderate Resolution Imaging Spectroradiometer (MODIS), Pathfinder Atmospheric Extended (PATMOS-X), and Microwave Surface and Precipitation Products System (MSPPS) indicates that SPARE-ICE appears to perform well even in difficult conditions. SPARE-ICE is available for public use.

Systematic and random errors between collocated satellite ice water path observations

There remains large disagreement between ice-water path (IWP) in observational data sets, largely because the sensors observe different parts of the ice particle size distribution. A detailed comparison of retrieved IWP from satellite observations in the Tropics (!30 " latitude) in 2007 was made using collocated measurements. The radio detection and ranging(radar)/light detection and ranging (lidar) (DARDAR) IWP data set, based on combined radar/lidar measurements, is used as a reference because it provides arguably the best estimate of the total column IWP. For each data set, usable IWP dynamic ranges are inferred from this comparison. IWP retrievals based on solar reflectance measurements, in the moderate resolution imaging spectroradiometer (MODIS), advanced very high resolution radiometer–based Climate Monitoring Satellite Applications Facility (CMSAF), and Pathfinder Atmospheres-Extended (PATMOS-x) datasets, were found to be correlated with DARDAR over a large IWP range (~20–7000 g m -2 ). The random errors of the collocated data sets have a close to lognormal distribution, and the combined random error of MODIS and DARDAR is less than a factor of 2, which also sets the upper limit for MODIS alone. In the same way, the upper limit for the random error of all considered data sets is determined. Data sets based on passive microwave measurements, microwave surface and precipitation products system (MSPPS), microwave integrated retrieval system (MiRS), and collocated microwave only (CMO), are largely correlated with DARDAR for IWP values larger than approximately 700 g m -2 . The combined uncertainty between these data sets and DARDAR in this range is slightly less MODIS-DARDAR, but the systematic bias is nearly an order of magnitude.

Monitoring scan asymmetry of microwave humidity sounding channels using simultaneous all angle collocations (SAACs)

Simultaneous all angle collocations (SAACs) of microwave humidity sounders (AMSU-B and MHS) on-board polar orbiting satellites are used to estimate scan-dependent biases. This method has distinct advantages over previous methods, such as that the estimated scan-dependent biases are not influenced by diurnal differences between the edges of the scan and the biases can be estimated for both sides of the scan. We find the results are robust in the sense that biases estimated for one satellite pair can be reproduced by double differencing biases of these satellites with a third satellite. Channel 1 of these instruments shows the least bias for all satellites. Channel 2 has biases greater than 5 K, thus needs to be corrected. Channel 3 has biases of about 2 K and more and they are time varying for some of the satellites. Channel 4 has the largest bias which is about 15 K when the data are averaged for 5 years, but biases of individual months can be as large as 30 K. Channel 5 also has large and time varying biases for two of the AMSU-Bs. NOAA-15 (N15) channels are found to be affected the most, mainly due to radio frequency interference (RFI) from onboard data transmitters. Channel 4 of N15 shows the largest and time varying biases, so data of this channel should only be used with caution for climate applications. The two MHS instruments show the best agreement for all channels. Our estimates may be used to correct for scan-dependent biases of these instruments, or at least used as a guideline for excluding channels with large scan asymmetries from scientific analyses.

Understanding inter-satellite biases of microwave humidity sounders using global simultaneous nadir overpasses

Simultaneous nadir overpasses (SNOs) of polar-orbiting satellites are most frequent in polar areas but can occur at any latitude when the equatorial crossing times of the satellites become close owing to orbital drift. We use global SNOs of polar orbiting satellites to evaluate the intercalibration of microwave humidity sounders from the more frequent high-latitude SNOs. We have found based on sensitivity analyses that optimal distance and time thresholds for defining collocations are pixel centers less than 5 km apart and time differences less than 300 s. These stringent collocation criteria reduce the impact of highly variable surface or atmospheric conditions on the estimated biases. Uncertainties in the estimated biases are dominated by the combined radiometric noise of the instrument pair. The effects of frequency changes between different versions of the humidity sounders depend on the amount of water vapor in the atmosphere. There are significant scene radiance and thus latitude dependencies in the estimated biases and this has to taken into account while intercalibrating microwave humidity sounders. Therefore the results obtained using polar SNOs will not be representative for moist regions, necessitating the use of global collocations for reliable intercalibration.

Revisiting trends in wetness and dryness in the presence of internal climate variability and water limitations over land

A theoretically expected consequence of the intensification of the hydrological cycle under global warming is that on average, wet regions get wetter and dry regions get drier (WWDD). Recent studies, however, have found significant discrepancies between the expected pattern of change and observed changes over land. We assess the WWDD theory in four climate models. We find that the reported discrepancy can be traced to two main issues: (1) unforced internal climate variability strongly affects local wetness and dryness trends and can obscure underlying agreement with WWDD, and (2) dry land regions are not constrained to become drier by enhanced moisture divergence since evaporation cannot exceed precipitation over multiannual time scales. Over land, where the available water does not limit evaporation, a “wet gets wetter” signal predominates. On seasonal time scales, where evaporation can exceed precipitation, trends in wet season becoming wetter and dry season becoming drier are also found.

Lightning as a space-weather hazard: UK thunderstorm activity modulated by the passage of the heliospheric current sheet

Lightning flash rates, RL, are modulated by corotating interaction regions (CIRs) and the polarity of the heliospheric magnetic field (HMF) in near-Earth space. As the HMF polarity reverses at the heliospheric current sheet (HCS), typically within a CIR, these phenomena are likely related. In this study, RL is found to be significantly enhanced at the HCS and at 27 days prior/after. The strength of the enhancement depends on the polarity of the HMF reversal at the HCS. Near-Earth solar and galactic energetic particle fluxes are also ordered by HMF polarity, though the variations qualitatively differ from RL, with the main increase occurring prior to the HCS crossing. Thus, the CIR effect on lightning is either the result of compression/amplification of the HMF (and its subsequent interaction with the terrestrial system) or that energetic particle preconditioning of the Earth system prior to the HMF polarity change is central to solar wind lightning coupling mechanism.

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.

The sensitivity of convective aggregation to diabatic processes in idealized radiative-convective equilibrium simulations

Idealized explicit convection simulations of the Met Office Unified Model exhibit spontaneous self-aggregation in radiative-convective equilibrium, as seen in other models in previous studies. This self-aggregation is linked to feedbacks between radiation, surface fluxes, and convection, and the organization is intimately related to the evolution of the column water vapor field. Analysis of the budget of the spatial variance of column-integrated frozen moist static energy (MSE), following Wing and Emanuel [2014], reveals that the direct radiative feedback (including significant cloud longwave effects) is dominant in both the initial development of self-aggregation and the maintenance of an aggregated state. A low-level circulation at intermediate stages of aggregation does appear to transport MSE from drier to moister regions, but this circulation is mostly balanced by other advective effects of opposite sign and is forced by horizontal anomalies of convective heating (not radiation). Sensitivity studies with either fixed prescribed radiative cooling, fixed prescribed surface fluxes, or both do not show full self-aggregation from homogeneous initial conditions, though fixed surface fluxes do not disaggregate an initialized aggregated state. A sensitivity study in which rain evaporation is turned off shows more rapid self-aggregation, while a run with this change plus fixed radiative cooling still shows strong self-aggregation, supporting a “moisture memory” effect found in Muller and Bony [2015]. Interestingly, self-aggregation occurs even in simulations with sea surface temperatures (SSTs) of 295 K and 290 K, with direct radiative feedbacks dominating the budget of MSE variance, in contrast to results in some previous studies.

Sahel decadal rainfall variability and the role of model horizontal resolution

Substantial low-frequency rainfall fluctuations occurred in the Sahel throughout the twentieth century, causing devastating drought. Modeling these low-frequency rainfall fluctuations has remained problematic for climate models for many years. Here we show using a combination of state-of-the-art rainfall observations and high-resolution global climate models that changes in organized heavy rainfall events carry most of the rainfall variability in the Sahel at multiannual to decadal time scales. Ability to produce intense, organized convection allows climate models to correctly simulate the magnitude of late-twentieth century rainfall change, underlining the importance of model resolution. Increasing model resolution allows a better coupling between large-scale circulation changes and regional rainfall processes over the Sahel. These results provide a strong basis for developing more reliable and skilful long-term predictions of rainfall (seasons to years) which could benefit many sectors in the region by allowing early adaptation to impending extremes.

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.

Climatology and variability of trace gases in extratropical double-tropopause regions from MLS, HIRDLS, and ACE-FTS measurements

Upper tropospheric and lower stratospheric measurements from the Aura Microwave Limb Sounder (MLS), the Aura High Resolution Dynamics Limb Sounder (HIRDLS), and the Atmospheric Chemistry Experiment-Fourier transform spectrometer (ACE-FTS) are used to present the first global climatological comparison of extratropical, nonpolar trace gas distributions in double-tropopause (DT) and single-tropopause (ST) regions. Stratospheric tracers, O3, HNO3, and HCl, have lower mixing ratios ∼2–8 km above the primary (lowermost) tropopause in DT than in ST regions in all seasons, with maximum Northern Hemisphere (NH) differences near 50% in winter and 30% in summer. Southern Hemisphere winter differences are somewhat smaller, but summer differences are similar in the two hemispheres. H2O in DT regions of both hemispheres shows strong negative anomalies in November through February and positive anomalies in July through October, reflecting the strong seasonal cycle in H2O near the tropical tropopause. CO and other tropospheric tracers examined have higher DT than ST values 2–7 km above the primary tropopause, with the largest differences in winter. Large DT-ST differences extend to high NH latitudes in fall and winter, with longitudinal maxima in regions associated with enhanced wave activity and subtropical jet variations. Results for O3 and HNO3 agree closely between MLS and HIRDLS, and differences from ACE-FTS are consistent with its sparse and irregular midlatitude sampling. Consistent signatures in climatological trace gas fields provide strong evidence that transport from the tropical upper troposphere into the layer between double tropopauses is an important pathway for stratosphere-troposphere exchange.