SEARCH · Engineering Papers

Results for “Cloud parameterizations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Exploring the Inner Edge of the Habitable Zone with Fully Coupled Oceans

The role of rotation in planetary atmospheres plays an important role in regulating atmospheric and oceanic heat flow, cloud formation and precipitation. Using the Goddard Institute for Space Studies (GISS) three dimension General Circulation Model (3D-GCM) we demonstrate how varying rotation rate and increasing the incident solar flux on a planet are related to each other and may allow the inner edge of the habitable zone to be much closer than many previous habitable zone studies have indicated. This is shown in particular for fully coupled ocean runs -- some of the first that have been utilized in this context. Results with a 100m mixed layer depth and our fully coupled ocean runs are compared with those of Yang et al. 2014, which demonstrates consistency across models. However, there are clear differences for rotations rates of 1-16x present earth day lengths between the mixed layer and fully couple ocean models, which points to the necessity of using fully coupled oceans whenever possible. The latter was recently demonstrated quite clearly by Hu & Yang 2014 in their aquaworld study with a fully coupled ocean when compared with similar mixed layer ocean studies and by Cullum et al. 2014. Atmospheric constituent amounts were also varied alongside adjustments to cloud parameterizations (results not shown here). While the latter have an effect on what a planet's global mean temperature is once the oceans reach equilibrium they do not qualitatively change the overall relationship between the globally averaged surface temperature and incident solar flux for rotation rates ranging from 1 to 256 times the present Earth day length. At the same time this study demonstrates that given the lack of knowledge about the atmospheric constituents and clouds on exoplanets there is still a large uncertainty as to where a planet will sit in a given star's habitable zone.

solar flux↗

NASA NTRS1990

A scheme for parameterizing ice-cloud water content in general circulation models

The optical properties of ice clouds are a primary issue for climate and climate change. Evaluating these optical properties in three-dimensional models for studying climate will require a method to calculate the ice water content of such clouds. A procedure is developed to parameterize ice water content as a function of large-scale meteorological characteristics for use in circulation models in which the ice water content is not calculated by means of a three-dimensional prognostic equation for condensed water. The technique identifies large-scale flows in which ice clouds exist and calculates their ice water content by reconstructing the trajectory associated with cloud formation. As the cloud forms, its ice content changes both by deposition of ice from water vapor and by ice removal by sedimentation. The sedimentation process is found to modify significantly the ice water content expected from deposition alone. Ice water contents predicted by the parameterization are compared with aircraft observations collected in the middle latitudes and the tropics, and show reasonable agreement over four orders-of-magnitude of ice water content. A parameterization for the sublimation of ice crystals settling into ice-subsaturated environments is also presented.

Heymsfield, Andrew J.↗

NASA NTRS1987

FIRE - The First ISCCP Regional Experiment

The First International Satellite Cloud Climatology Project Regional Experiment (FIRE) designed to study the roles of clouds, in particular marine stratocumulus and cirrus-cloud systems, in the global climate is discussed. The objectives of FIRE are: (1) to develop a cloud-classification scheme; (2) to validate and improve satellite cloud-retrieval techniques; (3) to improve cloud radiation models; (4) to collect cloud space/time statistics; (5) to improve cloud dynamics models; and (6) to validate and improve GCM cloud parameterizations. The methods used to acquire extended time data and intensive field observations are described. The extended time and intensive field data collected during the FIRE are to be archived in the NASA Pilot Climate Data System at Goddard Space Flight Center.

Cox, Stephen K.↗

NASA NTRS1989

Earth radiation budget and cloudiness simulations with a general circulation model

A GCM with new parameterizations of solar and terrestrial radiation, parameterized cloud optical properties, and a simple representation of the cloud liquid water feedback is used with several observational data sets to analyze the effects of cloudiness on the earth's radiation budget. The January and July results from the model are in reasonable agreement with data from Nimbus-7. It is found that the simulated cloudiness overpredicts subtropical and midlatitude cloudiness. The simulated atmospheric cloud radiative forcing is examined. The clear-sky radiation fields obtained by two methods of Cess and Potter (1987) are compared. Also, a numerical experiment was performed to determine the effects of the water vapor continuum on the model results.

HARSHVARDHAN↗

NASA NTRS1991

The effect of clouds on the earth's radiation budget

The radiative fluxes from the Earth Radiation Budget Experiment (ERBE) and the cloud properties from the International Satellite Cloud Climatology Project (ISCCP) over Indonesia for the months of June and July of 1985 and 1986 were analyzed to determine the cloud sensitivity coefficients. The method involved a linear least squares regression between co-incident flux and cloud coverage measurements. The calculated slope is identified as the cloud sensitivity. It was found that the correlations between the total cloud fraction and radiation parameters were modest. However, correlations between cloud fraction and IR flux were improved by separating clouds by height. Likewise, correlations between the visible flux and cloud fractions were improved by distinguishing clouds based on optical depth. Calculating correlations between the net fluxes and either height or optical depth segregated cloud fractions were somewhat improved. When clouds were classified in terms of their height and optical depth, correlations among all the radiation components were improved. Mean cloud sensitivities based on the regression of radiative fluxes against height and optical depth separated cloud types are presented. Results are compared to a one-dimensional radiation model with a simple cloud parameterization scheme.

Ziskin, Daniel↗

NASA NTRS1992

Models of atmosphere-ecosystem-hydrology interactions: Approaches and testing

Interactions among the atmosphere, terrestrial ecosystems, and the hydrological cycle have been the subject of investigation for many years, although most of the research has had a regional focus. The topic is broad, including the effects of climate and hydrology on vegetation, the effects of vegetation on hydrology, the effects of the hydrological cycle on the atmosphere, and interactions of the cycles via material flux such as solutes and trace gases. The intent of this paper is to identify areas of critical uncertainty, discuss modeling approaches to resolving those problems, and then propose techniques for testing. I consider several interactions specifically to illustrate the range of problems. These areas are as follows: (1) cloud parameterizations and the land surface, (2) soil moisture, and (3) the terrestrial carbon cycle.

Schimel, David S.↗

NASA NTRS1995

Comparison of simulated cloud cover with satellite obsrvations over the Western United States

Satellite imagery datasets and regional climate model results are intercompared for evaluation of model accuracy in the simulation of cloud cover. Both monthly average individual simulation times are analyzed. To provide a consistent comparison, satellite data are first mapped into the model's geographic projection, grid domain, and resolution. It is found that September 1988 monthly average cloud fraction results from the modeled simulations correspond to observations, in both spatial pattern and magnitude, with bias less than +/- 20% cloud fraction over the entire inland West. Agreement in the pattern of cloud fraction also is evident for monthly average cloud fraction in July, but there is no negative bias of 10%-30% cloud fraction in the model diagnosis of cloud cover. Correlations between the spatial distributions of model-derived and observed cloud fractions are found to exceed 0.80 for certain geographic regions of the West, and these correlations are largest over mountainous areas during summer. Case studies of a series of daily cloud cover demonstrate the ability of the model to simulate the effects of frontal passage on cloud distribution. The ability of the RegCM1 to simulate daily cloud fraction and diurnal cloud evolution is somewhat weak for the summer convective season. It is anticipated that a more recent version of the regional climate model may improve the simulation of summer season cloud cover, through changes in cloud parameterization and improvements in model resolution.

Wetzel, Melanie A.↗

NASA NTRS2001

The Cloudsat Mission and the EOS Constellation: A New Dimension of Space-Based Observation of Clouds and Precipitation

CloudSat is a satellite experiment designed to measure the vertical structure of clouds from space. The expected launch of CloudSat is planned for 2004 and, once launched, CloudSat will orbit in formation as part of a constellation of satellites including NASA's Aqua and Aura satellites, a NASA-CNES lidar satellite (P-C) and a CNES satellite carrying a polarimeter (PARASOL). A unique feature that CloudSat brings to this constellation is the ability to fly a precise orbit enabling the fields of view of the CloudSat radar to be overlapped with the P-C lidar footprint and the other measurements of the EOS constellation. The precision of this overlap creates a unique multi-satellite observing system for studying the atmospheric processes essential to the hydrological cycle. The vertical profile of cloud properties provided by CloudSat fills a critical gap in the investigation of feedback mechanisms linking clouds to climate. Measuring the vertical profile of cloud properties requires a combination of active and passive instruments, and this will be achieved by combining the radar data of CloudSat with active and passive data from other sensors of the constellation. This paper describes the underpinning science, and gives an overview of the mission, and provides some idea of the expected products and anticipated application of these products. Notably, the CloudSat mission is expected to provide new knowledge about global cloudiness, stimulating new areas of research on clouds including data assimilation and cloud parameterization. The mission also provides an important opportunity to demonstrate active sensor technology for future scientific and tactical applications. The CloudSat mission is a partnership between NASA/JPL, the Canadian Space Agency, Colorado State University, the US Air Force, and the US Department of Energy.

Stephens, Graeme L.↗

NASA NTRS2001

Sensitivity of the Tropical Atmospheric Energy Balance to ENSO-Related SST Changes: Comparison of Climate Model Simulations to Observed Responses

This paper focuses on how fresh water and radiative fluxes over the tropical oceans change during ENSO warm and cold events and how these changes affect the tropical energy balance. At present, ENSO remains the most prominent known mode of natural variability at interannual time scales. While this natural perturbation to climate is quite distinct from possible anthropogenic changes in climate, adjustments in the tropical water and energy budgets during ENSO may give insight into feedback processes involving water vapor and cloud feedbacks. Although great advances have been made in understanding this phenomenon and realizing prediction skill over the past decade, our ability to document the coupled water and energy changes observationally and to represent them in climate models seems far from settled (Soden, 2000 J Climate). In a companion paper we have presented observational analyses, based principally on space-based measurements which document systematic changes in rainfall, evaporation, and surface and top-of-atmosphere (TOA) radiative fluxes. Here we analyze several contemporary climate models run with observed SSTs over recent decades and compare SST-induced changes in radiation, precipitation, evaporation, and energy transport to observational results. Among these are the NASA / NCAR Finite Volume Model, the NCAR Community Climate Model, the NCEP Global Spectral Model, and the NASA NSIPP Model. Key disagreements between model and observational results noted in the recent literature are shown to be due predominantly to observational shortcomings. A reexamination of the Langley 8-Year Surface Radiation Budget data reveals errors in the SST surface longwave emission due to biased SSTs. Subsequent correction allows use of this data set along with ERBE TOA fluxes to infer net atmospheric radiative heating. Further analysis of recent rainfall algorithms provides new estimates for precipitation variability in line with interannual evaporation changes inferred from the da Silva, Young, Levitus COADS analysis. The overall results from our analysis suggest an increase (decrease) of the hydrologic cycle during ENSO warm (cold) events at the rate of about 5 W/sq m per K of SST change. Model results agree reasonably well with this estimate of sensitivity. This rate is slightly less than that which would be expected for constant relative humidity over the tropical oceans. There remain, however, significant quantitative uncertainties in cloud forcing changes in the models as compared to observations. These differences are examined in relationship to model convection and cloud parameterizations Analysis of the possible sampling and measurement errors compared to systematic model errors is also presented.

Robertson, Franklin R.↗

NASA NTRS1998

Making Intercontinental Burning Influences Quantitative: Lessons on Transport and Geochemistry From the 3-D Simulation of CO in PEM-T A and TRACE-A

We report on two aspects of the simulation of global transport of plumes originating from subtropical biomass fires. We use of meteorological assimilation (MM5) at 2-degree resolution with a Grail cloud parameterization and a Blackadar-based planetary boundary layer parameterization. Ames's GRACES model provides emissions, transport, and an appropriate level of simulated chemical transformation. We have worked with passive-tracer CO or linear chemistry. This is appropriate since we find major work to be done in evaluating CO source strengths and transport mechanisms before chemical integrations could be meaningful. First, we present mechanisms by which CO and other pollutants are introduced into the free troposphere, and are then transported with hate dilution from approx. 0 to approx. 180 degrees longitude. One principal conduit for these plumes is the vernal subtropical jet; however the plumes appear at various altitudes and latitudes as they influenced by frontal motions and (most likely) radiative processes. A common, repeated pattern of transport has pollutant plumes arriving in the distant Pacific Ocean from Africa and South America at 25 degrees south and 14 km altitude. Following this, there is then a general appearance of pollution at extending down to 5 kin at more equatorial (10 S) and polar latitudes (to 45 S). Second, we evaluate the quantitative success of our simulation. (Such success requires efforts considerably beyond trajectory analyses, and is necessary for our community to claim an understanding of the effects of biomass burning on global atmospheric chemistry and the planet's trend in oxidizing capacity.) We find that we simulate most pollution episodes sampled by Glenn Sachse's CO instrument and the Blake hydrocarbon analyses during PENT A. We will present our current ideas on why our general levels appear satisfactory when the observations are within 20 ppb of background levels, but substantially miss the variability associated with the most intense biomass burning plumes. As time allows, we will speculate on the origins of a class of plumes that appear in the mid-troposphere, but are missing from our simulations.

Chatfield, Robert↗

NASA NTRS2002

GEWEX Cloud System Study (GCSS) Working Group on Cirrus Cloud Systems (WG2)

Status, progress and plans will be given for current GCSS (GEWEX Cloud System Study) WG2 (Working Group on Cirrus Cloud Systems) projects, including: (a) the Idealized Cirrus Model Comparison Project, (b) the Cirrus Parcel Model Comparison Project (Phase 2), and (c) the developing Hurricane Nora extended outflow model case study project. Past results will be summarized and plans for the upcoming year described. Issues and strategies will be discussed. Prospects for developing improved cloud parameterizations derived from results of GCSS WG2 projects will be assessed. Plans for NASA's CRYSTAL-FACE (Cirrus Regional Study of Tropical Anvils and Layers - Florida Area Cirrus Experiment) potential opportunities for use of those data for WG2 model simulations (future projects) will be briefly described.

Starr, David↗

NASA NTRS2007

Low-Level Cloud Variability over the Equatorial Cold Tongue in Observations and Models

A fourth paper now in press is, Low-level cloud variability over the equatorial cold tongue in observations and models, by D. K. Mansbach and J. R. Norris (2007, J. Climate). This study examined cloud and meteorological observations from satellite, surface, and reanalysis datasets and fount that monthly anomalies in low-level cloud amount and near-surface temperature advection are strongly negatively correlated on the southern side of the equatorial Pacific cold tongue. This inverse correlation occurs independently of relationships between cloud amount and sea surface temperature (SST) or lower tropospheric static stability (LTS) and the combination of advection plus SST or LTS explains significantly more interannual cloud variability in a multilinear regression than does SST or LTS alone. Warm anomalous advection occurs when the equatorial cold tongue is well defined and the southeastern Pacific trade winds bring relatively warm air over colder water. Ship meteorological reports and soundings show that the atmospheric surface layer becomes stratified under these conditions, thus inhibiting the upward mixing of moisture needed to sustain cloudiness against subsidence and entrainment drying. Cold anomalous advection primarily occurs when the equatorial cold tongue is weak or absent and the air-sea temperature difference is substantially negative. These conditions favor a more convective atmospheric boundary layer, greater cloud amount, and less frequent occurrence of clear sky. Examination of output from global climate models developed by the Geophysical Fluid Dynamics Laboratory (GFDL) and the National Center for Atmospheric Research (NCAR) indicates that both models generally fail to simulate the cloud-advection relationships observed on the northern and southern sides of the equatorial cold tongue. Although the GFDL atmosphere model does reproduce the expected signs of cloud-advection correlations when forced with prescribed historical SST variations, it does not consistently do so when coupled to an ocean model. The NCAR model has difficulty reproducing the observed correlations in both atmosphere-only and coupled versions. This suggests that boundary layer cloud parameterizations could be improved through better representation of the effects of advection over varying SST.

Mansbach, David K.↗

NASA NTRS2008

Retrieved Latent Heating from TRMM

The global hydrological cycle is central to the Earth's climate system, with rainfall and the physics of precipitation formation acting as the key links in the cycle. Two-thirds of global rainfall occurs in the tropics with the associated latent heating (LH) accounting for three-fourths of the total heat energy available to the Earth's atmosphere. 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. In the last decade, it has been established that standard products of LH from satellite measurements, particularly TRMM measurements, would be a valuable resource for scientific research and applications. Such products would enable new insights and investigations concerning the complexities of convection system life cycles, the diabatic heating controls and feedbacks related to meso-synoptic circulations and their forecasting, the relationship of tropical patterns of LH to the global circulation and climate, and strategies for improving cloud parameterizations in environmental prediction models. The status of retrieved TRMM LH products, TRMM LH inter-comparison and validation project, current TRMM LH applications and critic issues/action items (based on previous five TRMM LH workshops) is presented in this article.

Tao, Wei-Kuo↗

NASA NTRS2013

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗

NASA NTRS2018

Climates of Warm Earth-Like Planets. I. 3D Model Simulations

We present a large ensemble of simulations of an Earth-like world with increasing insolation and rotation rate. Unlike previous work utilizing idealized aquaplanet congurations we focus our simulations on modern Earth-like topography. The orbital period is the same as modern Earth, but with zero obliquity and eccentricity. The atmosphere is 1 bar N2-dominated with CO2=400 ppmv and CH4=1 ppmv. The simulations include two types of oceans; one without ocean heat transport (OHT) between grid cells as has been commonly used in the exoplanet literature, while the other is a fully coupled dynamic bathtub type ocean. The dynamical regime transitions that occur as day length increases induce climate feedbacks producing cooler temperatures, rst via the reduction of water vapor with increasing rotation period despite decreasing shortwave cooling by clouds, and then via decreasing water vapor and increasing shortwave cloud cooling, except at the highest insolations. Simulations without OHT are more sensitive to insolation changes for fast rotations while slower rotations are relatively insensitive to ocean choice. OHT runs with faster rotations tend to be similar with gyres transporting heat poleward making them warmer than those without OHT. For slower rotations OHT is directed equator-ward and no high latitude gyres are apparent. Uncertainties in cloud parameterization preclude a precise determination of habitability but do not a affect robust aspects of exoplanet climate sensitivity. This is the first paper in a series that will investigate aspects of habitability in the simulations presented herein. The datasets from this study are open source and publicly available.

Planetary systems↗

NASA NTRS2020

AWARE: The Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment

The US Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) performed comprehensive meteorological and aerosol measurements, and ground-based atmospheric remote sensing at two Antarctic stations using the most advanced instrumentation available. A suite of cloud research radars, lidars, spectral and broadband radiometers, aerosol chemical and microphysical sampling equipment, and meteorological instrumentation was deployed at McMurdo Station on Ross Island from December 2015 through December 2016. A smaller suite of radiometers and meteorological equipment including radiosondes, optimized for surface energy budget measurement, was deployed on the West Antarctic Ice Sheet between 4 December 2015 and 17 January 2016. AWARE provided Antarctic atmospheric data comparable to several well-instrumented high Arctic sites that have operated for many years and that reveal numerous contrasts with the Arctic in aerosol and cloud microphysical properties. These include persistent differences in liquid cloud occurrence, cloud height and cloud thickness. Antarctic aerosol properties are also quite different from the Arctic in both seasonal cycle and composition, due to the continent’s isolation from lower latitudes by Southern Ocean storm tracks. Antarctic aerosol number and mass concentrations are not only non-negligible but perhaps play a more important role than previously recognized because of the higher sensitivities of clouds at the very low concentrations caused by the large-scale dynamical isolation. Antarctic aerosol chemical composition, particularly organic components, has implications for local cloud microphysics. The AWARE data set, fully available online in the ARM Program data archive, offers numerous case studies for unique and rigorous evaluation of mixed-phase cloud parameterization in climate models.

Atmospheric Radiation Measurement (ARM)↗

NASA NTRS2005

Evaluation of a New Cloud Droplet Activation Parameterization wtih in Situ Data from CRYSTAL-FACE and CSTRIPE

The accuracy of the 2003 prognostic, physically based aerosol activation parameterization of A. Nenes and J. H. Seinfeld (NS) with modification introduced by C. Fountoukis and A. Nenes in 2005 (modified NS) is evaluated against extensive microphysical data sets collected on board the Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft for cumuliform and stratiform clouds of marine and continental origin. The cumuliform cloud data were collected during NASA's Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE, Key West, Florida, July 2002), while the stratiform cloud data were gathered during Coastal Stratocumulus Imposed Perturbation Experiment (CSTRIPE, Monterey, California, July 2003). In situ data sets of aerosol size distribution, chemical composition, and updraft velocities are used as input for the NS parameterization, and the evaluation is carried out by comparing predicted cloud droplet number concentrations (CDNC) with observations. This is the first known study in which a prognostic cloud droplet activation parameterization has been evaluated against a wide range of observations. On average, predicted droplet concentration in adiabatic regions is within -20% of observations at the base of cumuliform clouds and -30% of observations at different altitudes throughout the stratiform clouds, all within experimental uncertainty. Furthermore, CDNC is well parameterized using either a single mean updraft velocity w or by weighting droplet nucleation rates with a Gaussian probability density function of w. This study suggests that for nonprecipitating warm clouds of variable microphysics, aerosol composition, and size distribution the modified NS parameterization can accurately predict cloud droplet activation and can be successfully implemented for describing the aerosol activation process in global climate models.

Meskhidze, Nicholas↗

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↗