Engineering topics

Tao, Wei-Kuo

Publications and source records attributed to Tao, Wei-Kuo.

At least 73 records · Page 4

NASA NTRS2010

The Impact of Microphysics and Model Resolution on Precipitation Associated with Typhoon Morakot 2009

In recent years, the heavy rainfall that was associated with severe weather events (e.g., typhoons, local heavy precipitation events) has caused significant damages in the economy and loss of human life throughout Taiwan. Especially, the extreme heavy rainfall (over 2500 mm over 24 hours) associated with Typhoon Morakot 2009 caused more than 600 human beings lost and more than $100 million US dollar damage. In this paper, we are using WRF to simulate the precipitation processes associated Typhoon Morakot 2009. The preliminary results indicated that the wrf model with using 2 km grid size and with utilizing the 310E scheme (cloud ice, snow and hail) can simulate more than 2500 mm rainfall over 24 hour integration. In this talk, we will evaluate the performance of the microphysical schemes for the Typhoon Morakot case. In addition, we will examine the impact of model resolution (in both horizontal and vertical) on the Typhoon Morakot case.

Lin, Pay-Liam↗

NASA NTRS2010

The Role of Atmospheric Aerosol Concentration on Deep Convective Precipitation: Cloud-Resolving Model Simulations

Aerosols and especially their effect on clouds are one of the key components of the climate system and the hydrological cycle [Ramanathan et al., 2001]. Yet, the aerosol effect on clouds remains largely unknown and the processes involved not well understood. A recent report published by the National Academy of Science states "The greatest uncertainty about the aerosol climate forcing - indeed, the largest of all the uncertainties about global climate forcing - is probably the indirect effect of aerosols on clouds NRC [2001]." The aerosol effect on Clouds is often categorized into the traditional "first indirect (i.e., Twomey)" effect on the cloud droplet sizes for a constant liquid water path and the "semi-direct" effect on cloud coverage. The aerosol effect on precipitation processes, also known as the second type of aerosol indirect effect, is even more complex, especially for mixed-phase convective clouds. In this paper, a cloud-resolving model (CRM) with detailed spectral-bin microphysics was used to examine the effect of aerosols on three different deep convective cloud systems that developed in different geographic locations: South Florida, Oklahoma and the Central Pacific, In all three cases, rain reaches the ground earlier for the low CCN (clean) case. Rain suppression is also evident in all three cases with high CCN (dirty) case. However, this suppression only occurs during the first hour of the simulations. During the mature stages of the simulations, the effects of increasing aerosol concentration range from rain suppression in the Oklahoma case, to almost no effect in the Florida case, to rain enhancement in the Pacific case. These results show the complexity of aerosol interactions with convection. The model results suggest that evaporative cooling is a key process in determining whether high CCN reduces or enhances precipitation. Stronger evaporative cooling can produce a stronger cold pool and thus stronger low-level convergence through interactions with the low-level wind shear. Consequently, precipitation processes can be more vigorous. For example,, the evaporative cooling is more than two times stronger in the lower troposphere with high CCN for the Pacific case. Sensitivity tests also suggest that ice processes are crucial for suppressing precipitation in the Oklahoma case with high CCN.

Tao, Wei-Kuo↗

NASA NTRS2010

Use the High-Resolution Numerical Model to Simulate Typhoon Morakot 2009

Typhoon Morakot struck Taiwan on the night of Friday August 7th, 2009 as a category 2 storm with sustained winds of 85 knots (92 mph). Although the center made landfall in Hualien county along the central cast coast of Taiwan and passed over the central northern part of the island, it was southern Taiwan that received the worst effects of the storm where locally as Much as 2200 mm (2.2 m) of rain were reported, resulting in the worst flooding there in 50 years. The result of the enormous amount of rain has been massive flooding and devastating mudslides. More than 600 people are confirmed dead. In this paper, we will present the results from high-resolution (2-km) WRF for this typhoon case. The results showed that the model captured both in terms of maximum rainfall area and intensity. The model results also showed that the heavy amounts of rain over the southern portion of the island is due to persistent southwesterly flow associated with Morakot and it's circulation was able to draw up copious amounts of moisture from the South China Sea into southern Taiwan where it was able to interact with the steep topography. In the paper, we will also present results from sensitivity test of terrain heights and SST on the precipitation processes (rainfall) associated with Typhoon Morakot (2009), In addition, we will present high-resolution visualization (36 second and 2-km) to show the evolution of Typhoon Morakot.

Tao, Wei-Kuo↗

NASA NTRS2010

Relating Convective and Stratiform Rain to Latent Heating

The relationship among surface rainfall, its intensity, and its associated stratiform amount is established by examining observed precipitation data from the Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR). The results show that for moderate-high stratiform fractions, rain probabilities are strongly skewed toward light rain intensities. For convective-type rain, the peak probability of occurrence shifts to higher intensities but is still significantly skewed toward weaker rain rates. The main differences between the distributions for oceanic and continental rain are for heavily convective rain. The peak occurrence, as well as the tail of the distribution containing the extreme events, is shifted to higher intensities for continental rain. For rainy areas sampled at 0.58 horizontal resolution, the occurrence of conditional rain rates over 100 mm/day is significantly higher over land. Distributions of rain intensity versus stratiform fraction for simulated precipitation data obtained from cloud-resolving model (CRM) simulations are quite similar to those from the satellite, providing a basis for mapping simulated cloud quantities to the satellite observations. An improved convective-stratiform heating (CSH) algorithm is developed based on two sources of information: gridded rainfall quantities (i.e., the conditional intensity and the stratiform fraction) observed from the TRMM PR and synthetic cloud process data (i.e., latent heating, eddy heat flux convergence, and radiative heating/cooling) obtained from CRM simulations of convective cloud systems. The new CSH algorithm-derived heating has a noticeably different heating structure over both ocean and land regions compared to the previous CSH algorithm. Major differences between the new and old algorithms include a significant increase in the amount of low- and midlevel heating, a downward emphasis in the level of maximum cloud heating by about 1 km, and a larger variance between land and ocean in the new CSH algorithm.

Tao, Wei-Kuo↗

NASA NTRS2010

Shallow and Deep Latent Heating Modes Over Tropical Oceans Observed with TRMM PR Spectral Latent Heating Data

The global hydrological cycle is central to the Earth's climate system, with rainfall and the physics of its formation acting as the key links in the cycle. Two-thirds of global rainfall occurs in the Tropics. Associated with this rainfall is a vast amount of heat, which is known as latent heat. It arises mainly due to the phase change of water vapor condensing into liquid droplets; three-fourths of the total heat energy available to the Earth's atmosphere comes from tropical rainfall. In addition, fresh water provided by tropical rainfall and its variability exerts a large impact upon the structure and motions of the upper ocean layer. Three-dimensional distributions of latent heating estimated from Tropical Rainfall Measuring Mission Precipitation Radar (TRMM PR)utilizing the Spectral Latent Heating (SLH) algorithm are analyzed. Mass-weighted and vertically integrated latent heating averaged over the tropical oceans is estimated as approx.72.6 J/s (approx.2.51 mm/day), and that over tropical land is approx.73.7 J/s (approx.2.55 mm/day), for 30degN-30degS. It is shown that non-drizzle precipitation over tropical and subtropical oceans consists of two dominant modes of rainfall systems, deep systems and congestus. A rough estimate of shallow mode contribution against the total heating is about 46.7 % for the average tropical oceans, which is substantially larger than 23.7 % over tropical land. While cumulus congestus heating linearly correlates with the SST, deep mode is dynamically bounded by large-scale subsidence. It is notable that substantial amount of rain, as large as 2.38 mm day-1 in average, is brought from congestus clouds under the large-scale subsiding circulation. It is also notable that even in the region with SST warmer than 28 oC, large-scale subsidence effectively suppresses the deep convection, remaining the heating by congestus clouds. Our results support that the entrainment of mid-to-lower-tropospheric dry air, which accompanies the large-scale subsidence is the major factor suppressing the deep convection. Therefore, representation of the realistic entrainment is very important for proper reproduction of precipitation distribution and resultant large-scale circulation.

Takayabu, Yukari N.↗

NASA NTRS2010

Using Multi-Scale Modeling Systems to Study the Precipitation Processes

In recent years, exponentially increasing computer power has extended Cloud Resolving Model (CRM) integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-km scales are resolved in horizontal domains as large as 10,000 km in two-dimensions, and 1,000 x 1,000 km2 in three-dimensions. Cloud resolving models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that NWP and mesoscale model can be run in grid size similar to cloud resolving model through nesting technique. Recently, a multi-scale modeling system with unified physics was developed at NASA Goddard. It consists of (1) a cloud-resolving model (Goddard Cumulus Ensemble model, GCE model), (2) a regional scale model (a NASA unified weather research and forecast, WRF), (3) a coupled CRM and global model (Goddard Multi-scale Modeling Framework, MMF), and (4) a land modeling system. The same microphysical processes, long and short wave radiative transfer and land processes and the explicit cloud-radiation, and cloud-land surface interactive processes are applied in this multi-scale modeling system. This modeling system has been coupled with a multi-satellite simulator to use NASA high-resolution satellite data to identify the strengths and weaknesses of cloud and precipitation processes simulated by the model. In this talk, a review of developments and applications of the multi-scale modeling system will be presented. In particular, the results from using multi-scale modeling system to study the interactions between clouds, precipitation, and aerosols will be presented. Also how to use of the multi-satellite simulator to improve precipitation processes will be discussed.

Tao, Wei-Kuo↗

NASA NTRS2010

Ten-Year Climatology of Summertime Diurnal Rainfall Rate Over the Conterminous U.S.

Diurnal cycles of summertime rainfall rates are examined over the conterminous United States, using radar-gauge assimilated hourly rainfall data. As in earlier studies, rainfall diurnal composites show a well-defined region of rainfall propagation over the Great Plains and an afternoon maximum area over the south and eastern portion of the United States. Zonal phase speeds of rainfall in three different small domains are estimated, and rainfall propagation speeds are compared with background zonal wind speeds. Unique rainfall propagation speeds in three different regions can be explained by the evolution of latent-heat theory linked to the convective available potential energy, than by gust-front induced or gravity wave propagation mechanisms.

Matsui, Toshihisa↗

NASA NTRS2010

Improving a Spectral Bin Microphysical Scheme Using TRMM Satellite Observations

Comparisons between cloud model simulations and observations are crucial in validating model performance and improving physical processes represented in the mod Tel.hese modeled physical processes are idealized representations and almost always have large rooms for improvements. In this study, we use data from two different sensors onboard TRMM (Tropical Rainfall Measurement Mission) satellite to improve the microphysical scheme in the Goddard Cumulus Ensemble (GCE) model. TRMM observed mature-stage squall lines during late spring, early summer in central US over a 9-year period are compiled and compared with a case simulation by GCE model. A unique aspect of the GCE model is that it has a state-of-the-art spectral bin microphysical scheme, which uses 33 different bins to represent particle size distribution of each of the seven hydrometeor species. A forward radiative transfer model calculates TRMM Precipitation Radar (PR) reflectivity and TRMM Microwave Imager (TMI) 85 GHz brightness temperatures from simulated particle size distributions. Comparisons between model outputs and observations reveal that the model overestimates sizes of snow/aggregates in the stratiform region of the squall line. After adjusting temperature-dependent collection coefficients among ice-phase particles, PR comparisons become good while TMI comparisons worsen. Further investigations show that the partitioning between graupel (a high-density form of aggregate), and snow (a low-density form of aggregate) needs to be adjusted in order to have good comparisons in both PR reflectivity and TMI brightness temperature. This study shows that long-term satellite observations, especially those with multiple sensors, can be very useful in constraining model microphysics. It is also the first study in validating and improving a sophisticated spectral bin microphysical scheme according to long-term satellite observations.

Li, Xiaowen↗

NASA NTRS2010

MJO Signals in Latent Heating: Results from TRMM Retrievals

The Madden-Julian Oscillation (MJO) is the dominant intraseasonal signal in the global tropical atmosphere. Almost all numerical climate models have difficulty to simulate realistic MJO. Four TRMM datasets of latent heating were diagnosed for signals in the MJO. In all four datasets, vertical structures of latent heating are dominated by two components, one deep with its peak above the melting level and one shallow with its peak below. Profiles of the two components are nearly ubiquitous in longitude, allowing a separation of the vertical and zonal/temporal variations when the latitudinal dependence is not considered. All four datasets exhibit robust MJO spectral signals in the deep component as eastward propagating spectral peaks centered at period of 50 days and zonal wavenumber 1, well distinguished from lower- and higher-frequency power and much stronger than the corresponding westward power. The shallow component shows similar but slightly less robust MJO spectral peaks. MJO signals were further extracted from a combination of band-pass (30 - 90 day) filtered deep and shallow components. Largest amplitudes of both deep and shallow components of the MJO are confined to the Indian and western Pacific Oceans. There is a local minimum in the deep components over the Maritime Continent. The shallow components of the MJO differ substantially among the four TRMM datasets in their detailed zonal distributions in the eastern hemisphere. In composites of the heating evolution through the life cycle of the MJO, the shallow components lead the deep ones in some datasets and at certain longitudes. In many respects, the four TRMM datasets agree well in their deep components, but not in their shallow components and the phase relations between the deep and shallow components. These results indicate that caution must be exercised in applications of these latent heating data.

Zhang, Chidong↗

NASA NTRS2009

Multiscale Cloud System Modeling

The central theme of this paper is to describe how cloud system resolving models (CRMs) of grid spacing approximately 1 km have been applied to various important problems in atmospheric science across a wide range of spatial and temporal scales and how these applications relate to other modeling approaches. A long-standing problem concerns the representation of organized precipitating convective cloud systems in weather and climate models. Since CRMs resolve the mesoscale to large scales of motion (i.e., 10 km to global) they explicitly address the cloud system problem. By explicitly representing organized convection, CRMs bypass restrictive assumptions associated with convective parameterization such as the scale gap between cumulus and large-scale motion. Dynamical models provide insight into the physical mechanisms involved with scale interaction and convective organization. Multiscale CRMs simulate convective cloud systems in computational domains up to global and have been applied in place of contemporary convective parameterizations in global models. Multiscale CRMs pose a new challenge for model validation, which is met in an integrated approach involving CRMs, operational prediction systems, observational measurements, and dynamical models in a new international project: the Year of Tropical Convection, which has an emphasis on organized tropical convection and its global effects.

Tao, Wei-Kuo↗

NASA NTRS2009

Coupling Advanced Modeling and Visualization to Improve High-Impact Tropical Weather Prediction

To meet the goals of extreme weather event warning, this approach couples a modeling and visualization system that integrates existing NASA technologies and improves the modeling system's parallel scalability to take advantage of petascale supercomputers. It also streamlines the data flow for fast processing and 3D visualizations, and develops visualization modules to fuse NASA satellite data.

Shen, Bo-Wen↗

NASA NTRS2009

Diagnosing the Ice Crystal Enhancement Factor in the Tropics

Recent modeling studies have revealed that ice crystal number concentration is one of the dominant factors in the effect of clouds on radiation. Since the ice crystal enhancement factor and ice nuclei concentration determine the concentration, they are both important in quantifying the contribution of increased ice nuclei to global warming. In this study, long-term cloud-resolving model (CRM) simulations are compared with field observations to estimate the ice crystal enhancement factor in tropical and midlatitudinal clouds, respectively. It is found that the factor in tropical clouds is ~10 3-104 times larger than that of mid-latitudinal ones, which makes physical sense because entrainment and detrainment in the Tropics are much stronger than in middle latitudes. The effect of entrainment/detrainment on the enhancement factor, especially in tropical clouds, suggests that cloud microphysical parameterizations should be coupled with subgrid turbulence parameterizations within CRMs to obtain a more accurate depiction of cloud-radiative forcing.

Zeng, Xiping↗

NASA NTRS2009

A Multiscale Modeling System: Developments, Applications, and Critical Issues

The foremost challenge in parameterizing convective clouds and cloud systems in large-scale models are the many coupled dynamical and physical processes that interact over a wide range of scales, from microphysical scales to the synoptic and planetary scales. This makes the comprehension and representation of convective clouds and cloud systems one of the most complex scientific problems in Earth science. During the past decade, the Global Energy and Water Cycle Experiment (GEWEX) Cloud System Study (GCSS) has pioneered the use of single-column models (SCMs) and cloud-resolving models (CRMs) for the evaluation of the cloud and radiation parameterizations in general circulation models (GCMs; e.g., GEWEX Cloud System Science Team 1993). These activities have uncovered many systematic biases in the radiation, cloud and convection parameterizations of GCMs and have led to the development of new schemes (e.g., Zhang 2002; Pincus et al, 2003; Zhang and Wu 2003; Wu et al. 2003; Liang and Wu 2005; Wu and Liang 2005, and others). Comparisons between SCMs and CRMs using the same large-scale forcing derived from field campaigns have demonstrated that CRMs are superior to SCMs in the prediction of temperature and moisture tendencies (e.g., Das et al. 1999; Randall et al 2003b; Xie et al. 2005).

Tao, Wei-Kuo↗

NASA NTRS2009

Combining Satellite Microwave Radiometer and Radar Observations to Estimate Atmospheric Latent Heating Profiles

In this study, satellite passive microwave sensor observations from the TRMM Microwave Imager (TMI) are utilized to make estimates of latent + eddy sensible heating rates (Q1-QR) in regions of precipitation. The TMI heating algorithm (TRAIN) is calibrated, or "trained" using relatively accurate estimates of heating based upon spaceborne Precipitation Radar (PR) observations collocated with the TMI observations over a one-month period. The heating estimation technique is based upon a previously described Bayesian methodology, but with improvements in supporting cloud-resolving model simulations, an adjustment of precipitation echo tops to compensate for model biases, and a separate scaling of convective and stratiform heating components that leads to an approximate balance between estimated vertically-integrated condensation and surface precipitation. Estimates of Q1-QR from TMI compare favorably with the PR training estimates and show only modest sensitivity to the cloud-resolving model simulations of heating used to construct the training data. Moreover, the net condensation in the corresponding annual mean satellite latent heating profile is within a few percent of the annual mean surface precipitation rate over the tropical and subtropical oceans where the algorithm is applied. Comparisons of Q1 produced by combining TMI Q1-QR with independently derived estimates of QR show reasonable agreement with rawinsonde-based analyses of Q1 from two field campaigns, although the satellite estimates exhibit heating profile structure with sharper and more intense heating peaks than the rawinsonde estimates. 2

Grecu, Mircea↗

NASA NTRS2009

Production of Lightning NO(x) and its Vertical Distribution Calculated from 3-D Cloud-scale Chemical Transport Model Simulations

A 3-D cloud scale chemical transport model that includes a parameterized source of lightning NO(x), based on observed flash rates has been used to simulate six midlatitude and subtropical thunderstorms observed during four field projects. Production per intracloud (P(sub IC) and cloud-to-ground (P(sub CG)) flash is estimated by assuming various values of P(sub IC) and P(sub CG) for each storm and determining which production scenario yields NO(x) mixing ratios that compare most favorably with in-cloud aircraft observations. We obtain a mean P(sub CG) value of 500 moles NO (7 kg N) per flash. The results of this analysis also suggest that on average, P(sub IC) may be nearly equal to P(sub CG), which is contrary to the common assumption that intracloud flashes are significantly less productive of NO than are cloud-to-ground flashes. This study also presents vertical profiles of the mass of lightning NO(x), after convection based on 3-D cloud-scale model simulations. The results suggest that following convection, a large percentage of lightning NO(x), remains in the middle and upper troposphere where it originated, while only a small percentage is found near the surface. The results of this work differ from profiles calculated from 2-D cloud-scale model simulations with a simpler lightning parameterization that were peaked near the surface and in the upper troposphere (referred to as a "C-shaped" profile). The new model results (a backward C-shaped profile) suggest that chemical transport models that assume a C-shaped vertical profile of lightning NO(x) mass may place too much mass neat the surface and too little in the middle troposphere.

Ott, Lesley↗

NASA NTRS2009

Modeling and Observational Framework for Diagnosing Local Land-Atmosphere Coupling on Diurnal Time Scales

Land-atmosphere interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture states. The degree of coupling between the land surface and PBL in numerical weather prediction and climate models remains largely unexplored and undiagnosed due to the complex interactions and feedbacks present across a range of scales. Further, uncoupled systems or experiments (e.g., the Project for Intercomparison of Land Parameterization Schemes, PILPS) may lead to inaccurate water and energy cycle process understanding by neglecting feedback processes such as PBL-top entrainment. In this study, a framework for diagnosing local land-atmosphere coupling is presented using a coupled mesoscale model with a suite of PBL and land surface model (LSM) options along with observations during field experiments in the U. S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to the Land Information System (LIS), which provides a flexible and high-resolution representation and initialization of land surface physics and states. Within this framework, the coupling established by each pairing of the available PBL schemes in WRF with the LSMs in LIS is evaluated in terms of the diurnal temperature and humidity evolution in the mixed layer. The co-evolution of these variables and the convective PBL is sensitive to and, in fact, integrative of the dominant processes that govern the PBL budget, which are synthesized through the use of mixing diagrams. Results show how the sensitivity of land-atmosphere interactions to the specific choice of PBL scheme and LSM varies across surface moisture regimes and can be quantified and evaluated against observations. As such, this methodology provides a potential pathway to study factors controlling local land-atmosphere coupling (LoCo) using the LIS-WRF system, which will serve as a testbed for future experiments to evaluate coupling diagnostics within the community.

Santanello, Joseph A., Jr.↗

NASA NTRS2009

Convective and Stratiform Precipitation Processes and their Relationship to Latent Heating

The global hydrological cycle is central to the Earth's climate system, with rainfall and the physics of its formation acting as the key links in the cycle. Two-thirds of global rainfall occurs in the Tropics. Associated with this rainfall is a vast amount of heat, which is known as latent heat. It arises mainly due to the phase change of water vapor condensing into liquid droplets; three-fourths of the total heat energy available to the Earth's atmosphere comes from tropical rainfall. In addition, fresh water provided by tropical rainfall and its variability exerts a large impact upon the structure and motions of the upper ocean layer. An improved convective -stratiform heating (CSH) algorithm has been developed to obtain the 3D structure of cloud heating over the Tropics based on two sources of information: 1) rainfall information, namely its amount and the fraction due to light rain intensity, observed directly from the Precipitation Radar (PR) on board the TRMM satellite and 2) synthetic cloud physics information obtained from cloud-resolving model (CRM) simulations of cloud systems. The cloud simulations provide details on cloud processes, specifically latent heating, eddy heat flux convergence and radiative heating/cooling, that. are not directly observable by satellite. The new CSH algorithm-derived heating has a noticeably different heating structure over both ocean and land regions compared to the previous CSH algorithm. One of the major differences between new and old algorithms is that the level of maximum cloud heating occurs 1 to 1.5 km lower in the atmosphere in the new algorithm. This can effect the structure of the implied air currents associated with the general circulation of the atmosphere in the Tropics. The new CSH algorithm will be used provide retrieved heating data to other heating algorithms to supplement their performance.

Tao, Wei-Kuo↗

NASA NTRS2009

A Contribution by Ice Nuclei to Global Warming

Ice nuclei (IN) significantly affect clouds via supercooled droplets, that in turn modulate atmospheric radiation and thus climate change. Since the IN effect is relatively strong in stratiform clouds but weak in convective ones, the overall effect depends on the ratio of stratiform to convective cloud amount. In this paper, 10 years of TRMM (Tropical Rainfall Measuring Mission) satellite data are analyzed to confirm that stratiform precipitation fraction increases with increasing latitude, which implies that the IN effect is stronger at higher latitudes. To quantitatively evaluate the IN effect versus latitude, large-scale forcing data from ten field campaigns are used to drive a CRM (cloud-resolving model) to generate longterm cloud simulations. As revealed in the simulations, the increase in the net downward radiative flux at the TOA (top of the atmosphere) from doubling the current IN concentrations is larger at higher latitude, which is attributed to the meridional tendency in the stratiform precipitation fraction. Surface warming from doubling the IN concentrations, based on the radiative balance of the globe, is compared with that from anthropogenic COZ . It is found that the former effect is stronger than the latter in middle and high latitudes but not in the Tropics. With regard to the impact of IN on global warming, there are two factors to consider: the radiative effect from increasing the IN concentration and the increase in IN concentration itself. The former relies on cloud ensembles and thus varies mainly with latitude. In contrast, the latter relies on IN sources (e.g., the land surface distribution) and thus varies not only with latitude but also longitude. Global desertification and industrialization provide clues on the geographic variation of the increase in IN concentration since pre-industrial times. Thus, their effect on global warming can be inferred and then be compared with observations. A general match in geographic and seasonal variations between the inferred and observed warming suggests that IN may have contributed positively to global warming over the past decades, especially in middle and high latitudes.

Zeng, Xiping↗