משוואת זריחת השמש

ממוזער|שמאל|350px|הסקיצה הצבעונית שבאיור ממחישה את מספר השעות של אור היום כפונקציה של קו הרוחב (הציר האנכי) והיום בשנה (הציר האופקי).
באסטרונומיה, **משוואת זריחת השמש** (באנגלית: **Sunrise equation**) מאפשרת לחשב את שעות הזריחה והשקיעה בעבור כל נטייה שמשית וקו רוחב:

:\cos \omega_\circ = -\tan \phi \times \tan \delta 

כאשר:
: \omega_\circ היא זווית השעה בזריחה (כאשר נלקח הערך השלילי) או בשקיעה (כאשר נלקח הערך החיובי);
: \phi היא קו הרוחב של הצופה על כדור הארץ;
: \delta היא זווית הנטייה השמשית.

## תאוריה של המשוואה
### נטייה שמשית
ממוזער|250px|שמאל|הזווית <math>\beta</math> היא הזווית בין הכיוון ממנו מגיעות קרני השמש לקו ייחוס מסוים המונח במישור המילקה.
הנטייה השמשית מוגדרת כזווית המינימלית שיוצר קו הראייה מן הצופה למרכז השמש עם ישר הזנית ועבור צופה בקו המשווה; כיוון שכיוון ציר הסיבוב של כדור הארץ קבוע ביחס למישור המילקה, במהלך מסלולו השנתי סיבוב הרדיוס וקטור המחבר בינו לבין השמש גורם לשינוי הזווית המינימלית הזאת, החל מ-0° (כלומר כאשר השמש יכולה להיות בזנית) ביום השוויון ועד לזוויות המרביות של \pm 23.44^\circ שמתקבלות ביום ההיפוך החורפי ויום ההיפוך הקיצי, בהמיספירה הצפונית והדרומית בהתאמה. אפקט זה של שינוי הנטייה השמשית הוא שגורם לשינויים בהספק שבו מחממת השמש את הקרקע, מה שמחולל את עונות השנה.

ביממה נתונה, הזווית המרבית של השמש ביחס לאופק מתקבלת כאשר הרדיוס וקטור (המחבר בין השמש לארץ), כיוון ציר הסיבוב וכיוון ישר הזנית נחים כולם על מישור אחד. לכן, כדי לקבל את הביטוי לזווית הנטייה השמשית כפונקציה של הזווית \beta של הרדיוס וקטור במישור המילקה וזווית נטיית כדור הארץ (ממנו ניתן גם לחשב את הנטייה השמשית כתלות בתאריך) יש להיעזר במכפלה סקלרית של וקטור הכיוון של ציר כדור הארץ והרדיוס וקטור שבמישור המילקה, מה שנותן:

: \cos(\pi/2 - \delta) = \cos\beta \cdot \cos (\pi/2 - 23.44^\circ) 

או:

:  \sin \delta = \sin\beta \cdot \sin ( 23.44^\circ) 

כאשר \beta מכוילת כך שתהיה שווה לאפס ביום השוויון.

### גזירת המשוואה
זווית הנטייה השמשית הרגעית \delta מאפשרת לחלק את כדור הארץ לשלושה אזורים: קווי הרוחב \pi/2 - \delta \le \phi \le \pi/2 (בהמיספירה הצפונית), קווי הרוחב "האמצעיים"  -(\pi/2 - \delta) , וקווי הרוחב הדרומיים: -\pi/2 \le \phi \le -(\pi/2 - \delta).
בהתאם לעונת השנה, כל אחד מתחומי קווי הרוחב "הקיצוניים" מתאימים לאזורים על כדור הארץ בהם יש לילה לאורך כל היממה (השמש לא זורחת) או יש יום לאורך כל היממה (השמש לא שוקעת), כאשר המצב מתהפך בין שני אזורים קיצוניים אלו.
שמאל|200px|ממוזער|באיור מוראה החתך המעגלי המתאים לקו הרוחב <math>\phi</math> (המעגל הגדול), בעוד המעגל הקטן שמוכל בו הוא חיתוך החרוט המתאים לקו הרוחב <math>\pi/2 - \delta</math> עם המישור היוצר את חתך קו הרוחב <math>\phi</math>. הקשת של האזור הלבן (האזור שמימין לקו הישר ובתוך המעגל הגדול) היא הקשת שמתאימה לחצי הכדור המואר (היכן שיש יום), בעוד הקשת של שאר המעגל הגדול (המיוצג בהצללה קלה) מתאימה לחצי הכדור החשוך.
כדי לחשב את השעות המדויקת בהן מתרחשת השקיעה או הזריחה בקו רוחב \phi, יש להסתייע בדמיון מרחבי כדי להבין איך נראה המצב ההדדי של המעגל הגדול שמפריד בין חצי הכדור המואר לחצי הכדור החשוך והמעגל הקטן המתאים לקו הרוחב \phi. שני המעגלים הללו נחתכים בשתי נקודות, אשר מחלקות את המעגל הקטן לשתי קשתות: האחת מתאימה ליום והשנייה ללילה - כלומר המיקום הזוויתי של שתי הנקודות הללו הוא בדיוק זווית השעה המתאימה לזריחה ולשקיעה, בהתאמה. הזווית המרכזית של הקשת המתאימה ליום תהיה שווה ל-2\omega_{o}, ולפיכך כדי למצאה יש להסתייע בטריגונומטריה כדי לחשב את אורך הקטע *CP* שבאיור.

ברור שרדיוס המעגל הגדול שבאיור, המתאים לקו הרוחב \phi, הוא  r_{large} = R\cdot \cos\phi, כאשר *R* הוא רדיוס כדור הארץ. לעומת זאת, מן הגאומטריה של הטלות גרפיות מקבלים, שרדיוס המעגל הקטן שבאיור שווה לרדיוס החיתוך של החרוט המתאים לקו הרוחב \pi/2 - \delta עם המישור היוצר את חתך קו הרוחב \phi. הרדיוס של קו הרוחב \pi/2 - \delta הוא: R\cdot \sin \delta, ולכן, מדמיון משולשים מקבלים, שרדיוס המעגל הפנימי הוא:  CP = r_{small} = \frac .

מכאן נקבל, שקוסינוס זווית השעה המתאימה לשקיעה או לזריחה (כלומר \omega_{0}) הוא:

\cos \omega_{0} = \frac } = \frac  = \tan\phi \cdot \tan\delta

כאשר סימן המינוס במשוואה שבתחילת הערך הוא עניין של קונבנציה בלבד. נוסחה זאת ישימה רק לקווי הרוחב שבין החוג הארקטי והחוג האנטארקטי, שכן מצפון לראשון ומדרום לשני יש תמיד לפחות יום אחד בשנה ללא שקיעה או ללא זריחה.

באופן כללי ניתן להראות שהגובה הזוויתי \theta של השמש מעל לאופק כפונקציה של זווית השעה ביממה \omega מקיים:

\sin\theta = -\sin\delta\cdot \sin\phi + \cos\phi\cdot \cos\omega\cdot \cos \delta

להסבר מורחב על יישומי נוסחה כללית יותר זו ראו גם זווית הזניט השמשית . נוסחה זאת משמשת בין היתר כדי לקבל את שעות הזריחה והשקיעה המתוקנות (ראו פרק המשוואה המוכללת למטה).

### מסקנות מן הנוסחה
* מן הנוסחה עולה, למשל, שבקו המשווה השקיעה והזריחה חלות כל השנה באותה שעה.

## המשוואה המוכללת
משוואת זריחת השמש הסטנדרטית מזניחה את ההשפעה של שבירה אטמוספירית (אשר מרימה את הדיסקה השמשית בערך ב-0.6° כאשר היא על האופק) ואת קיומו של גודל זוויתי לדיסקה הזו (קוטר זוויתי של בערך 0.5°). זמני הזריחה והשקיעה המתוקנים באלמנכים אסטרונומיים ניתנים לחישוב מהמשוואה הכללית יותר:

:\cos \omega_\circ = \dfrac{\sin a - \sin \phi \times \sin \delta}{\cos \phi \times \cos \delta }

כאשר (a) מייצג את הגובה הזוויתי של מרכז הדיסקה השמשית מתחת לאופק ברגע האחרון של השקיעה, ויש להציב בו את הערך 0.83°- (או 50- דקות קשת).

## שימושים
ממוזער|שמאל|300px|מדידות שונות של השקיעה והזריחה אפשרו לימאים להתמצא באוקיינוס הפתוח.
כמו נוסחאות שימושיות רבות אחרות באסטרונומיה, המשוואה יכולה לשמש למשל, לצורכי ניווט ימי. עקרונית ניתן להיעזר בנתוני מיקום ונסיקה כלליים של השמש ולאו דווקא בערכים שלהם כאשר השמש נמצאת על האופק (זאת בעזרת שימוש בנוסחה הכללית יותר למיקום השמש לפי הזמן) אולם כאן ניצמד לתהליך המתמטי שעושה שימוש ב"משוואת זריחת השמש".

בעבר הלא רחוק, כאשר מלחים נזדקקו לקבוע את מיקומם הגאוגרפי בים, הם יכלו להשתמש בקריאתו של שעון מדויק (כרונומטר) כדי למדוד מדדים זמניים שונים הקשורים בשקיעה ובזריחת השמש, וכך לדעת, בעזרתן של טבלאות נתונים מדויקים שחושבו ישירות מנוסחאות אסטרונומיות, את קו הרוחב בו הם נמצאים. עם זאת, משוואת זריחת השמש בצורתה הישירה **אינה** יכולה לשמש כדי לקבוע את קו הרוחב הגאוגרפי, שכן השעון שעל הספינה מראה את השעה T בקו האורך המקורי אליו כוון, וכיוון שקו האורך הנוכחי גם הוא אינו ידוע, לא ניתן לדעת את השעה המקומית.

בפועל, כדי למצוא את קו הרוחב יש להשתמש בשעון כדי למדוד את **קצב ** שינוי הגובה הזוויתי של השמש מעל לאופק, ולא כדי למדוד את השעה. מתמטית, צריך להשתמש במשוואה הכללית יותר לגובה הזוויתי של השמש מעל לאופק כפונקציה של השעה, ולגזור אותה לפי הזמן. כך מקבלים משוואה של קצב הנסיקה הזוויתי \dot \theta(t) של השמש בזריחה, אשר שילובה עם משוואת זריחת השמש מאפשר להרכיב שלוש משוואות עם נעלמים שהם קו הרוחב \phi, הנטייה השמשית \delta ושעת הזריחה \omega_0 בשעון המקומי. כאשר השמש על האופק הזווית \theta  קטנה ולכן \sin\theta\approx \theta , ומכאן מקבלים שקצב הנסיקה הזוויתי בזריחה הוא:

: \cos\phi \cdot \sin\omega_0 = -\frac 

ושתי המשוואות האחרות הן:

: \cos\omega_0\cdot \cot\phi = \tan\delta 
: \delta = \phi - \theta_{max} + 90

כאן \alpha = 15^\circ/hour הוא קצב הסיבוב הזוויתי של כדור הארץ סביב צירו, ו-\theta_{max} הוא הגובה הזוויתי המרבי של השמש בשמיים (שנמדד בעזרת סקסטנט). לסיכום, מדידת קצב הנסיקה הזוויתי של השמש בזריחה מאפשרת להסיק את שלושת הפרמטרים הנעלמים, בעוד השוואת שעת הזריחה בשעון המקומי עם קריאת הכרונומטר מאפשרת לקבוע גם את קו האורך המקומי l, בהתאם לקשר: l = l_0 + \alpha \cdot(\omega_0 - T) .

## ראו גם
* אנלמה

## קישורים חיצוניים

## הערות שוליים

קטגוריה:אסטרונומיה