Coupled atmosphere–ocean data assimilation experiments with a low-order model and CMIP5 model data

Climate Dynamics - Coupled atmosphere–ocean data assimilation (DA) experiments are performed for estimating the Atlantic meridional overturning circulation (AMOC). Recovery of the AMOC with...     

The UHF wind profiler at Vienna airport – data quality control and comparisons to rawinsonde data

Summary A 1290MHz wind profiler (Radian Lap-3000), at present one of three operational wind profilers in Austria, is operated at Vienna airport. In spite of quality assurance procedures as consensus averaging included in the data evaluation process from profiler raw data, some spurious peaks of wind speed and unrealistic changes of the wind vector in time or height occur in the wind measurements. This is especially true for sampling intervals of only 5 minutes which are used to resolve the temporal evolution of summer thunderstorms and frontal passages. Averaging periods of only a few minutes are rather the lower limit apt for wind profiler observations and result in a low data availability of 28%, whereas about 55% of data (relative to the maximum height range according to the parameter setting) are available for 10 to 30 minutes profiles.Approaches to a posteriori quality control using checks for automatic error detection are proposed and tested on a one and a half year data-set: Flagging data when the three-dimensional wind divergence exceeds a predefined limit (0.5s^–1) is in most cases successful in combination with thresholds for wind speed (2 times the median of the daily data-set) or wind shear (0.2s^–1).The wind profiler data is compared to wind profiles from the next radiosonde station where soundings are launched 4 times a day at Hohe Warte, approx. 20km northwest, at the hill-side of the Viennese Woods. Deviations of about 1m s^–1 in wind speed are found between the observations of the two systems. Differences between the wind profiles within the boundary layer can be explained by local differences in the wind regime observed at the airport and the radiosounding – blocking effects of the Viennese Woods during south-easterly flow. Comparing the profiler data to radiosoundings on a monthly basis gives a tool to monitor the profiler performance.     

The KNMI HEXMAX stress data – a reanalysis

An earlier discussion of the wind stress measurements of the 1986 HEXMAX experiment off the Dutch coast left some uncertainty about the actual relationship between the stress and the wave field. In this paper we try to find a more definitive answer by looking at the consequences for gradients of wind speed, stress and phase speed of the waves if we assume the Charnock coefficient to be either a constant or to have an inverse wave age dependence. It turns out that both assumptions have a wind speed range where they produce coherent results, but that they break down when the wave length exceeds a value that is related to the water depth.     

A data assimilation technique based on the π-algorithm

A practical implementation of the data assimilation algorithm based on the Kalman filter in its complete formulation is impossible due to high dimension of the associated equation sets and to nonlinearity of the predicted processes. The main direction in the implementation of the Kalman filter is an ensemble approach. Under the assumption of ergodicity of random forecast errors, an alternative algorithm with respect to the ensemble Kalman filter can be considered, in which probability averaging is replaced by time averaging. The proposes algorithm is based this assumption. The algorithm is easy to implement; however, its convergence, applicability to the data assimilation problems, and connection to the Kalman filter have not been studied. In the paper, applicability of the π-algorithm to data assimilation is considered on an example of a simple one-dimensional advection equation. Use of this simple equation allows comparing the classical Kalman filter algorithm with various practical approaches to its implementation.     

An observational‐data based evapotranspiration function for general circulation models

Abstract

In this paper Nappo’ s (1975) formulation of the moisture availability function was used to derive the ß function of the form employed in the evapotranspiration calculations of various GCMs (Carson, 1981). An inverse calculation using the planetary boundary‐layer parameterizations of the GLAS general circulation model was made to derive this function. For this purpose two ground temperatures, namely those of saturated and naturally dry ground, were prognostically carried in a 47‐day integration with the model. The form of the calculated ß function is different from those reviewed by Carson. An example of global evapotranspiration obtained with the derived ß function is shown. Two separate five‐day simulations, one using ß as derived here and the other using Nappo's (1975) M function as a ß function, are compared. Large differences in the calculated evapotranspiration occur in dry regions.     

Coupled atmosphere–ocean data assimilation experiments with a low-order climate model

A simple idealized atmosphere–ocean climate model and an ensemble Kalman filter are used to explore different coupled ensemble data assimilation strategies. The model is a low-dimensional analogue of the North Atlantic climate system, involving interactions between large-scale atmospheric circulation and ocean states driven by the variability of the Atlantic meridional overturning circulation (MOC). Initialization of the MOC is assessed in a range of experiments, from the simplest configuration consisting of forcing the ocean with a known atmosphere to performing fully coupled ensemble data assimilation. “Daily” assimilation (that is, at the temporal frequency of the atmospheric observations) is contrasted with less frequent assimilation of time-averaged observations. Performance is also evaluated under scenarios in which ocean observations are limited to the upper ocean or are non-existent. Results show that forcing the idealized ocean model with atmospheric analyses is inefficient at recovering the slowly evolving MOC. On the other hand, daily assimilation rapidly leads to accurate MOC analyses, provided a comprehensive set of oceanic observations is available for assimilation. In the absence of sufficient observations in the ocean, the assimilation of time-averaged atmospheric observations proves to be more effective for MOC initialization, including the case where only atmospheric observations are available.     

Ice‐floe collisions interpreted from acceleration data during LIMEX ‘89

Abstract

Wave‐induced ice motions measured during the Labrador Ice Margin Experiment (LIMEX ‘89) were interpreted to determine the cause and the frequency of collisions between floes. The LIMEX acceleration data were acquired with an optimal resolution near the predominant wave frequency and did not contain information above 0.5 Hz. It was therefore possible to establish the frequency of collisions, but not the magnitude of the events. Events were defined by any contact between floes in a wave cycle, and the distribution of times between events indicates that floes are more likely to collide in adjacent wave cycles than if the events were independent. Periods of continuous and intermittent collisions were related to the wave characteristics, and the frequency of events increased with a decrease in air temperature and an increase in local wind speed. Contrary to expectations, there was not a positive relation between collision frequency and wave amplitude.     

The global climatology of an interannually varying air–sea flux data set

The air–sea fluxes of momentum, heat, freshwater and their components have been computed globally from 1948 at frequencies ranging from 6-hourly to monthly. All fluxes are computed over the 23 years from 1984 to 2006, but radiation prior to 1984 and precipitation before 1979 are given only as climatological mean annual cycles. The input data are based on NCEP reanalysis only for the near surface vector wind, temperature, specific humidity and density, and on a variety of satellite based radiation, sea surface temperature, sea-ice concentration and precipitation products. Some of these data are adjusted to agree in the mean with a variety of more reliable satellite and in situ measurements, that themselves are either too short a duration, or too regional in coverage. The major adjustments are a general increase in wind speed, decrease in humidity and reduction in tropical solar radiation. The climatological global mean air–sea heat and freshwater fluxes (1984–2006) then become 2 W/m² and ?0.1 mg/m² per second, respectively, down from 30 W/m² and 3.4 mg/m² per second for the unaltered data. However, decadal means vary from 7.3 W/m² (1977–1986) to ?0.3 W/m² (1997–2006). The spatial distributions of climatological fluxes display all the expected features. A comparison of zonally averaged wind stress components across ocean sub-basins reveals large differences between available products due both to winds and to the stress calculation. Regional comparisons of the heat and freshwater fluxes reveal an alarming range among alternatives; typically 40 W/m² and 10 mg/m² per second, respectively. The implied ocean heat transports are within the uncertainty of estimates from ocean observations in both the Atlantic and Indo-Pacific basins. They show about 2.4 PW of tropical heating, of which 80% is transported to the north, mostly in the Atlantic. There is similar good agreement in freshwater transport at many latitudes in both basins, but neither in the South Atlantic, nor at 35°N.     

Cloudiness over the Black Sea region in 1985–2009 from satellite data

The seasonal and interannual variability ofcloud fraction over the Black Sea region for the period of1985-2009 is analyzed using the CM SAF dataset obtained from the satellite measurements of a high-resolution AVHRR instrument. The features of geographic distribution and seasonal variations in cloudiness are investigated. The causes for its spatial inhomogeneity in different months are analyzed. It is demonstrated using the long-term dataset that the dramatic decrease in the amount of cloudiness occurred over the Black Sea region from 67% in 1985 to 54% in 2008. The value of the trend is -0.4% per year. Both the trends and the features of interannual variability of cloudiness, in particular, strongly pronounced four-year cycles, are in antiphase with variations in sea surface temperature. The cloudiness reduction accompanied by the increase in the influx of short-wave radiation may be the basic reason for the warming and sea surface temperature variations in the Black Sea region.     

Mediterranean drought fluctuation during the last 500 years based on tree-ring data

A 2.5 × 2.5° gridded summer (April–September) drought reconstruction over the larger Mediterranean land area (32.5°/47.5°N, 10°W/50°E; 152 grid points) is described, based on a network of 165 tree-ring series. The drought index used is the self-calibrated Palmer Drought Severity Index, and the period considered is 1500–2000. The reconstruction technique combines an analogue technique for the estimation of missing tree-ring data with an artificial neural network for optimal non-linear calibration, including a bootstrap error assessment. Tests were carried out on the various sources of error in the reconstructions. Errors related to the temporal variations of the number of proxies were tested by comparing four reconstructions calibrated with four different sized regressor datasets, representing the decrease in the number of available proxies over time. Errors related to the heterogeneous spatial density of predictors were tested using pseudo-proxies, provided by the global climate model ECHO-G. Finally the errors related to the imperfect climate signal recorded by tree-ring series were tested by adding white noise to the pseudo-proxies. Reconstructions pass standard cross-validation tests. Nevertheless tests using pseudo-proxies show that the reconstructions are less good in areas where proxies are rare, but that the average reconstruction curve is robust. Finally, the noise added to proxies, which is by definition a high frequency component, has a major effect on the low frequency signal, but not on the medium frequencies. The comparison of the low frequency trends of our mean reconstruction and the GCM simulation indicates that the detrending method used is able to preserve the long-term variations of reconstructed PDSI. The results also highlight similar multi-decadal PDSI variations in the central and western parts of the Mediterranean basin and less clear low frequency changes in the east. The sixteenth and the first part of the seventeenth centuries are characterized by marked dry episodes in the west similar to those observed in the end of the twentieth century. In contrast, the eighteenth and nineteenth centuries (Little Ice Age) are characterized by dominant wet periods. In the eastern part of the Mediterranean basin the observed strong drought period of the end of the twentieth century seems to be the strongest of the last 500 years.     

Changes in the Brewer-Dobson circulation for 1980–2009 revealed in MERRA reanalysis data

Changes in the Brewer-Dobson circulation (BDC) during the 30 years 1980–2009 are investigated using Modern Era Retrospective-analysis for Research and Applications (MERRA) reanalysis data. The mass streamfunction that is induced by wave forcings in the transformed Eulerian-mean (TEM) equation through the downward-control principle is used as a proxy for the BDC. The changes in the BDC are investigated using two aspects: the wave propagation conditions in the stratosphere and the wave activity in the upper troposphere. They are compared in the first (P1) and second (P2) 15-year periods. The resolved wave forcing, expressed by the Eliassen-Palm (EP) flux divergence (EPD), is significantly enhanced during the December-January-February (DJF) season in P2 in both the Northern Hemisphere (NH) high latitudes and the Southern Hemisphere (SH) mid- and high latitudes. The increased zonal mean zonal wind at high latitudes in the SH, caused by ozone depletion, leads to an upward shift of the Rossby-wave critical layer and this allows more transient planetary waves to propagate into the stratosphere. In the NH, the enhanced EPD in DJF leads to an increase in the frequency of Sudden Stratospheric Warming (SSW) events. The gravity wave drag (GWD) is smaller than the EPD and the change in it between the two time periods is insignificant. The residual term in the TEM equation is similar to the GWD in the two periods, but its change between the two periods is as large as the change in the EPD. Among the four components of the EP flux at 250 hPa, the meridional heat flux played a dominant role in the enhancement of the BDC in P2.     

MODIS brightness temperature data assimilation under cloudy conditions Ⅱ: Impacts on rainstorm forecasting

Although satellite observations provide large amount of information of clouds and precipitation and play an important role in the forecast of heavy rainfall, they have not been fully taken advantage of in data assimilation of numerical weather predictions, especially those in infrared channels. Assimilating radiances is common under clear-sky conditions since it is extremely difficult to simulate infrared transmittance in cloudy sky. Based on the Global and Regional Assimilation and Prediction Enhanced System 3-dimensional variance (GRAPES-3DVar), cloud liquid water content, ice-water content and cloud cover are employed as governing variables in the assimilation system. This scheme can improve the simulation of infrared transmittance by a fast radiative transfer model for TOVS (RTTOV) and adjust the atmospheric and cloud parameters based on infrared radiance observations. In this paper, we investigate a heavy rainfall over Guangdong province on May 26, 2007, which is right after the onset of a South China Sea monsoon. In this case, channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) for observing water vapor (Channel 27) and cloud top altitude (Channel 36) are selected for the assimilation. The process of heavy rainfall is simulated by the Weather Research and Forecasting (WRF) model. Our results show that the assimilated MODIS data can improve the distribution of water vapor and temperature in the first guess field and indirectly adjust the upper-level wind field. The tendency of adjustment agrees well with the satellite observations. The assimilation scheme has positive impacts on the short-range forecasting of rainstorm.     

Atmospheric transparency over central Yakutia from the sun photometer data in 2004—2014

The transparency of the atmosphere over the central part of Yakutia is classified using ground-based observations of aerosol optical depth (AOD) for the period of 2004-2014. The annual variations in monthly mean values of AOD are compared with satellite monitoring data. It is revealed that every year the days with the values of AOD corresponding to the class III of atmospheric transparency (turbid atmosphere) made up 25-30% of the total number of AOD measurement days.     

AMOC variations in 1979–2008 simulated by NCEP operational ocean data assimilation system

Variations in the Atlantic meridional overturning circulation (AMOC) between 1979 and 2008 are documented using the operational ocean analysis, the Global Ocean Data Assimilation System (GODAS), at the National Centers for Climate Prediction (NCEP). The maximum AMOC at 40°N is about 16?Sv in average with peak-to-peak variability of 3–4?Sv. The AMOC variations are dominated by an upward trend from 1980 to 1995, and a downward trend from 1995 to 2008. The maximum AMOC at 26.5°N is slightly weaker than hydrographic estimates and observations from mooring array. The dominant variability of the AMOC in 20°–65°N (the first EOF, 51% variance) is highly correlated with that in the subsurface temperature (the first EOF, 33% variance), and therefore, with density (the first EOF, 25% variance) in the North Atlantic, and is consistent with the observational estimates based on the World Ocean Database 2005. The dominant variabilities of AMOC and subsurface temperature are also analyzed in the context of possible links with the net surface heat flux, deep convection, western boundary current, and subpolar gyre. Variation in the net surface heat flux is further linked to the North Atlantic Oscillation (NAO) index which is found to lead AMOC variations by about 5?years. Our results indicate that AMOC variations can be documented based on an ocean analysis system such as GODAS.     

Analysis of rainfall characteristics of the Madden–Julian oscillation using TRMM satellite data

Rainfall characteristics of the Madden–Julian oscillation (MJO) are analyzed primarily using tropical rainfall measuring mission (TRMM) precipitation radar (PR), TRMM microwave imager (TMI) and lighting imaging sensor (LIS) data. Latent heating structure is also examined using latent heating data estimated with the spectral latent heating (SLH) algorithm.The zonal structure, time evolution, and characteristic stages of the MJO precipitation system are described. Stratiform rain fraction increases with the cloud activity, and the amplitude of stratiform rain variation associated with the MJO is larger than that of convective rain by a factor of 1.7. Maximum peaks of both convective rain and stratiform rain precede the minimum peak of the outgoing longwave radiation (OLR) anomaly which is often used as a proxy for the MJO convection. Stratiform rain remains longer than convective rain until ∼4000 km behind the peak of the mature phase. The stratiform rain contribution results in the top-heavy heating profile of the MJO.Associated with the MJO, there are tri-pole convective rain top heights (RTH) at 10–11, ∼7 and ∼3 km, corresponding to the dominance of afternoon showers, organized systems, and shallow convections, respectively. The stratiform rain is basically organized with convective rain, having similar but slightly lower RTH and slightly lags the convective rain maximum. It is notable that relatively moderate (∼7 km) RTH is dominant in the mature phase of the MJO, while very tall rainfall with RTH over 10 km and lightning frequency increase in the suppressed phase. The rain-yield-per flash (RPF) varies about 20–100% of the mean value of ∼2–10 × 10⁹ kg fl⁻¹ over the tropical warm ocean and that of ∼2–5 × 10⁹ kg fl⁻¹ over the equatorial Islands, between the convectively suppressed phase and the active phase of MJO, in the manner that RPF is smaller in the suppressed phase and larger in the active phase.     

Spatial response of two European atmospheric circulation classifications (data 1901–2010)

Air pressure field and circulation pattern frequencies were investigated to (1) locate and compare positions of the underlying pressure fields, (2) analyse the spatial dimension of affected areas, (3) create schematic maps of important circulation types and (4) compare the classification types in their response to the data. Two manual classifications were used, selected for the length of their time series and their applicability to a larger region: the Grosswetterlagen classification (GWLc) and the Vangengeim–Girs classification (VGc). Their time series were correlated with a global set of gridded monthly sea-level pressure data. Results show the different conceptual orientation of VGc (hemispheric) and GWLc (continental). The highest correlation values and the largest affected areas are visible in winter, where patterns frequently extended into northern Africa and western Asia. Schematic maps, illustrating the average location of main pressure centres, are provided for basic classes of both classifications. Re-arranging GWLc subtypes increases the classifications comparability with the VGc. Analysis of moving correlation coefficients reveals high fluctuations in the relation of both classifications over time.     

Time changes of atmospheric precipitation in Russia from the corrected data during 1936–2000

Homogeneous time series of atmospheric precipitation with corrected systematic errors of measurements at 100 stations in Russia for the period of 1936–2000 are obtained. Combined effects are considered of all kinds of systematic errors of standard network precipitation-measuring instruments (the raingauge with the Nifer shield and the Tret’yakov raingauge) on the measured precipitation totals. Comparative analysis is carried out of the measured and corrected long-term mean characteristics of precipitation amounts (annual totals, warm and cold season totals, and different types of precipitation). On the basis of the obtained archives of precipitation homogeneous time series, linear trends are estimated for the period under consideration with estimation of their statistical significance. Schematic charts are plotted and analyzed of time changes in the annual precipitation amounts and in the amounts of different types of precipitation.     

The effect of flood experience on household mitigation—Evidence from longitudinal and insurance data

Proactive flood damage mitigation on the household level is seen as a crucial element of comprehensive flood risk management. The ongoing socio-economic development and climate change will further increase the relevance of flood risks in the future. This paper analyses the causal effect of flood experience on mitigation decisions of households. It contributes to the emerging literature on the determinants of private flood mitigation and adopts data which has rarely been used in the field: Instead of cross-sectional data, longitudinal data of household surveys before and after a major flood event in Germany is used (N > 7,400). Flood experience is elicited by three different data sources: reports of heads of households, issuance of flood alerts, and damage reports of the insurance industry. The results show a causal effect of insured flood damage on private flood mitigation and a correlation of mitigation with self-reported flood experience. The estimated effects differ for mitigation types (higher for behavioural options, such as “adapted use”) and household types (higher for better educated heads of households).