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.     

Variability of the Somali Current and eddies during the southwest monsoon regimes

The meso-scale eddies and currents in the Arabian Sea are analyzed using different satellite observations, Simple Oceanic Data Assimilation (SODA) reanalysis, and Ocean Reanalysis System 4 (ORAS4) from 1993 to 2016 to investigate the impacts of Southwest (SW) Monsoon strength on Somali Current (SC) mesoscale circulations such as the Great Whirl (GW), the Socotra Eddy (SE), the Southern Gyre (SG), and smaller eddies. Increased Ekman pumping during stronger SW monsoons strengthens coastal upwelling along the Somali coast. The Arabian Sea basin-wide anticyclonic circulation and presence of the GW form mesoscale circulation patterns favourable to advection of upwelled waters eastward into the central Arabian Sea. In September, after the SW monsoon winds reach peak strength in July and August, a higher number of discrete anticyclonic eddies with higher ( > 20 cm) sea surface height anomalies develop in strong and normal intensity SW monsoon seasons than weaker SW monsoon seasons.     

The seasonal variation of the upper layers of the South China Sea (SCS) circulation and the Indonesian through flow (ITF): An ocean model study

The upper layer, wind-driven circulation of the South China Sea (SCS), its through-flow (SCSTF) and the Indonesian through flow (ITF) are simulated using a high resolution model, FVCOM (finite volume coastal ocean model) in a regional domain comprising the Maritime Continent. The regional model is embedded in the MIT global ocean general circulation model (ogcm) which provides surface forcing and boundary conditions of all the oceanographic variables at the lateral open boundaries in the Pacific and Indian oceans. A five decade long simulation is available from the MITgcm and we choose to investigate and compare the climatologies of two decades, 1960–1969 and 1990–1999.The seasonal variability of the wind-driven circulation produced by the monsoon system is realistically simulated. In the SCS the dominant driving force is the monsoon wind and the surface circulation reverses accordingly, with a net cyclonic tendency in winter and anticyclonic in summer. The SCS circulation in the 90s is weaker than in the 60s because of the weaker monsoon system in the 90s. In the upper 50 m the interaction between the SCSTF and ITF is very important. The southward ITF can be blocked by the SCSTF at the Makassar Strait during winter. In summer, part of the ITF feeds the SCSTF flowing into the SCS through the Karimata Strait. Differently from the SCS, the ITF is primarily controlled by the sea level difference between the western Pacific and eastern Indian Ocean. The ITF flow, consistently southwestward below the surface layer, is stronger in the 90s.The volume transports for winter, summer and yearly are estimated from the simulation through all the interocean straits. On the annual average, there is a ∼5.6 Sv of western Pacific water entering the SCS through the Luzon Strait and ∼1.4 Sv exiting through the Karimata Strait into the Java Sea. Also, ∼2 Sv of SCS water enters the Sulu Sea through the Mindoro Strait, while ∼2.9 Sv flow southwards through the Sibutu Strait merging into the ITF. The ITF inflow occurs through the Makassar Strait (up to ∼62%) and the Lifamatola Strait (∼38%). The annual average volume transport of the ITF inflow from the simulation is ∼15 Sv in the 60s and ∼16.6 Sv in the 90s, very close to the long term observations. The ITF outflow through the Lombok, Ombai and Timor straits is ∼16.8 Sv in the 60s and 18.9 Sv in the 90s, with the outflow greater by 1.7 Sv and 2.3 Sv respectively. The transport estimates of the simulation at all the straits are in rather good agreement with the observational estimates.We analyze the thermal structure of the domain in the 60s and 90s and assess the simulated temperature patterns against the SODA reanalysis product, with special focus on the shallow region of the SCS. The SODA dataset clearly shows that the yearly averaged temperatures of the 90s are overall warmer than those of the 60s in the surface, intermediate and some of the deep layers and the decadal differences (90s − 60s) indicate that the overall warming of the SCS interior is a local effect. In the simulation the warm trend from the 60s to the 90s in well reproduced in the surface layer. In particular, the simulated temperature profiles at two shallow sites at midway in the SCSTF agree rather well with the SODA profiles. However, the warming trend in the intermediate (deep) layers is not reproduced in the simulation. We find that this deficiency is mostly due to a deficiency in the initial temperature fields provide by the MITgcm.     

The East Asian Monsoon Simulation with IAP AGCMs-A Composite Study

Ｔｈｅ　Ｅａｓｔ　Ａｓｉａｎ　Ｍｏｎｓｏｏｎ　Ｓｉｍｕｌａｔｉｏｎ　ｗｉｔｈ　ＩＡＰ　ＡＧＣＭｓ－Ａ　Ｃｏｍｐｏｓｉｔｅ　ＳｔｕｄｙＷａｎｇＨｕｉｊｕｎａｎｄＢｉＸｕｎｑｉａｎｇ（ＩｎｓｔｉｔｕｔｅｏｆＡｔｍｏｓｐｈｅｒｉｃＰｈｙｓｉｃｓ（ＩＡＰ），...     

An analysis of extreme intraseasonal rainfall events during January–March 2010 over eastern China

The precipitation over eastern China during January–March 2010 exhibited a marked intraseasonal oscillation (ISO) and a dominant period of 10-60 days. There were two active intraseasonal rainfall periods. The physical mechanisms responsible for the onset of the two rainfall events were investigated using ERA-interim data. In the first ISO event, anomalous ascending motion was triggered by vertically integrated (1000–300 hPa) warm temperature advection. In addition to southerly anomalies on the intraseasonal (10–60-day) timescale, synoptic-scale southeasterly winds helped advect warm air from the South China Sea and western Pacific into the rainfall region. In the second ISO event, anomalous convection was triggered by a convectively unstable stratification, which was caused primarily by anomalous moisture advection in the lower troposphere (1000–850 hPa) from the Bay of Bengal and the Indo-China Peninsula. Both the intraseasonal and the synoptic winds contributed to the anomalous moisture advection. Therefore, the winter intraseasonal rainfall events over East Asia in winter could be affected not only by intraseasonal activities but also by higher frequency disturbances.     

Variability of rainfall and effective onset and length of the monsoon season over a sub-humid climatic environment

The study has analyzed the variability and trends in monthly, seasonal and annual rainfall and rainy days of four locations over different agro-ecological zones of Bihar, namely Samastipur (zone-I), Madhepura (zone-II), Sabour (zone-IIIA) and Patna (zone-IIIB). The Mann–Kendall nonparametric test was employed for detection of statistical significance and slopes of the trend lines were determined using the method of least square linear fitting. The variability and trends of onset of effective monsoon and length of monsoon period were also analyzed using the same method. The mean annual rainfall varies from 1137 mm at Patna to 1219 mm at Sabour. July is the rainiest month in all the zones followed by August. Maximum increase in annual rainfall was found at Sabour (40.1% of mean/30 years at 95% confidence level) and minimum for Patna (10.1% of mean/30 years). Significant increasing trend of rainfall during July, August and September at rates of 41.9, 83.2, and 112.7% of the mean/30 years, respectively has been noticed at Madhepura. Analysis of rainy days indicates that rainy days increased during winter and annually for all the sites. The mean effective onset of monsoon varies from 18th June at Sabour to 28th June at Patna. The trends in the date of effective onset of monsoon indicate that the date tends to be early in all the sites except Madhepura. But a significant delayed trend in the onset at a rate of 2.8% of the mean/30 years has been observed for Madhepura. The early trend of the effective onset of monsoon and increasing trends of length of monsoon season have been observed for Samastipur, Sabour and Patna.     

Analysis of quantitative precipitation forecasts using the Dynamic State Index

The German Weather Service (DWD) has two non-hydrostatic operational weather prediction models with different spatial resolution and precipitation parametrisations. The coarser COSMO-EU model has a spatial resolution of 7 km, whereas the higher-resolution COSMO-DE model has a gridspace of 2.8 km and explicitly resolves deep convection. To improve the numerical weather prediction (NWP) models it is necessary to understand precipitation processes. A central goal is the statistical evaluation of precipitation forecasts with dynamic parameters. Here, the Dynamic State Index (DSI) is used as a dynamic threshold parameter. The DSI theoretically describes the change of atmospheric flow fields as deviations from a stationary adiabatic solution of the primitive equations (Névir, 2004). For seasonal area means the DSI shows a remarkably high correlation with the precipitation forecasts provided by the COSMO-DE model. This is especially the case for the summer of 2007. The same analysis has been performed with the COSMO-EU forecast data and the results were compared with those from the COSMO-DE model. Moreover, an independent precipitation analysis, with a resolution corresponding to 7 km and 2.8 km, has been compared with respect to modelled precipitation and the DSI. In addition, correlations between the DSI and modelled as well as observed precipitation as a function of the forecast time for the different grid resolutions are also presented. The results show, that after 12 h, the correlation of the persistence forecast with the DSI reaches two thirds of the initial value. Thus, the DSI offers itself as a new dynamic forecast tool for precipitation events.     

Using the mean pressure gradient and NCEP/NCAR reanalysis to estimate the strength of the South Atlantic Anticyclone

A new methodology is proposed to estimate the strength of the South Atlantic Anticyclone (SAA), using the gridded sea level pressure (SLP) of the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) Reanalysis data. The top quartile (1017.3 hPa) of the SLP data was found a reasonable criterion to delimit the SAA area. Consequently, we defined the SAA area as the quadrangle containing 80% of the observations with pressure >1017.3 hPa. In this quadrangle, an area weighted pressure gradient (AWPG) was computed for the whole area and for the north–south and west–east halves. When compared with maximum pressure, the AWPG showed a better correlation with the significant wave height (SWH) and wind speed (WS) derived from altimetry. The mean value of the AWPG was 8 × 10⁻⁴ Pa/m, with representative values of 9.1 × 10⁻⁴ Pa/m and 7.4 × 10⁻⁴ Pa/m for austral winter and summer, respectively. The phase difference between the monthly AWPG in the north and south sub-quadrangles accounts for the evolution of the spatial pattern of the anticyclone throughout a year. This quantitative approach proved to be a useful estimate of the strength of South Atlantic Anticyclone. Further improvements of this approach are discussed.     

The 183-WSL fast rain rate retrieval algorithm: Part I: Retrieval design

The Water vapour Strong Lines at 183 GHz (183-WSL) fast retrieval method retrieves rain rates and classifies precipitation types for applications in nowcasting and weather monitoring. The retrieval scheme consists of two fast algorithms, over land and over ocean, that use the water vapour absorption lines at 183.31 GHz corresponding to the channels 3 (183.31 ± 1 GHz), 4 (183.31 ± 3 GHz) and 5 (183.31 ± 7 GHz) of the Advanced Microwave Sounding Unit module B (AMSU-B) and of the Microwave Humidity Sounder (MHS) flying on NOAA-15-18 and Metop-A satellite series, respectively.The method retrieves rain rates by exploiting the extinction of radiation due to rain drops following four subsequent steps. After ingesting the satellite data stream, the window channels at 89 and 150 GHz are used to compute scattering-based thresholds and the 183-WSLW module for rainfall area discrimination and precipitation type classification as stratiform or convective on the basis of the thresholds calculated for land/mixed and sea surfaces. The thresholds are based on the brightness temperature difference Δ_win = T_B89 ? T_B150 and are different over land (L) and over sea (S): cloud droplets and water vapour (Δ_win < 3 K L; Δ_win < 0 K S), stratiform rain (3 K < Δ_win < 10 K L; 0 K < Δ_win < 10 K S), and convective rain (Δ_win > 10 K L and S). The thresholds, initially empirically derived from observations, are corroborated by the simulations of the RTTOV radiative transfer model applied to 20000 ECMWF atmospheric profiles at midlatitudes and the use of data from the Nimrod radar network. A snow cover mask and a digital elevation model are used to eliminate false rain area attribution, especially over elevated terrain. A probability of detection logistic function is also applied in the transition region from no-rain to rain adjacent to the clouds to ensure continuity of the rainfall field. Finally, the last step is dedicated to the rain rate retrieval with the modules 183-WSLS (stratiform) and 183WSLC (convective), and the module 183-WSL for total rainfall intensity derivation.A comparison with rainfall retrievals from the Goddard Profiling (GPROF) TRMM 2A12 algorithm is done with good results on a stratiform and hurricane case studies. A comparison is also conducted with the MSG-based Precipitation Index (PI) and the Scattering Index (SI) for a convective-stratiform event showing good agreement with the 183-WSLC retrieval. A complete validation of the product is the subject of Part II of the paper.     

Mixing states of cloud interstitial particles between water-soluble and insoluble materials at Mt. Tateyama,Japan: Effects of meteorological conditions

Mixing states of cloud interstitial particles between water-soluble and insoluble materials apparently differ under various cloud-forming conditions. To study the mixing states of cloud interstitial particles, we made observations at Mt. Tateyama, Japan (2300 m a.s.l.) during June 2007 using fog (> 10 μm)-cut inlets. Number concentrations of dried particles (0.3–0.5 μm diameter) selected for less-grown (LG) particles (particles smaller than 0.56 μm diameter at 88% relative humidity) were used to quantify tendencies of the growth characteristics of cloud interstitial particles. Size-segregated soot mass concentrations (< 0.4 and < 1.1 μm) were also measured for cloud interstitial particles. Three samples of cloud interstitial LG particles at 88% RH were investigated for water-soluble and insoluble components using dialysis (extraction) of water-soluble materials with transmission electron microscopy (TEM). For one TEM sample with high fractions of the LG particles and high soot mass concentrations under high precipitation (2–6 mm/h), most particles (0.1–0.5 μm) were found to be water insoluble. More than half of the water-insoluble particles were considered to be soot particles showing chain aggregations of electron-opaque spherules. Regarding the other two TEM samples with low fractions of the LG particles under less intense precipitation (ca. 1 mm/h), most particles were partly water soluble. The scavenging process in the precipitating cloud can change the population of particles left behind, preferentially leaving insoluble particles according to cloud formation conditions.     

Severe weather reports and proximity to deep convection over Northern Argentina

Conventional surface data and quantitative estimations of precipitation are used to document the occurrence and spatial distribution of severe weather phenomena associated with deep moist convection over southeastern South America.Data used in this paper are 24-hour rainfall, maximum hourly gusts and present weather reports from the surface station network for Argentina to the north of 40°S and cover the period 2000–2005. Hourly rainfall estimated with the CMORPH technique (CPC MORPHing technique, R. J. Joyce et al., 2004) is included in the analysis in order to increase the density of the precipitation database from January 2003 to December 2005. Extreme events are detected by means of a 95th-percentile analysis of the 24-hour rainfall and wind; values greater than 30 mm and 25 m s^?1 respectively are considered extreme in the study area. These results are related to the presence of deep convection by considering the 235 K and 218 K cloud shield evolution in Geostationary Operational Environmental Satellite-12 Infrared (GOES-IR) imagery evaluated by the Forecasting and Tracking of Cloud Cluster (FORTRACC) technique. Rainfall above 30 mm day^?1 and present convection-related weather events tend to occur in the northeast of the country.Finally, an analysis is made of the relationship between severe phenomena and the location and lifecycle of Mesoscale Convective Systems (MCSs) defined by the 218 K or 235 K levels. According to the reports, favorable locations for severe weather concentrate to the northeast of the cloud shield anvil centroid although most of the cases are found in the northwest. This feature can be seen in systems with anvil areas larger than 250,000 km² in association to the predominant mid-level wind shear direction from the northwest over the area. Moreover, systems with centers located north of 30°S present a more circular shape while those to the south are more elongated with a NW–SE main axis clearly related to the presence and interaction with frontal zones over the area. Most of the events occur previous to the moment when the systems reach their maximum extension, between 2 and 10 h after the initiation of the system depending on the size of the MCSs.     

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.     

Circulation analysis in the northwest Indian Ocean using ARGO floats and surface drifter observations,and SODA reanalysis output

This study incorporates observations from Array of Real-time Geostrophic Oceanography (ARGO) floats and surface drifters to identify seasonal circulation patterns at the surface, 1000 m, 1500 m, and 2000 m in the northwest Indian Ocean, and quantify velocities associated with them. A skill comparison of the Simple Ocean Data Assimilation (SODA) reanalysis output was also performed to contribute to the understanding of the circulation dynamics in this region.Subsurface currents were quantified and validated using the ARGO float data. Surface currents were identified using surface drifter data and compared to the subsurface observations to enhance our previous understanding of surface circulations. Quantified Southwest Monsoon surface currents include the Somali Current (v_max = 179.5 cm/s), the East Arabian Current (v_max = 52.3 cm/s), and the Southwest Monsoon Current (v_max = 51.2 cm/s). Northeastward flow along the Somali coast is also observed at 1000 m (v_max = 26.1 cm/s) and 1500 m (v_max = 12.7 cm/s). Currents associated with the Great Whirl are observed at the surface (v_max = 161.4 cm/s) and at 1000 m (v_max = 16.2 cm/s). In contrast to previous studies, both ARGO and surface drifter data show the Great Whirl can form as early as the boreal Spring intermonsoon, lasting until the boreal Fall intermonsoon. The Arabian Sea exhibits eastward/southeastward flow at the surface, 1000 m, 1500 m, and 2000 m. Quantified Northeast Monsoon surface currents include the Somali Current (v_max = 97.3 cm/s), Northeast Monsoon Current (v_max = 30.0 cm/s), and the North Equatorial Current (v_max = 28.5 cm/s). Southwestward flow along the Somali coast extends as deep as 1500 m.Point-by-point vector and scalar correlations of SODA output to ARGO and surface drifter data showed that surface SODA output and surface drifter data generally produced a strong correlation attributed to surface currents strongly controlled by the monsoons, while subsurface correlations of SODA output and ARGO were mostly insignificant due to variability associated with intermonsoonal transitions. SODA output produced overall smaller velocities than both observational datasets. Assimilating ARGO velocities into the SODA reanalysis could improve subsurface velocity assimilation, especially during the boreal fall and spring when ARGO observations suggest that flow is highly variable.     

Effect of Horizontal Resolution on the Representation of the Global Monsoon Annual Cycle in AGCMs

The sensitivity of the representation of the global monsoon annual cycle to horizontal resolution is compared in three AGCMs:the Met Office Unified Model-Global Atmosphere 3.0;the Meteorological Research Institute AGCM3;and the Global High Resolution AGCM from the Geophysical Fluid Dynamics Laboratory.For each model,we use two horizontal resolution configurations for the period 1998–2008.Increasing resolution consistently improves simulated precipitation and low-level circulation of the annual mean and the first two annual cycle modes,as measured by the pattern correlation coefficient and equitable threat score.Improvements in simulating the summer monsoon onset and withdrawal are region-dependent.No consistent response to resolution is found in simulating summer monsoon retreat.Regionally,increased resolution reduces the positive bias in simulated annual mean precipitation,the two annual-cycle modes over the West African monsoon and Northwestern Pacific monsoon.An overestimation of the solstitial mode and an underestimation of the equinoctial asymmetric mode of the East Asian monsoon are reduced in all high-resolution configurations.Systematic errors exist in lower-resolution models for simulating the onset and withdrawal of the summer monsoon.Higher resolution models consistently improve the early summer monsoon onset over East Asia and West Africa,but substantial differences exist in the responses over the Indian monsoon region,where biases differ across the three low-resolution AGCMs.This study demonstrates the importance of a multi-model comparison when examining the added value of resolution and the importance of model physical parameterizations for simulation of the Indian monsoon.     

北半球平流层环流与亚洲夏季风关系的初步分析

季风是大气环流季节变化的某些特性在局部地区的特殊表现。人们逐渐认识到,大气环流的季节变化是全球性的现象。但长期以来关于季风的研究一般只限于对流层。而对流层与平流层大气之间存在着一定的联系,平流层的季节变化往往早于对流层。本文对本世纪70年代亚洲夏季风强年和弱年的平流层环流及其季节变化作了对比分析,指出:亚洲夏季风强弱与平流层环流有密切联系。弱季风年前冬(11—4月)50 hPa绕极地的西风环流偏强,季节变化开始得晚,但完成得早。50与     

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.     

A numerical comparison study of cloud seeding by silver iodide and liquid carbon dioxide

A comparison study on dynamic and microphysical effects of cloud seeding by silver iodide (AgI) and liquid carbon dioxide (liquid CO₂) was made using a 3D cloud model with seeding processes. The model was initialized based on the rawinsonde sounding taken from Pinliang station located in the western China on 20 April 2001. The sounding air reflects moist and stable characteristics at middle and low layers. The model results show that the seeding by liquid CO₂ and AgI at − 15 to − 20 °C levels of cloud has almost the same dynamic effect on the simulated clouds. The seeding is able to induce the formation of weak convective cells in both seeded and unseeded regions due to latent heat released by the transformation from liquid saturation to ice saturation. However, the initial seeding conducted by liquid CO₂ in the region of maximum supercooled water with temperature of 0 to − 5 °C enable to produce much stronger dynamic effect and precipitation by forming many convective new cells at low levels in the later stage of seeded clouds. The accumulated precipitation at the surface can be increased and redistributed, and more concentrated in the downstream region of seeded clouds.     

Interannual spring Wyrtki jet variability and its regional impacts

The role of spring Wyrtki jets in modulating the equatorial Indian Ocean and the regional climate is an unexplored problem. The source of interannual variability in the spring Wyrtki jets is explored in this study. The relationship between intraseasonal and interannual variability from 1958 to 2008 and its relation with Indian Summer Monsoon is further addressed. Analysis reveals that the interannual variability in spring Wyrtki jets is controlled significantly by their intraseasonal variations. These are mostly defined by a single intraseasonal event of duration 20 days or more which either strengthens or weakens the seasonal mean jet depending on its phase. The strong spring jets are driven by such intraseasonal westerly wind bursts lasting for 20-days or more, whereas the weak jets are driven by weaker intraseasonal westerlies. During the years of strong jets, the conventional westward phase propagation of Wyrtki jets is absent and instead there is an eastward phase propagation indicating the possible role of Madden Julian Oscillation (MJO) in strengthening the spring Wyrtki jets. These strong intraseasonal westerly wind bursts with eastward phase propagation during strong years are observed mainly in late spring and have implications on June precipitation over the Indian and adjoining land mass. Anomalously strong eastward jets accumulate warm water in the eastern equatorial Indian Ocean (EIO), leading to anomalous positive upper ocean heat content and supporting more local convection in the east. This induces subsidence over the Indian landmass and alters monsoon rainfall by modulating monsoon Hadley circulation. In case of weak current years such warm anomalies are absent over the eastern EIO. Variations in the jet strength are found to have strong impact on sea level anomalies, heat content, salinity and sea surface temperature over the equatorial and north Indian Ocean making it a potentially important player in the north Indian Ocean climate variability.     

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.     

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.