使用suncalc包getSunlightTimes函数时，北海区域早段日期的night/nightEnd字段返回NA的问题求助

各位好，我在使用R语言的suncalc包处理日照时间数据时遇到了个棘手的问题，想请教下大家：

我有一个名为diel的数据框，日期范围覆盖2023年5月25日到2023年9月16日，里面的经纬度坐标对应北海区域内的5个点位。

我用这段代码来获取相关日照时间：

sun <- getSunlightTimes(data = diel,
                        keep = c("nightEnd", "sunrise", "sunset", "night"), 
                        tz = "UTC")

结果发现个奇怪的情况：2023年5月25日到7月22日这段时间里，nightEnd和night字段的结果全是NA，但sunrise和sunset却能正常返回数值。我还特意测试了nauticalDawn和nauticalDusk参数，这两个并没有出现NA的问题。

我已经试过调整时区，也仔细验证过日期列的格式，但问题始终没解决。有没有朋友遇到过类似情况，或者能帮我分析下原因呢？

备注：内容来源于stack exchange，提问作者AtlasA