Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “clear-sky”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

A New ENSO Index Derived from Satellite Measurements of Column Ozone

Column Ozone measured in tropical latitudes from Nimbus 7 total ozone mapping spectrometer (TOMS), Earth Probe TOMS, solar backscatter ultraviolet (SBUV), and Aura ozone monitoring instrument (OMI) are used to derive an El Nino-Southern Oscillation (ENSO) index. This index, which covers a time period from 1979 to the present, is defined as the Ozone ENSO Index (OEI) and is the first developed from atmospheric trace gas measurements. The OEI is constructed by first averaging monthly mean column ozone over two broad regions in the western and eastern Pacific and then taking their difference. This differencing yields a self-calibrating ENSO index which is independent of individual instrument calibration offsets and drifts in measurements over the long record. The combined Aura OMI and MLS ozone data confirm that zonal variability in total column ozone in the tropics caused by ENSO events lies almost entirely in the troposphere. As a result, the OEI can be derived directly from total column ozone instead of tropospheric column ozone. For clear-sky ozone measurements a +1K change in Nino 3.4 index corresponds to +2.9 Dobson Unit (DU) change in the OEI, while a +1 hPa change in SOI coincides with a -1.7DU change in the OEI. For ozone measurements under all cloud conditions these numbers are +2.4DU and -1.4 DU, respectively. As an ENSO index based upon ozone, it is potentially useful in evaluating climate models predicting long term changes in ozone and other trace gases.

Ziemke, J. R.↗
NASA NTRS2011

Use and Limitations of a Climate-Quality Data Record to Study Temperature Trends on the Greenland Ice Sheet

Enhanced melting of the Greenland Ice Sheet has been documented in recent literature along with surface-temperature increases measured using infrared satellite data since 1981. Using a recently-developed climate-quality data record, 11- and 12-year trends in the clear-sky ice-surface temperature (IST) of the Greenland Ice Sheet have been studied using the Moderate-Resolution Imaging Spectroradiometer (MODIS) IST product. Daily and monthly MODIS ISTs of the Greenland Ice Sheet beginning on 1 March 2000 and continuing through 31 December 2010 are now available at 6.25-km spatial resolution on a polar stereographic grid as described in Hall et al. (submitted). This record will be elevated in status to a climate-data record (CDR) when more years of data become available either from the MODIS on the Terra or Aqua satellites, or from the Visible Infrared Imager Radiometer Suite (VIIRS) to be launched in October 2011. Maps showing the maximum extent of melt for the entire ice sheet and for the six major drainage basins have been developed from the MODIS IST dataset. Twelve-year trends of the duration of the melt season on the ice sheet vary in different drainage basins with some basins melting progressively earlier over the course of the study period. Some (but not all) of the basins also show a progressively-longer duration of melt. IST 12-year trends are compared with in-situ data, and climate data from the Modern Era Retrospective-Analysis for Research and Applications (MERRA) Reanalysis.

Hall, Dorothy K.↗
NASA NTRS2011

Discrepancies Between MODIS and ISCCP Land Surface Temperature Products Analyzed with Microwave Measurements

This paper compares land surface temperature (LST) products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the International Satellite Cloud Climatology Project (ISCCP). With both sources, the LST data are derived from infrared measurements. For ISCCP, LST is a secondary product in support of the primary cloud analyses, but the LST data have been used for several other purposes. The MODIS measurements from the Aqua spacecraft are taken at about 01:30 and 13:30 local time, and the ISCCP three-hourly data, based on several geostationary and polar orbiting satellites. were interpolated to the MODIS measurement times. For July 2003 monthly averages over all clear-sky locations, the ISCCP-MODIS differences were +5.0 K and +2.5 K for day and night, respectively, and there were areas with differences as large as 25 K. The day-night differences were as much as approximately 10 K higher for ISCCP than for MODIS. The MODIS measurements were more consistent with independent microwave measurements from AMSR-E, by several measures, with respect to day-night differences and day-to-day variations.

Moncet, Jean-Luc↗
NASA NTRS2011

Analysis of Co-Located MODIS and CALIPSO Observations Near Clouds

The purpose of this paper is to help researchers combine data from different satellites and thus gain new insights into two critical yet poorly understood aspects of anthropogenic climate change, aerosol-cloud interactions and aerosol radiative effects, For this, the paper explores whether cloud information from the Aqua satellite's MODIS instrument can help characterize systematic aerosol changes near clouds by refining earlier perceptions of these changes that were based on the CALIPSO satellite's CALIOP instrument. Similar to a radar but using visible and ncar-infrared light, CALIOP sends out laser pulses and provides aerosol and cloud information along a single line that tracks the satellite orbit by measuring the reflection of its pulses. In contrast, MODIS takes images of reflected sunlight and emitted infrared radiation at several wavelengths, and covers wide areas around the satellite track. This paper analyzes a year-long global dataset covering all ice-free oceans, and finds that MODIS can greatly help the interpretation of CALIOP observations, especially by detecting clouds that lie outside the line observed by CALlPSO. The paper also finds that complications such as differences in view direction or clouds drifting in the 72 seconds that elapse between MODIS and CALIOP observations have only a minor impact. The study also finds that MODIS data helps refine but does not qualitatively alter perceptions of the systematic aerosol changes that were detected in earlier studies using only CALIOP data. It then proposes a statistical approach to account for clouds lying outside the CALIOP track even when MODIS cannot as reliably detect low clouds, for example at night or over ice. Finally, the paper finds that, because of variations in cloud amount and type, the typical distance to clouds in maritime clear areas varies with season and location. The overall median distance to clouds in maritime clear areas around 4-5 km. The fact that half of all clear areas is closer than 5 km to clouds implies that pronounced near-cloud changes in aerosol properties have significant implications for overall clear-sky characteristics, including the radiative impact of aerosols.

Varnai, Tamas↗
NASA NTRS2012

Application of Artificial Neural Networks to the Development of Improved Multi-Sensor Retrievals of Near-Surface Air Temperature and Humidity Over Ocean

Improved estimates of near-surface air temperature and air humidity are critical to the development of more accurate turbulent surface heat fluxes over the ocean. Recent progress in retrieving these parameters has been made through the application of artificial neural networks (ANN) and the use of multi-sensor passive microwave observations. Details are provided on the development of an improved retrieval algorithm that applies the nonlinear statistical ANN methodology to a set of observations from the Advanced Microwave Scanning Radiometer (AMSR-E) and the Advanced Microwave Sounding Unit (AMSU-A) that are currently available from the NASA AQUA satellite platform. Statistical inversion techniques require an adequate training dataset to properly capture embedded physical relationships. The development of multiple training datasets containing only in-situ observations, only synthetic observations produced using the Community Radiative Transfer Model (CRTM), or a mixture of each is discussed. An intercomparison of results using each training dataset is provided to highlight the relative advantages and disadvantages of each methodology. Particular emphasis will be placed on the development of retrievals in cloudy versus clear-sky conditions. Near-surface air temperature and humidity retrievals using the multi-sensor ANN algorithms are compared to previous linear and non-linear retrieval schemes.

Roberts, J. Brent↗
NASA NTRS2011

Enhancing a Simple MODIS Cloud Mask Algorithm for the Landsat Data Continuity Mission

The presence of clouds in images acquired by the Landsat series of satellites is usually an undesirable, but generally unavoidable fact. With the emphasis of the program being on land imaging, the suspended liquid/ice particles of which clouds are made of fully or partially obscure the desired observational target. Knowing the amount and location of clouds in a Landsat scene is therefore valuable information for scene selection, for making clear-sky composites from multiple scenes, and for scheduling future acquisitions. The two instruments in the upcoming Landsat Data Continuity Mission (LDCM) will include new channels that will enhance our ability to detect high clouds which are often also thin in the sense that a large fraction of solar radiation can pass through them. This work studies the potential impact of these new channels on enhancing LDCM's cloud detection capabilities compared to previous Landsat missions. We revisit a previously published scheme for cloud detection and add new tests to capture more of the thin clouds that are harder to detect with the more limited arsenal channels. Since there are no Landsat data yet that include the new LDCM channels, we resort to data from another instrument, MODIS, which has these bands, as well as the other bands of LDCM, to test the capabilities of our new algorithm. By comparing our revised scheme's performance against the performance of the official MODIS cloud detection scheme, we conclude that the new scheme performs better than the earlier scheme which was not very good at thin cloud detection.

Wilson, Michael J.↗
NASA NTRS2012

A Study of Global Cirrus Cloud Morphology with AIRS Cloud-clear Radiances (CCRs)

Version 6 (V6) AIRS cloud-clear radiances (CCR) are used to derive cloud-induced radiance (Tcir=Tb-CCR) at the infrared frequencies of weighting functions peaked in the middle troposphere. The significantly improved V 6 CCR product allows a more accurate estimation of the expected clear-sky radiance as if clouds are absent. In the case where strong cloud scattering is present, the CCR becomes unreliable, which is reflected by its estimated uncertainty, and interpolation is employed to replace this CCR value. We find that Tcir derived from this CCR method are much better than other methods and detect more clouds in the upper and lower troposphere as well as in the polar regions where cloud detection is particularly challenging. The cloud morphology derived from the V6 test month, as well as some artifacts, will be shown.

Wu, Dong L.↗
NASA NTRS2012

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.↗
NASA NTRS2011

Comparative Analysis of Aerosol Retrievals from MODIS, OMI and MISR Over Sahara Region

MODIS is a wide field-of-view sensor providing daily global observations of the Earth. Currently, global MODIS aerosol retrievals over land are performed with the main Dark Target algorithm complimented with the Deep Blue (DB) Algorithm over bright deserts. The Dark Target algorithm relies on surface parameterization which relates reflectance in MODIS visible bands with the 2.1 micrometer region, whereas the Deep Blue algorithm uses an ancillary angular distribution model of surface reflectance developed from the time series of clear-sky MODIS observations. Recently, a new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm has been developed for MODIS. MAIAC uses a time series and an image based processing to perform simultaneous retrievals of aerosol properties and surface bidirectional reflectance. It is a generic algorithm which works over both dark vegetative surfaces and bright deserts and performs retrievals at 1 km resolution. In this work, we will provide a comparative analysis of DB, MAIAC, MISR and OMI aerosol products over bright deserts of northern Africa.

Lyapustin, A.↗
NASA NTRS2012

Using Satellites to Investigate the Sensitivity of Longwave Downward Radiation to Water Vapor at High Elevations

Many studies suggest that high-elevation regions may be among the most sensitive to future climate change. However, in situ observations in these often remote locations are too sparse to determine the feedbacks responsible for enhanced warming rates. One of these feedbacks is associated with the sensitivity of longwave downward radiation (LDR) to changes in water vapor, with the sensitivity being particularly large in many high-elevation regions where the average water vapor is often low. We show that satellite retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Clouds and the Earth's Radiant Energy System (CERES) can be used to expand the current ground-based observational database and that the monthly averaged clear-sky satellite estimates of humidity and LDR are in good agreement with the well-instrumented Center for Snow and Avalanche Studies ground-based site in the southwestern Colorado Rocky Mountains. The relationship between MODIS-retrieved precipitable water vapor and surface specific humidity across the contiguous United States was found to be similar to that previously found for the Alps. More important, we show that satellites capture the nonlinear relationship between LDR and water vapor and confirm that LDR is especially sensitive to changes in water vapor at high elevations in several midlatitude mountain ranges. Because the global population depends on adequate fresh water, much of which has its source in high mountains, it is critically important to understand how climate will change there. We demonstrate that satellites can be used to investigate these feedbacks in high-elevation regions where the coverage of surface-based observations is insufficient to do so.

Naud, Catherine M.↗
NASA NTRS2012

Large-Scale Covariability Between Aerosol and Precipitation Over the 7-SEAS Region: Observations and Simulations

One of the seven scientific areas of interests of the 7-SEAS field campaign is to evaluate the impact of aerosol on cloud and precipitation (http://7-seas.gsfc.nasa.gov). However, large-scale covariability between aerosol, cloud and precipitation is complicated not only by ambient environment and a variety of aerosol effects, but also by effects from rain washout and climate factors. This study characterizes large-scale aerosol-cloud-precipitation covariability through synergy of long-term multi ]sensor satellite observations with model simulations over the 7-SEAS region [10S-30N, 95E-130E]. Results show that climate factors such as ENSO significantly modulate aerosol and precipitation over the region simultaneously. After removal of climate factor effects, aerosol and precipitation are significantly anti-correlated over the southern part of the region, where high aerosols loading is associated with overall reduced total precipitation with intensified rain rates and decreased rain frequency, decreased tropospheric latent heating, suppressed cloud top height and increased outgoing longwave radiation, enhanced clear-sky shortwave TOA flux but reduced all-sky shortwave TOA flux in deep convective regimes; but such covariability becomes less notable over the northern counterpart of the region where low ]level stratus are found. Using CO as a proxy of biomass burning aerosols to minimize the washout effect, large-scale covariability between CO and precipitation was also investigated and similar large-scale covariability observed. Model simulations with NCAR CAM5 were found to show similar effects to observations in the spatio-temporal patterns. Results from both observations and simulations are valuable for improving our understanding of this region's meteorological system and the roles of aerosol within it. Key words: aerosol; precipitation; large-scale covariability; aerosol effects; washout; climate factors; 7- SEAS; CO; CAM5

Huang, Jingfeng↗
NASA NTRS2011

Estimation of Surface Air Temperature Over Central and Eastern Eurasia from MODIS Land Surface Temperature

Surface air temperature (T(sub a)) is a critical variable in the energy and water cycle of the Earth.atmosphere system and is a key input element for hydrology and land surface models. This is a preliminary study to evaluate estimation of T(sub a) from satellite remotely sensed land surface temperature (T(sub s)) by using MODIS-Terra data over two Eurasia regions: northern China and fUSSR. High correlations are observed in both regions between station-measured T(sub a) and MODIS T(sub s). The relationships between the maximum T(sub a) and daytime T(sub s) depend significantly on land cover types, but the minimum T(sub a) and nighttime T(sub s) have little dependence on the land cover types. The largest difference between maximum T(sub a) and daytime T(sub s) appears over the barren and sparsely vegetated area during the summer time. Using a linear regression method, the daily maximum T(sub a) were estimated from 1 km resolution MODIS T(sub s) under clear-sky conditions with coefficients calculated based on land cover types, while the minimum T(sub a) were estimated without considering land cover types. The uncertainty, mean absolute error (MAE), of the estimated maximum T(sub a) varies from 2.4 C over closed shrublands to 3.2 C over grasslands, and the MAE of the estimated minimum Ta is about 3.0 C.

Shen, Suhung↗
NASA NTRS2012

Measuring the Surface Temperature of the Cryosphere using Remote Sensing

A general description of the remote sensing of cryosphere surface temperatures from satellites will be provided. This will give historical information on surface-temperature measurements from space. There will also be a detailed description of measuring the surface temperature of the Greenland Ice Sheet using Moderate-Resolution Imaging Spectroradiometer (MODIS) data which will be the focus of the presentation. Enhanced melting of the Greenland Ice Sheet has been documented in recent literature along with surface-temperature increases measured using infrared satellite data since 1981. Using a recently-developed climate data record, trends in the clear-sky ice-surface temperature (IST) of the Greenland Ice Sheet have been studied using the MODIS IST product. Daily and monthly MODIS ISTs of the Greenland Ice Sheet beginning on 1 March 2000 and continuing through 31 December 2010 are now freely available to download at 6.25-km spatial resolution on a polar stereographic grid. Maps showing the maximum extent of melt for the entire ice sheet and for the six major drainage basins have been developed from the MODIS IST dataset. Twelve-year trends of the duration of the melt season on the ice sheet vary in different drainage basins with some basins melting progressively earlier over the course of the study period. Some (but not all) of the basins also show a progressively-longer duration of melt. The consistency of this IST record, with temperature and melt records from other sources will be discussed.

Hall, Dorothy K.↗
NASA NTRS2012

Cross-calibration of Imaging Sensors using Model-based, SI-traceable Predictions of At-sensor Radiance

Many inter-consistency efforts force empirical agreement between sensors viewing a source nearly coincident in time and geometry that ensures consistency between sensors rather than obtain an SI-traceable calibration with documented error budgets. The method described here provides interconsistency via absolute radiometric calibration with defensible error budget avoiding systematic errors through prediction of at-sensor radiance for a site viewed by multiple sensors but not necessarily viewed at coincident times. The method predicts spectral radiance over a given surface site for arbitrary view and illumination angles and for any date dominated by clear-sky conditions. The foundation is a model-based, SI-traceable prediction of at-sensor radiance over selected sites based on physical understanding of the surface and atmosphere. The calibration of the ground site will include spatial, spectral, and sunview geometric effects based on satellite and ground-based data. The result is an interconsistency of hyperspectral and multispectral sensors spanning spatial resolutions from meters to kilometers all relative to the surface site rather than a single sensor. The source-centric philosophy of calibrating the site inherently accounts for footprint size mismatch, spectral band mismatch, and temporal and spatial sampling effects. The method for characterizing the test site allows its use for SI-traceable calibration of any sensor that can view the calibrated test site. Interconsistency is obtained through the traceability and error budget rather than coincident views. Such an approach to inter-consistency provides better understanding of biases between sensors as well producing more accurate results with documented SI-traceability that reduces the need for overlapping data sets.

Thome, Kurtis↗
NASA NTRS2012

Applications for Near-Real Time Satellite Cloud and Radiation Products

At NASA Langley Research Center, a variety of cloud, clear-sky, and radiation products are being derived at different scales from regional to global using geostationary satellite (GEOSat) and lower Earth-orbiting (LEOSat) imager data. With growing availability, these products are becoming increasingly valuable for weather forecasting and nowcasting. These products include, but are not limited to, cloud-top and base heights, cloud water path and particle size, cloud temperature and phase, surface skin temperature and albedo, and top-of-atmosphere radiation budget. Some of these data products are currently assimilated operationally in a numerical weather prediction model. Others are used unofficially for nowcasting, while testing is underway for other applications. These applications include the use of cloud water path in an NWP model, cloud optical depth for detecting convective initiation in cirrus-filled skies, and aircraft icing condition diagnoses among others. This paper briefly describes a currently operating system that analyzes data from GEOSats around the globe (GOES, Meteosat, MTSAT, FY-2) and LEOSats (AVHRR and MODIS) and makes the products available in near-real time through a variety of media. Current potential future use of these products is discussed.

Minnis, Patrick↗
NASA NTRS2009

Case Studies of Aerosol Remote Sensing in the Vicinity of Clouds

Studying the spatial variability of aerosol properties in the vicinity of clouds is essential to our ability to determine aerosol direct and indirect effects on climate. In this paper, we describe aerosol observations collected near cloud edges by an airborne Sun photometer over dark ocean waters. Focusing on case studies of aerosol measurements near eight cloud edges within a dissipating stratiform cloud deck, we compare the airborne Sun photometer observations to retrievals of aerosol properties using the standard Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol algorithm applied to 500-m-resolution MODIS spectral reflectances. We find a persistent, spectrally neutral increase in the Sun photometer-derived aerosol optical depth (AOD) of up to 10% (0.015) in the 2-km distances closest to the edges of several distinct clouds. At midvisible wavelengths, the MODIS AOD retrievals show similar increases toward cloud edges, although a larger increase in AOD is found in the MODIS along-scan direction. At shortwave infrared (SWIR) wavelengths (1240-2130 nm), the MODIS-derived AOD increases near cloud edges are of the order of 0.03 and as such three times as large as the Sun photometer-derived values. Hence, in contrast to recently discussed bluing of aerosols near cloud edges, i.e., a preferential apparent increase in the visible reflectances of clear-sky pixels due to 3-D radiative transfer effects in the vicinity of clouds, we find a reddening of aerosols in the MODIS 500-m-resolution aerosol retrievals near clouds. This reddening in our study can be traced to larger absolute increases in SWIR reflectances when compared to visible reflectances near clouds, which in turn seem to stem from larger electronic cross talk in the MODIS SWIR bands (5-7). We note that the lack of bluing in our MODIS observations is likely due to the small geometric and optical thicknesses of the clouds considered.

Redemann, Jens↗
NASA NTRS2013

Sensitivity of Downward Longwave Surface Radiation to Moisture and Cloud Changes in a High-elevation Region

Several studies have suggested enhanced rates of warming in high-elevation regions since the latter half of the twentieth century. One of the potential reasons why enhanced rates of warming might occur at high elevations is the nonlinear relationship between downward longwave radiation (DLR) and specific humidity (q). Using ground-based observations at a high-elevation site in southwestern Colorado and coincident satellite-borne cloud retrievals, the sensitivity of DLR to changes in q and cloud properties is examined and quantified using a neural network method. It is also used to explore how the sensitivity of DLR to q (dDLR/dq) is affected by cloud properties. When binned by season, dDLR/dq is maximum in winter and minimum in summer for both clear and cloudy skies. However, the cloudy-sky sensitivities are smaller, primarily because (1) for both clear and cloudy skies dDLR/dq is proportional to 1/q, for q>0.5 g/kg, and (2) the seasonal values of q are on average larger in the cloudy-sky cases than in clear-sky cases. For a given value of q, dDLR/dq is slightly reduced in the presence of clouds and this reduction increases as q increases. In addition, DLR is found to be more sensitive to changes in cloud fraction when cloud fraction is large. In the limit of overcast skies, DLR sensitivity to optical thickness decreases as clouds become more opaque. These results are based on only one high-elevation site, so the conclusions here need to be tested at other high-elevation locations.

cloud physics↗
NASA NTRSDate not supplied

TPSAS-NF1676L-11084-DND

The solar radiation must be considered as the primary factor when we study the climatic system because it provides almost all energy that drives the climatic system and most atmospheric phenomena. The solar radiation at the TOA is partially transferred, transformed into other forms of energy that are eventually dissipated by the general circulation of the atmosphere and oceans, and partially used in chemical and biological processes. The absorption of solar energy in the atmosphere (Acol) plays a key role in the climate system, not only because it affects the amount of solar energy absorbed by the Earth system, but also because it changes the vertical distribution of absorbed energy, which in turn impacts surface evaporation, atmospheric convection and precipitation. To investigate Acol, we have collected the collocated satellite-surface observations over the 35 selected BSRN surface sites during the period 2000-2004. The surface radiation budgets are averaged from the BSRN Shortwave (SW) observations over a 1-hour interval centered at the time of the NASA Terra and Aqua satellites overpass, and the TOA radiation budgets are averaged from the closest FOV satellite data to the BSRN stations. Then the atmospheric column SW absorption will be inferred from the satellite TOA albedo and BSRN surface absorption (through the estimated surface albedo). Monthly means of Acol, TOA albedo, and surface absorption have been calculated under both clear- and all-sky conditions over the 35 selected BSRN sites which represent different climate regimes. To evaluate the GCM simulations, we have collected the ECHAM5-HAM simulated surface, TOA and atmospheric radiation budgets by researchers at ETH Zurich and compared those simulations with observations. The preliminary comparisons have shown that the ECHAM5 simulated clear-sky surface absorption, TOA albedo, and ACOL agree very well (1-3%) with observations. Under all-sky conditions, they have strong correlations with CERES-derived cloud fraction. The simulated surface absorption is lower and TOA albedo is higher than those under clear skies, but ACOL does not increase too much. However, there are relatively large differences over some regions and months. A further study is needed.

X. Dong↗