Search NASA⌕ Search

SEARCH · Search NASA

Results for “height of different cloud layers”

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.

101 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

An airborne study of boundary layer heat and moisture fluxes for project FIFE

A summary of work accomplished during the first International Satellite Cloud Climatology Project (ISLSCP) Field Experiment is presented. The Wyoming King Air participated in Intensive Field Campaigns (IFCs) 3 and 4 in 1987, which were committed to a combination of airborne sensible and latent heat flux measurements and soil moisture mappings, with the University of Kansas' X-band side-looking airborne radar (SLAR). A total of 9 flux missions were flown in 1987 using several different flight designs. A critical decision in the first round of aircraft flux analyses, agreed to by all the flux aircraft investigators, was to pass the aircraft data through a high-pass filter prior to the eddy-correlation flux calculations. Several conclusions were drawn from this study: (1) boundary layer (BL) profiles of heat fluxes were usually linear; (2) high-pass filtering applied to the aircraft data did not add to the disagreement between the profile; (3) undersampling at high frequencies could have accounted for as much as a 15 percent underestimate of surface fluxes; (4) disagreement between the aircraft and surface latent heat fluxes changed signs between summer and fall; (5) magnitude and sign of disagreement between aircraft and surface heat fluxes varied systematically with BL depth or height of aircraft profile; and (6) Bowen ratios from both detrended and filtered aircraft data agreed with surface values better for moist, summer cases than for dry cases.

Kelly, Robert D.↗
NASA NTRS2000

Sensing Water Vapon via Spacecraft Radio Occultation Observations

The radio occultation technique has been used to characterize planetary atmospheres since the 1960's spanning atmospheric pressures from 16 microbars to several bars. In 1988, the use of GPS signals to make occultation observations of Earth's atmosphere was realized by Tom Yunck and Gunnar Lindal at JPL. In the GPS to low-Earth-orbiter limb- viewing occultation geometry, Fresnel diffraction yield a unique combination of high vertical resolution of 100 m to 1 km at long wavelengths (approx. 20 cm) insensitive to particulate scattering which allows routine limb sounding from the lower mesosphere through the troposphere. A single orbiting GPS/GLONASS receiver can observe - 1000 to 1400 daily occultations providing as many daily, high vertical resolution soundings as the present global radiosonde network, but with far more evenly distributed, global coverage. The occultations yield profiles of refractivity as a function of height. In the cold, dry conditions of the upper troposphere and above (T less than 240 K), profiles of density, pressure (geopotential), and temperature can be derived. Given additional temperature information, water vapor can be derived in the midddle and lower troposphere with a unique combination of vertical resolution, global distribution and insensitivity to clouds and precipitation to an accuracy of approx. 0.2 g/kg. At low latitudes, moisture profiles will be accurate to 1-5% within the convective boundary layer and better than 20% below 6 to 7 km. Accuracies of climatological averages should be approx. 0. 1 g/kg limited by the biases in the temperature estimates. To use refractivity to constrain water vapor, knowledge of temperature is required. The simplest approach is to use the temperature field from an analysis such as the 6 hour ECMWF global analysis interpolated to the locations of each occultation. A better approach is to combine the temperature and moisture fields from such an analysis with the occultation refractivity in a weighting scheme based on the errors in each data field. A ID variational combinational approach has been developed at the UKMO. We win present results from both approaches from GPS/MET data taken in June and July 1995 and compare them with the ECMWF global 6 hour moisture analyses which are derived largely from TOVS and radiosonde data. Overall, the atmosphere below the 500 mb level appears somewhat drier in general than the ECNIWF humidity field. A 2-D (latitude vs. height) climatological snapshot derived from a 2-week span of GPS/MET data will be compared to the humidity climatology of Peixoto and Oort derived from radiosonde data from 1963-1973. Differences between the GPS results and Peixoto and Oort may be the signature of a climate trend over the past 30 years.

Kursinski, E. Robert↗
NASA NTRS2012

Application of the CALIOP Layer Product to Evaluate the Vertical Distribution of Aerosols Estimated by Global Models: AeroCom Phase I Results

The CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) layer product is used for a multimodel evaluation of the vertical distribution of aerosols. Annual and seasonal aerosol extinction profiles are analyzed over 13 sub-continental regions representative of industrial, dust, and biomass burning pollution, from CALIOP 2007-2009 observations and from AeroCom (Aerosol Comparisons between Observations and Models) 2000 simulations. An extinction mean height diagnostic (Z-alpha) is defined to quantitatively assess the models' performance. It is calculated over the 0-6 km and 0-10 km altitude ranges by weighting the altitude of each 100 m altitude layer by its aerosol extinction coefficient. The mean extinction profiles derived from CALIOP layer products provide consistent regional and seasonal specificities and a low inter-annual variability. While the outputs from most models are significantly correlated with the observed Z-alpha climatologies, some do better than others, and 2 of the 12 models perform particularly well in all seasons. Over industrial and maritime regions, most models show higher Z-alpha than observed by CALIOP, whereas over the African and Chinese dust source regions, Z-alpha is underestimated during Northern Hemisphere Spring and Summer. The positive model bias in Z-alpha is mainly due to an overestimate of the extinction above 6 km. Potential CALIOP and model limitations, and methodological factors that might contribute to the differences are discussed.

Aerosols and particles↗
NASA NTRS2001

The Vertical Error Characteristics of GOES-derived Winds: Description and Impact on Numerical Weather Prediction

Errors in the height assignment of some satellite-derived winds exist because the satellites sense radiation emitted from a finite layer of the atmosphere rather than a specific level. Potential problems in data assimilation may arise because the motion of a measured layer is often represented by a single-level value. In this research, cloud and water vapor motion winds that are derived from the Geostationary Operational Environmental Satellites (GOES winds) are compared to collocated rawinsonde observations (RAOBs). An important aspect of this work is that in addition to comparisons at each assigned height, the GOES winds are compared to the entire profile of the collocated RAOB data to determine the vertical error characteristics of the GOES winds. The impact of these results on numerical weather prediction is then investigated. The comparisons at individual vector height assignments indicate that the error of the GOES winds range from approx. 3 to 10 m/s and generally increase with height. However, if taken as a percentage of the total wind speed, accuracy is better at upper levels. As expected, comparisons with the entire profile of the collocated RAOBs indicate that clear-air water vapor winds represent deeper layers than do either infrared or water vapor cloud-tracked winds. This is because in cloud-free regions the signal from water vapor features may result from emittance over a thicker layer. To further investigate characteristics of the clear-air water vapor winds, they are stratified into two categories that are dependent on the depth of the layer represented by the vector. It is found that if the vertical gradient of moisture is smooth and uniform from near the height assignment upwards, the clear-air water vapor wind tends to represent a relatively deep layer. The information from the comparisons is then used in numerical model simulations of two separate events to determine the forecast impacts. Four simulations are performed for each case: 1) A control simulation that assimilates no satellite wind data, 2) assimilation of all GOES winds according to their assigned single level height, 3) assimilation of all GOES winds spread over multiple levels, and 4) assimilation of all GOES winds spread over multiple levels, but with variations in the vertical influence of clear-air water vapor winds based on the moisture profile in the model. In the first case, a strong mid-latitude cyclone is present and the use of the satellite data results in improved storm tracks during the initial approx. 36 h forecast period. This is because the satellite data improves the analysis of the environment into which the storm progresses. Statistics for mean wind vector and height differences show that, with the exception of the height field at later times in the first case, the use of GOES winds improves the simulation with time. The simulation results suggest that it is beneficial to spread the GOES wind information over multiple levels, particularly when the moisture profile is used to define the vertical influence.

Rao, P. Anil↗
NASA NTRS2020

The Evolution of Lidar Networks: A US Perspective

Atmospheric aerosols and clouds play an important role in climate by directly scattering and absorbing sunlight. Aerosol-cloud interactions modify the particle properties causing indirect effects, alter aerosol deposition and rainfall, and contribute substantially to the uncertainties in predicting climate effects. Aerosols also affect air quality and both constituents modulate boundary layer dynamics to a certain extent, in turn impacting aerosol transport and aerosol-cloud interactions. The spatiotemporal distribution of aerosol and cloud layers is thus important, and lidar remains the primary instrument for determining the vertical distribution of aerosols and thin clouds. Lidar still provides important information for opaque cloud layers by determining base heights. Combined lidar and radar data provide a comprehensive coverage of aerosol and cloud vertical distributions, and thus locations where and when they interact. Lidar data provide a means of constraining column aerosol loading observations (e.g. AOD) to vertical extents. In addition, lidar has proven effective at determining a proxy for boundary layer height by examining aerosol gradients and/or cloud base heights in lidar profiles. Polarized lidars provide additional information on particle shape, allowing estimates of cloud phase and separation of dust or smoke from sulfate and sea salt aerosols. Given the importance of lidar observations and a critical maturation in technology and retrieval techniques, several organizations began to build operational lidar networks around 2000. The European Aerosol Research Lidar Network (EARLINET) was started as a research focused network of advanced lidars across Europe. The Asian Dust Lidar Network (ADNET) was also developed as a regional network providing lidar profiles of dust and pollution across Eastern Asia. In the US, the NASA Micro Pulse Lidar Network (MPLNET) was created to provide global lidar profiling at key sites in the NASA Aerosol Robotic Network (AERONET). The Network for the Detection of Atmospheric Composition Change (NDACC) pre-dates these networks and many sites have lidar, but it is not strictly a lidar network and at the time focused less on the lower troposphere. Finally, there were already existing ceilometer networks operated by various meteorological agencies, but in particular here in the US the profile data has not been available. The ceilometer networks were used to provide only clouds base heights and estimates of PBL height. Thus, the distinguishing feature between lidar and ceilometer networks was historically the ability to actually provide profile information (in addition to differences in wavelength, and advanced retrievals such as the raman technique). As time progressed, each lidar network matured and developed more operational capabilities and data sets, coupled with viable data centers providing DAAC services and access to near-real-time (NRT) data. In 2008 under WMO guidance, the Global Atmospheric Watch (GAW) Aerosol Lidar Observation Network (GALION) was formed as a global network made up of the existing lidar networks. The goal was to share information, best practices, and develop frameworks and techniques for quality data. GALION grew to include several other regional lidar networks and this has led to a vast increase in quality lidar sites worldwide.

Welton, Ellsworth J.↗
NASA NTRS2021

A Dark Target research aerosol algorithm for MODIS observations over eastern China: increasing coverage while maintaining accuracy at high aerosol loading

Satellite aerosol products such as the Dark Target (DT) produced from the MODerate resolution Imaging Spectroradiometer (MODIS) are useful for monitoring the progress of air pollution. Unfortunately, the DT often fails to retrieve during the heaviest aerosol events as well as the more moderate events in winter. Some of the literature at-tributes this lack of retrieval to the cloud mask. However, we found this lack of retrieval is mainly traced to thresholds used for masking of inland water and snow. Modifications to these two masks greatly increase 50 % of the retrievals of aerosol optical depth at 0.55 μm (AOD) greater than 1.0. The “extra”-high-AOD retrievals tend to be biased when com-pared with a ground-based sun photometer (AErosol RObotic NETwork, AERONET). Reducing bias in new retrievals re-quires two additional steps. One is an update to the assumed aerosol optical properties (aerosol model); the haze in this region is both less absorbing and lower in altitude than what is assumed in the global algorithm. The second is account-ing for the scale height of the aerosol, specifically that the heavy-aerosol events in the region are much closer to the surface than what is assumed by the global DT algorithm. The resulting combination of modified masking thresholds, new aerosol model, and lower aerosol layer scale height was applied to 3 months of MODIS observations (January–March2013) over eastern China. After these two additional steps are implemented, the significant increase in new retrievals introduces no overall bias at a high-AOD regime but does degrade other overall validation statistics. We also find that the research algorithm is able to identify additional pollution events that AERONET instruments may not due to different spatial sampling. Mean AOD retrieved from the re-search algorithm increases from 0.11 to 0.18 compared to values calculated from the operational DT algorithm during January to March of 2013 over the study area. But near Beijing, where the severe pollution occurs, the new algorithm increases AOD by as much as 3.0 for each 0.5°grid box over the previous operational-algorithm values.

Dark Target↗
NASA NTRS2023

Earth-System-Model Evaluation of Cloud and Precipitation Occurrence for Supercooled and Warm Clouds Over the Southern Ocean's Macquarie Island

Over the remote Southern Ocean (SO), cloud feedbacks contribute substantially to Earth system model (ESM) radiative biases. The evolution of low Southern Ocean clouds (cloud-top heights < ∼ 3 km) is strongly modulated by precipitation and/or evaporation, which act as the primary sink of cloud condensate. Constraining precipitation processes in ESMs requires robust observations suitable for process-level evaluations. A year-long subset (April 2016–March 2017) of ground-based profiling instrumentation deployed during the Macquarie Island Cloud and Radiation Experiment (MICRE) field campaign (54.5∘ S, 158.9∘ E) combines a 95 GHz (W-band) Doppler cloud radar, two lidar ceilometers, and balloon-borne soundings to quantify the occurrence frequency of precipitation from the liquid-phase cloud base. Liquid-based clouds at Macquarie Island precipitate ∼ 70 % of the time, with deeper and colder clouds precipitating more frequently and at a higher intensity compared to thinner and warmer clouds. Supercooled cloud layers precipitate more readily than layers with cloud-top temperatures > 0 ∘C, regardless of the geometric thickness of the layer, and also evaporate more frequently. We further demonstrate an approach to employ these observational constraints for evaluation of a 9-year GISS-ModelE3 ESM simulation. Model output is processed through the Earth Model Column Collaboratory (EMC2) radar and lidar instrument simulator with the same instrument specifications as those deployed during MICRE, therefore accounting for instrument sensitivities and ensuring a coherent comparison. Relative to MICRE observations, the ESM produces a smaller cloud occurrence frequency, smaller precipitation occurrence frequency, and greater sub-cloud evaporation. The lower precipitation occurrence frequency by the ESM relative to MICRE contrasts with numerous studies that suggest a ubiquitous bias by ESMs to precipitate too frequently over the SO when compared with satellite-based observations, likely owing to sensitivity limitations of spaceborne instrumentation and different sampling methodologies for ground- versus space-based observations. Despite these deficiencies, the ESM reproduces the observed tendency for deeper and colder clouds to precipitate more frequently and at a higher intensity. The ESM also reproduces specific cloud regimes, including near-surface clouds that account for ∼ 25 % of liquid-based clouds during MICRE and optically thin, non-precipitating clouds that account for ∼ 27 % of clouds with bases higher than 250 m. We suggest that the demonstrated framework, which merges observations with appropriately constrained model output, is a valuable approach to evaluate processes responsible for cloud radiative feedbacks in ESMs.

cloud feedbacks↗
NASA NTRS1993

A satellite-borne radar wind sensor (RAWS)

Modeling global atmospheric circulations and forecasting the weather would improve if worldwide information on winds aloft were available. Accurate prediction of weather is important to agriculture, shipping, air traffic, and many other fields. Global system models of climate are of great importance. Current global atmospheric models use pressure measurements and thermodynamic properties to calculate the effects of wind for use in Numerical Weather Prediction (NWP) models. Inputs to the NWP models are temperature, pressure and wind velocities at different heights. Clearly direct wind measurements could significantly improve the NWP model performance. The RAdar Wind Sounder (RAWS) program at the University of Kansas is a study of the feasibility and the trade-offs in the design of a space-based radar system to measure wind vectors. This can be done by measuring the Doppler shift of cloud and rain returns from three or more points and calculating the components of the wind vector. The RAWS study to date uses the candidate system selected after preliminary study of frequencies and sensitivities. Two frequencies chosen, 10 and 35 GHz, allow higher sensitivity for clouds and more penetration for rain. The past year was devoted to modeling the signal-to-noise ratio (SNR) achievable for the two frequencies. The determination of SNR versus cloud penetration depth used a cloud backscattering and attenuation model in the appropriate radar equation. Calculations assumed reasonable losses in reception and transmission, in addition to the atmospheric attenuation. We discovered that ice clouds provide a higher SNR than previously calculated, but some water clouds give lower SNRs than we calculated before. One of the primary issues in the SNR calculation was the choice of the drop size distribution. Although Xin used several distributions (e.g., log normal, Khrigian and Mazin), this year we used the Deirmendjian cloud model. SNR versus cloud penetration plots were generated to validate the candidate system. Rain, which appears in the cloud models at the lower altitudes, provides ample SNR, as do the higher clouds composed of ice particles. However, in some cloud situations we found the sensitivity for the clouds was marginal or inadequate. At 35 GHz, two of the cloud models characterized by 1 to 2 g/cu m of water content at altitudes extending from 150 to 1500 meters, produced a sufficient SNR. Other models, however, with water contents ranging from 0.5 to 4 g/cu m and altitudes up to 4000 meters, exhibit SNR of -3 to -23 dB, largely because of attenuation in the upper cloud layers. These results coupled with the lower SNR at 10 GHz, led to an investigation of alternate frequencies. The rain present beneath these clouds provides adequate SNR at 10 GHz, and in most cases, at GHz.

Moore, Richard K.↗
NASA NTRS2007

Electrification in Hurricanes over the Tropical Americas: Implication for Stratospheric Water Vapor

This study explores the relation between lightning activity and water vapor in the Tropical Tropopause Layer (TTL) over hurricane systems in the Tropical Americas. The hypothesis herein is that hurricanes that exhibit enhanced lightning activity are associated with stronger updrafts that can transport more moisture directly into the TTL (and subsequently into the tropical stratosphere) or even directly into the tropical stratosphere over this region. The TTL over the Tropical Americas, which includes the Caribbean and Gulf of Mexico, is of particular interest, because summertime cold point tropopause is the lowest in height and thus the warmest in temperature over the tropics. The latter condition implies higher saturation values and thus potential for more water vapor to enter the stratosphere. Climate forecast is very sensitive to stratospheric water vapor abundance, because of the key role that water vapor plays in regulating the chemical and radiative properties of the stratosphere. Given the potential for increases in hurricane intensity and frequency under predicted warmer conditions, it becomes essential to understand the effect of hurricanes on stratospheric water vapor. In this study, we use a combination of ground and space-borne observations as well as trajectory calculations. The observations include: cloud-to-ground (CG) lightning data from the U.S. National Lightning Detection Network (NLDN), geostationary infrared observations from the National Climatic Data Center Hurricane Satellite (HURSAT) data set, cloud properties from Aqua-MODIS, and water vapor from Aura-MLS. We analyze hurricanes from the 2005 season when Aura-MLS data are available, namely: Dennis, Emily, Katrina, Rita, and Wilma. Our analysis consists of examining CG lightning, cloud-top properties, and TTL water vapor (i.e., 100 and 147 mb) over the hurricane while it remains over water in the Tropical Americas region. We investigate daily as well as diurnal statistical properties. The hurricanes analyzed in this study showed that lightning activity is negatively correlated with minimum infrared brightness temperature and positively correlated with 100-mb water vapor. An examination of the maxima in water vapor observed over the hurricane not only shows larger magnitudes, but also larger differences between water vapor averages and water vapor maxima over the hurricane as lightning activity increases. Trajectory calculations are performed using the Flextra model in order to investigate the fate of the moister air masses found in the TTL.

Pittman, Jasna V.↗
NASA NTRS2012

Evaluating Nighttime CALIOP 0.532 micron Aerosol Optical Depth and Extinction Coefficient Retrievals

NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) Version 3.01 5-km nighttime 0.532 micron aerosol optical depth (AOD) datasets from 2007 are screened, averaged and evaluated at 1 deg X 1 deg resolution versus corresponding/co-incident 0.550 micron AOD derived using the US Navy Aerosol Analysis and Prediction System (NAAPS), featuring two-dimensional variational assimilation of quality-assured NASA Moderate Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) AOD. In the absence of sunlight, since passive radiometric AOD retrievals rely overwhelmingly on scattered radiances, the model represents one of the few practical global estimates available from which to attempt such a validation. Daytime comparisons, though, provide useful context. Regional-mean CALIOP vertical profiles of night/day 0.532 micron extinction coefficient are compared with 0.523/0.532 micron ground-based lidar measurements to investigate representativeness and diurnal variability. In this analysis, mean nighttime CALIOP AOD are mostly lower than daytime (0.121 vs. 0.126 for all aggregated data points, and 0.099 vs. 0.102 when averaged globally per normalised 1 deg. X 1 deg. bin), though the relationship is reversed over land and coastal regions when the data are averaged per normalised bin (0.134/0.108 vs. 0140/0.112, respectively). Offsets assessed within single bins alone approach +/- 20 %. CALIOP AOD, both day and night, are higher than NAAPS over land (0.137 vs. 0.124) and equal over water (0.082 vs. 0.083) when averaged globally per normalised bin. However, for all data points inclusive, NAAPS exceeds CALIOP over land, coast and ocean, both day and night. Again, differences assessed within single bins approach 50% in extreme cases. Correlation between CALIOP and NAAPS AOD is comparable during both day and night. Higher correlation is found nearest the equator, both as a function of sample size and relative signal magnitudes inherent at these latitudes. Root mean square deviation between CALIOP and NAAPS varies between 0.1 and 0.3 globally during both day/night. Averaging of CALIOP along-track AOD data points within a single NAAPS grid bin improves correlation and RMSD, though day/night and land/ocean biases persist and are believed systematic. Vertical profiles of extinction coefficient derived in the Caribbean compare well with ground-based lidar observations, though potentially anomalous selection of a priori lidar ratios for CALIOP retrievals is likely inducing some discrepancies. Mean effective aerosol layer top heights are stable between day and night, indicating consistent layer-identification diurnally, which is noteworthy considering the potential limiting effects of ambient solar noise during day.

aerosol optical depth↗
NASA NTRS2004

Florida Thunderstorms: A Faucet of Reactive Nitrogen to the Upper Troposphere

During the NASA Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) enhanced mixing ratios of nitric oxide were measured in the anvils of thunderstorms and in clear air downwind of storm systems on flights of a Wl3-57F high-altitude aircraft. Mixing ratios greater than l0 - 20 times background were readily observed over distances of 25-120 km due to lightning activity. In many of the Florida storms deposition of NO occurred up to near the tropopause but major deposition usually occurred 1 - 2 km below the tropopause, or mostly within the visible anvil volume formed prior to storm decay. Observations from two storms of very different anvil size and electrical activity allowed estimates of the total mass of NO, vented to the middle and upper troposphere. Using the cloud-to ground (CG) flash accumulations from the National Lightning Detection Network, climatological intra-cloud (IC) to CG ratios, and assuming that CG and IC flashes were of equivalent efficiency for NO production, the ranges of production per flash for a moderate-sized and a large storm were (0.51 - 1.0) x l0(exp 26) and (2.3 - 3.1) x 10(exp 26) molecules NO/flash, respectively. Using the recently determined average global flash rate of 44 8, a gross extrapolation of these two storms to represent possible global annual production rates yield 1.6 - 3.2 and 7.3 - 9.9 Tg(N)/yr, respectively. If the more usual assumption is made that IC efficiency is l/l0th that of CG activity, the ranges of production for the moderate-sized and large storm were (1.3 - 2.7) x l0(exp 26) and (6.0 - 8.1) x l0(exp 26) molecules NO/CG flash, respectively. The estimates from the large storm may be high because there is indirect evidence that the IC/CG ratio was larger than would be derived from climatology. These two storms and others studied did not have flash rates that scaled as approx. H(sup 5) where H is the cloud top altitude. The observed CG flash accumulations and NO(x) mass production estimate for the month of July over the Florida area were compared with a representative 3D global Chemistry-Transport Model (CTMJ that uses the Price et al. lightning parameterization. For two land grid points representing the Florida peninsula the model compared well with the observations: CG flash rates were low by only a factor of approx. 2. When the model grid points included the coastal regions of Florida the flash accumulations were lower than observed by a factor of 3.4 - 4.6. It is recommended that models using the Price et al. parameterization allow any global coastal grid point to maintain the land rather than the marine flash rate parameterization. The convection in this CTM underestimated the actual cloud top heights over Florida by 1 - 2 km and thus the total lightning flash rates and the altitude range of reactive nitrogen deposition. Broad scale (20 - 120 km) median mixing ratios of NO within anvils over Florida were significantly larger than in storms previously investigated over Colorado and New Mexico.

Ridley, B.↗