ALBERT

Description: In this study, the authors discuss two different parameterizations for the effect of mixed layer eddies, one based on ageostrophic linear stability analysis (ALS) and the other one based on a scaling of the potential energy release by eddies (PER). Both parameterizations contradict each other in two aspects. First, they predict different functional relationships between the magnitude of the eddy fluxes and the Richardson number (Ri) related to the background state. Second, they also predict different vertical structure functions for the horizontal eddy fluxes. Numerical simulations for two different configurations and for a large range of different background conditions are used to evaluate the parameterizations. It turns out that PER is better suited to capture the Ri dependency of the magnitude of the eddy fluxes. On the other hand, the vertical structure of the meridional eddy fluxes predicted by ALS is more accurate than that of PER, while the vertical structure of the vertical eddy fluxes is well predicted by both parameterizations. Therefore, this study suggests the use of the magnitude of PER and the vertical structure functions of ALS for an improved parameterization of mixed layer eddy fluxes.

Description: Northern Hemisphere stratospheric variability is investigated with respect to chaotic behavior using time series from three different variables extracted from four different reanalysis products and two numerical model runs with different forcing. The time series show red spectra at all frequencies and the probability distribution functions show persistent deviations from a Gaussian distribution. An exception is given by the numerical model forced with perpetual winter conditions—a case that shows more variability and follows a Gaussian distribution, suggesting that the deviation from Gaussianity found in the observations is due to the transition between summer and winter variability. To search for the presence of a chaotic attractor the correlation dimension and entropy, the Lyapunov spectrum, and the associated Kaplan–Yorke dimension are estimated. A finite value of the dimensions can be computed for each variable and data product, with the correlation dimension ranging between 3.0 and 4.0 and the Kaplan–Yorke dimension between 3.3 and 5.5. The correlation entropy varies between 0.6 and 1.1. The model runs show similar values for the correlation and Lyapunov dimensions for both the seasonally forced run and the perpetual-winter run, suggesting that the structure of a possible chaotic attractor is not determined by the seasonality in the forcing, but must be given by other mechanisms.

Description: Southern Hemisphere (SH) stratospheric variability is investigated with respect to chaotic behavior using time series from three different variables extracted from four different reanalysis products. The results are compared with the same analysis applied to the Northern Hemisphere (NH). The probability density functions (PDFs) for the SH show persistent deviations from a Gaussian distribution. The variability is given by white spectra for low frequencies, a slope of −1 for intermediate frequencies, and −3 slopes for high frequencies. Considering the time series for winter and summer separately, PDFs show a Gaussian distribution and the variability spectra change their slopes, indicating the role of the transition between winter and summer variability in shaping the time series. The correlation (D2) and the Kaplan–Yorke (DKY) dimensions are estimated. A finite value of the dimensions can be computed for each variable and data product, except for the NCEP zonal-mean zonal wind and temperature data, which violate the requirement D2 ≤ DKY, possibly owing to the presence of spurious trends and inconsistencies in the data. The value of D2 ranges between 2.6 and 3.9, while DKY ranges between 3.0 and 4.5. The results show that both D2 and DKY display large variability in their values both for different datasets and for different variables within the same dataset. The variability of the values of D2 and DKY thus leaves open the question about the existence of a low-dimensional attractor or if the finite dimensions of the system are the result of the projection of a larger attractor in a low-dimensional embedding space.

Description: The sensitivities of the Brewer–Dobson circulation (BDC) to different distributions of tropical SST heating are investigated in an idealized aquaplanet model. It is found that an increase in tropical SSTs generally leads to an acceleration of tropical upwelling and an associated reduction in the age of air (AOA) in the polar stratosphere and that the AOA near the subtropical tropopause is correlated with local isentropic mixing of tropospheric air with stratospheric air. The zonal distribution of SST perturbations has a major impact on the vertical and meridional structure of the BDC as compared with other SST characteristics. Zonally localized SST heatings tend to generate a shallow acceleration of the stratospheric residual circulation, enhanced isentropic mixing associated with a weakened stratospheric jet, and a reduction in AOA mostly within the polar vortex. In contrast, SST heatings with a zonally symmetric structure tend to produce a deep strengthening of the stratospheric residual circulation, suppressed isentropic mixing associated with a stronger stratospheric jet, and a decrease of AOA in the entire stratosphere. The shallow versus deep strengthening of the stratospheric residual circulation change has been linked to wave propagation and dissipation in the subtropical lower stratosphere rather than wave generation in the troposphere, and the former can be strongly affected by the vertical position of the subtropical jet. These results suggest that, while the longitudinally localized SST trends under climate change may contribute to the change in the shallow branch of the BDC, the upward shift of the subtropical jet associated with the zonal SST heating can impact the deep branch of the BDC.

Description: The authors test the hypothesis that recent observed trends in surface westerlies in the Southern Hemisphere are directly consequent on observed trends in the timing of stratospheric final warming events. The analysis begins by verifying that final warming events have an impact on tropospheric circulation in a simplified GCM driven by specified equilibrium temperature distributions. Seasonal variations are imposed in the stratosphere only. The model produces qualitatively realistic final warming events whose influence extends down to the surface, much like what has been reported in observational analyses. The authors then go on to study observed trends in surface westerlies composited with respect to the date of final warming events. If the considered hypothesis were correct, these trends would appear to be much weaker when composited with respect to the date of the final warming events. The authors find that this is not the case, and accordingly they conclude that the observed surface changes cannot be attributed simply to this shift toward later final warming events.

Description: Precipitation is highly variable in space and time; hence, rain gauge time series generally exhibit additional random small-scale variability compared to area averages. Therefore, differences between daily precipitation statistics simulated by climate models and gauge observations are generally not only caused by model biases, but also by the corresponding scale gap. Classical bias correction methods, in general, cannot bridge this gap; they do not account for small-scale random variability and may produce artifacts. Here, stochastic model output statistics is proposed as a bias correction framework to explicitly account for random small-scale variability. Daily precipitation simulated by a regional climate model (RCM) is employed to predict the probability distribution of local precipitation. The pairwise correspondence between predictor and predictand required for calibration is ensured by driving the RCM with perfect boundary conditions. Wet day probabilities are described by a logistic regression, and precipitation intensities are described by a mixture model consisting of a gamma distribution for moderate precipitation and a generalized Pareto distribution for extremes. The dependence of the model parameters on simulated precipitation is modeled by a vector generalized linear model. The proposed model effectively corrects systematic biases and correctly represents local-scale random variability for most gauges. Additionally, a simplified model is considered that disregards the separate tail model. This computationally efficient model proves to be a feasible alternative for precipitation up to moderately extreme intensities. The approach sets a new framework for bias correction that combines the advantages of weather generators and RCMs.

Description: There is an ongoing discussion in the community concerning the wave-averaged momentum equations in the hybrid vertically Lagrangian and horizontally Eulerian (VL) framework and, in particular, the form stress term (representing the residual effect of pressure perturbations) which is thought to restrict the handling of higher order waves in terms of a perturbation expansion. The present study shows that the traditional pressure-based form stress term can be transformed into a set of terms that do not contain any pressure quantities but do contain the time derivative of a wave-induced velocity. This wave-induced velocity is referred to as the pseudomomentum in the VL framework, as it is analogous to the generalized pseudomomentum in Andrews and McIntyre. This enables the second expression for the wave-averaged momentum equations in the VL framework (this time for the development of the total transport velocity minus the VL pseudomomentum) to be derived together with the vortex force. The velocity-based expression of the form stress term also contains the residual effect of the turbulent viscosity, which is useful for understanding the dissipation of wave energy leading to transfer of momentum from waves to circulation. It is found that the concept of the virtual wave stress of Longuet-Higgins is applicable to quite general situations: it does not matter whether there is wind forcing or not, the waves can have slow variations, and the viscosity coefficient can vary in the vertical. These results provide a basis for revisiting the surface boundary condition used in numerical circulation models.

Description: Mooring observations and model simulations point to an instability of the Labrador Current (LC) during winter, with enhanced eddy kinetic energy (EKE) at periods between 2 to 5 days, and much less EKE during other seasons. Linear stability analysis using vertical shear and stratification from the model reveals three dominant modes of instability in the LC: - a balanced interior mode with along-flow wavelengths of about 30–45 km, phase velocities of 0.3 m/s, maximal growth rates of 1 d−1 and surface intensified, but deep reaching amplitudes, - a balanced shallow mode with along-flow wavelengths of about 0.3–1.5 km, about three times larger phase speeds and growth rates, but amplitudes confined to the mixed layer (ML), - and an unbalanced symmetric mode with largest growth rates, vanishing phase speeds and along-flow structure, and very small cross-flow wavelengths, also confined to the ML. Both balanced modes are akin to baroclinic instability, but operate at moderate to small Richardson numbers Ri with much larger growth rates as for the quasi-geostrophic limit of Ri ≫ 1. The interior mode is found to be responsible for the instability of the LC during winter. Weak stratification and enhanced vertical shear due to local buoyancy loss and the advection of convective water masses from the interior result in small Ri within the LC, and to three times larger growth rates of the interior mode in March compared to summer and fall conditions. Both the shallow and the symmetric mode are not resolved by the model, but it is suggested that they might also play an important role for the instability in the LC and for lateral mixing.

Description: Sensible and latent heat fluxes were estimated from turbulence measurements gathered during several Atlantic transects of the R/V Polarstern. The inertial dissipation method was used to analyze the data. Resulting bulk transfer coefficients were then applied to the data from the ship’s meteorological system to get continuous time series of the heat fluxes. Combined to the measured downward solar and longwave radiation fluxes allows for an estimate of the total energy budget at the air-sea interface. Comparing these parameterized energy fluxes to ones based on the COARE 3.0 bulk flux algorithm show very strong agreement.

Description: The dynamical origin of the spectral and autocorrelation structure of annular variability in the troposphere is investigated by a deductive approach. Specifically, the structure of the power spectrum and autocorrelation function of the zonal-mean geopotential is analyzed for the case of a quasigeostrophic spherical atmosphere subject to a white noise mechanical forcing applied in a single Hough mode and concentrated at a particular level in the vertical, with vertically uniform Newtonian cooling and Rayleigh drag concentrated at a rigid lower boundary. Analytic expressions for the power spectrum are presented together with expressions for an approximate red noise (i.e., a Lorentzian-shaped) power spectrum. It is found that for an infinitely deep atmosphere the power spectrum can be well approximated by a red noise process for the first few Hough modes (associated with large Rossby heights), provided the distance from the forcing is not larger than about one Rossby height. When a frictional rigid lower boundary is included, however, the approximation is generally bad. The high-frequency part of the power spectrum exhibits near-exponential behavior and the autocorrelation function shows a transition from a rapid decay at short lags to a much slower decay at longer lags, if the thermal and mechanical damping time scales are sufficiently well separated. Since observed annular variability exhibits the same characteristics, the above results lead to the hypothesis that these characteristics may, to some extent, be intrinsic to the linear zonal-mean response problem—although the need for an additional contribution from eddy feedbacks is also implied by the results.

Description: A representation of an equatorial basin mode excited in a shallow water model for a single high order baroclinic vertical normal mode is used as a simple model for the equatorial deep jets. The model is linearized about both a state of rest and a barotropic mean flow corresponding to the observed Atlantic Equatorial Intermediate Current System. We found that the eastward mean flow associated with the North and South Intermediate Counter Currents (NICC and SICC, respectively) effectively shields the Equator from off-equatorial Rossby waves. The westward propagation of these waves is blocked and focusing on the Equator due to beta dispersion is prevented. This leads to less energetic jets along the Equator. On the other hand, the westward barotropic mean flow along the Equator reduces the gradient of absolute vorticity and hence widens the cross-equatorial structure of the basin mode. Increasing lateral viscosity predominantly affects the width of the basin modes’ Kelvin wave component in the presence of the mean flow while the Rossby wave is confined by the flanking NICC and SICC. Independent of the presence of the mean flow, the application of sufficient lateral mixing also hinders the focusing of off-equatorial Rossby waves, which is hence an unlikely feature of a low-frequency basin mode in the real ocean.

Description: A decadal change in the character of ENSO was observed around year 2000 toward weaker-amplitude, higher-frequency events with an increased occurrence of central Pacific El Niños. Here these changes are assessed in terms of the Bjerknes stability index (BJ index), which is a measure of the growth rate of ENSO-related SST anomalies. The individual terms of the index are calculated from ocean reanalysis products separately for the time periods 1980–99 and 2000–10. The spread between the products is large, but they show a robust weakening of the thermocline feedback due to a reduced thermocline slope response to anomalous zonal wind stress as well as a weakened wind stress response to eastern equatorial Pacific SST anomalies. These changes are consistent with changes in the background state of the tropical Pacific: cooler mean SST in the eastern and central equatorial Pacific results in reduced convection there together with a westward shift in the ascending branch of the Walker circulation. This shift leads to a weakening in the relationship between eastern Pacific SST and longitudinally averaged equatorial zonal wind stress. Also, despite a steeper mean thermocline slope in the more recent period, the thermocline slope response to wind stress anomalies weakened due to a smaller zonal wind fetch that results from ENSO-related wind anomalies being more confined to the western basin. As a result, the total BJ index is more negative, corresponding to a more strongly damped system in the past decade compared to the 1980s and 1990s.

Description: In his comment, G. Bürger criticizes the conclusion that inflation of trends by quantile mapping is an adverse effect.He assumes that the argument would be ‘‘based on the belief that long-term trends and along with them future climate signals are to be large scale.’’ His line of argument reverts to the so-called inflated regression. Here it is shown, by referring to previous critiques of inflation and standard literature in statistical modeling as well as weather forecasting, that inflation is built upon a wrong understanding of explained versus unexplained variability and prediction versus simulation. It is argued that a sound regressionbased downscaling can in principle introduce systematic local variability in long-term trends, but inflation systematically deteriorates the representation of trends. Furthermore, it is demonstrated that inflation by construction deteriorates weather forecasts and is not able to correctly simulate small-scale spatiotemporal structure.

Description: The relationship between the Agulhas Current and the Agulhas leakage is not well understood. Here, this is investigated using two basin-scale and two global ocean models, of incrementally increasing resolution. The response of the Agulhas Current is evaluated under a series of sensitivity experiments, in which idealised anomalies, designed to geometrically modulate zonal trade wind stress, are applied across the Indian Ocean basin. The imposed wind stress changes exceed ±2 standard deviations from the annual mean trade winds and, in the case of intensification, are partially representative of recently observed trends. The Agulhas leakage is quantified using complimentary techniques based on Lagrangian virtual floats and Eulerian passive tracer flux. As resolution increases, model behavior converges and the sensitivity of the leakage to Agulhas Current transport anomalies is reduced. In the two eddy-resolving configurations tested, the leakage is insensitive to changes in Agulhas Current transport at 32°S, though substantial eddy kinetic energy anomalies are evident. Consistent with observations, the position of the retroflection remains stable. The decoupling of Agulhas Current variability from the Agulhas leakage suggests that, while correlations between the two may exist, they may not have a clear dynamical basis. It is suggested that present and future Agulhas leakage proxies be considered in the context of potentially transient forcing regimes.

Description: A surface diurnal warm layer is diagnosed from Seaglider observations, and develops on half the days in the CINDY/DYNAMO Indian Ocean experiment. The diurnal warm layer occurs on days of high solar radiation flux (〉 80 W m−2) and low wind speed (〈 6 m s−1), and preferentially in the inactive stage of the Madden–Julian Oscillation. Its diurnal harmonic has an exponential vertical structure with a depth scale of 4–5 m (dependent on chlorophyll concentration), consistent with forcing by absorption of solar radiation. The effective sea surface temperature (SST) anomaly due to the diurnal warm layer often reaches 0.8°C in the afternoon, with a daily mean of 0.2°C, rectifying the diurnal cycle onto longer time scales. This SST anomaly drives an anomalous flux of 4 W m−2 that cools the ocean. Alternatively, in a climate model where this process is unresolved, this represents an erroneous flux that warms the ocean. A simple model predicts a diurnal warm layer to occur on 30–50% of days across the tropical warm pool. On the remaining days, with low solar radiation and high wind speeds, a residual diurnal cycle is observed by the Seaglider, with a diurnal harmonic of temperature that decreases linearly with depth. As wind speed increases, this already weak temperature gradient decreases further, tending towards isothermal conditions.

Description: Sea surface temperature (SST) anomalies in the eastern equatorial Atlantic are connected to modulations in the strength of the South Atlantic subtropical high-pressure system, referred to as the South Atlantic Anticyclone (SAA). Using ocean and atmosphere reanalysis products we show here that the strength of the SAA from February to May impacts the timing of the cold tongue onset and the intensity of its development in the eastern equatorial Atlantic (EEA) via anomalous tropical wind power. This modulation of the timing and amplitude of the seasonal cold tongue development manifests as anomalous SST events peaking between June and August. The timing and impact of this connection is not completely symmetric for warm and cold events. For cold events, an anomalously strong SAA in February and March leads to positive wind power anomalies from February to June resulting in an early cold tongue onset and subsequent cold SST anomalies in June and July. For warm events the anomalously weak SAA persists until May, generating negative wind power anomalies that lead to a late cold tongue onset as well as a suppression of the cold tongue development and associated warm SST anomalies. Mechanisms by which SAA induced wind power variations south of the equator influence EEA SST are discussed, including ocean adjustment via Rossby and Kelvin wave propagation, meridional advection, and local intraseasonal wind variations

Description: We present a detailed quality assessment of a novel underwater sensor for the measurement of CO2 partial pressure (pCO2) based on surface water field deployments carried out between 2008 and 2011. The commercially available sensor, which is based on membrane equilibration and NDIR spectrometry is small and can be integrated into mobile platforms. It is calibrated in water against a proven flow-through pCO2 instrument within a custom-built calibration setup. The aspect of highest concern with respect to achievable data quality of the sensor is the compensation for signal drift inevitably connected to absorption measurements. We use three means to correct for drift effects: (i) a filter correlation or dual-beam setup, (ii) regular zero gas measurements realized automatically within the sensor and (iii) a zero-based transformation of two sensor calibrations flanking the time of sensor deployment. Three sensors were tested against an underway pCO2 system during two major research cruises providing an in situ temperature range from 7.4 to 30.1°C and pCO2 values between 289 and 445 μatm. The average difference between sensor and reference pCO2 was found to be -0.6 ± 3 μatm with a RMSE of 3.7 μatm.

Description: Ammassalik in southeast Greenland is known for strong wind events that can reach hurricane intensity and cause severe destruction in the local town. Yet, these winds and their impact on the nearby fjord and shelf region have not been studied in detail. Here, data from two meteorological stations and the European Centre for Medium-Range Weather Forecasts Interim Re-Analysis (ERA-Interim) are used to identify and characterize these strong downslope wind events, which are especially pronounced at a major east Greenland fjord, Sermilik Fjord, within Ammassalik. Their local and regional characteristics, their dynamics and their impacts on the regional sea ice cover, and air–sea fluxes are described. Based on a composite of the events it is concluded that wind events last for approximately a day, and seven to eight events occur each winter. Downslope wind events are associated with a deep synoptic-scale cyclone between Iceland and Greenland. During the events, cold dry air is advected down the ice sheet. The downslope flow is accelerated by gravitational acceleration, flow convergence inside the Ammassalik valley, and near the coast by an additional thermal and synoptic-scale pressure gradient acceleration. Wind events are associated with a large buoyancy loss over the Irminger Sea, and it is estimated that they drive one-fifth of the net wintertime loss. Also, the extreme winds drive sea ice out of the fjord and away from the shelf.

Description: In 2013, the vast majority of the monitored climate variables reported here maintained trends established in recent decades. ENSO was in a neutral state during the entire year, remaining mostly on the cool side of neutral with modest impacts on regional weather patterns' around the world. This follows several years dominated by the effects of either La Nina. or El Nino events. According to several independent analyses, 2013 was again among the 10 warmest years on record at the global scale, both at the Earth's surface and through the troposphere. Some regions in the Southern Hemispherehad record or near-record high temperatures for the year. Australia observed its hottest year on record, while Argentina and New Zealand reported their second and third hottest years, respectively. In Antarctica, Amundsen-Scott South Pole Station reported its highest annual temperature since records began in 1957 At the opposite pole, the Arctic observed its seventh warmest year since record's began in the early 20th century. At 20-m depth, record high temperatures were measured at some permafrost stations on the North Slope Of Alaska and in the Brooks Range. In the Northern Hemisphere extratropics, anomalous meridional atmospheric circulation occurred throughout much of the year leading to marked regional extremes of both temperature and precipitation. Cold temperature anomalies during Winter across Eurasia were followed by warm Spring temperature anomalies, which, were linked to a new record Eurasian snow cover extent in the May. Minimum sea ice extent in the :Arctic was the sixth lowest since satellite Observations began in 1979. Including 2013, all seven lowest extents on record have occurred in the past seven years Antarctica, on the other hand, had above average sea ice extent throughout 2013, with 116 days Of new daily high extent records, inclding a new daily maximum sea ice area of 19.57 million km(2) reached on 1 October. ENSO-neutral conditions in the eastern central Pacific Ocean and a negative Pacific decadal oscillation pattern in the North Pacific had the largest impacts on the global sea surface temperature in 2013: The North Pacific reached a historic high temperature in 2013 and on balance the globally-averaged sea surface temperature was among the 10 highest on record. Overall, the salt Content in near-surface ocean waters increased while in intermediate waters it decreased. Global mean sea level continued to rise during 2013, on pace with a trend of 3,2 mm yr(-1) over the past two decades. A portion of this trend (0.5 mm yr(-1)) has been attributed to natural variability associated with the Pacific decadal oscillation as Well as to ongoing contribution from the melting of glaciers and ice sheets and ocean warming. Global tropical cyclone frequency during 2013 was slightly above average with a total of 94 storms, although the North Atlantic Basin had its quietest hurricane season since 1994. In the Western North Pacific Basin, Super Typhoon Haiyan, the deadliest tropical Cyclone of 2013, had 1-minute sustained winds estimated to be 170 kt (OS m s(-1)) on 7 November, the highest Wind speed ever assigned to a tropical cyclone. High storm surge was also associated with Haiyan at it Made landfall over the central Philippines, an area where sea level is currently at historic highs increasing by 200 mm since 1970. In the atmosphere, carbon dioxide, methane, and nitrous oxide all Continued to increase in 2013. As in previous our years, each of these Major greenhouse gases once again reached historic high Concentrations. In the Arctic carbon dioxide and methane:, Increased at the same rate as:the global increase. These increases are likely due to export from lower latitudes rather than a consequence Of increases in Arctic sources such as thawing permafrost., At Mauna Loa, Hawaii, for the first time since measurement began in 1958, the daily average Mixing ratio Of carbon dioxide exceeded 400 ppm on 9 May. The state of these variables, along with dozens of others, and the 2013 climate conditions of regions around the world are discussed in further detail in this 24th, edition of the State of the Climate series