Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convective storms”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Nighttime observations of thunderstorm electrical activity from a high altitude airplane

Nocturnal thunderstorms were observed from above and features of cloud structure and lightning which are not generally visible from the ground are discussed. Most, lightning activity seems to be associated with clouds with strong convective cauliflower tops. In both of the storms lightning channels were visible in the clear air above the cloud. It is shown that substances produced by thunderstorm electrical discharges can be introduced directly into the stratosphere. The cause and nature of the discharges above the cloud are not clear. They may be produced by accumulations of space charge in the clear air above the cloud. The discharges may arise solely because of the intense electric fields produced by charges within the cloud. In the latter case the ions introduced by these discharges will increase the electrical conductivity of the air above the cloud and increase the conduction current that flows from the cloud to the electrosphere. More quantitative data at higher resolution may show significant spectral differences between cloud to ground and intracloud strokes. It is shown that electric field change data taken with an electric field change meter mounted in an airplane provide data on lightning discharges from above that are quite similar to those obtained from the ground in the past. The optical signals from dart leaders, from return strokes, and from continuing currents are recognizable, can be used to provide information on the fine structure of lightning, and can be used to distinguish between cloud to ground and intracloud flashes.

Brook, M.↗
NASA NTRS1989

Goes I-M and beyond: Science requirements and technology challenges

Measurement of Earth observation parameters from geosynchronous satellites; severe local storm observational guidelines; major instrument requirements; rapid convective cell growth; geosynchronous satellite requirements; and major technological challenges for the next 25 years are discussed.

Shenk, William E.↗
NASA NTRS1995

Observations of the Global Characteristics and Regional Radiative Effects of Marine Cloud Liquid Water

The large-scale spatial distribution and temporal variability of cloud liquid water path (LWP) over the world's oceans and the relationship of cloud LWP to temperature and the radiation budget are investigated using recent satellite measurements from the Special Sensor Microwave/Imager (SSM/I), the Earth Radiation Budget Experiment (ERBE), and the International Satellite Cloud Climatology Project (ISCCP). Observations of cloud liquid water on a 2.5 deg x 2.5 deg and are used over a 53-month period beginning July 1987 and ending in December 1991. The highest values of cloud liquid water (greater than 0.13 kg/sq m) occur largely along principal routes of northern midlatitude storms and in areas dominated by tropical convection. The zonally averaged structure is distinctly trimodal, where maxima appear in the midlatitudes and near the equator. The average marine cloud LWP over the globe is estimated to be about 0.113 kg/sq m. Its highest seasonal variability is typically between 15% and 25% of the annual mean but in certain locations can exceed 30%. Comparisons of cloud LWP to temperature for low clouds during JJA and DJF of 1990 show significant positive correlations at colder temperatures and negative correlations at warmer temperatures. The correlations also exhibit strong seasonal and regional variation. Coincident and collocated observations of cloud LWP from the SSM/I and albedo measurements from the Earth Radiation Budget Satellite (ERBS) and the NOAA-10 satellite are compared for low clouds in the North Pacific and North Atlantic. The observed albedo-LWP relationships correspond reasonably well with theory, where the average cloud effective radius (r(sub e)) is 11.1 microns and the standard deviation is 5.2 microns. The large variability in the inferred values of r(sub e) suggests that other factors may be important in the albedo-LWP relationships. In terms of the effect of the LWP on the net cloud forcing, the authors find that a 0.05 kg/sq m increase in LWP (for LWP less than 0.2 kg/sq m) results in a -25 W/sq m change in the net cloud forcing at a solar zenith angle of 75 deg.

Greenwald, Thomas J.↗
NASA NTRS2017

Utilizing the Precessing Orbit of TRMM to Produce Hourly Corrections of Geostationary Infrared Imager Data with the VIRS Sensor

Accurate characterization of the Earth's radiant energy is critical for many climate monitoring and weather forecasting applications. For example, groups at the NASA Langley Research Center rely on stable visible- and infrared-channel calibrations in order to understand the temporal/spatial distribution of hazardous storms, as determined from an automated overshooting convective top detection algorithm. Therefore, in order to facilitate reliable, climate-quality retrievals, it is important that consistent calibration coefficients across satellite platforms are made available to the remote sensing community, and that calibration anomalies are recognized and mitigated. One such anomaly is the infrared imager brightness temperature (BT) drift that occurs for some Geostationary Earth Orbit satellite (GEOsat) instruments near local midnight. Currently the Global Space-Based Inter-Calibration System (GSICS) community uses the hyperspectral Infrared Atmospheric Sounding Interferometer (IASI) sensor as a common reference to uniformly calibrate GEOsat IR imagers. However, the combination of IASI, which has a 21:30 local equator crossing time (LECT), and hyperspectral Atmospheric Infrared Sounder (AIRS; 01:30 LECT) observations are unable to completely resolve the GEOsat midnight BT bias. The precessing orbit of the Tropical Rainfall Measuring Mission (TRMM) Visible and Infrared Scanner (VIRS), however, allows sampling of all local hours every 46 days. Thus, VIRS has the capability to quantify the BT midnight effect observed in concurrent GEOsat imagers. First, the VIRS IR measurements are evaluated for long-term temporal stability between 2002 and 2012 by inter-calibrating with Aqua-MODIS. Second, the VIRS IR measurements are assessed for diurnal stability by inter-calibrating with Meteosat-9 (Met-9), a spin-stabilized GEOsat imager that does not manifest any diurnal dependency. In this case, the Met-9 IR imager is first adjusted with the official GSICS calibration coefficients. Then VIRS is used as a diurnal calibration reference transfer to produce hourly corrections of GEOsat IR imager BT. For the 9 three-axis stabilized GEO imagers concurrent with VIRS, the midnight effect increased the BT on average by 0.5 K (11 microns) and 0.4 K (12 microns), with a peak at approx.01:00 local time. As expected, the spin-stabilized GEOsats revealed a smaller diurnal temperature cycle (mostly < 0.2 K) with inconsistent peak hours.

Scarino, Benjamin↗
DOE OSTI2023

A semi-Lagrangian method for detecting and tracking deep convective clouds in geostationary satellite observations

Automated methods for the detection and tracking of deep convective clouds in geostationary satellite imagery have a vital role in both the forecasting of severe storms and research into their behaviour. Studying the interactions and feedbacks between multiple deep convective clouds (DCC), however, poses a challenge for existing algorithms due to the necessary compromise between false detection and missed detection errors. We utilise an optical flow method to determine the motion of deep convective clouds in GOES-16 ABI imagery in order to construct a semi-Lagrangian framework for the motion of the cloud field, independently of the detection and tracking of cloud objects. The semi-Lagrangian framework allows severe storms to be simultaneously detected and tracked in both spatial and temporal dimensions. For the purpose of this framework we have developed a novel Lagrangian convolution method and a number of novel implementations of morphological image operations that account for the motion of observed objects. These novel methods allow the accurate extension of computer vision techniques to the temporal domain for moving objects such as DCCs. By combining this framework with existing methods for detecting DCCs (including detection of growing cores through cloud top cooling and detection of anvil clouds using brightness temperature), we show that the novel framework enables reductions in errors due to both false and missed detections compared to any of the individual methods, reducing the need to compromise when compared with existing frameworks. The novel framework enables the continuous tracking of anvil clouds associated with detected deep convection after convective activity has stopped, enabling the study of the entire life cycle of DCCs and their associated anvils. Furthermore, we expect this framework to be applicable to a wide range of cases including the detection and tracking of low-level clouds and other atmospheric phenomena. In addition, this framework may be used to combine observations from multiple sources, including satellite observations, weather radar and reanalysis model data.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

Roles of divergent and rotational winds in the kinetic energy balance during intense convective activity

Contributions of divergent and rotational wind components to the synoptic-scale kinetic energy balance are described using rawinsonde data at 3 and 6 h intervals from NASA's fourth Atmospheric Variability experiment. Two intense thunderstorm complexes occurred during the period. Energy budgets are described for the entire computational region and for limited volumes that enclosed storm-induced, upper level wind maxima located poleward of convection. Although small in magnitude, the divergent wind component played an important role in the cross-contour generation and horizontal flux divergence of kinetic energy. The importance of V(D) appears directly related to the presence and intensity of convection. Although K(D) usually comprised less than 10 percent of the total kinetic energy content, generation of kinetic energy by V(D) was a major factor in the creation of upper-level wind maxima to the north of the storm complexes. Omission of the divergent wind apparently would lead to serious misrepresentations of the energy balance. A random error analysis is presented to assess confidence limits in the various energy parameters.

Fuelberg, H. E.↗
NASA NTRS2015

Tropical Cyclone Diurnal Cycle as Observed by TRMM

Using infrared satellite data, previous work has shown a consistent diurnal cycle in the pattern of cold cloud tops around mature tropical cyclones. In particular, an increase in the coverage by cold cloud tops often occurs in the inner core of the storm around the time of sunset and subsequently propagates outward to several hundred kilometers over the course of the following day. This consistent cycle may have important implications for structure and intensity changes of tropical cyclones and the forecasting of such changes. Because infrared satellite measurements are primarily sensitive to cloud top, the goal of this study is to use passive and active microwave measurements from the Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) and Precipitation Radar (PR), respectively, to examine and better understand the tropical cyclone diurnal cycle throughout a larger depth of the storm's clouds. The National Hurricane Center's best track dataset was used to extract all PR and TMI pixels within 1000 km of each tropical cyclone that occurred in the Atlantic basin between 1998-2011. Then the data was composited according to radius (100-km bins from 0-1000 km) and local standard time (LST; 3-hr bins). Specifically, PR composites involved finding the percentage of pixels with reflectivity greater than or equal to 20 dBZ at various heights (i.e., 2-14 km in increments of 2 km) as a function of radius and time. The 37- and 85- GHz TMI channels are especially sensitive to scattering by precipitation-sized ice in the mid to upper portions of clouds. Hence, the percentage of 37- and 85-GHz polarization corrected temperatures less than various thresholds were calculated using data from all storms as a function of radius and time. For 37 GHz, thresholds of 260 K, 265 K, 270 K, and 275 K were used, and for 85 GHz, thresholds of 200-270 K in increments of 10 K were utilized. Note that convection forced by the interactions of a tropical cyclone with land (e.g., due to frictional convergence) may disrupt the natural convective cycle of a cyclone. Hence, only data pertaining to storms whose centers were greater than 300 km from land were included in the composites. Early results suggest the presence of a diurnal cycle in the PR composites of all Atlantic basin tropical cyclones from a height of 2-12 km from approximately 0-400 km radius, but the cycle is most apparent above 6 km. At a height of 8 km, there is a peak (minimum) in the percentage of PR pixels greater than or equal to 20 dBZ near 0 (21) LST in the inner core with some indication that this signal propagates outward with time. In contrast, the 37- and 85-GHz composites show little indication of a diurnal cycle at any radii, regardless of the threshold used. Ongoing work with this project will involve sub-setting the composites according to storm intensity to see if the diurnal cycle varies with storm strength. Moderate to strong vertical wind shear often leads to asymmetries in tropical cyclone convection and may disrupt the cyclone's natural diurnal cycle. Therefore, wind shear thresholds will be applied to the composites to determine if the diurnal cycle becomes more apparent in a low shear environment. Finally, other work to be completed will involve developing composites for other tropical cyclone basins, including the East Pacific, Northwest Pacific, South Pacific, and Indian Ocean.

tropical cyclones↗
NASA NTRS1984

A kinetic energy study of the meso beta-scale storm environment during AVE-SESAME 5 (20-21 May 1979)

Kinetic energy of the near storm environment was analyzed by meso beta scale data. It was found that horizontal winds in the 400 to 150 mb layer strengthen rapidly north of the developing convection. Peak values then decrease such that the maximum disappears 6 h later. Southeast of the storms, wind speeds above 300 mb decrease nearly 50% during the 3 h period of most intense thunderstorm activity. When the convection dissipates, wind patterns return to prestorm conditions. The mesoscale storm environment of AVE-SESAME 5 is characterized by large values of cross contour generation of kinetic energy, transfers of energy to nonresolvable scales of motion, and horizontal flux divergence. These processes are maximized within the upper troposphere and are greatest during times of strongest convection. It is shown that patterns agree with observed weather features. The southeast area of the network is examined to determine causes for vertical wind variations.

Printy, M. F.↗
NASA NTRS1995

Temperature measurements of a Martian local dust storm

A technique for estimating the ground and near-ground atmospheric temperatures within a Martian local dust storm is presented. It is applied to soundings taken by the Viking orbiter infrared thermal mapper (IRTM) instrument at four times-of-day for one storm. Essentially, a comparison is made between infrared radiances emerging from the storm interior and those from the region surrounding the storm. Particle extinction properties are assumed to be independent of position in the storm region, and scattering properties must be selected arbitrarily. For the storm studied here, the ground temperature in the interior is at least 6 K cooler, whereas the near-ground atmospheric temperature may be less than or comparable to, those of the surroundings. The thermal structure of the storm interior did not change measurably between 11.5 and 16.6 hours local time. These observations favor the theories of dust storm development in which regional winds rather than local, dust driven convection initiate the mobilization of dust from the surface. It also concluded that the optical properties of dust particles in this local storm differ from those observed by Mariner 9 during the 1971-1972 global dust storm.

Kahn, Ralph↗
NASA NTRS1999

Global Lightning Variations Caused by Changes in Flash Rate and by Changes in Number of Thunderstorms

Observations from the Optical Transient Detector in space are used to examine the response of the number of storms and the total flash rate to thermodynamic forcing of deep convection on the diurnal and annual time scales. On the diurnal time scale, the flash rate appears to track with cloud buoyancy, and storm numbers dominate later in the cycle. On the annual time scale, almost all of the total lightning response is due to changes in the number of storms, with a negligible contribution from flash rate.

Williams, E.↗
DOE OSTI2026

Differences in Lightning—Storm Microphysical Relationships in Subtropical South America Versus the Continental U.S. and Their Impact on Model Lightning Parameterizations

Satellite studies have shown that some of the deepest convection with high lightning flash rates occurs downstream of the Andes in subtropical South America (SSA). Using high resolution ground‐based observations in SSA, we explore the relationship between lightning and radar‐inferred storm microphysics. Results show these storms have significantly larger 35‐dBZ volumes and more graupel per lightning flash compared to storms in the U.S. We then test microphysics‐based lightning parameterizations for storms in SSA on a 1‐km convection‐permitting simulation of a mesoscale convective system. Results from the simulation match well with observed flash rates. If parameterizations based on U.S. storms had been applied, flash rates would have been overestimated by a factor of eight. These results suggest that storms in SSA are fundamentally different compared to storms in the U.S., and caution should be exercised when applying microphysics‐based lightning parameterizations in different regions of the world.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

Observations of Two Sprite-Producing Storms in Colorado

Two sprite-producing thunderstorms were observed on 8 and 25 June 2012 in northeastern Colorado by a combination of low-light cameras, a lightning mapping array, polarimetric and Doppler radars, the National Lightning Detection Network, and charge moment change measurements. The 8 June event evolved from a tornadic hailstorm to a larger multicellular system that produced 21 observed positive sprites in 2 h. The majority of sprites occurred during a lull in convective strength, as measured by total flash rate, flash energy, and radar echo volume. Mean flash area spiked multiple times during this period; however, total flash rates still exceeded 60 min(sup 1), and portions of the storm featured a complex anomalous charge structure, with midlevel positive charge near 20degC. The storm produced predominantly positive cloud-to-ground lightning. All sprite-parent flashes occurred on the northeastern flank of the storm, where strong westerly upper level flow was consistent with advection of charged precipitation away from convection, providing a pathway for stratiform lightning. The 25 June event was another multicellular hailstorm with an anomalous charge structure that produced 26 positive sprites in less than 1 h. The sprites again occurred during a convective lull, with relatively weaker reflectivity and lower total flash rate but relatively larger mean flash area. However, all sprite parents occurred in or near convection and tapped charge layers in adjacent anvil cloud. The results demonstrate the sprite production by convective ground strokes in anomalously charged storms and also indicate that sprite production and convective vigor are inversely related in mature storms.

Radar↗
NASA NTRS1991

Spatial distribution of visible lightning on Jupiter

An analysis of Voyager images of Jupiter is used to evaluate the horizontal spatial distribution of visible lightning over most of one hemisphere. Essentially all the detectable activity is confined to very narrow latitude bands at 13.5 N and 49 N. The active regions at 49 N are the brightest, most numerous, and periodic in longitude. Activity at this latitude is long-lived and is most likely associated with moist convective regions deep in Jupiter's atmosphere. The longitudinal periodicity of the lightning storms may represent the effects of a planetary scale atmospheric wave trapped at the depth of the moist convection.

Magalhaes, J. A.↗
NASA NTRS2004

The GOES-R Lightning Mapper Sensor

The Lightning Mapper Sensor on GOES-R builds on previous measurements of lightning from low earth orbit by the OTD (Optical Transient Detector) and LIS (Lightning Imaging Sensor) sensors. Unlike observations from low earth orbit, the GOES-R platform will allow continuous monitoring of lightning activity over the Continental United States and southern Canada, Central and South America, and portions of the Atlantic and Pacific Oceans. The LMS will detect total (cloud-to-ground and intracloud) lightning at storm scale resolution (approx. 8 km) using a highly sensitive Charge Coupled Device (CCD) detector array. Discrimination between lightning optical transients and a bright sunlit background scene is accomplished by employing spectral, spatial, and temporal filtering along with a background subtraction technique. The result is 24 hour detection capability of total lightning. These total lightning observations can be made available to users within about 20 seconds. Research indicates a number of ways that total lightning observations from LMS could benefit operational activities, including 1) potential increases in lead times and reduced false alarms for severe thunderstorm and tornado Warnings, 2) improved routing of &rail around thunderstorms, 3) support for spacecraft launches and landings, 4) improved ability to monitor tropical cyclone intensity, 5) ability to monitor thunderstorm intensification/weakening during radar outages or where radar coverage is poor, 6) better identification of deep convection for the initialization of numerical prediction models, 7) improved forest fire forecasts, 8) identification of convective initiation, 9) identification of heavy convective snowfall, and 10) enhanced temporal resolution of storm evolution (1 minute) than is available from radar observations. Total lightning data has been used in an operational environment since July 2003 at the Huntsville, Alabama National Weather Service office. Total lightning measurements are obtained by the North Alabama Lightning Mapping Array (LMA) and have successfully been used in warning decisions. Every 2 minutes, total lightning counts in 2 km by 2 km horizontal, 1 km vertical grids are available to forecasters on an AWIPS (Advanced Weather Interactive Processing System) workstation. Experience with the LMA total lightning data is used to illustrate the potential use of LMS data that would be available to forecasters across the US. This abstract is for submission as a presentation to the National Weather Association Annual Meeting to be held 16-21 October 2004 in Portland, OR. This abstract will be published in the conference proceedings.

Buechler, Dennis↗
NASA NTRS2000

Kinematic and Microphysical Structures of Hurricane Bob (1991) Determined from a 1.3-km-Resolution Numerical Simulation

Hurricane Bob (1991) is simulated using the Penn State/NCAR mesoscale model MM5. The simulation is conducted for a 24-h period at 4-km resolution and for a 6-h period at 1.3-km resolution. The 4-km simulation is able to fairly realistically capture the intensity and structure of the storm. The 1.3-km simulation depicts very small-scale convective structures and produces a convective band outside of the eye wall that did not occur in the 4-km simulation. The 1.3-km results are used to characterize several kinematic and cloud microphysical structures in the storm. Characteristics of air parcels flowing into and rising within the eye wall will be examined through trajectory calculations.

Braun, Scott A.↗
DOE OSTI2024

Parameterization of the Elevated Convection With a Unified Convection Scheme (UNICON) and Its Impacts on the Diurnal Cycle of Precipitation

To improve the simulation of nocturnal precipitation, we develop a parameterization for the elevated convection that is launched from the level of maximum grid-mean moist static energy, when grid-mean vertical flow is upward at the launching interface. The parameterized elevated convection is forced by both cold pool-driven and partially resolved external mesoscale organized flows. Properties of the external mesoscale flow are estimated from grid-mean vertical velocity and three dimensional advection tendencies of temperature and moisture. The new parameterization is implemented into a unified convection scheme (UNICON) and tested for both the single-column case at the Southern Great Plain (SGP) site in US and in global simulations. At the SGP site, the parameterized elevated convection strengthens nocturnal convection, increases nocturnal precipitation, and better simulates the observed diurnal cycle of precipitation. It appears that the elevated and surface-based convections interact with each other by stabilizing the atmospheric column, which affects subsequent convection. Global simulation shows that the elevated convection mostly occurs over the continents during the night, and also over the oceanic mid-latitude warm air advection and storm track regions during summer. Without degrading global mean climate, the elevated convection improves the simulation of nocturnal precipitation in the mid-US but simulates somewhat strong nocturnal precipitation in northeastern Asia. It seems that a key to simulate observed nocturnal precipitation is to appropriately parameterize the impacts of external organized flow on the elevated convection.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

A Study of Precipitation Mean Particle Size with the GPM Combined Radar-Radiometer Algorithm

The Global Precipitation Measurement (GPM) mission provides quantitative estimates of the mass-weighted mean diameter, Dm, a characteristic parameter of the precipitation particle size distribution (PSD). This study investigates the three-dimensional distribution of Dm retrieved by the GPM Combined Radar-Radiometer Algorithm (CORRA) using observations from GPM’s dual-frequency precipitation radar and microwave imager. We examine Dm in individual precipitation systems in different climate regimes and investigate a six-year global climatology. The vertical structure of Dm is demonstrated with cases of deep convection, frontal rain and snow, and stratocumulus light rain. The Dm values, detectable by GPM, range from appx 0.7 mm in stratocumulus precipitation to >3.5 mm in the ice layers of intense convection. Within the constraint of the 12-dBZ detectability threshold, the smallest climatological mean Dm (~ 0.8 mm) are found in the eastern oceans, and the largest values (~ 2 mm) occur above the melting levels in convection over land in summer. The standard deviation is generally < 0.45 mm below 6 km. Climate regimes are characterized with Dm annual/seasonal variations, its convective/stratiform components, and vertical variabilities (2-10 km). The US Central Plains and Argentina are associated with the largest Dm in a deep layer. Tropical Africa has larger Dm and standard deviation than Amazon. Large convective Dm occurs at high latitudes of Eurasia and North America in summer; the southern hemisphere high latitudes have shallower systems with smaller Dm. Oceanic storm tracks in both hemispheres have relatively large Dm, particularly for convective Dm in winter.

Mei Han↗
NASA NTRS2021

Understanding the Global Three-Dimensional Distribution of Precipitation Mean Particle Size with the Global Precipitation Measurement Mission

This study addresses the global distribution of precipitation mean particle size using data from the Global Precipitation Measurement (GPM) mission. The mass-weighted mean diameter Dm is a characteristic parameter of the precipitation particle size distribution (PSD), estimated from the GPM Combined Radar–Radiometer Algorithm (CORRA) using data from GPM’s dual-frequency precipitation radar and microwave imager. We examine Dm in individual precipitation systems in different climate regimes and investigate a 6-yr (2014–20) global climatology within 70°N–70°S. The vertical structure of Dm is demonstrated with cases of deep convection, frontal rain and snow, and stratocumulus light rain. The Dm values, detectable by GPM, range from ~0.7 mm in stratocumulus precipitation to >3.5 mm in the ice layers of intense convection. Within the constraint of the 12-dBZ detectability threshold, the smallest annual mean Dm (~0.8 mm) are found in the eastern oceans, and the largest values (~2 mm) occur above the melting levels in convection over land in summer. The standard deviation of the annual mean is generally <0.45 mm below 6 km. Climate regimes are characterized with Dm annual/seasonal variations, its convective/stratiform components, and vertical variabilities (2–10 km). The U.S. Central Plains and Argentina are associated with the largest Dm in a deep layer. Tropical Africa has larger Dm and standard deviation than Amazon. Large convective Dm occurs at high latitudes of Eurasia and North America in summer; the Southern Hemisphere high latitudes have shallower systems with smaller Dm. Oceanic storm tracks in both hemispheres have relatively large Dm, particularly for convective Dm in winter. Relatively small Dm occurs over tropical oceans, including ITCZ, requiring further investigation.

Mei Han↗