Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble simulation”

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.

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

At least 235 records · Page 13

Similarity Assessment of Land Surface Model Outputs in the North American Land Data Assimilation System

Multimodel ensembles are often used to produce ensemble mean estimates that tend to have increased simulation skill over any individual model output. If multimodel outputs are too similar, an individual LSM would add little additional information to the multimodel ensemble, whereas if the models are too dissimilar, it may be indicative of systematic errors in their formulations or configurations. The article presents a formal similarity assessment of the North American Land Data Assimilation System (NLDAS) multimodel ensemble outputs to assess their utility to the ensemble, using a confirmatory factor analysis. Outputs from four NLDAS Phase 2 models currently running in operations at NOAA/NCEP and four new/ upgraded models that are under consideration for the next phase of NLDAS are employed in this study. The results show that the runoff estimates from the LSMs were most dissimilar whereas the models showed greater similarity for root zone soil moisture, snow water equivalent, and terrestrial water storage. Generally, the NLDAS operational models showed weaker association with the common factor of the ensemble and the newer versions of the LSMs showed stronger association with the common factor, with the model similarity increasing at longer time scales. Trade-offs between the similarity metrics and accuracy measures indicated that the NLDAS operational models demonstrate a larger span in the similarity-accuracy space compared to the new LSMs. The results of the article indicate that simultaneous consideration of model similarity and accuracy at the relevant time scales is necessary in the development of multimodel ensemble.

land surface models↗

Seasonal Predictability of Cloud Droplet Number Concentration

Aerosol emissions modify the properties of clouds hence impacting climate. The aerosol indirect effect may have offset part of the global warming caused by anthropogenic greenhouse gas emissions during the industrial era. It however remains unclear whether the same effect is significant over time scales relevant for seasonal and weather climate prediction. Answering such a question has been difficult since most weather prediction systems lack a proper representation of the aerosol evolution and transport and their interaction with clouds. Even in advanced systems it is not clear to what extent cloud microphysical properties are predictable over subseasonal to seasonal time scales. Such an issue is addressed in this study. We use a set of 30 year, four ensemble member, 9 month lead hindcast simulations of the NASA GEOS seasonal prediction system (GEOS-S2S) to study the predictability of cloud droplet number concentration in warm stratocumulus clouds. The latest version GEOS-S2S system implements interactive aerosol as well as a two moment cloud microphysics scheme therefore it is suitable for studying the aerosol indirect effect on climate. Long term retrievals from the MODIS (Moderate Resolution Imaging Spectroradiometer) are used to validate the model predictions and assess its skill in predicting cloud droplet number concentration.

Barahona, Donifan↗

Mid-21st Century Ozone Air Quality and Health Burden in China Under Emissions Scenarios and Climate Change

Despite modest emissions reductions of air pollutants in recent years, China still suffers from poor air quality, and the outlook for future air quality in China is uncertain. We explore the impact of two disparate 2050 emissions scenarios relative to 2015 in the context of a changing climate with the Geophysical Fluid Dynamics Laboratory Atmospheric Model version 3 (GFDL-AM3) chemistry-climate model. We impose the same near-term climate change for both emission scenarios by setting global sea surface temperature (SST) and sea ice cover (SIC) to the average over 2010–2019 and 2046–2055, respectively, from a three-member ensemble of GFDL coupled climate model simulations under the RCP8.5 (Representative Concentration Pathway) scenario. By the 2050s, annual mean surface ozone increases throughout China by up to 8 ppbv from climate change alone (estimated by holding air pollutants at 2015 levels while setting SIC and SST to 2050 conditions in the model) and by 8–12 ppbv in a scenario in which emissions of ozone precursors nitrogen oxides (NO (sub x) ) and anthropogenic volatile organic compounds (VOCs) increase by ~10%. In a scenario in which NO (sub x) and anthropogenic VOC emissions decline by 60%, annual mean surface ozone over China decreases by 16–20 ppbv in the 2050s relative to the 2010s. The ozone increase from climate change alone results in an additional 62 000 premature deaths in China as compared to 330 000 fewer premature deaths by the 2050s under a strong emissions mitigation scenario. In springtime over Southwestern China in the 2050s, the model projects 9–12 ppbv enhancements to surface ozone from the stratosphere (diagnosed with a model tracer) and from international anthropogenic emissions (diagnosed by differencing AM3 simulations with the same emissions within China but higher versus lower emissions in the rest of the world). Our findings highlight the effectiveness of emissions controls in reducing the health burden in China due to air pollution, and also the potential for climate change and rising global emissions to offset, at least partially, some of the ozone decreases attained with regional emission reductions in China.

Westervelt, D. M.↗

Seasonal Variation of the Quasi-Biennial Oscillation Descent

This study examines variations through the year in the rate of descent of the Quasi-Biennial Oscillation (QBO) wind in the tropical lower stratosphere. The work is based on the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) multi-year global reanalysis. A 10-member ensemble of free-running atmospheric simulations using the Global Earth Observing System (GEOS) atmospheric general circulation model is also examined. The QBO winds are decomposed into two orthogonal modes that describe amplitude and phase propagation. This reveals a semi-annual signal in descent rate, with slow descent at times of large middle-latitude planetary wave activity, suggesting that this wave activity influences equatorial vertical velocity. This result also explains the suppression of QBO wind transitions during active wave times. Results also show a variation in descent speed when plotted against QBO phase. These results provide a baseline of QBO behavior for future QBO modeling.

Quasi-biennial oscillation↗

Impact of a Regional Drought on Land and Atmospheric Carbon over the North America

The spatiotemporal extent of a regional drought’s impact on land-atmosphere carbon and their mechanistic linkages are still poorly understood. Allowing carbon-water-energy feedback in an ensemble of coupled land-atmosphere simulations, this study quantifies induced changes in meteorology and carbon caused by an idealized 3-month meteorological drought over to the lower Mississippi River Valley (~500,000 km2). The imposed drought leads to a 23% Gross Primary Production (GPP) reduction in the drought area and the GPP reductions in some remote areas that are adjacent to the imposed drought through the changes in remote meteorology, particularly by the induced precipitation. The effect of the meteorological anomalies on GPP is greater, by at least an order of magnitude, than that of the atmospheric CO2 anomalies. While the magnitude of induced CO2 anomalies near the land surface can be as large as 3.57ppm, the column-averaged CO2 increases up to 0.78ppm and thus is below the precision limit (~1ppm) of current greenhouse gas (GHG) observing satellites. Spatially, the CO2 anomalies cover an area up to three times of that of the imposed drought, which suggests that the atmospheric transport needs to be considered in the interpretation of carbon anomalies in the atmosphere.

Eunjee Lee↗

Excessive Subsidence of the Southeast Atlantic Biomass Burning Aerosol Plume in the NASA Goddard Earth Observing System (GEOS) Model

Global models struggle to accurately represent the meteorological environment over the Southeast Atlantic Ocean, as evidenced by the poor representation of persistent marine stratocumulus decks off the west coast of southern Africa. This complex regime also experiences significant seasonal biomass burning aerosol plumes that are transported from southern Africa over the Atlantic Ocean, which gradually descend via large scale subsidence and interact radiatively and microphysically with the cloudy environment of the region. In comparisons to space borne lidar and aircraft observations we find that the NASA Goddard Earth Observing System (GEOS) model’s MERRA-2 aerosol and meteorology reanalysis shows that the aerosol plume descends too rapidly due to exaggerated large-scale downward vertical motion. This complicates the evaluation of aerosol properties utilizing aircraft observations collected throughout the ORACLES campaign as well as the model’s representation of the heating rate profiles due to the aerosol itself and its radiative interaction with the semi-permanent cloud deck typically located lower in altitude than the smoke. Two types of modeling experiments using GEOS will be presented to diagnose deficiencies in the model related to the excessive subsidence. Ensembles of free-running model simulations with and without radiatively active biomass burning aerosol will be used to quantify the impact of smoke on the atmospheric circulation. It will be demonstrated that heating within the aerosol layer results in self lofting through a reduction of the subsidence but did not alter the vertical placement or speed of the African Easterly jet. The amount of aerosol absorption required to eliminate the bias in the height of the aerosol plume through self-lofting will then be quantified. Finally, a second type of modeling experiment will be discussed in which the meteorology is constrained to observations. This will allow for biases in the height of the aerosol plume due to the general circulation model to be distinguished from biases that may be induced through the data assimilation process.

Allison Collow↗

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in A Mach 4.5 Boundary Layer

We use deep learning, an ensemble variationaltechnique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (Deep-ONets), which have the known ability to learn com-plex nonlinear operators within dynamical systems,are used for machine learning. For the baseline config-uration of a smooth flat plate, the second-mode wavesat the DNS inflow cause a quick nonlinear breakdownof the high-speed boundary layer within the computa-tional domain. Results reported in the present studyvalidate the ability of DeepONets to model the tran-sition delay via a given topography of the roughnesselement. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substan-tially lowered by the DeepONets-based reduced-ordermodel. In comparison to the baseline method of EnVaroptimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition pastthe outflow boundary of the computational domainwhile utilizing almost 5–6 times fewer DNS.

Machine Learning↗

Multiyear Dry Periods in Southern Africa

Characteristics and physical features related to low precipitation across many years in Southern Africa that lead to societal disruptions are diagnosed using observed analyses and an ensemble of historical coupled climate model simulations during 1921 to 2014. Four regions are evaluated, as identified through a hierarchical clustering algorithm applied to the Standardized Precipitation Index (SPI) during the October–April precipitation season. Although dryness spanning many October–April occurs periodically in each region, they seldom occur simultaneously, consistent with largely insignificant SPI cross-correlations between them. However, characteristics relevant to low precipitation across many years are generalizable between the four regions, including the serial persistence of October–April precipitation, the likelihood of consecutive dry October–April, and the likelihood of dry October–April in temporal extents of up to 10 consecutive such 7-month seasons. Systematic precipitation persistence is not a feature in any of the four Southern Africa regions, as serial correlations of October–April SPI are not statistically significant at any time lags. It follows that there is an exponential-folding decay in the likelihood of consecutive October–April for various SPI thresholds and that there is a large spread in the likelihood of low October–April SPI across many years. In terms of physical features, low October–April SPI in each Southern Africa region is closely related to local atmospheric circulations; however, they are not as closely related to sea surface temperatures (SSTs). These results suggest that dryness spanning many years is determined primarily by persistent local circulations related to atmospheric variability and to a lesser extent variability related to SST anomalies, including the El Niño–Southern Oscillation.

subtropical Indian Ocean dipole↗

Hierarchical Statistical 3D ' Atomistic' Simulation of Decanano MOSFETs: Drift-Diffusion, Hydrodynamic and Quantum Mechanical Approaches

When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain voltage only single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. A pilot version of a hydrodynamic 'atomistic' simulator has been developed in order to study the effect of the nonequilibrium, non local transport in decanano MOSFETs on the random dopant induced current fluctuations. For the first time we have also applied the density gradient approach in 3D to investigate the effect of the quantum confinement on the threshold voltage fluctuations. The developed 'atomistic' simulation techniques have been applied to study various fluctuation resistant MOSFET architectures including epitaxial and delta doped devices.

Asenov, Asen↗

Predictability of the 1997 and 1998 South Asian Summer Monsoons

The predictability of the 1997 and 1998 south Asian summer monsoon winds is examined from an ensemble of 10 Atmospheric General Circulation Model (AGCM) simulations with prescribed sea surface temperatures (SSTs) and soil moisture, The simulations are started in September 1996 so that they have lost all memory of the atmospheric initial conditions for the periods of interest. The model simulations show that the 1998 monsoon is considerably more predictable than the 1997 monsoon. During May and June of 1998 the predictability of the low-level wind anomalies is largely associated with a local response to anomalously warm Indian Ocean SSTs. Predictability increases late in the season (July and August) as a result of the strengthening of the anomalous Walker circulation and the associated development of easterly low level wind anomalies that extend westward across India and the Arabian Sea. During these months the model is also the most skillful with the observations showing a similar late-season westward extension of the easterly CD wind anomalies. The model shows little predictability or skill in the low level winds over southeast Asia during, 1997. Predictable wind anomalies do occur over the western Indian Ocean and Indonesia, however, over the Indian Ocean they are a response to SST anomalies that were wind driven and they show no skill. The reduced predictability in the low level winds during 1997 appears to be the result of a weaker (compared with 1998) simulated anomalous Walker circulation, while the reduced skill is associated with pronounced intraseasonal activity that is not well captured by the model. Remarkably, the model does produce an ensemble mean Madden-Julian Oscillation (MJO) response that is approximately in phase with (though weaker than) the observed MJ0 anomalies. This is consistent with the idea that SST coupling may play an important role in the MJO.

Schubert, Siegfred D.↗

Precipitation processes developed during TOGA COARE (1992), GATE (1974), SCSMEX (1998), and KWAJEX (1999): 3D Cloud Resolving Model Simulation

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research (NCAR), NOAA GFDL, the U.K. Met. Office, Colorado State University and NASA Goddard Space Flight Center. An improved 3D Goddard Cumulus Ensemble (GCE) model was recently used to simulate periods during TOGA COARE (December 19-27, 1992), GATE (september 1-7, 1974), SCSMEX (May 18-26, June 2-11, 1998) and KWAJEX (August 7-13, August 18-21, and August 29-September 12, 1999) using a 512 by 512 km domain and 41 vertical layers. The major objectives of this paper are: (1) to identify the differences and similarities in the simulated precipitation processes and their associated surface and water energy budgets in TOGA COARE, GATE, KWAJEX, and SCSMEX, and (2) to asses the impact of microphysics, radiation budget and surface fluxes on the organization of convection in tropics.

Tao, W.-K.↗

On the Influence of North Pacific Sea Surface Temperature on the Arctic Winter Climate

Differences between two ensembles of Goddard Earth Observing System Chemistry-Climate Model simulations isolate the impact of North Pacific sea surface temperatures (SSTs) on the Arctic winter climate. One ensemble of extended winter season forecasts is forced by unusually high SSTs in the North Pacific, while in the second ensemble SSTs in the North Pacific are unusually low. High Low differences are consistent with a weakened Western Pacific atmospheric teleconnection pattern, and in particular, a weakening of the Aleutian low. This relative change in tropospheric circulation inhibits planetary wave propagation into the stratosphere, in turn reducing polar stratospheric temperature in mid- and late winter. The number of winters with sudden stratospheric warmings is approximately tripled in the Low ensemble as compared with the High ensemble. Enhanced North Pacific SSTs, and thus a more stable and persistent Arctic vortex, lead to a relative decrease in lower stratospheric ozone in late winter, affecting the April clear-sky UV index at Northern Hemisphere mid-latitudes.

Hurwitz, Margaret M.↗

The PMIP4 Last Glacial Maximum experiments: preliminary results and comparison with the PMIP3 simulations

The Last Glacial Maximum (LGM, ∼ 21 000 years ago) has been a major focus for evaluating how well state-of-the-art climate models simulate climate changes as large as those expected in the future using paleoclimate reconstructions. A new generation of climate models has been used to generate LGM simulations as part of the Paleoclimate Modelling Intercomparison Project (PMIP) contribution to the Coupled Model Intercomparison Project (CMIP). Here, we provide a preliminary analysis and evaluation of the results of these LGM experiments (PMIP4, most of which are PMIP4-CMIP6) and compare them with the previous generation of simulations (PMIP3, most of which are PMIP3-CMIP5). We show that the global averages of the PMIP4 simulations span a larger range in terms of mean annual surface air temperature and mean annual precipitation compared to the PMIP3-CMIP5 simulations, with some PMIP4 simulations reaching a globally colder and drier state. However, the multi-model global cooling average is similar for the PMIP4 and PMIP3 ensembles, while the multi-model PMIP4 mean annual precipitation average is drier than the PMIP3 one. There are important differences in both atmospheric and oceanic circulations between the two sets of experiments, with the northern and southern jet streams being more poleward and the changes in the Atlantic Meridional Overturning Circulation being less pronounced in the PMIP4-CMIP6 simulations than in the PMIP3-CMIP5 simulations. Changes in simulated precipitation patterns are influenced by both temperature and circulation changes. Differences in simulated climate between individual models remain large. Therefore, although there are differences in the average behaviour across the two ensembles, the new simulation results are not fundamentally different from the PMIP3-CMIP5 results. Evaluation of large-scale climate features, such as land–sea contrast and polar amplification, confirms that the models capture these well and within the uncertainty of the paleoclimate reconstructions. Nevertheless, regional climate changes are less well simulated: the models underestimate extratropical cooling, particularly in winter, and precipitation changes. These results point to the utility of using paleoclimate simulations to understand the mechanisms of climate change and evaluate model performance.

Last Glacial Maximum↗

Crop Model Improvement Reduces the Uncertainty of the Response to Temperature of Multi-Model Ensembles

To improve climate change impact estimates and to quantify their uncertainty, multi-model ensembles (MMEs) have been suggested. Model improvements can improve the accuracy of simulations and reduce the uncertainty of climate change impact assessments. Furthermore, they can reduce the number of models needed in a MME. Herein, 15 wheat growth models of a larger MME were improved through re-parameterization and/or incorporating or modifying heat stress effects on phenology, leaf growth and senescence, biomass growth, and grain number and size using detailed field experimental data from the USDA Hot Serial Cereal experiment (calibration data set). Simulation results from before and after model improvement were then evaluated with independent field experiments from a CIMMYT worldwide field trial network (evaluation data set). Model improvements decreased the variation (10th to 90th model ensemble percentile range) of grain yields simulated by the MME on average by 39% in the calibration data set and by 26% in the independent evaluation data set for crops grown in mean seasonal temperatures greater than 24 C. MME mean squared error in simulating grain yield decreased by 37%. A reduction in MME uncertainty range by 27% increased MME prediction skills by 47%. Results suggest that the mean level of variation observed in field experiments and used as a benchmark can be reached with half the number of models in the MME. Improving crop models is therefore important to increase the certainty of model-based impact assessments and allow more practical, i.e. smaller MMEs to be used effectively.

damage assessment↗

Simulated KWAJEX Convective Systems Using a 2D and 3D Cloud Resolving Model and Their Comparisons with Radar Observations

The 1999 Kwajalein Atoll field experiment (KWAJEX), one of several major TRMM (Tropical Rainfall Measuring Mission) field experiments, has successfully obtained a wealth of information and observation data on tropical convective systems over the western Central Pacific region. In this paper, clouds and convective systems that developed during three active periods (Aug 7-12, Aug 17-21, and Aug 29-Sep 13) around Kwajalein Atoll site are simulated using both 2D and 3D Goddard Cumulus Ensemble (GCE) models. Based on numerical results, the clouds and cloud systems are generally unorganized and short lived. These features are validated by radar observations that support the model results. Both the 2D and 3D simulated rainfall amounts and their stratiform contribution as well as the heat, water vapor, and moist static energy budgets are examined for the three convective episodes. Rainfall amounts are quantitatively similar between the two simulations, but the stratiform contribution is considerably larger in the 2D simulation. Regardless of dimension, fo all three cases, the large-scale forcing and net condensation are the two major physical processes that account for the evolution of the budgets with surface latent heat flux and net radiation solar and long-wave radiation)being secondary processes. Quantitative budget differences between 2D and 3D as well as between various episodes will be detailed.Morover, simulated radar signatures and Q1/Q2 fields from the three simulations are compared to each other and with radar and sounding observations.

Shie, Chung-Lin↗

Predictability of Zonal Means During Boreal Summer

This study examines the predictability of seasonal means during boreal summer. The results are based on ensembles of June-July-August (JJA) simulations (started in mid May) carried out with the NASA Seasonal-to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed sea surface temperatures (SSTS) and sea ice for the years 1980-1999. We find that the predictability of the JJA extra-tropical height field is primarily in the zonal mean component of the response to the SST anomalies. This contrasts with the cold season (January-February-March) when the predictability of seasonal means in the boreal extratropics is primarily in the wave component of the El Nino/Southern Oscillation (ENSO) response. Two patterns dominate the interannual variability of the ensemble mean JJA zonal mean height field. One has maximum variance in the tropical/subtropical upper troposphere, while the other has substantial variance in middle latitudes of both hemispheres. Both are symmetric with respect to the equator. A regression analysis suggests that the tropical/subtropical pattern is associated with SST anomalies in the far eastern tropical Pacific and the Indian Ocean, while the middle latitude pattern is forced by SST anomalies in the tropical Pacific just east of the dateline. The two leading zonal height patterns are reproduced in model runs forced with the two leading JJA SST patterns of variability. A comparison with observations shows a signature of the middle latitude pattern that is consistent with the occurrence of dry and wet summers over the United States. We hypothesize that both patterns, while imposing only weak constraints on extratropical warm season continental-scale climates, may play a role in the predilection for drought or pluvial conditions.

Schubert, Siegfried↗

Numerical Simulations of TRMM LBA, TOGA, COARE, GATE, ARM and PRESTORM Convective Systems: Sensitivity tests on Microphysical Processes

The 3D Goddard Cumulus Ensemble (GCE) model was utilized to examine the behavior and response of simulated deep tropical cloud systems that occurred over the west Pacific warm pool region, the Atlantic ocean and the central United States. The periods chosen for simulation were convectively active periods during TOGA-COARE (February 22 1993, December 11-17, 1992; December 19-28, February 9-13, 1993), GATE (September 4, 1974), LBA (January 26 and February 23, 1998), ARM (1997 IOP) and PRESTORM (June 11, 1985). We will examine differences in the microphysics for both warm rain and ice processes (evaporation /sublimation and condensation/ deposition), Q1 (Temperature), Q2 (Water vapor) and Q3 (momentum both U and V) budgets for these three convective events from different large-scale environments. The contribution of stratiform precipitation and its relationship to the vertical shear of the large-scale horizontal wind will also be examined. New improvements to the GCE model (i.e., microphysics: 4ICE two moments and 3ICE one moment; advection schemes) as well as their sensitivity to the model results will be discussed. Preliminary results indicated that various microphysical schemes could have a major impact on stratiform formation as well as the size of convective systems. However, they do not change the major characteristics of the convective systems, such as: arc shape, strong rotational circulation on both ends of system, heavy precipitation along the leading edge of systems.

Tao, W.-K.↗