The melting of the Greenland ice sheet (GrIS) in summer has rapidly and significantly increased in recent decades, especially for the northern GrIS. Circulation related to GrIS melting is important for understanding the contribution of the GrIS to the global sea level. In this paper, we used the SOM method to obtain three spatial patterns of GrIS melting based on model output data: overall melting, northern melting, and southern melting patterns. We also examined their linkages to the observed atmospheric circulation. GrIS melting is primarily related to Greenland blocking (GB), while different types of GB are responsible for different melting patterns. The overall and northern melting patterns are both associated with high-latitude GB, which is associated with the decay and westward movement of mid-latitude and high-latitude European blocking (EB), respectively. It is found that the generation of mid-latitude and high-latitude EBs are related to mid-latitude eastern North Atlantic warming and Greenland–Norskehavet–North Sea warming, respectively, while the movement of EB is related t upstream weakened zonal wind anomalies. Furthermore, the southern melting pattern is linked to mid-latitude GB, which is associated with the wave train from the East Pacific to Southern Greenland through North America and further related to the mid-latitude of East Pacific warming.

In atmospheric modeling, sub-grid parameterization is an important method for studying the topographic effects of solar radiation using high-resolution Digital Elevation Model (DEM) data. For reducing the amount of computation, some approximate methods that can lead to errors are used in existing sub-grid parameterization schemes for clear-sky direct solar radiation (SPS-CSDSR). The lack of a high-precision fast terrain occlusion algorithm (HPFTOA) remains one of the biggest constraints in this field. This study proposed an HPFTOA. It mainly uses two kinds of acceleration algorithms. One method is to use a dynamic, lossless, and fast occlusion search radius. Another way is to use the rectangular grid for calculations within the accuracy of DEM data to avoid coordinate projection conversions. The test results indicate that the HPFTOA can carry out large-scale computation based on DEM data with a resolution of 90 m. Because it rarely uses approximation algorithms and considers the curvature of the Earth, SPS-CSDSR can achieve unprecedented precision. The HPFTOA can also be used in the fields of mountain solar energy assessment, remote sensing, and telemetry, including terrain-obscuring the probe. As computer performance improves and algorithms and execution code are optimized, the application prospects will be very broad.

Existing sub-grid parameterization schemes for clear-sky direct solar radiation (SPS-CSDSR) assume that the sub-grid cells have the same atmospheric transparency. This study shows that in undulating terrain, significant errors can occur when the sub-grid is in turbid weather or partly above the cloud top. A correction factor was proposed. It can effectively eliminate errors under a cloudless sky and can reduce some errors when part of the sub-grid is above the cloud or fog top. For atmospheric models with high horizontal resolution, example test verification shows that the cast shadowless coverage method can lead to large errors. It should no longer be used based on current computing power. These improvements and the high-precision fast terrain occlusion algorithm in Part I will allow SPS-CSDSR to achieve unprecedented high accuracy. Based on the proposed daily interpolation method, the high-precision SPS-CSDSR is also feasible for regional climate simulation. The analysis pointed out that the sub-grid terrain radiative effect (STRE) is distributed over inclined surfaces with larger areas and at different heights. Existing methods of coupling STRE on one flat surface have certain physical drawbacks. This paper suggests introducing parameterization of STRE at different altitudes and improving the coupling of land–air.