The origin of sub-surface source waters define the sea–air flux of methane in the Mauritanian Upwelling,NW Africa

Concentrations and flux densities of methane were determined during a Lagrangian study of an advective filament in the permanent upwelling region off western Mauritania. Newly upwelled waters were dominated by the presence of North Atlantic Central Water and surface CH₄ concentrations of 2.2 ± 0.3 nmol L⁻¹ were largely in equilibrium with atmospheric values, with surface saturations of 101.7 ± 14%. As the upwelling filament aged and was advected offshore, CH₄ enriched South Atlantic Central Water from intermediate depths of 100–350 m was entrained into the surface mixed layer of the filament following intense mixing associated with the shelf break. Surface saturations increased to 198.9 ± 15% and flux densities increased from a mean value over the shelf of 2.0 ± 1.1 μmol m⁻² d⁻¹ to a maximum of 22.6 μmol m⁻² d⁻¹. Annual CH₄ emissions for this persistent filament were estimated at 0.77 ± 0.64 Gg which equates to a maximum of 0.35% of the global oceanic budget. This raises the known outgassing intensity of this area and highlights the importance of advecting filaments from upwelling waters as efficient vehicles for air-sea exchange.     

Exploring the vertical extent of breaking internal wave turbulence above deep-sea topography

A mooring equipped with 200 high-resolution temperature sensors between 6 and 404 m above the bottom was moored in 1890 m water depth above a steep, about 10° slope of Mount Josephine, NE-Atlantic. The sensors have a precision of less than 0.5 mK. They are synchronized via induction every 4 h so that the 400 m range is measured to within 0.02 s, every 1 s. Thin cables and elliptical buoyancy assured vertical mooring motions to be smaller than 0.1 m under maximum 0.2 m s⁻¹ current speeds. The local bottom slope is supercritical for semidiurnal internal tides by a factor of two. Exploring a one-month record in detail, the observations show: 1/semidiurnal tidal dominance in variations of dissipation rate ε, eddy diffusivity K_z and temperature, but no significant correlation between the records of ε and total kinetic energy, 2/a variation with time over four orders of magnitude of 100-m vertically averaged ε, 3/a local minimum in density stratification between 50 and 100 m above the bottom, 4/a gradual decrease in daily or longer averaged ε and K_z by one order of magnitude over a vertical distance of 250 m, upwards from 150 m above the bottom, 5/monthly mean values of <[ε]> = 2 ± 0.5 × 10⁻⁷ m² s⁻³, <[K_z]> = 8 ± 3 × 10⁻³ m² s⁻¹ averaged over the lower 150 m above the bottom.     

Climate change impacts on ecosystems and the terrestrial carbon sink: a new assessment

Climate output from the UK Hadley Centre's HadCM2 and HadCM3 experiments for the period 1860 to 2100, with IS92a greenhouse gas forcing, together with predicted patterns of N deposition and increasing CO₂, were input (offline) to the dynamic vegetation model, Hybrid v4.1 (Friend et al., 1997; Friend and White, 1999). This model represents biogeochemical, biophysical and biogeographical processes, coupling the carbon, nitrogen and water cycles on a sub-daily timestep, simulating potential vegetation and transient changes in annual growth and competition between eight generalized plant types in response to climate.Global vegetation carbon was predicted to rise from about 600 to 800 PgC (or to 650 PgC for HadCM3) while the soil carbon pool of about 1100 PgC decreased by about 8%. By the 2080s, climate change caused a partial loss of Amazonian rainforest, C₄ grasslands and temperate forest in areas of southern Europe and eastern USA, but an expansion in the boreal forest area. These changes were accompanied by a decrease in net primary productivity (NPP) of vegetation in many tropical areas, southern Europe and eastern USA (in response to warming and a decrease in rainfall), but an increase in NPP of boreal forests. Global NPP increased from 45 to 50 PgC y⁻¹ in the 1990s to about 65 PgC y⁻¹ in the 2080s (about 58 PgC y⁻¹ for HadCM3). Global net ecosystem productivity (NEP) increased from about 1.3 PgC y⁻¹ in the 1990s to about 3.6 PgC y⁻¹ in the 2030s and then declined to zero by 2100 owing to a loss of carbon from declining forests in the tropics and at warm temperate latitudes — despite strengthening of the carbon sink at northern high latitudes. HadCM3 gave a more erratic temporal evolution of NEP than HadCM2, with a dramatic collapse in NEP in the 2050s.     

Kinematic metrics of the Loop Current in the Gulf of Mexico from satellite altimetry

We analyzed a 20-year time series (January 1st, 1993 through December 31st, 2012) of Loop Current (LC) surface area derived from satellite altimetry in the eastern Gulf of Mexico to estimate kinematical metrics of this potent flow. On average the LC intrudes to its maximum northward position about 216 ± 126 days after the previous eddy separation; and ∼30 ± 31 days later sheds a large anticyclonic eddy. When the northern extent of the LC intrusion following the previous eddy separation is greater than 27°N, the current retreats very quickly until it sheds another eddy with the entire separation process occurring on the order of 30 days. To first order the change in areal extent of the LC during intrusion into the Gulf occurs at an average rate of 225 km² day⁻¹, which corresponds to an intrusion velocity of 1.7 cm s⁻¹ of the LC front, and adds Caribbean water to the Gulf at a rate of 2.6 ± 0.7 Sv.     

Intraseasonal variability of winter precipitation over central asia and the western tibetan plateau from 1979 to 2013 and its relationship with the North Atlantic Oscillation

Winter precipitation over Central Asia and the western Tibetan Plateau (CAWTP) is mainly a result of the interaction between the westerly circulation and the high mountains around the plateau. Empirical Orthogonal Functions (EOFs), Singular Value Decomposition (SVD), linear regression and composite analysis were used to analyze winter daily precipitation and other meteorological elements in this region from 1979 to 2013, in order to understand how interactions between the regional circulation and topography affect the intraseasonal variability in precipitation. The SVD analysis shows that the winter daily precipitation variability distribution is characterized by a dipole pattern with opposite signs over the northern Pamir Plateau and over the Karakoram Himalaya, similar to the second mode of EOF analysis. This dipole pattern of precipitation anomaly is associated with local anomalies in both the 700 hPa moisture transport and the 500 hPa geopotential height and is probably caused by oscillations in the regional and large-scale circulations, which can influence the westerly disturbance tracks and water vapor transport. The linear regression shows that the anomalous mid-tropospheric circulation over CAWTP corresponds to an anti-phase variation of the 500 hPa geopotential height anomalies over the southern and northern North Atlantic 10 days earlier (at 95% significance level), that bears a similarity to the North Atlantic Oscillation (NAO). The composite analysis reveals that the NAO impacts the downstream regions including CAWTP by controlling south-north two branches of the middle latitude westerly circulation around the Eurasian border. During the positive phases of the NAO, the northern branch of the westerly circulation goes around the northwest Tibetan Plateau, whereas the southern branch encounters the southwest Tibetan Plateau, which leads to reduced precipitation over the northern Pamir Plateau and increased precipitation over the Karakoram Himalaya, and vice versa.     

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.     

Observations and operational model simulations reveal the impact of Hurricane Matthew (2016) on the Gulf Stream and coastal sea level

In October 7–9, 2016, Hurricane Matthew moved along the southeastern coast of the U.S., causing major flooding and significant damage, even to locations farther north well away from the storm’s winds. Various observations, such as tide gauge data, cable measurements of the Florida Current (FC) transport, satellite altimeter data and high-frequency radar data, were analyzed to evaluate the impact of the storm. The data show a dramatic decline in the FC flow and increased coastal sea level along the U.S. coast. Weakening of the Gulf Stream (GS) downstream from the storm’s area contributed to high coastal sea levels farther north. Analyses of simulations of an operational hurricane-ocean coupled model reveal the disruption that the hurricane caused to the GS flow, including a decline in transport of ∼20 Sv (1 Sv = 10⁶ m³ s⁻¹). In comparison, the observed FC reached a maximum transport of ∼40 Sv before the storm on September 10 and a minimum of ∼20 Sv after the storm on October 12. The hurricane impacts both the geostrophic part of the GS and the wind-driven currents, generating inertial oscillations with velocities of up to ±1 m s⁻¹. Analysis of the observed FC transport since 1982 indicated that the magnitude of the current weakening in October 2016 was quite rare (outside 3 standard deviations from the mean). Such a large FC weakening in the past occurred more often in October and November, but is extremely rare in June-August. Similar impacts on the FC from past tropical storms and hurricanes suggest that storms may contribute to seasonal and interannual variations in the FC. The results also demonstrated the extended range of coastal impacts that remote storms can cause through their influence on ocean currents.     

One year measurements of aerosol optical properties over an urban coastal site: Effect on local direct radiative forcing

We present results of direct aerosol radiative forcing over a French Mediterranean coastal zone based on one year of continuous observations of aerosol optical properties during 2005–2006. Monthly-mean aerosol optical depth at 440 nm ranged between 0.1 and 0.34, with high Angstrom coefficient (α > 1.2). The single scattering albedo (at 525 nm) estimated at the surface ranged between 0.7 and 0.8, indicating significant absorption. The presence of aerosols over the Mediterranean zone during summer decreases the shortwave radiation reaching the surface by as much as 26 ± 3.9 W m^− 2, and increases the top of the atmosphere reflected radiation by as much as 5.2 ± 1.0 W m^− 2. The shortwave atmospheric absorption translates to an atmospheric heating of 2.5 to 4.6 K day^− 1. Concerted efforts are needed for investigating the possible impact of the increase in heating rate on the maintenance of heat-waves frequently occurring over this coastal region during summer time.     

Water on an urban planet: Urbanization and the reach of urban water infrastructure

Urban growth is increasing the demand for freshwater resources, yet surprisingly the water sources of the world's large cities have never been globally assessed, hampering efforts to assess the distribution and causes of urban water stress. We conducted the first global survey of the large cities’ water sources, and show that previous global hydrologic models that ignored urban water infrastructure significantly overestimated urban water stress. Large cities obtain 78 ± 3% of their water from surface sources, some of which are far away: cumulatively, large cities moved 504 billion liters a day (184 km³ yr⁻¹) a distance of 27,000 ± 3800 km, and the upstream contributing area of urban water sources is 41% of the global land surface. Despite this infrastructure, one in four cities, containing $4.8 ± 0.7 trillion in economic activity, remain water stressed due to geographical and financial limitations. The strategic management of these cities’ water sources is therefore important for the future of the global economy.     

Carbonaceous materials in size-segregated atmospheric aerosols from urban and coastal-rural areas at the Western European Coast

A high-volume cascade impactor, equipped with a PM10 inlet, was used to collect size-segregated aerosol samples during the summer of 2004 at two Portuguese locations: a coastal-rural area (Moitinhos) and an urban area (Oporto). Concentrations of airborne particulate matter (PM), total carbon (TC), organic carbon (OC), elemental carbon (EC), and water-soluble organic carbon (WSOC) were determined for the following particle size ranges: < 0.49, 0.49–0.95, 0.95–3.0, and 3.0–10 µm. The total PM mass concentrations at the urban and coastal-rural sites ranged from 22.8 to 79.6 μg m^− 3 and 19.9 to 28.2 μg m^− 3, respectively, and more than 56% of the total aerosol mass was found in the fractions below 3.0 μm. At both locations the highest concentrations of OC and EC were found in the submicrometer size range. The regional variability for the OC and EC concentrations, with the highest concentrations being found in the urban area, was related to the contribution of local primary sources (mostly traffic emissions). It was also verified an enrichment of the small size particles in WSOC, representing on average 37.3(± 12.4)% and 59.7(± 18.0)% of OC in the very fine aerosol at the coastal-rural and urban areas, respectively. The amount of secondary OC calculated by the minimum OC/EC ratio method indicates that secondary organic aerosol formation was important throughout the study at both sites. The obtained results suggest that long-range transport and favourable summer conditions for photochemical oxidation are key factors determining secondary OC formation in the coastal-rural and urban areas. The ultraviolet absorption properties of the chromophoric constituents of the WSOC fractions were also different among the different particle size ranges and also between the two sampling locations, thus suggesting the strong impact of the diverse emission sources into the composition of the size-segregated organic aerosol.     

Investigations on meteorological conditions for elevated PM2.5 in Fairbanks,Alaska

The relationships between meteorological conditions (temperature, wind-speed and direction, relative humidity, surface-inversion depth and strength, and stability) and PM_2.5 concentrations in Fairbanks, Alaska were investigated using ten years of observational data. The results show that during wintertime (November through February) PM_2.5 concentrations exceeding the 24 h National Air Quality Standard (35 μg/m³) occurred under calm wind, extremely low temperature (≤20 °C) and moisture (water-vapor pressure < 2 hPa) multiday surface-inversion conditions that trap the pollutants in the breathing level and inhibit transport of polluted air out of Fairbanks. PM_2.5 concentrations tend to be higher under stable than other conditions, but are not sensitive to the degree of stability. The presence of a surface inversion and calm wind are necessary, but in combination with low temperatures and humidity, the conditions are sufficient for high PM_2.5 concentrations. The low temperatures are required because they lead to increased emission rates from domestic heating and power production. During multiday inversions with temperatures above ? 20 °C, high relative humidity (> 75%) partly caused by water-vapor emission reduces PM_2.5 concentrations.     

Size-resolved aerosol chemical concentrations at rural and urban sites in Central California,USA

Aerosol size distributions were measured with Micro Orifice Uniform Deposit Impactor (MOUDI) cascade impactors at the rural Angiola and urban Fresno Supersites in California's San Joaquin Valley during the California Regional PM₁₀/PM_2.5 Air Quality Study (CRPAQS) winter campaign from December 15, 2000 to February 3, 2001. PM_2.5 filter samples were collected concurrently at both sites with Sequential Filter Samplers (SFS). MOUDI nitrate (NO₃⁻) concentrations reached 66 μg/m³ on January 6, 2001 during the 1000–1600 PST (GMT-8) period. Pair-wise comparisons between PM_2.5 MOUDI and SFS concentrations revealed high correlations at the Angiola site (r > 0.93) but more variability (r < 0.85) at the Fresno site for NO₃⁻, sulfate (SO₄⁼), and ammonium (NH₄⁺). Correlations were higher at Fresno (r > 0.87) than at Angiola (r < 0.7) for organic carbon (OC), elemental carbon (EC), and total carbon (TC). NO₃⁻ and SO₄⁼ size distributions in Fresno were multi-modal and wider than the uni-modal distributions observed at Angiola. Geometric mean diameters (GMD) were smaller for OC and EC than for NO₃⁻ and SO₄⁼ at both sites. OC and EC were more concentrated on the lowest MOUDI stage (0.056 µm) at Angiola than at Fresno. The NO₃⁻ GMD increased from 0.97 to 1.02 µm as the NO₃⁻ concentration at Angiola increased from 43 to 66 µg m^− 3 during a PM_2.5 episode from January 4–7, 2001. There was a direct relationship between GMD and NO₃⁻ and SO₄⁼ concentrations at Angiola but no such relationships for OC or EC. This demonstrates that secondary aerosol formation increases both concentration and particle size for the rural California environment.     

Distribution of CO2 parameters in the Western Tropical Atlantic Ocean

The variability of sea surface Total Alkalinity (TA) and sea surface Total Inorganic Carbon (C_T) is examined using all available data in the western tropical Atlantic (WTA: 20°S-20°N, 60°W-20°W). Lowest TA and C_T are observed for the region located between 0°N-15°N/60°W-50°W and are explained by the influence of the Amazon plume during boreal summer. In the southern part of the area, 20°S-10°S/40°W-60°W, the highest values of TA and C_T are linked to the CO₂–rich waters due to the equatorial upwelling, which are transported by the South Equatorial Current (SEC) flowing from the African coast to the Brazilian shore. An increase of C_T of 0.9 ± 0.3 μmol kg⁻¹yr⁻¹ has been observed in the SEC region and is consistent with previous published estimates. A revised C_T-Sea Surface Salinity (SSS) relationship is proposed for the WTA to take into account the variability of C_T at low salinities. This new C_T-SSS relationship together with a published TA-SSS relationship allow to calculate pCO₂ values that compare well with observed pCO₂ (R² = 0.90).     

A method for determining thermophysical properties of organic material in aqueous solutions: Succinic acid

A method for determining evaporation rates and thermodynamic properties of aqueous solution droplets is introduced. The method combines evaporation rate measurements using modified TDMA technique with data evaluation using an accurate evaporation model. The first set of data has been collected and evaluated for succinic acid aqueous solution droplets.Evaporation rates of succinic acid solution droplets have been measured using a TDMA system at controlled relative humidity (65%) and temperature (298 K). A temperature-dependent expression for the saturation vapour pressure of pure liquid phase succinic acid at atmospheric temperatures has been derived by analysing the evaporation rate data with a numerical model. The obtained saturation vapour pressure of liquid phase succinic acid is ln(p) = 118.41 − 16204.8/T − 12.452ln(T). The vapour pressure is in unit of Pascal and the temperature in Kelvin. A linear expression for the enthalpy of vaporization for liquid state succinic acid is also presented.According to the results presented in the following, a literature expression for the vapour pressure of liquid phase succinic acid defined for temperatures higher than 461 K [Yaws, C.L., 2003. Yaws' Handbook of Thermodynamic and Physical Properties of Chemical Compounds, Knovel] can be extrapolated to atmospheric temperatures with very good accuracy. The results also suggest that at 298 K the mass accommodation coefficient of succinic acid is unity or very close to unity.     

The possible role of Brazilian promontory in Little Ice Age

The Gulf Stream, one of the strongest currents in the world, transports approximately 31 Sv of water (Kelly and Gille, 1990, Baringer and Larsen, 2001, Leaman et al., 1995) and 1.3 × 10¹⁵ W (Larsen, 1992) of heat into the Atlantic Ocean, and warms the vast European continent. Thus any change of the Gulf Stream could lead to the climate change in the European continent, and even worldwide (Bryden et al., 2005). Past studies have revealed a diminished Gulf Stream and oceanic heat transport that was possibly associated with a southward migration of intertropical convergence zone (ITCZ) and may have contributed to Little Ice Age (AD ∼1200 to 1850) in the North Atlantic (Lund et al., 2006). However, the causations of the Gulf Stream weakening due to the southward migration of the ITCZ remain uncertain. Here we use satellite observation data and employ a model (oceanic general circulation model – OGCM) to demonstrate that the Brazilian promontory in the east coast of South America may have played a crucial role in allocating the equatorial currents, while the mean position of the equatorial currents migrates between northern and southern hemisphere in the Atlantic Ocean. Northward migrations of the equatorial currents in the Atlantic Ocean have little influence on the Gulf Stream. Nevertheless, southward migrations, especially abrupt large southward migrations of the equatorial currents, can lead to the increase of the Brazil Current and the significant decrease of the North Brazil Current, in turn the weakening of the Gulf Stream. The results from the model simulations suggest the mean position of the equatorial currents in the Atlantic Ocean shifted at least 180–260 km southwards of its present-day position during the Little Ice Age based on the calculations of simple linear equations and the OGCM simulations.     

Impressions of the turbulence variability in a weakly stratified,flat-bottom deep-sea ‘boundary layer’

The character of turbulent overturns in a weakly stratified deep-sea is investigated in some detail using 144 high-resolution temperature sensors at 0.7 m intervals, starting 5 m above the bottom. A 9-day, 1 Hz sampled record from the 912 m depth flat-bottom (<0.5% bottom-slope) mooring site in the central-north Alboran Sea (W-Mediterranean) demonstrates an overall conservative temperature range of only 0.05 °C, a typical mean buoyancy period as large as 3 h and a 1 Hz-profile-vertically-averaged turbulence dissipation rate maximum of only 10⁻⁸ m² s⁻³. Nonetheless, this ‘boundary layer’ varies in height between <6 and >104 m above the bottom and is thus not homogeneous throughout; the temperature variations are seldom quiescent and are generally turbulent in appearance, well exceeding noise levels. The turbulence character is associated with small-scale internal waves; examples are found of both shear- and convection-driven turbulence; particular association, although not phase-locked, is found between turbulence variations and tidal rather than with inertial motions; the mean buoyancy frequency of a few times the inertial frequency implies the importance of ‘slantwise convection’ in the direction of the earth rotational vector rather than in the direction of gravity. Such convection is observed both in near-homogeneous and weakly stratified form.     

Globalisation of water resources: international virtual water flows in relation to crop trade

The water that is used in the production process of a commodity is called the ‘virtual water’ contained in the commodity. International trade of commodities brings along international flows of virtual water. The objective of this paper is to quantify the volumes of virtual water flows between nations in the period 1995–1999 insofar related to international crop trade and to analyse national virtual water balances in relation to national water needs and water availability. The basic approach is to multiply international crop trade flows (ton/yr) by their associated virtual water content (m³ ton⁻¹). The calculations show that the global volume of crop-related international virtual water flows between nations was 695 Gm³ yr⁻¹ in average over the period 1995–1999. For comparison: the total water use by crops in the world has been estimated at 5400 Gm³ yr⁻¹. This means that 13% of the water used for crop production in the world is not used for domestic consumption but for export (in virtual form). This is a conservative estimate because only a limited number of crops––although the most important ones––have been taken into account and because crop products (such as cotton clothes) have been excluded from the study. The countries with the largest net virtual water export are United States, Canada, Thailand, Argentina and India. The largest net import appears to be in Japan, the Netherlands, the Republic of Korea, China and Indonesia.     

Distribution of inorganic nitrogen-containing species in atmospheric particles from an island in the Yellow Sea

In this study, total suspended particles (TSP) and size-segregated atmospheric aerosol samples were measured on Qianliyan Island in the Yellow Sea in spring (April–May), summer (July–August) and fall (October–November) of 2006 and in water (January–February) of 2007. The mass concentration of the TSP varied from 75.6 to 132.0 μg/m³. The average concentration were 9.37 ± 7.56 μg/m³ and 5.32 ± 4.25 μg/m³ for nitrate and ammonium in the TSP, respectively. TSP concentration showed a significant correlation with those of nitrate (n = 27, r = 0.73) and ammonium (n = 27, r = 0.60). The mass-size distribution of atmospheric particles exhibited two modes with an accumulation mode at 0.43–1.1 μm and a coarse mode at 3.3–4.7 μm throughout the sampling months. A bi-modal size distribution of nitrate in concentration occurred in the April–May, October–November and January–February, but a uni-modal size distribution occurred in the August. The uni-modal size distribution of ammonium at 0.43–0.65 μm was observed throughout the sampling months. The average of inorganic nitrogen in mass concentration accounted for 4.0% of the total mass of aerosol particles while ammonium-N was the dominant fraction of TIN (Total Inorganic Nitrogen), contributing to 62–71% of the TIN.     

Influence of spatial and temporal scales on statistical analyses of rainfall variability in the River Nile basin

In this study, empirical orthogonal function was applied to analyze rainfall variability in the Nile basin based on various spatio-temporal scales. The co-occurrence of rainfall variability and the variation in selected climate indices was analyzed based on various spatio-temporal scales. From the highest to the lowest, the cumulative amount of variance explained by the first two principal components (PCs) for any selected size of the spatial domain was obtained for the annual, seasonal, and monthly rainfall series respectively. The variability in the annual rainfall of 1° × 1° spatial coverage explained by only the first PC was about 55% on average. However, this percentage reduced to about 40% on average across the study area when the size of the spatial domain was increased from 1° × 1° to 10° × 10°. The variation in climate indices was shown to explain rainfall variability more suitably at a regional than location-specific spatial scale. The magnitudes and sometimes signs of the correlation between rainfall variability and the variation in climate indices tended to vary from one time scale to another. These findings are vital in the selection of spatial and temporal scales for more considered attribution of rainfall variability across the study area.     

Spatio-temporal variability of atmospheric rivers and associated atmospheric parameters in the Euro-Atlantic region

We study the spatio-temporal variability of Atmospheric Rivers (ARs) and associated integrated water vapor and atmospheric parameters over the Euro-Atlantic region using long-term reanalysis datasets. Winds, temperature, and specific humidity at different pressure levels during 1979–2018 are used to study the water vapor transport integrated between 1000 and 300 hPa (IVT300) in mapping ARs. The intensity of ARs in the North Atlantic has been increasing in recent times (2009–2018) with large decadal variability and poleward shift (~ 5° towards the North) in landfall during 1999–2018. Though different reanalysis datasets show similar spatial patterns of IVT300 in mapping ARs, bias in specific humidity and wind components led to IVT300 mean bias of 50 kg m⁻¹ s⁻¹ in different reanalysis products compared to ERA5. The magnitude of winds and specific humidity in the lower atmosphere (below 750 hPa) dominates the total column water vapor and intensity of ARs in the North Atlantic. Reanalysis datasets in the central North Atlantic show an IVT300 standard deviation of 200 kg m⁻¹ s⁻¹ which is around 33% of the ARs climatology (~ 600 kg m⁻¹ s⁻¹). Though ARs have a higher frequency of landfalling over Western Europe in winter half-year, the intensity of IVT300 in winter ARs is 3% lower than the annual mean. The lower frequency of ARs in the summer half-year shows 3% higher IVT300 than the annual mean. While ARs in the North Atlantic show a strong decadal change in frequency and path, the impact of the North Atlantic Oscillation (NAO) and Scandinavian blocking on the location of landfall of ARs are significant. Furthermore, there is a strong latitudinal dependence of the source of moisture flux in the open ocean, contributing to the formation and strengthening ARs.