1. Trang chủ
  2. » Khoa Học Tự Nhiên

AIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGION

135 635 0
Tài liệu đã được kiểm tra trùng lặp

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

THÔNG TIN TÀI LIỆU

Thông tin cơ bản

Định dạng
Số trang 135
Dung lượng 9,26 MB

Nội dung

AIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGIONAIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGIONAIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGIONAIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGIONAIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGIONAIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGION

Đại học Quốc Gia Hà Nội VIET NAM NATIONAL UNIVERSITY, HANOI Đinh Văn Ưu Dinh Van Uu Tương tác biển-khí Khu vực nhiệt đới ấn độ Thái bình dương AIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGION Nhµ xuÊt Đại học Quốc Gia Hà Nội - 2003 VIET NAM NATIONAL UNIVERSITY PRESS - 2003 VIET NAM NATIONAL UNIVERSITY, HANOI Dinh Van Uu Table of contents INTRODUCTION Chapter I SMALL SCALE AIR-SEA INTERACTION 10 AIR- SEA INTERACTION IN THE INDIAN-PACIFIC TROPICAL REGION 1.1 Introduction 10 1.2 Surface Processes 10 1.2.1 Equation of motion with Viscosity 12 1.2.2 Turbulence in boundary layer 13 1.2.3 Turbulent Stresses: The Reynolds Stress 14 1.2.4 Thickness of the lower atmospheric boundary layer 15 1.2.5 Vertical profile of wind speed within the surface layer 16 1.2.6 Characteristics of the lower atmospheric boundary layer on the sea 18 1.3 Influence of the stratification and wind wave on the turbulent characteristics of the lower atmospheric boundary layer 21 1.3.1 Stratification and its influence on the turbulent characteristics of the lower atmospheric boundary layer 21 1.3.2 Dependence of the sea surface roughness on wind wave characteristics 23 1.4 VIET NAM NATIONAL UNIVERSITY PRESS - 2003 Calculation of turbulent fluxes on sea surface 25 1.4.1 Direct Calculation of Fluxes 25 1.4.2 Gradient method 26 1.4.3 Indirect Calculation of Fluxes: Bulk Formulas 28 1.4.4 Field and model estimates of the sea surface drag 30 1.4.5 Climatic estimation of the turbulent fluxes 32 1.4.6 Calculation fluxes in the storm condition 35 2.5.1 Introduction 59 Chapter II OCEAN’S RESPONSE TO THE ATMOSPHERE 36 2.5.2 The Oceanic Mixed Layer and Thermocline 59 Dynamic interaction and Ekman layer 37 2.5.3 Geographical Distribution of Surface Temperature and Salinity 63 2.1.1 Response of the upper ocean to winds and inertial motion 37 2.5.4 Thermal transformation of the lower atmospheric boundary layer on 2.1.2 Wind circulation and Ekman Layer at the Sea Surface 38 2.1.3 Langmuir Circulation 41 2.1.4 Influence of stability in the Ekman Layer 42 2.1.5 Ekman Mass Transports 42 2.1 2.4.4 2.5 Meridional Heat Transport 57 Thermodynamic interaction of air-sea boundary layers and the upper ocean active layer 59 the sea -67 3Chapter III AIR SEA INTERACION IN THE INDIAN-PACIFIC TROPICAL ZONE 69 3.1 Introduction 69 3.2 Thermal Land-Ocean-Atmosphere interaction and general atmospheric and ocean circulation 70 2.1.6 2.2 Application of Ekman Theory 43 2.2.1 Boundary layers model with constant turbulent viscosity 45 2.2.2 Boundary layers model with variable turbulent viscosity 47 2.3 3.2.1 Importance of the Ocean in Earth’s Heat Budget 71 3.2.2 The development of a thermal circulation 71 3.2.3 The Coriolis effect (deflection caused by Earth’s rotation) 74 3.2.4 Sea level pressure and wind 75 3.2.5 Intertropical Convergence Zone (ITCZ) 76 3.2.6 Atmospheric circulation in tropical zone 78 3.2.7 General ocean circulation 81 Hydrodynamic interaction model of air-sea boundary layers 45 Wind Waves as the ocean’s response to atmosphere 48 2.3.1 Introduction 48 2.3.2 Linear Theory of Ocean Surface Waves 49 2.3.3 Waves and the Concept of a Wave Spectrum 51 2.3.4 Generation of Waves by Wind 52 2.3.5 Influence of atmospheric stratification on wave development 52 2.3.6 Wave Forecasting 53 3.3 2.4 Scales of variability in the ocean-atmosphere system 91 3.3.1 Diurnal and synoptic scales 92 3.3.2 Seasonal to interannual time scales 95 3.3.3 Interannual cycle 96 The Oceanic Heat Budget 54 3.4 2.4.1 Introduction 54 2.4.2 Factors influencing Heat-Budget Terms 54 2.4.3 Geographic Distribution of Terms in the Heat Budget 56 Variability of the monsoon system 98 3.4.1 Monsoon activity in the tropical ocean-atmosphere system 98 3.4.2 Variability of the Asian-Australian monsoon system 101 3.5 1051 Interannual climate variability in the Indian-Pacific tropical zone -100 3.5.1 ENSO activity and variability 105 3.5.2 The Quasi-Biennial zonal wind Oscillation (QBO) 112 3.5.3 Decadal variability in the Pacific and Indian Oceans 116 3.5.4 The variability in the Bien Dong (South China) Sea 121 INTRODUCTION REFERENCES 128 Air and water are the two fluids we know most about They are both essential to the maintenance of our lives, providing a hospitable environment for all living things We have a direct, lifelong experience in observing how they behave and how we utilize them In the natural environment we can see, feel or hear examples of almost all the kinds of fluid flows that we will study in this book The atmosphere is a layer of gas held to the surface of the earth by gravitational attraction Most of the mass of the atmosphere is confined to the first 15 kilometers above sea level, yet the small amount above this level is responsible for filtering out the deadly high energy radiation from the sun which would otherwise destroy life The interaction of the atmosphere with sunlight helps to maintain the earth's surface temperature above that of an airless planet, like the moon (This increase in temperature, called the greenhouse effect, threatens to grow in the future because of anthropogenic emissions of heat absorbing gases.) The motion of the atmosphere that we observe around us is driven by the diurnal pattern of heating by the sun and cooling by radiation to outer space Over the year, these heating/cooling patterns shift to different latitudes, giving rise to annual climate variations An intimate part of this process is the evaporation of water from the earth's surface, the formation of clouds from this water at high elevations and its precipitation back to the surface This distillation of ocean water, moved to land by winds, provides the sweet water that maintains terrestrial life Local weather provides a variety of wind motions Sometimes the wind speed is quite small, especially at night when radiative cooling stabilizes the atmosphere But storms driven by precipitation of water vapor, such as thunderstorms and hurricanes, can have very high wind speeds Cold air is more dense than warm air, so that a cold air mass tends to flow under a warm air mass, forming a cold front Large scale weather patterns drift past our locality, bringing changes that are not greatly affected by local conditions As we can see from watching the daily television weathercast, the main features of the weather pattern extend over many thousands of kilometers, a distance hundreds of times larger the atmospheric height Yet the changes of pressure, temperature and humidity are much greater in the vertical direction than in the horizontal direction, despite the much greater horizontal size of a weather pattern The pull of gravity is so strong over large distances that it forces the atmosphere to flow mostly in the horizontal direction Because we are so small compared to the vertical and horizontal dimensions of the atmosphere, we can observe how the wind blows in a tiny portion of the atmosphere that is nearby us It is noticeable that the wind speed and direction are somewhat variable, especially over time intervals of less than a minute These changes are much more rapid than the changes accompanying a weather pattern, which may take a day to pass us by Fluid flows that exhibit variability over time and length scales which are small compared to that of the overall flow are called turbulent flows The atmospheric wind is a turbulent flow The atmospheric motion is responsible for diluting air pollutants, such as those emitted by power plants and automobiles When these pollutant streams are marked by smoke, we can readily observe how the smoke intensity decreases as the wind turbulence mixes the pollutant stream with clean air, diluting the strength of the pollutant within the plume (or, if you prefer, dirtying more and more of the atmosphere) Most of these pollutants mix no higher than a few kilometers in the atmosphere , and are eventually carried far downwind and deposited back to the earth's surface Some, however, not soon return to the earth and instead mix gradually throughout the entire atmosphere, including the stratosphere Some of these gases lead to the destruction of stratospheric ozone and to increased average surface temperature The mixing properties of the atmosphere are extremely important in determining the degree of atmospheric pollution in urban areas The water in the ocean, lakes and rivers, as well as that in underground aquifers, is called the hydrosphere The volume of fresh water on the continents is small compared to the oceanic volume, but it is the part of the hydrosphere that is most important to the maintenance of terrestrial life The management and use of this water for agricultural and other purposes form an important branch of engineering But the ocean is important too It provides the source of precipitation over the land and, in its surface layer, an environment for the 10 growth of microscopic plants and animals that form the base of the oceanic food chain The ocean tends to make the climate more uniform in the latitudinal direction by moving warm tropical waters towards the poles and displacing cold polar water towards the equator We are all familiar with the downhill motion of streams and rivers, flowing toward the sea The energy in this flow can be tapped by building dams that force the river flow through turbines to generate electric power Sometimes this energy is dramatically dissipated as the river plunges over a precipice to form a waterfall The violently turbulent motion at the base of the falls converts the river's directed energy into heat When the river reaches the sea, its fresh water, being lighter than the sea water, floats on top of the sea, sometimes many miles beyond the river's mouth before it is diluted with sea water Most of the fresh water on the continents is out of sight, below ground It exists in the pores between mineral deposits and is fed by precipitation that percolates through the ground under the influence of gravity The fluid velocity in the underground aquifers is much lower than it is in rivers, the water being impeded by the frictional force of the porous medium through which it flows Underground water is often the source of potable water Locating underground water and pumping it from the ground for human use is limited by the characteristics of the underground aquifer Preserving the purity of this water from the contamination by toxic fluids buried or dumped on the surface and subsequently traveling down into the aquifer is a major problem worldwide At the edge of the ocean we see the ocean waves crashing on the shore The waves carry to the shore energy generated by 11 the wind blowing over the ocean surface Of course, the ocean surface doesn't move (on the average) in the direction of the waves, but it oscillates as the wave passes by Ocean waves are called gravity waves because the pull of gravity at the air-sea interface is responsible for the propagation of these waves, which not penetrate far below the ocean surface The other oceanic motion we notice at the sea shore is the tidal rise and fall of the sea surface This twice-a-day cycle is caused by the difference in gravitational pull of the moon (and to a lesser extent, the sun) on opposite sides of the earth The differential gravity force gives rise to a bulging of the ocean surface in the direction of the moon, which passes a given location twice in the lunar day of 25 hours The tidal motion, consisting of both a vertical and horizontal oscillation, may be amplified greatly along the continental coastline, sometimes by a factor of ten above the general oceanic amplitudes Oceans may contain localized currents, like the Gulf stream, that are mighty rivers flowing across a nearly stationary ocean Earthquakes can generate tsunami waves that travel many thousands of kilometers before crashing ashore, sometimes wreaking devastation on lowlying coastlines Even hurricanes can generate storm surge waves that cause coastal flooding The rise and fall of the tide can be utilized to produce mechanical power, but at the present time this is seldom economical compared to river power Mechanisms have been devised to extract power from ocean waves, but this has also proved to be uneconomical But the forces exerted on ships and wave barriers by the ocean waves can be very substantial, and protecting them against such forces very expensive Knowledge 12 of the dynamics of the ocean is important to many of mankind's pursuits Air-sea interaction provides a comprehensive account of how the atmosphere and the ocean interact to control the global climate, what physical laws govern this interaction, and its prominent mechanisms In recent years air-sea interaction has emerged as a subject in its own right, encompassing small-scale and large-scale processes in both air and sea By developing its subject from basic physical (thermodynamic) principles, the book is broadly accessible to a wide audience It is mainly directed towards graduate students and research scientists in meteorology, oceanography, and environmental engineering The book will be of value on entry level courses in meteorology and oceanography, and also to the broader physics community interested in the treatment of transfer laws, and thermodynamics of the atmosphere and ocean Current research in air-sea interaction is attempting to apply our understanding of molecular and turbulence scale processes to global scale phenomena There are many unknowns in the fundamental processes that control the global climate For example, what processes control the amount and distribution of water in the atmosphere; what is the effect of cloud variability on the sea surface temperature? How changes in the ocean circulation affect the atmospheric circulation and, hence, the distribution of wind stress, temperature and precipitation, and how does this feed back to the ocean? These, and many other related questions, encurage some of the studies of the interaction of the atmosphere and the ocean Major uncertainties remain in our understanding of the fundamental processes of air-sea interaction, particularly, in heterogeneous and nonequilibrium conditions; for example, we 13 not know enough about the relationships between the directional wave spectra, surface fluxes and the properties of the oceanic and atmospheric boundary layers to develop satisfactory predictive models The ocean and atmosphere are interdependent, or coupled, because of the dependence of the atmosphere on heat and moisture at the sea surface and the dependence of the ocean circulation on the wind Studies rarely combine investigations of both environments to determine the extent of the fedbacks between the ocean and the atmosphere Understanding the coupled ocean-atmosphere system depends largely on the scales of interaction between the two fluids and the processes that provide the strongest feedbacks Advances are most likely through multidisciplinary process studies that connect the upper ocean and lower atmosphere For global studies more comprehensive parameterizations of the surface processes are required as well as improvements in satellite retrievals and assimilation in numerical models The high spatial and temporal variability of surface processes needs to be properly considered In situ measurements are revealing very complex horizontally and vertically heterogeneous fields that cannot be resolved by current remote sensing techniques High resolution models, which include the physics of the processes that contribute to this variability, combined with satellite data seems to be the best tool for global analysis and prediction The large-scale dynamics of the ocean and the atmosphere are closely related Energy is transferred from the atmosphere to the ocean surface mixed layer driving the circulation of the upper ocean In turn, energy from the ocean is fed back to the atmosphere affecting the atmospheric circulation, the weather 14 and the climate The concept is deceptively simple, but as we explore the coupled earth system, we are frequently limited by our lack of understanding of the interchanges between the atmosphere and ocean Kraus and Businger (1994) highlighted several areas that continue to require more attention, for example, the interaction of the wind and surface waves, the parameterization of subgrid scale processes in large-scale circulation models, and the transfer of gases across the air-sea interface It has been convenient to divide air-sea interaction studies into two categories: small- and large- scale ocean-atmosphere interactions However, this often belies the fundamental precept that the basis of the interaction of the atmosphere and the ocean is the exchange of matter and energy across a material interface - the sea surface An exchange that, for the most part, occurs on molecular scales, involving both turbulent and laminar processes modified by wave breaking, surface tension, the structure of the planetary boundary layer and the ocean mixed layer and other effects A satisfactory understanding of these processes remains elusive, but is essential if we are to address adequately the larger scale oceanatmosphere problems In the past many proponents of large-scale studies, such as global climate and ocean circulation relied heavily on the veracity of the parameterization of small-scale air-sea exchange processes, often overlooking the uncertainties in the basic measurements and their interpretation More recently, we have recognized the importance of connecting small-scale process studies, investigating the exchange of heat, moisture, momentum and trace constituents across the air-sea interface, 15 with the large-scale problems of global climate change and ocean circulation that rely heavily on numerical models and highly averaged fields Studies, such as the Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment (TOGA COARE) and the World Ocean Circulation Experiment (WOCE), are leading the way by highlighting the importance of process studies for a satisfactory understanding of global climate and ocean general circulation problems The small-scale exchange processes are generally related to the global-scale problems via parameterizations of the fluxes that use mean quantities obtained by measuring on various platforms such as buoys, ships and satellites These parameterizations are also used extensively in operational meteorological models as well as many research general circulation models of the coupled ocean-atmosphere system Large uncertainties exist in the derivations of the bulk parameterizations due to the difficulty of measuring surface fluxes directly, and the difficulty of applying these measurements to scales greater than a few hundred kilometers and several hours The problems are particularly acute in the tropics where low wind speeds and very high sea surface temperatures result in primarily buoyancy-driven fluxes that are not well parameterized by most prevailing methods, and in coastal regions where fetch, topography and water depth vary considerably While there is resistance to the establishment of canonical parameterizations of the fluxes, the TOGA COARE research community recognized that such an approach has the advantage of focusing the attention of a large group of researchers to exchange and compare data and rapidly transfer 16 this information to a broader community whose interests require surface fluxes Despite these efforts, it is unlikely that in this decade of large scale climate and ocean circulation studies we will know the surface fluxes as well as enough we would like A continued effort will be required to improve our knowledge of heat, moisture, momentum and trace constituent fluxes, to increase our understanding of the uncertainties in the fluxes we measure, and to improve our knowledge of the relationship between boundary layer processes on each side of the interface to the surface exchange mechanisms research to the global and Indian-Pacific regional issues of atmosphereocean climate and its variations The Indian- Pacific tropical zone is most important heat storage of the World Ocean with the highest value of the mean sea surface temperature (SST); there is the Western Pacific/ Eastern Indian or AsianAustralian warm pool The formation of the warm pool and its variation is the favourable condition for typhoon formation and development in the western part of the Pacific Ocean The variability of the global-scale ocean and atmosphere circulation including Trade’s wind, the Walker circulation, the Asian-Australian monsoon system, the El Niño/Southern Oscillation (ENSO) and sea surface temperature (SST) is related to the oscillations of the ocean-atmosphere system This course will acquentance physical oceanography and meteorology students with the one-dimensional theories of turbulent boundary layers, to observations of the planetary boundary layer of the atmosphere, and to observations of the surface mixed layer of the ocean Coupled one-dimensional models of air-sea interaction will be studied with a view towards understanding the importance of interactions of the turbulent boundary layers with each other, and with the interior of their respective fluids Finally, we will review current progress in our understanding of the surface processes, and the application of this 17 18 2.2 Chaptre SMALL SCALE AIR-SEA INTERACTION 2.1 Introduction This part of the course will introduce physical oceanography and meteorology students to the one-dimensional theories of turbulent boundary layers, to observations of the planetary boundary layer of the atmosphere, and to observations of the surface mixed layer of the ocean Necessary concepts of turbulence theory will be studied, along with the probability and statistical tools that are appropriate Coupled one-dimensional models of air-sea interaction will be studied with a view towards understanding the importance of interactions of the turbulent boundary layers with each other, and with the interior of their respective fluids Finally, the subject of the influence of horizontal variability will be opened, with the objective of appreciating the role that ocean and atmosphere dynamics play in modulating the thermodynamics of the turbulent boundary layers, and as an introduction to the associated part of the course on Large-Scale Ocean-Atmosphere Interactions: Air-sea interaction and tropical meteorology and oceanography 19 Surface Processes Interactions between the ocean and atmosphere occur at the air-sea interface The ocean surface is a material interface that is a barrier to the exchange of heat, moisture, momentum and trace constituents Away from the surface both fluids are usually in turbulent motion, but near the interface turbulence is suppressed and the transport is controlled primarily by molecular processes To quantify the exchanges at the interface it is necessary to understand how the turbulent layers of the ocean and the atmosphere are connected via the molecular sublayers in either side of the sea surface In turn, we need to understand how the turbulent layers transport the properties of the interface into the interior of these fluids, the extent to which changes in the interior feed back to the interface, and how processes at the surface affect the structure of the deep ocean and free atmosphere The fundamental processes that connect the atmosphere and the ocean are the energy input to the ocean by the wind, the net freshwater flux, expressed primarily as precipitation minus evaporation, and the net surface heat flux As Charnock (1951) pointed out the energy transmitted by the wind to the ocean is a tiny fraction of the radiation received at the surface, yet wind-driven currents largely determine the regions where the ocean energy is fed back into the atmosphere that sets the pattern of cloudiness, which in turn determines the radiation input The ocean-atmosphere system is intrinsically coupled, although feedbacks across the air-sea interface are often masked by temporal and spatial differences As in the example above, we may understand the processes that connect the ocean and the atmosphere, but we not understand well enough the distribution of energy within the 20 This schema illustrated two equations for ice covrage in the Northern Ice Ocean (I) and heat storage of warm current branch in north Atlantic (Q): dI   m(Q  Q0 ) dt (3.9) dQ  n( I  I ) dt Figure 3.34 Schema of circulation system for Atlantic Northern Ice Oceans (from Doronin, 1980): The period of this cycle could be estimated according 2 following formula:   , its approximated value is about 3.5 mn years  - Florida warm water,  Gulf Stream warm current branch, H – New Foundland,  Interaction of the tropics and subtropics is another potentially important means by which the upper ocean processes influence climate variability, possibly including the decadal modulation of ENSO - north Atlantic warm current,  - Labrador cold current,  deep water circulation in the Northern Ice Ocean The acceleration of water circulation is conduced to decrease the remaining time of water in the warm tropical region before reaching Gulf Stream thus the water temperature is not reached its normal value In this case the temperature difference between west and east regions of the ocean will be decreased too, wind decreases and the remaining time of water in the warm tropical region begins increase The water reaching Gulf Stream has temperature warmer than its normal value then the east-west temperature gradient and its anomaly increase that increases meridional wind and new cycle will be developed again The period of this cycle is changed from years to years 241 4.5.4 The variability in the Bien Dong (South China) Sea The western Pacific/South East Asia region has an important influence on global climate because the heating of atmosphere associate with the heavy rains of the AsianAustralian monsoon system drives mean and variable globalscale Ocean and Atmosphere Circulation The Bien Dong (South China) Sea with the maximum depth greater than 5000 m is the largest marginal sea in Western Pacific This region is most important heat storage of the World Ocean with the highest value of the Mean Sea Surface Temperature (SST); there is the Western Pacific/ Eastern Indian or AsianAustralian Warm pool (figure 3.35) The mechanism of formation of the warm pool and its variation is very important in order to understand processes 242 air-sea interaction in the South East Asian monsoon region and global climate variability The high heat storage of the SCS and of the adjacent Indonesian Seas is the favourable condition for typhoon formation and development in the western part of the Pacific Ocean The convergence of moisture supplied by the surrounding oceans in the InterTropical Convergence Zone (ITCZ) gives the diabatic heat sources, migrating seasonally across the equator with reversing monsoon The Asian-Australian monsoon system is one of the most important circulation systems in the general circulation of the global atmosphere Interannual variation in the intensity and location of these heat sources are responsible for the quasibiennial oscillation and ENSO This variation is associated with the intensity and location of atmospheric cells of the Walker and Hadley circulation in the South East Asian region The water circulation in marginal seas such as the Bien Dong Sea and the complex straits and seas of the Maritime Continent allow a strongly variable flow from the Pacific to Indian Ocean: Pacific Indian Throughflow The water circulation of the South China Sea (SCS) is principally determined by reversing monsoon winds The seasonal circulation pattern in the SCS and its adjacent seas has been investigated by several researchers such as Wyrtki, 1961, Pohlmann T., 1987, Ping-Tung Shaw and Shenn-Yu Chao, 1994, etc They have pointed out that the SCSC circulation is related to water exchanges between the SCS and the Kurioshio through the Louzon Strait, between the SCS and Indonesian seas through the Mindoro Strait (to Sulu Sea) and Karimata Strait (to Java Sea) and between the SCS and the East China Sea through Taiwan strait 243 The variability of the marginal circulation has not, in general, been adequately observed with respect to interannual variability, and in some areas with respect to the seasonal cycle The air-sea flux fields of sufficient accuracy are not available and the ocean thermohaline structure observations are too sparse to provide adequate spiritual coverage for mapping without appeal to models for either statistical or dynamical interpolation and such models may not be accurate enough The bifurcation latitude of the North Equatorial Current (NEC) at the Philippine coast is affected by the amount of transport around the Philippines For the ocean-scale circulation model the absence of the analysis thermohaline structure information and the absence of the Pacific to Indian Ocean throughflow may seriously distort the model boundary currents The marginal sea circulation models have become more realistic; they useful perspectives on certain aspects of the real world Figure 3.35 Distribution of the mean sea surface temperature in the South-eastern Asian seas (Youhai He and Cuihua Guan, 1998) 244 The variation of the wind direction is very important not only by reversing monsoon wind from one season to another (fig.3.36) but there is also strong variation from one region to another during one season (Helleman & Rosenstein 1983) During the SW monsoon the variation of the wind directions and velocities is very important because the position of the ITCZ is more variable and the wind come from Indochina continent is influenced by orographic transformation due to local conditions (relief, coastal line ) The SW monsoon season is also season of typhoon and tropical cyclone with strong perturbation of the wind field in the sea and in the adjacent regions There are two seasons of the precipitation in this region with about 70 % of the total amount accumulated in the SW monsoon season Along the Vietnamese coastal zone the precipitation maximum and its time have a high correlation with the position of the ITCZ Beside of the wind and the precipitation-evaporation variation, the air- sea flux of the heat is strong variable in the space and in the time The seasonal variation of the air temperature for north-west region is very strong in compare with the southern region of the sea and another tropical region During the NE monsoon the air temperature difference between the north and the south region of the sea is very big, its difference my reaches about 15-20C Due to the low air temperature and strong wind, the total amount of the heat fluxes in the north-west coastal region is about 150 - 200 wt/m2 from sea surface to the atmosphere At the same time, in the south-east region, the sea water always receives the heat fluxes from atmosphere The Variational Inverse Method (VIM) was developed in the GeoHydrodynamic and Environment (GHER), Liege University and applied to make Mediterranean Ocean Data Base (MODB) For the South China Sea we used VIM with the grid size 15’x15’ longitude- latitude and 30 depth levels In the figures 3.37 showed the analyzed surface temperature and salinity distributions for two monsoon seasons using VIM Figure 3.36 Wind stress on the Bien Dong Sea during winter (top) and summer (bottom) from Uu, 2001 245 246 The Asian-Australian monsoon is one of the most remarkable phenomena of the ocean-atmosphere-land system in response to a seasonal march of solar insolation with complex feedback mechanisms The exchange fluxes in the ocean-atmosphere-land interfaces are dependent of atmospheric and oceanic variables such as land and sea surface temperature (SST), winds, etc In particular, the SST is one of the most important variables in the prediction of monsoon and marine system Interannual variation in the intensity and location of the heat sources are responsible for the quasibiennial oscillation (QBO) and ENSO This variation is associated with the intensity and location of atmospheric cells of the Walker and Hadley circulation in the equatorial zone, especially in the South East Asian region The variations of the SST in the region of Indonesian seas may influence in interannual variability of the monsoon Capturing the physics of the interactions between the monsoon diabatic heat sources and the warm pool heat sources in numerical models of the coupled oceanatmosphere-land system is critical for improve ENSO and monsoon forecasts Figures 3.37 Analysed surface layer temperature fields during winter (top) and summer (bottom) seasons from Uu D.V and Brankart, 1997 247 The shape of the warm pool along the boundary is influenced by the surface western Pacific boundary currents In particular, the bifurcation of the North Equatorial Current near 14N into the Mindanao Current and the Kurioshio appears to be subject to considerable interannual variability related to the QBO and ENSO (Lukas et al., 1996) The simulation results (Metzger and Hurburt, 1996) show that when the Sulu archipelago is opened, a net cyclonic flow develops around the Philippines, which is essentially an extension of the northern Pacific tropical gyre The bifurcation latitude of the north Equatorial Current at the Philippines coast is also affected by the amount of transport through the Sulu archipelago 248 The water circulation of the South China (SCS) and adjacent seas is principally determined by reversing monsoon winds (fig 3.38) The seasonal circulation pattern in the SCS and its adjacent seas has been investigated by several researchers They have pointed out that the SCS circulation is related to water exchanges between the SCS and the Kurioshio through the Luzon Strait, between the SCS and Indonesian seas through the Mindoro Strait (to Sulu Sea) and Karimata Strait (to Java Sea) and between the SCS and the East China Sea through Taiwan strait Intrusion of the Kurioshio to South China Sea (Faris and Wimbush, 1996) and its variation makes the complexity of the basin-scale circulation in the SCS and adjacent seas The results of the three-dimensional thermohaline and circulation modelling of the SCS show the possibility to estimate the water fluxes exchanged between Pacific and Indian Oceans via Indonesian seas and the interaction between SCS and Kuroshio Current (KC) system There are two principal branches of the circulation in the SCS: one is permanently passing by eastern part of the sea from Luzon strait to Sulu Sea and another passing western part of the sea in direction to Java Sea in winter and from Java Sea in summer The simulated permanent loop current in the north-eastern region of the SCS shows the high intensity of the interaction between SCS and KC system (fig 3.39, 3.40) Figure 3.38 the water circulation in the South China Sea in winter (top) and summer (bottom) after Wyrtki (1961) 249 The Kuroshio penetrates through the Luzon Strait in to the northern South China Sea forming a loop current When this happens part or all the Kuroshio flows to the west through the Luzon strait just north of 20N The water may travel as far west about 117E before turning north and then east In the passage to Java Sea the current during southwest monsoon is conducted to warm water transport from Indian Ocean to SCS 250 (c) (d) Figure 3.39 The water circulation in the South China Sea in winter after Fang and Fang 1998 (a) and model simulation (c) after Uu, 2001 251 Figure 3.40 the water circulation in the South China Sea in summer after Fang and Fang 1998 (b) and model simulation (d) after Uu, 2001 252 The simulated salinity fields (Uu, 2001) shows that the north Pacific subtropical maximum salinity water mass enters into north part of the South China sea following the Kuroshio loop current through Luzon Strait The occurrence of the Kuroshio loop current was modelled during the winter and summer monsoon season The simulated results show that the Kurioshio, the western boundary current of the North Pacific, has been found frequently to intrude into the northern SCS through the Luzon strait, which is the only deep and wide passage connecting the semi-enclosed SCS to western North Pacific The intruded Kurioshio moves westward, impinges upon the continental slope near the Dongsha Islands, and then splits The major part turns right to move eastward and leaves the SCS through the northern portion of the Luzon Strait, while the splinter turns left and moves southward, forming the SCS Branch (SCSB) The north-eastward current is the South China Sea Warm Current (SCSWC) The SCSWC is just onshore of the southwestward SCSB, but moves in opposite direction Their climatological appearance is shown in figure 3.41 for two seasons Figure 3.41 Climatological patches of SCSWC and SCSB in winter (top) and summers (bottom) 253 254 Darnell W.L., Staylor, F., Gupta, S.K & Denn, F.M 1988 Estimation of surface insolation using sun-synchronous satellite data J Climate (8): 820–835 REFERENCES Bortkovskii N.Z Ariel’, E.D Nadyozhina, V.V Simonov and A.V Sternzat, 1993: Over-Ocean atmospheric boundary layer dynamics under varying wind-wave development stage, Journal of Marine Systems, 4, 155160 Bjerknes J 1966 The possible response of the atmospheric Hadley circulation to equatorial anomalies of ocean temperature Tellus 4: 820-929 Cess R.D., M.H Zhang, P.Minnis, L.Corsetti, E.G Dutton, B.W Forgan, D.P Garber,W.L Gates, J.J Hack, E.F Harrison, X Jing, J.T Kiehl, C.N Long, J.-J Morcrette,G.L Potter, V Ramanathan, B Subasilar, C.H Whitlock, D.F Young, and Y Zhou 1995 Absorption of solar radiation by clouds: Observations versus models Science 267 (5197): 496–499 Charnock H., 1951.Meteorology-Energy transfer between the atmosphere and the ocean, Science Progress, 153, 80-95 Charnock H 1955 Wind stress on a water surface Quarterly Journal Royal Meteorological Society 81: 639–640 Charnock H 1987 Ocean currents and meridional transfers Quarterly J the Royal Meteorological Society Q.J.R Meteorol Soc 113 (475): 3–18 Cutler A N and J C Swallow (1984) Surface currents of the Indian Ocean Institute of Oceanographic Sciences Report 187, 8pp & 36 charts 255 Darnell W.L., W.F Staylor, S.K Gupta, N.A Ritchey, and A.C Wilbur 1992 Seasonal variation of surface radiation budget derived from International Satellite Cloud Climatology Project C1 data Journal of Geophysical Research 97: 15,741–15,760 10 DaSilva A., C.C Young, and S Levitus 1995 Atlas of surface marine data 1994 Vol.1: Algorithms and procedures National Oceanic and Atmospheric Administration Report 11 Dietrich G., K Kalle, W Krauss, and G Siedler 1980 General Oceanography 2nd ed Translated by Susanne and Hans Ulrich Roll New York: John Wiley and Sons (Wiley- Interscience) 12 Doronin Iu P., 1981, Atmosphere Ocean Interaction, Hydromet Publ House, Leningrad, 288 pp (in Russian) 13 Duncan C P., S G Schladow and W G Williams (1982) Surface currents near the Greater and Lesser Antilles International Hydrographic Review 59(2), 67-78 14 Duvanin A.I., 1968 Interaction models of largscale processes in the atmosphere and ocean, Journal Oceanology, VIII-4 (in Russian) 15 Ekman V.W 1905 On the in.uence of the Earth’s rotation on ocean currents Arkiv for Matematik, Astronomi, och Fysik: (11) 16 Fang W D and G H Fang, 1998: The recent progress in the study of the southern South China Sea circulation Advance in Earth Sciences, 13(2), 166-172 256 17 Farris A and M Wimbush, 1996: Wind-induced Kuroshio into the South China Sea, J of Oceanography, vol 52 Pp 771-784 26 Hinze J.O 1975 Turbulence 2nd ed New York: McGraw-Hill 18 Frouin R., Gautier, C & Morcrette, J-J 1988 Downward long wave irradiance at the ocean surface from satellite data: Methodology and in situ validation J Geophysical Res 93 (C1): 597–619 28 Houghton J.T., L.G.M Filho, B.A Callander, N Harris, A Kattenberg, and K Maskell 1996 Climate Change 1995: The Science of Climate Change Cambridge: Cambridge University Press 19 Gill A.E 1982 Atmosphere-Ocean Dynamics New York: Academic Press 29 Kawamura, R., 1994: A rotated EOF analysis of global sea surface temperature variability with interannual and interdecadal time scales J Phys Oceanogr., 24, 707-715 20 Godfrey J.S., A.C Hirst and J Wilkin, 1993: Why does the Indonesian through-flow appear to originate in the North Pacific? J Phys Oceanogr 23, 1087-1098 21 Goldstein S 1965 Modern Developments in Fluid Dynamics: Two Volumes New York: Dover Publications, Inc 22 Gu D and S.G.H Philander, 1997: Interdecadal climate fluctuations that depend on exchanges between the tropics and subtropics Science, 275, 805-807 23 Hasselmann K 1970 Wind–driven inertial Geophysical Fluid Dynamics 1: 1463–502 oscillations 24 Hasselmann K., T.P Barnett, E Bouws, H Carlson, D.E Cartwright, K Enke, J.A Ewing, H Gienapp, D.E Hasselmann, P Kruseman, A Meerburg, P Mller, D.J Olbers, K Richter, W Sell, and H Walden 1973 Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP) Ergnzungsheft zur Deutschen Hydrographischen Zeitschrift Reihe A(8) (Nr 12): 95 25 Helleman S., and Rosenstein M., 1983 Normal monthly windstress over the world ocean with error estimates, J Phys Oceanogr., 13, pp 1093-1104 257 27 Houghton J.T 1977 The Physics of Atmospheres Cambridge: Cambridge University Press 30 Kitaigorotxkii C.A., 1970: Physics of atmosphere Ocean Interaction, Hydromet Publ House, Leningrad, 284 pp (in Russian) 31 Komen G.J., L Cavaleri, M Donelan, K Hasselmann, S Hasselmann, and P.A.E.M Janssen.1996 Dynamics and Modelling of Ocean Waves 1st paperback ed Cambridge:Cambridge University Press 32 Kraus E B., and J A Businger, 1994: AtmosphereOcean Interaction, 2nd Edition, Oxford University Press, 362 pp 33 Kundu P.K 1990 Fluid Mechanics San Diego: Academic Press, Inc 34 Kunze E and J.M Toole 1997 Tidally driven vorticity, diurnal shear, and turbulenceatop Fieberling Seamount Journal of Physical Oceanography 27 (2): 2,663–2,693 35 Lagerloef G.S.E., G Mitchum, R Lukas, and P Niiler 1999 Tropical Paci.c near-surfacecurrent estimates from altimeter, 258 including the molecular constraints at the interface Journal of Atmospheric Science 36: 1722–1735 wind and drifter data Journal of Geophysical Research 104 (C10): 23,313–23,326 36 Lamb H 1945 Hydrodynamics 6th ed New York: Dover 37 Latif M and T.P Barnett, 1994: Causes of decadal climate variability over the North Pacific and North America Science, 266, 634- 637 46 Liu W.T 1988 Moisture and latent heat ux variabilities in the tropical Paci.c derived from satellite data J Geophysical Res 93 (C6): 6749–6760 38 Latif M and T.P Barnett, 1996: Decadal climate variability over the North Pacific and North America: Dynamics and predictability J Climate, 9, 2407-2423 47 Lowell T.V., C.J Heusser, B.G Andersen, P.I Moreno, A Houser, L.E Heusser, C Schlăuchter,D.R Marchant, and G.H Denton 1995 Interhemispheric correlations of late Pleistocene Glacial Events Science 269 (5230): 1541–1549 39 Lau K.-M., 1992: East Asian summer monsoon variability and climate teleconnection J Meteor Soc Japan, 70, 211-242 48 Lukas R (1986) The termination region of the Equatorial Undercurrent in the eastern Pacific Progress in Oceanography 16, 63-90 40 Lau N.-C, 1997: Interactions between global SST anomalies and the mid-latitude atmospheric circulation Bull Amer Meteor Soc., 78, 21-33 49 Lukas R., T Yamagata, and J McCreary, 1996: Pacific low-latitude western boundary currents and the Indonesian throughflow J Geophys Res., 101, 1220912216 41 Lau K.-M and W Bua, 1998: Mechanism of monsoonSouthern Oscillation coupling: insights from GCM experiments Climate Dynamics 42 Lau K.-M., and C.-H Sui, 1997: Mechanisms of ShortTerm Sea Surface Temperature Regulation: Observations during TOGA COARE J Climate , 10, 465-472 50 Madden R A., and P R Julian, 1972: Description of global-scale circulation cells in the tropics with a 40-50 day period J Atmos Sci., 29, 1109-1123 51 McPhaden M J (1982) Variability in the central equatorial Indian Ocean part I: Ocean dynamics Journal of Marine Research 40, 157-176 43 Lau K.-M and S Yang, 1996: The Asian monsoon and the predictability of the tropical ocean-atmosphere system Quart J Roy Meteor Soc., 122, 945-957 52 Metzger E.J and H.E Hurlburt, 1996 Coupled dynamics of South China Sea, the Sulu Sea, and the Pacific Ocean, J Geophysical Research, C5, 15 44 Levitus S 1982 Climatological Atlas of the World Ocean NOAA, Professional Paper 13 53 Miles J.W 1957 On the generation of surface waves by shear ows J Fluid Mechanics 3(2) 185–204 45 Liu W T., Katsaros, K B., & Businger, J A (1979) Bulk parameterization of air-sea exchanges of heat and water vapor, 259 260 54 Minobe S., 1996: Climate variability with periodicity of 50-70 years over the North Pacific and North America Geophys Res Lett., 24, 683-686 65 Ping-Tung Shaw and Shenn-Yu Chao, 1994 Surface circulation in the South China Sea, J Deep-Sea Research, 1, 41 55 Moskowitz L., 1964 Estimates of the power spectrums for fully developed seas for wind speeds of 20 to 40 knots J Geophysical Research 69 (24): 5161–5179 66 Plumb R.A 1977: The Interaction of two internal waves with the mean flow: implications for the theory of the quasi-biennial oscillation.Journal of Atmospheric Sciences Vol 34, pp 1847-1858 56 Open University, 1989 Ocean Circulation Oxford: Pergamon Press 57 Pedlosky J., 1996 Ocean Circulation Theory Berlin: Springer– Verlag 67 Pokudov V.V and V.A Golovastov (ed.), 1990 Ocean Atmosphere Interaction and Monsoon Dynamic, Hydromet Publ House, Leningrad, 336 pp (in Russian) 58 Petty G.W., 1995 Frequencies and characteristics of global oceanic precipitation from shipboard present-weather reports Bulletin of the American Meteorological Society 76 (9): 1593– 1616 68 Pohlmann T., 1987 A three-dimensional circulation model of the South China Sea, In Three-dimensional model of Marine and Estuarine Dynamics (Ed by J.C.J Nihoul and B.M.Jamart) Amsterdam, 245-268 59 Peterson R G and L Stramma (1991) Upper-level circulation in the South Atlantic Ocean Progress in Oceanography 26, 1-73 69 Ralph E.A., and P.P Niiler 2000 Wind-driven currents in the tropical Pacific Journal of Physical Oceanography 29 (9): 2121– 2129 60 Philander S G., 1990 El Niño, La Niña, and the Southern Oscillation Academic Press, San Diego 61 Phillips O.M 1957 On the generation of waves by turbulent wind J Fluid Mechanics (5): 417–445 62 Phillips O.M 1977 The Dynamics of the upper ocean Second ed Cambridge: Cambridge University Press 63 Perry A.H and J.M Walker, 1977 The ocean-atmosphere system, Longman Group Ltd., 190 pp 64 Pierson W.J., and L Moskowitz 1964 A proposed spectral form for fully developed wind seas based on the similarity theory of S.A Kitaigordskii J Geophysical Research 69: 5181–5190 261 70 Ramanathan V., B Subasilar, G.J Zhang, W Conant, R.D Cess, J.T.Kiehl, H Grassl, and L Shi 1995 Warm pool heat budget shortwave cloud forcing: A missing physics? Science 267 (5197): 499–503 71 Reynolds O 1883 An experimental investigation of the circumstances which determine whether the motion of water will be direct or sinuous, and the law of resistance in parallel channels Philosophical Transactions, Royal Society London 174: 935 72 Reynolds R.W., and T.M Smith 1995 A high-resolution global sea surface climatology J Climate (6): 1571–1583 262 73 Rossow W.B & Schiffer, R.A 1991 isccp Cloud Data Products Bull American Meteorology Society, 72 (1): 2–20 82 Stewart R., 2002 Introduction to Physical Oceanography, Texas A&M University 74 Sarachik E.S., M Winton, and F.L Yin, 1996: Mechanisms for decadal-to-centennial climate variability In: NATO ASI Vol 144, Decadal Climate Variability: Dynamics and Predictability, Eds.: D.L.T Anderson and J Willebrand, Springer Verlag, Heidelberg, Germany, 157-210 83 Stramma, L., Y Ikeda and R G Peterson, 1990 Geostrophic transport in the Brazil Current region north of 20°S Deep-Sea Research 37, 1875-1886 75 Schmetz J 1989 Towards a surface radiation climatology: Retrieval of downward irradiances from satellites Atmos Res 23: 287–321 76 Schubert S.D., R.B Rood, and J Pfaendtner 1993 An assimilated dataset for Earth science applications Bulletin American Meteorological Society 74 (12): 2331–2342 77 Schulz J., J Meywerk, S Ewald, and P Schluessel 1997 Evaluation of satellite-derived latent heat uxes Journal of Climate 10 (11): 2782–2795 78 Shenn-Yu Chao, Ping-Tung Shaw and Chau-Ron Wu, 1998 Seasonal and interannual variations in velocity field of the South China Sea, Journal of Oceanography, vol 54 1998, p 125-142 79 Smith S D., 1988 Coefficients for sea surface wind stress, J Geophys Res., 93, 15467-15472 80 Smith R.D., J.K Dukowitz, and R.C Malone 1992 Parallel ocean general circulation modeling Physica D 60: 38–45 81 Smith W.H.F., and D.T Sandwell, 1994 Bathymetric prediction from dense satellite altimetry and sparse shipboard bathymetry J Geophysical Research 99 (B11): 21,803– 21,824 263 84 SWAMP Group 1985 Ocean Wave Modeling New York: Plenum Press 85 Talley L.D 1984 Meridional heat transport in the Paci.c Ocean Journal of Physical Oceanography 14 (2): 231–241 100: 114– 121 86 Taylor V.R 1990 Implementation of reflectance models in operational AVHRR radiation budget processing Washington, D.C.: NOAA NESDIS Technical Report 49 87 Tennekes H and J.L Lumley 1990 A First Course in Turbulence Boston: MIT Press 88 The Open University, 1989 Ocean Circulation Milton Keynes, England: The Open University 89 Tomczak M and J.S Godfrey 1994 Regional Oceanography: An Introduction London: Pergamon 90 Tomita, T and T Yasunari, 1996: Role of the northeast winter monsoon on the biennial oscillation of the ENSO/monsoon system J Meteor Soc Japan, 74, 399413 91 Trenberth K.E., 1990: Recent observed interdecadal climate changes in the Northern Hemisphere Bull Amer Meteor Soc., 71, 988-993 92 Trenberth K.E and J.W Hurrell, 1994: Decadal atmosphere-ocean variations in the Pacific Climate Dynamics, 9, 303-319 264 93 Trenberth K.E., and A Solomon 1994 The global heat balance: heat transports in the atmosphere and ocean Climate Dynamics 10 (3): 107–134 94 Trenberth K E., and D P Stepaniak, Indices of El Niño evolution, J Clim., 2001 95 Tsuchiya M., R Lukas, R A Fine, E Firing and E Lindstrom (1989) Source waters of the Pacific Equatorial Undercurrent Progress in Oceanography 23, 101-147 96 Uu D V., 2001, Investigation of the Bien Dong threedimensional thermohaline and circulation structure and its aplications Final report of the Vietnamese Marine Research Project KHCN 06-02, Hanoi (in Vietnamese) 97 Uu D V and J-M Brankart, 1997 Seasonal variation of temperature and salinity fields and water masses in the Bien Dong (South China Sea), J Mathematical and Computer Modelling , No 12 98 Uu D.V., 1998: The seasonal variability of the circulation and thermohaline structure of the Bien Dong (South China) Sea in the condition of reversing monsoon winds: preliminary result of a three-dimensional model for its analysis and simulation Proceedings of The Fourth International Scientific Symposium UNESCO/IOC/WESTPAC, February, 1998, Okinawa, Japan, 100-109 99 WAMDI Group S Hasselmann, K Hasselmann, E Bauer, P.A.E.M Janssen, G.J Komen, L Bertotti, P Lionello, A Guillaume, V.C Cardone, J.A Greenwood, M Reistad, L Zambresky, and J.A Ewing 1988 The WAM model—A third 265 generation wave prediction model J Physical Oceanography 18: 1775–1810 100 Webster F 1968 Observations of inertial-period motions in the deep sea Reviews of Geophysics (4): 473–490 101 Weller R.A., J.P Dean, J Marra, J.F Price, E.A Francis, and D.C Boardman 1985.Three-dimensional flow in the upper ocean Science 227: 1552–1556 102 Wiegel R.L 1964 Oceanographical Engineering Englewood Cliffs, New Jersey: Prentice Hall 103 Wolter K., 1987: The Southern Oscillation in surface circulation and climate over the tropical Atlantic, Eastern Pacific, and Indian Oceans as captured by cluster analysis J Climate Appl Meteor., 26, 540558 104 Wolter K., and M.S Timlin, 1993: Monitoring ENSO in COADS with a seasonally adjusted principal component index Proc of the 17th Climate Diagnostics Workshop, Norman, OK, NOAA/N MC/CAC, NSSL, Oklahoma Clim Survey, CIMMS and the School of Meteor., Univ of Oklahoma, 52-57 105 Wyrtki K., 1962 Physical Oceanography of Southeast Asian Waters Scientific, Result of marine Investigation of South China Sea and Gulf of Thai Land 1959-1961 NAGA Report 2, California 106 Wyrtki K (1973) An equatorial jet in the Indian Ocean Science 181, 262-264 107 Wyrtki K 1975 El Ni˜no-The dynamic response of the equatorial Pacific Ocean to atmospheric forcing Journal of Physical Oceanography (4): 572–584 266 108 Xie P and P.A Arkin 1997 Global precipitation: A 17year monthly analysis based on gauge observations, satellite estimate, and numerical model outputs Bulletin of the American Meteorological Society 78 (1): 2539–2558 109 Youhai He and Cuihua Guan, 1998 South China Sea Warm Pool, Proceeding of The IV IOC/WESTPAC Scientific Symposium, Okinawa, 89-99 110 Zhang, Y., J.M Wallace, and D.S Battisti, 1997: ENSO-like interdecadal variability: 1900-93 J Climate, 10, 1004-1020 267 268 269 270

Ngày đăng: 06/09/2016, 14:02

TỪ KHÓA LIÊN QUAN

TÀI LIỆU CÙNG NGƯỜI DÙNG

TÀI LIỆU LIÊN QUAN

w