Glacier length is a key morphological element that has many glaciological applications; however, it is often difficult to determine, especially for glaciers that cover larger spatial areas or those that exhibit frequent temporal change. In this paper, we describe a new Arc GIS-based method that can derive glacier flow lines for determining glacier length based on digital elevation model and glacier outlines. This method involves(1) extraction of the highest and lowest points on a glacier,(2) calculation of 10-m contour lines on the glacier from 10 m to 100 m height, and(3) connection of the midpoints of each contour line with the highest and the lowest points in order to create a flow line, which is subsequently smoothed. In order to assess the reliability of this method, we tested the algorithm’s results against flow lines calculated using field measurements, analysing data from the Chinese Glacier Inventory, and manual interpretation. These data showed that the new automated method is effective in deriving glacier flow lines when contour lines are relatively large; in particular, when they are between 70 m and 100 m. Nonetheless, a key limitation of the algorithm is the requirement to automatically delete repeated and closed curves in the pre-treatment processes. In addition to calculating glacier flow lines for derivation of glacier length, this method also can be used to effectively determine glacier terminus change. Glacier length is a key morphological element that has many glaciological applications; however, it is often difficult to determine, especially for glaciers that cover larger spatial areas or those that exhibit appear frequent temporal change. In this paper, we describe a new Arc GIS-based method that can derive glacier flow lines for determining glacier length based on digital elevation model and glacier outlines. This method includes (1) extraction of the highest and lowest points on a glacier, (2) calculation of 10-m contour lines on the glacier from 10 m to 100 m height, and (3) connection of the midpoints of each contour line with the highest and the lowest points in order to create a flow line, which is subsequently smoothed. In order to assess the reliability of this method, we tested the algorithm’s results against flow lines calculated using field measurements, analyzing data from the Chinese Glacier Inventory, and manual interpretation. These data showed that the new automated method i s effective in deriving glacier flow lines when contour lines are relatively large; in particular, when they are between 70 m and 100 m. Nonetheless, a key limitation of the algorithm is the requirement to automatically delete repeated and closed curves in the pre-treatment In addition to calculating glacier flow lines for derivation of glacier length, this method also can be used to be more determined glacier terminus change.