Changes in HIRS Detection of Cloud over Australia from 1985 to 2001

: A long-term archive of cloud properties (cloud top pressure, CTP; and cloud effective emissivity, ε ) determined from High-resolution Infrared Radiation Sounder (HIRS) data is investigated for evidence of regional cloud cover change. In the 17 years between 1985 and 2001, different cloud types are analysed over the Australian region (10 ◦ S–45 ◦ S, 105 ◦ E–160 ◦ E) and areas of change in total cloud frequency examined. Total cloud frequency change over the Australian region between two adjacent eight-year time periods (1994 to 2001 minus 1985 to 1992) shows the largest increases (ranges between 6% and 12%) of average HIRS total cloud cover occurring over the offshore regions to the northwest and northeast of the continent. Over land, the largest reduction of average HIRS total cloud frequency is in the southwestern region of Australia (between 2% and 8%). Through central Australia, there is a 2% to 7% increase in average HIRS total cloud frequency when comparing these eight-year periods. This paper examines the regional cloud changes in 17 years over Australia that are embedded in global cloud statistics. Examining total HIRS cloud cover frequency over Australia and comparing two different eight-year time periods, has highlighted notable areas of average change. Preliminary reporting of satellite-derived HIRS cloud products and Global Precipitation Climatology Project (GPCP) rainfall products during La Niña seasons between 1985 and 2001 has also been undertaken.


Introduction
Since the first cloud-observing satellite was launched in 1959, numerous generations of satellites and instruments have contributed to the creation of large archives of physical and microphysical cloud data. As satellites have performed beyond lifespan specifications and the continuity of subsequently launched platforms, research into clouds and weather systems has expanded from comprehensive short-term analysis to include long-term climate trends.
A number of long-term cloud data archives exist within the Advanced Very High-Resolution Radiometer (AVHRR) Pathfinder Atmospheres-Extended (PATMOS-x) program [1], the International Satellite Cloud Climatology Project (ISCCP) [2], and the University of Wisconsin/National Oceanic and Atmospheric Administration (UW/NOAA) Pathfinder High-resolution Infrared Radiation Sounder (HIRS) data archive [3].
As the spectral range of the respective instruments and the methodologies used to interpret the sensor information can be different, a thorough comparison among different climatologies was reported by the Global Energy and Water Cycle Experiment (GEWEX) Cloud Assessment Project [4]. Along with AVHRR, ISCCP, and HIRS data, results from numerous other data sets such as MODIS (Moderate-resolution Imaging Spectroradiometer), AIRS (Atmospheric Infrared Sounder), CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations), and CloudSat were reported.

Materials and Methods
The HIRS data archive consists of radiometric measurements recorded from four different generations of the HIRS instrument [9] beginning in the late 1970s through to present day. The following cloud analysis focusses on data retrieved from the HIRS/2 instrument that was aboard satellites NOAA-5 through to NOAA-14 (operational during the period 1979 to 2001).
The measurements recorded by the HIRS/2 instrument exist for 20 different spectral bands, 19 within the infrared spectral region (ranging between 3.7 µm and 15 µm) and one in the visible spectrum (at 700 nm) [10].
Five of the narrow sensor spectral bands selected along the transmission/absorption boundary of the 15 µm CO 2 absorption band are able to detect different depths of the atmosphere. The data collected from these infrared bands are then applied through the CO 2 slicing method to derive cloud parameters.
The CO 2 slicing algorithm was theorised by Chahine [11] and Smith et al. [12] with some initial successful validation results of cloud height derived using this method reported by Smith and Platt [13]. Further validation was conducted by Menzel et al. [14] and Frey et al. [15]. The technique has been outlined in detail in Wylie et al. [5] and can be briefly summarised as being based on the radiative transfer equation for a single cloud layer in the atmosphere that aims to obtain measurements of cloud top pressure (CTP) and effective emissivity (ε). When considering observed radiance incident on a satellite sensor, the contributions from a clear field-of-view (FOV) would be divided between radiance from the Earth's surface and radiance from the atmosphere. If a cloud is present, the radiation would come from the Earth's surface, from the atmosphere below the cloud and above the cloud, and also from the cloud itself. By finding the difference between the clear FOV radiance and the cloud observed radiance, the cloud contribution can be determined. In order to reveal clouds at different atmospheric levels, specific combinations of spectrally close bands (between 13 µm and 15 µm and adjacent to the CO 2 absorption band) are used. The HIRS/2 data provided by Wylie [16], generated by the CO 2 slicing algorithm, were the uncorrected average counts of monthly cloud observations for each satellite with a 1 • latitude/longitude resolution. Apart from the total number of cloud observations and the total number of clear sky observations observed during each month, cloud observations were also arranged into combinations of nine different divisions of cloud height and six different divisions of cloud emissivity. Data were initially examined from the ascending orbit satellites (NOAA-5, NOAA-7, NOAA-9, NOAA-11, and NOAA-14) with a local equator crossing time in the early afternoon that offered a continuous data record.
Unexplained performance issues with NOAA-7 during the early 1980s were noted in Wylie et al. [6] and this has led to the HIRS analysis conducted here to focus only on the consecutive data collected between 1985 and 2001 from three satellites, NOAA-9, NOAA-11, and NOAA-14.
When assessing changes in cloud frequency observations and changes in rainfall, the Global Precipitation Climatology Project (GPCP) data set was selected as a precipitation data source with similar spatial and temporal attributes to the available HIRS cloud products. GPCP development is described in Adler et al. [17], Huffman et al. [18], Huffman and Bolvin [19], and Adler et al. [20]. The Version 2.3 monthly precipitation product merges rain-gauge data with satellite-derived rainfall estimated from several different microwave and infrared sensor observations and was sourced online from the NOAA Earth System Research Laboratory [21].
When assessing seasons with above-and below-average rainfall, different cloud types and classifications from El Niño and La Niña seasonal months are investigated to provide a clearer understanding of the type of clouds associated with rainfall.
The El Niño Southern Oscillation (ENSO) is described as a coupled ocean-atmosphere system with two distinct phases that can affect weather patterns in tropical and subtropical climates, the warm phase called El Niño and the cold phase, La Niña [22].
The oceanic aspect of the ENSO system involves sea surface temperature (SST) anomalies. The Oceanic Niño Index (ONI) calculates three-month running means of sea surface temperature (SST) anomalies determined for a particular oceanic region; the Niño 3.4 region (5 • N-5 • S, 120 • W-170 • W) [23]. If the mean anomaly exceeds 0.5 • C for more than five consecutive index values, the contributing months are deemed a warm episode. A cold episode is determined by five continuous index values registering less than −0.5 • C.
The atmospheric part of the ENSO system involves mean sea level pressure (MSLP) differences. The Southern Oscillation Index (SOI) is determined from the mean sea level pressure (MSLP) difference between Tahiti in the central Pacific Ocean and Darwin, Australia (MSLP difference = Tahiti − Darwin) [24]. During an El Niño event, the SOI experiences negative values (higher surface pressure over Darwin) and higher than normal ocean temperatures in the eastern Pacific Ocean. La Niña events have positive values of the SOI (high surface pressure over Tahiti) combined with cooler eastern Pacific Ocean temperatures.
An El Niño event within Australia often means lower than average levels of rainfall while a La Niña event would usually be associated with higher-than-average levels of rainfall. In order to determine the El Niño and La Niña events that occurred between 1985 and 2001, a combination of SOI values and variations in the NCEP/NOAA Oceanic Niño Index (ONI) were assessed. The defined events used in this study are outlined in Table 1. While low Earth orbit platforms target a specific local overpass time upon launch, after many years of operation, the equator crossing times (EXT) of satellites can change. After launch, the early NOAA satellites EXT shifted from their original orbits. Adjusting the satellite radiometric data for orbit drift ensures cloud frequency observations can be compared at the same time of day over extended periods. The HIRS data in this analysis were collected between 1985 and 2001 aboard three afternoon satellites, NOAA-9, NOAA-11, and NOAA-14. Using the platform parameters tabled in Ignatov et al. [25], the drift of each of these satellites was calculated and is shown in Figure 1. compared at the same time of day over extended periods. The HIRS data in this analysis 151 were collected between 1985 and 2001 aboard three afternoon satellites, NOAA-9, NOAA-152 11, and NOAA-14. Using the platform parameters tabled in Ignatov et al. [25], the drift of 153 each of these satellites was calculated and is shown in Figure 1. Similar to the procedure in Wylie et al. [6], scatter plots of cloud frequency and EXT 160 were generated for all combinations of the HIRS cloud categories (high, mid and low; thin, 161 thick and opaque) over the Australian region (10° S -45° S, 105° E -160° E) for the lifetime 162 of each platform. 163 The linear cloud frequency changes per hour of satellite orbit drift were determined 164 from each plot then used to linearly-correct the cloud frequency data from each platform 165 to a 2 pm EXT. The orbit drift adjustment was applied to the HIRS cloud data prior to any 166 analysis being conducted. 167 Overall, the orbit drift corrections had little influence on the summarised Australian 168 cloud frequencies. High cloud frequency observations remained unaltered while there 169 were minor increases calculated in mid-and low-level cloud observations. The main ef-170 fects of the orbit drift corrections were seen in the low opaque cloud observation category 171 where an overall increase of 1% resulted from the orbit drift corrections. A slight increase 172 Similar to the procedure in Wylie et al. [6], scatter plots of cloud frequency and EXT were generated for all combinations of the HIRS cloud categories (high, mid and low; thin, thick and opaque) over the Australian region (10 • S-45 • S, 105 • E-160 • E) for the lifetime of each platform.
The linear cloud frequency changes per hour of satellite orbit drift were determined from each plot then used to linearly-correct the cloud frequency data from each platform to a 2 pm EXT. The orbit drift adjustment was applied to the HIRS cloud data prior to any analysis being conducted.
Overall, the orbit drift corrections had little influence on the summarised Australian cloud frequencies. High cloud frequency observations remained unaltered while there were minor increases calculated in mid-and low-level cloud observations. The main effects of the orbit drift corrections were seen in the low opaque cloud observation category where an overall increase of 1% resulted from the orbit drift corrections. A slight increase was calculated for mid-level cloud leading to a 1% increase due to rounding. This has led to the total cloud frequency observations over Australia increasing by 2% (by rounding) from 56% to 58% after the data was corrected for orbit drift (see Table 2). The HIRS cloud parameters are calculated (via the CO 2 slicing algorithm) using a static CO 2 value of 380 ppm. In reality, CO 2 levels within the atmosphere have changed over time. At the beginning of the HIRS/2 data archive in 1979 the CO 2 value was 335 ppm and by the end of 2001, the CO 2 value had increased to 375 ppm. While the CO 2 difference wasn't corrected for in this study, Wylie et al. [6] discuss the effect that the difference in CO 2 has on the algorithm performance. Wylie suggests that the clouds detected in earlier years are slightly too high (by 13 hPa-50 hPa). Accordingly, any high cloud frequency differences calculated using this data, between two adjacent eight-year periods will therefore be slightly underestimated.

Results
The Wylie data [16] reported 75% cloud coverage globally [6]. When reprocessed for the Southern region of the globe (20 • S-60 • S) and separated into land and ocean statistics, all cloud over land was observed 64% of the time and over ocean, 80% of the time.
Over the Australian region (10 • S-45 • S, 105 • E-160 • E), all cloud was observed over land 58% of the time (see Table 2). This is a reduction of 17% from the global cloud coverage observed over land and ocean, and a reduction of 6% from the Southern region cloud coverage observed over land.
High-level cloud was observed at 25% (43% of the total) over Australia with mid-level cloud being observed 10% of the time (17% of the total) and low-level cloud observed 23% of the time (40% of the total). The high-, mid-, and low-level percentage ratios of total cloud observed over Australia are shown to align more with the GEWEX results (45-15-40) than the elevated mid-level and reduced low-level ratios seen in the HIRS global results (44-24-32).
During the summer months (December, January, and February) there exists an onset and retreat of monsoonal cloud cover over northern Australia. A region of diminished continental cloud cover occurs during the winter months (June, July, and August) over the northern half of the continent and over the north-west ocean region. The area affected contracts during the spring months (September, October, and November) to the north-west coast of Western Australia.
This latter extended period of infrequent cloud cover is consistent with the Gascoyne region of Western Australia being the most cloud-free region in Australia. Research conducted by Udelhofen et al. [26] reported the mid-west region of Western Australia (bounded by the latitudes of 20 • S and 30 • S which incorporates the Gascoyne and Pilbara regions) to have the lowest annual average percentage of cloud cover (33.7% ± 3.7%) when compared to 13 other Australian regions. These results were calculated from 14 years of TOMS (NASA Total Ozone Mapping Spectrometer) data between 1979 and 1992.  Figure 3 illus-228 trates the differences in HIRS total cloud frequency between these two periods by sub-229 tracting the averaged data for 1985 to 1992 from the more recently averaged data for the 230 period 1994 to 2001. An overlaid reference of state boundaries and regional locations 231 across Australia are also included in Figure 3. 232 Figure 3 shows the largest increases in average total cloud cover, from the earlier 233 eight-year period to the more recent eight-year period, occurring over the waters to the 234 northwest and northeast of the Australian continent (between 6% and 12%). Smaller in-235 creases (between 2% and 7%) in average total cloud frequency occur over the central re-236 gion of Australia. 237 Small reductions of average total cloud frequency over time, up to 2% and 3%, are 238 observed over southern South Australia, Tasmania, the southeast corner of Australia and 239 smaller zones within the Northern Territory and Northern Queensland. 240 In Western Australia, there are two distinct areas of total cloud cover change. The 241 largest reduction of average total cloud frequency between the two temporal periods in 242 Figure 3 (shown in blue) is seen to cover a substantial area of southwestern and mid-west-243 ern Australia. The average cloud frequency of total cloud during the latter part of the 244 1990s (between 1994 and 2001) was between 2% and 8% less than the average cloud fre-245 quency observed in the late 1980s (1985 to 1992). The two primarily offshore regions (cen-246 tred on 17° S, 118° E) off the NW coast and off the NE coast (centred on 20° S, 150° E) show 247 an increase in cloud cover frequency of up to 12%. 248 Changes in rainfall across Australia can also be depicted by differences between the 249 same two adjacent eight-year time periods. An attempt was made to determine the types 250 of clouds best associated with geographic differences in rainfall during El Niño and La 251 Niña seasons. Analysing an extended period of HIRS data allows for a comparison between two separate time periods where any systematic trends in cloud frequency might be highlighted. For the 17 years between 1985 and 2001, an initial eight-year period consists of averaged HIRS cloud frequency data between the years of 1985 and 1992. A second eight-year period averages HIRS cloud frequency data between 1994 and 2001. Figure 3 illustrates the differences in HIRS total cloud frequency between these two periods by subtracting the averaged data for 1985 to 1992 from the more recently averaged data for the period 1994 to 2001. An overlaid reference of state boundaries and regional locations across Australia are also included in Figure 3.    Figure 3 shows the largest increases in average total cloud cover, from the earlier eight-year period to the more recent eight-year period, occurring over the waters to the northwest and northeast of the Australian continent (between 6% and 12%). Smaller increases (between 2% and 7%) in average total cloud frequency occur over the central region of Australia.
Small reductions of average total cloud frequency over time, up to 2% and 3%, are observed over southern South Australia, Tasmania, the southeast corner of Australia and smaller zones within the Northern Territory and Northern Queensland.
In Western Australia, there are two distinct areas of total cloud cover change. The largest reduction of average total cloud frequency between the two temporal periods in Figure 3 (shown in blue) is seen to cover a substantial area of southwestern and midwestern Australia. The average cloud frequency of total cloud during the latter part of the 1990s (between 1994 and 2001) was between 2% and 8% less than the average cloud frequency observed in the late 1980s (1985 to 1992). The two primarily offshore regions (centred on 17 • S, 118 • E) off the NW coast and off the NE coast (centred on 20 • S, 150 • E) show an increase in cloud cover frequency of up to 12%.
Changes in rainfall across Australia can also be depicted by differences between the same two adjacent eight-year time periods. An attempt was made to determine the types of clouds best associated with geographic differences in rainfall during El Niño and La Niña seasons. Table 3 examines the observational frequencies of total cloud as well as the frequency of high clouds for the Australian region. Summaries are detailed for three different subsets of data; the entire 17-year archive, El Niño years and La Niña years. Over the entire 17-year archive, cloud of any type is observed 58% of the time over the Australian continent. During El Niño seasons the cloud observation percentage remains the same as the archival value, 58%, while the La Niña season frequency increases to 60%. After expanding the total cloud frequencies into different height classifications, the leading influence behind the increased La Niña total cloud frequency values is seen in the high cloud (CTP < 440 hPa) frequency summaries. Listed in Table 3, the archive high cloud observation percentage is 25%. The La Niña high cloud observation percentage is 27%; 2% higher than the archival and the El Niño high cloud observation percentages.
High clouds are then categorised as thin, thick, or opaque and the calculated frequencies from the HIRS data for El Niño and La Niña seasons are summarised in Table 4. La Niña season high thick clouds are observed 15% of the time, a frequency 3% higher than in the El Niño seasons and 2% higher than the 17-year archive. The differences in high thick clouds between common months within the La Niña and El Niño seasons are illustrated in Figure 4. While El Niño events are associated with lower−than−average rainfall and La Niña events are associated with higher−than−average rainfall, the influence of these ENSO phases, within Australia, centre around the eastern and northern regions of the continent. Higher than average levels of rainfall would be expected in these areas of the country during La Niña events and lower than average rainfall expected during El Niño events. If increasing levels of high thick clouds are associated with La Niña rainfall then a response from the La Niña/El Niño difference on the East coast of Australia should be seen.    Figure 4, Spring and Autumn show the largest spatial areas where this is the case. The western half of Australia also shows positive high thick cloud/La Niña differences during Summer.
One notable negative cloud frequency difference in the seasonal images of Figure 4 is in Autumn (March, April, and May), in the southern half of Western Australia. The ENSO season classification highlighted negative cloud frequency differences in this area (indicating El Niño classified months reported a higher presence of high thick clouds than La Niña classified months). Rainfall in this area is not greatly affected by the ENSO processes [27]. Figures 5 and 6 show the investigation into the changes in average high thick cloud between two selected eight−year periods (average between 1994 and 2001 minus the average between 1985 and 1992) and the differences in average rainfall for the same periods. The rainfall data were sourced from the Global Precipitation Climatology Project (GPCP) Version 2.3 data set.
Remote Sens. 2021, 13, x FOR PEER REVIEW 9 of 12 areas where this is the case. The western half of Australia also shows positive high thick 302 cloud/La Niña differences during Summer. 303 One notable negative cloud frequency difference in the seasonal images of Figure 4 304 is in Autumn (March, April, and May), in the southern half of Western Australia. The 305 ENSO season classification highlighted negative cloud frequency differences in this area 306 (indicating El Niño classified months reported a higher presence of high thick clouds than 307 La Niña classified months). Rainfall in this area is not greatly affected by the ENSO pro-308 cesses [27]. 309 Figures 5 and 6 show the investigation into the changes in average high thick cloud 310 between two selected eight−year periods (average between 1994 and 2001 minus the av-311 erage between 1985 and 1992) and the differences in average rainfall for the same periods. 312 The rainfall data were sourced from the Global Precipitation Climatology Project (GPCP) 313 Version 2.3 data set. Average frequency of high thick cloud cover over the northern parts of the Northern 319 Territory and Queensland and over the northern waters of Australia are 2% to 5% higher 320 in the second temporal period (1994 − 2001) when compared to the first temporal period 321 (1985 − 1992) ( Figure 5). There is also an increase in rainfall in these areas of between 18 322 mm and 30 − 35 mm ( Figure 6). Across the centre of Australia, there is again the E−W band 323 that experiences a 2% to 3% increase in high thick clouds along with between 4 mm and 324 20 mm increase in rainfall. 325 Decreases in the average frequency of high thick clouds (shown in blue in Figure 5,326 where there are more clouds present between 1985 and 1992 than between 1994 and 2001) 327 are found in two main zones in Figure 5. The main zone is on the west coast of the country 328 while lesser decreases are distributed along the full extent of the east coast. 329 The area spanning between the Gascoyne region and the southern coastline of West-330 ern Australia (20° S to 35° S) experiences up to a 3% reduction in high thick cloud. 331 Looking at the same zone in Figure 6, the GPCP rainfall shows a small decrease (ap-332 proximately 3 mm − 4 mm) in rainfall for the southwest corner but a slight increase (ap-333 proximately 6 mm) for the remaining part of the reduced high thick cloud zone. Average frequency of high thick cloud cover over the northern parts of the Northern Territory and Queensland and over the northern waters of Australia are 2% to 5% higher in the second temporal period (1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)) when compared to the first temporal period (1985-1992) ( Figure 5). There is also an increase in rainfall in these areas of between 18 mm and 30-35 mm ( Figure 6). Across the centre of Australia, there is again the E-W band that experiences a 2% to 3% increase in high thick clouds along with between 4 mm and 20 mm increase in rainfall.
Decreases in the average frequency of high thick clouds (shown in blue in Figure 5, where there are more clouds present between 1985 and 1992 than between 1994 and 2001) are found in two main zones in Figure 5. The main zone is on the west coast of the country while lesser decreases are distributed along the full extent of the east coast.
The area spanning between the Gascoyne region and the southern coastline of Western Australia (20 • S to 35 • S) experiences up to a 3% reduction in high thick cloud.
Looking at the same zone in Figure 6, the GPCP rainfall shows a small decrease (approximately 3 mm-4 mm) in rainfall for the southwest corner but a slight increase (approximately 6 mm) for the remaining part of the reduced high thick cloud zone. In Figure 5, the eastern side of Australia contains areas (Tasmania, Victoria, and the 338 southeast corner of Queensland) that show up to a 2% decrease in high thick clouds. Fig-339 ure 6 shows decreases of GPCP rainfall all along the eastern coastline, in some areas up to 340 18 mm. 341 A question remains -how significant are the cloud changes that are reported here? 342 Frey and Menzel [28] reported the mean absolute monthly changes of HIRS−derived, high 343 cloud frequency are "−0.2% (with scatter about those means from 15.5% to 18.4%)… 344 Clearly there is considerable variability in the monthly means from year to year". This 345 implies that when you average over eight years the scatter gets reduced by an order of 346 magnitude, so changes greater than 2% can be considered significant and less than 2% 347 near random. 348

349
The initial motivation for this research was to undertake a study of regional changes 350 in cloud climatology across Australia. 351 Figure 3 in this work reveals that the regional scale processes that we have observed 352 suggest there are significant processes at work that are seriously impacting the regional 353 cloud climatology of Australia and that the direction of trends being observed are sus-354 tained over significant timescales. The magnitudes of these changes, as shown in Figure 355 3, are up to 12%; surprisingly large over this relatively short period of observation. The 356 increasing cloud cover is shown to occur predominantly over oceanic regions in latitudes 357 from ~10° S to ~22° S with the trends for the NW and NE regions well correlated. 358 We find that the cloud cover over the bulk of continental Australia shows (in Figure 359 3) large regions of cloud cover reduction over the 17−year observing period. These areas 360 primarily correspond to Australia's agriculturally productive regions which is a very sig-361 nificant concern if these trends continue. It is, in a sense, anomalous that the only regions 362 of increased continental cloud cover occur over the extensive desert and low productivity 363 Australian rangelands. The extensive SW region of the Australian continent (from 22° S 364 to 35° S by 112° E to 125° E), of area ~1.5 x 10 6 km 2 , shows the largest reduction in cloud 365 cover in the order of ~2% to 8%. The GPCP precipitation data ( Figure 6) shows a small 366 average decrease (~3 mm -4 mm) in rainfall for the southwest corner of this region but a 367 slight increase (~6 mm) for the remaining area of the SW zone. Again, this is a concerning 368 In Figure 5, the eastern side of Australia contains areas (Tasmania, Victoria, and the southeast corner of Queensland) that show up to a 2% decrease in high thick clouds. Figure 6 shows decreases of GPCP rainfall all along the eastern coastline, in some areas up to 18 mm.
A question remains-how significant are the cloud changes that are reported here? Frey and Menzel [28] reported the mean absolute monthly changes of HIRS−derived, high cloud frequency are "−0.2% (with scatter about those means from 15.5% to 18.4%) . . . Clearly there is considerable variability in the monthly means from year to year". This implies that when you average over eight years the scatter gets reduced by an order of magnitude, so changes greater than 2% can be considered significant and less than 2% near random.

Discussion and Conclusions
The initial motivation for this research was to undertake a study of regional changes in cloud climatology across Australia. Figure 3 in this work reveals that the regional scale processes that we have observed suggest there are significant processes at work that are seriously impacting the regional cloud climatology of Australia and that the direction of trends being observed are sustained over significant timescales. The magnitudes of these changes, as shown in Figure 3, are up to 12%; surprisingly large over this relatively short period of observation. The increasing cloud cover is shown to occur predominantly over oceanic regions in latitudes from~10 • S to~22 • S with the trends for the NW and NE regions well correlated.
We find that the cloud cover over the bulk of continental Australia shows (in Figure 3) large regions of cloud cover reduction over the 17-year observing period. These areas primarily correspond to Australia's agriculturally productive regions which is a very significant concern if these trends continue. It is, in a sense, anomalous that the only regions of increased continental cloud cover occur over the extensive desert and low productivity Australian rangelands. The extensive SW region of the Australian continent (from 22 • S to 35 • S by 112 • E to 125 • E), of area~1.5 × 10 6 km 2 , shows the largest reduction in cloud cover in the order of~2% to 8%. The GPCP precipitation data ( Figure 6) shows a small average decrease (~3 mm-4 mm) in rainfall for the southwest corner of this region but a slight increase (~6 mm) for the remaining area of the SW zone. Again, this is a concerning trend since this southwest corner is a major agricultural area typically producing 30% to 40% of the national cereal crop.
The preliminary examination of HIRS high thick cloud cover frequency ( Figure 5) and the GPCP precipitation pattern (Figure 6) shows that a spatial relationship is evident. There are some regions where the links between the selected cloud parameters and rainfall are not clear, and it would suggest that further investigation into other cloud parameters is necessary. The addition of agricultural productivity data would also provide useful information.
In conclusion, we find that cloud cover changes over these 17 years are regionally defined and argue that continuing regional investigations such as that reported here are important given the timescale of change, the magnitudes of change, and the spatial scales of impacted areas.

Future Outlook
There is the prospect of improved cloud diagnostic tools from future on−orbit instruments that should aid an improved understanding of cloud microphysics and associated precipitation. The modelling of clouds and their role in climate change into the future remain a challenging and productive area for research that may identify the sources of the observed changes in cloud cover and give guidance as to what we should expect in the future.  Data Availability Statement: The HIRS data are available from the corresponding author upon request. The GPCP precipitation data are available from the NOAA/OAR/ESRL PSL, Boulder, Colorado, U.S.A. from their website at http://www.psl.noaa.gov/data/gridded/data.gpcp.html.